On the Use of Entropy Issues to Evaluate and Control the Transients in Some Epidemic Models

Abstract

This paper studies the representation of a general epidemic model by means of a first-order differential equation with a time-varying log-normal type coefficient. Then the generalization of the first-order differential system to epidemic models with more subpopulations is focused on by introducing the inter-subpopulations dynamics couplings and the control interventions information through the mentioned time-varying coefficient which drives the basic differential equation model. It is considered a relevant tool the control intervention of the infection along its transient to fight more efficiently against a potential initial exploding transmission. The study is based on the fact that the disease-free and endemic equilibrium points and their stability properties depend on the concrete parameterization while they admit a certain design monitoring by the choice of the control and treatment gains and the use of feedback information in the corresponding control interventions. Therefore, special attention is paid to the evolution transients of the infection curve, rather than to the equilibrium points, in terms of the time instants of its first relative maximum towards its previous inflection time instant. Such relevant time instants are evaluated via the calculation of an “ad hoc” Shannon’s entropy. Analytical and numerical examples are included in the study in order to evaluate the study and its conclusions.

1. Introduction

Some classical works by Boltzmann, Gibbs and Maxwell have defined entropy under a statistical framework. A useful entropy concept is the Shannon entropy since it is a basic tool to quantify the amount of uncertainty in many kinds of physical or biological processes [1,2,3,4,5,6]. It may be interpreted as a quantification of information loss [1,2,3,7,8,9]. On the other hand, entropy-based tools have been also proposed to evaluate the propagation of epidemics and related public control interventions (see, for instance, [10,11,12,13,14,15,16,17] and some of the references therein). There are also models whose basic framework relies on the use of entropy tools, as for instance [13,14,15,16]. It can be also pointed out that the control designs might be incorporated to some epidemic propagation and other biological problems, see, for instance, [18,19,20,21,22,23,24,25,26,27], and, in particular, for the synthesis of decentralized control in patchy (or network node-based) interlaced environments [24,27]. A typical situation is that of several towns each with its own health center, whose susceptible and infectious populations, apart from their coupled self-dynamics among their integrating subpopulations, might also mutually interact with the subpopulations of the neighboring nodes through in-coming and out-coming fluxes.

It can be pointed out that the knowledge or estimation of the transient behavior of the infection is very relevant for the hospital management of the disease since it is necessary to manage the availability of beds and other sanitary utensils and sanitary means, in general. The work by Wang et al. in [11] pays mainly attention to the description of the transient behavior of the evolution of epidemics rather than to the equilibrium states. The main purpose in that paper was to formulate the time interval occurring between the time instant of the maximum of the infection, which gives a relative maximum of the infection evolution through time (and which zeroes the first time-derivative of the infection function), and the time instant giving its previous inflection time instant. It turns out that the knowledge of the first part of the transient evolution is very relevant to fight against the initial exploding of the illness since any eventual control intervention is typically much more efficient as far as it is taken as quickly as possible. The model proposed in [11] is a time-varying differential equation of first-order describing the infectious population which is the unique explicit one in the model. It is also pointed out in that paper that the time-varying coefficient might potentially contain the supplementary environment information to make such an equation well-posed to practically describe a concrete disease evolution. An interesting point of that work is that the infection evolution is identified with a log-normal distribution whose parameterization is selected in such a way that the entropy production rate is maximized. The above proposed theoretical first-order model has been proved to be very efficient to describe the data of SARS 2003. Alternative interpretations of the entropy in terms of maximum entropy or maximum entropy rate are given, for instance, in [12,13,14] and some references therein.

This paper studies how to link the extension of the first-order differential system proposed in [11] for the study of infection propagations to epidemic models with more integrated coupled subpopulations (such as susceptible, immune, vaccinated etc.) by introducing the coupling and control information through the time-varying coefficient which drives the basic differential equation model. It is considered relevant the control of the infection along its transient to fight more efficiently against a potential initial exploding transmission. Note that the disease-free and endemic equilibrium points and their stability properties depend on the concrete parameterization while they admit a certain design monitoring by the choice of the control and treatment gains and the use of feedback information in the corresponding controls. See, for instance [19,27]. Therefore, special attention is paid to the transients of the infection curve evolution in terms of the time instants of its first relative maximum towards its previous inflection time instant since there is a certain gap in the background literature concerning the study of such transients. The ratio of such time instants is later on considered subject to some worst-case uncertainty relations via the calculation and analysis of an “ad hoc” Shannon’s entropy. Note that entropy issues have been considered in the study of biological, evolution and epidemic models by incorporating techniques of information theory. See, for instance [11,12,13,28,29,30,31,32]. It is well-known that the entropy production theorems might be classified according to a generalized sequence of stable thermodynamic states. Also, the thermodynamic equilibrium, which is characterized by the absence of gradients of state or kinematic variables, is in a state of maximum entropy and zero entropy production [33,34]. Furthermore, linear non-equilibrium processes are associated with entropy production so that the entropy concept may be also invoked in transient processes [35]. On the other hand, it may be pointed out that uncertainties can appear in the characterization of the infection evolution through time, even in deterministic models, due to parameterization uncertainties, fluxes of populations or existing uncertainties in the initial conditions. Other mathematical techniques of interest which combine analytical and numerical issues have been also been applied to the analysis and discussion of epidemic models with eventual support of mathematical techniques on homotopy analysis and distribution functions as, forinstance, the log-normal distribution [36,37]. For instance, in [38], the SIR and SIS epidemic models are solved through the homotopy analysis method. A one-parameter family of series solutions is obtained which gives a method to ensure convergent series solutions for those kinds of models. On the other hand, in [39], the analytic solutions of an SIR epidemic model are investigated in parametric form. It is also found that the generalization of a SIR model including births and mortality with vital dynamics might be reduced to an Abel-type which greatly simplify the analysis.

The paper is organized as follows: Section 2 gives an extension of the basic model of [11] to be then compared in subsequent sections with some existing models with several subpopulations. Such a model only considers the infection evolution through time and it is based on the action of two auxiliary non-negative functions which define appropriately the time-varying coefficient which defines the first-order differential equation of the infection evolution. The model includes, as particular case, that of the abovementioned reference where both such auxiliary functions are identical to the time argument. Particular choices of those functions make it possible to consider alternative effects linked to the basic model like, for instance, the influence on the infectious subpopulation of other coupled subpopulations in more general models like, for instance, the susceptible, exposed, recovered or vaccinated ones. It is also possible to include the control effects through such a varying coefficient, if any, like for instance, the vaccination and treatment controls. Some basic formal results are stated and proved mainly concerning with the first relative maxima and inflection time instants of the infection curve through time. The above two time instants are relevant to take appropriate control interventions to fight against an initially exploding infectious disease.

Section 3 links the basic model of Section 2 with some known epidemic models which integrate more subpopulations than just the infectious one, like for instance, the susceptible and recovered subpopulations, The time-varying coefficient driving the infection evolution is defined explicitly for each of the discussed epidemic models. Basically, it is taken in mind that some relevant information of higher-order differential epidemic models concerning the transient trajectory solution can be captured by a parameter-dependent and time-varying coefficient which drives a first-order differential equation to the light of the basic model of Section 2. So, the time-varying coefficient describing the infection evolution depends in those cases of the remaining subpopulations integrated in the model. The maximum and inflection time instants are characterized for some given examples involving epidemic models of several subpopulations. In particular, the last one of the discussed theoretical examples includes the effects of vaccination and treatment intervention controls generated by linear feedback of the susceptible and infectious subpopulations, respectively. Later on, Section 4 investigates the entropy associated with the infection accordingly to the generalizations of Section 2 concerning the specific structure of the time-varying coefficient describing the infection dynamics and its links with the theoretical examples discussed in Section 3. The error of the entropy related to the reference one associated with the log-normal distribution is estimated. In practice, that property can be interpreted in terms of public medical and social interventions which control the disease propagation when introducing the controls of the last example discussed in Section 3. The second part of Section 4 is devoted to linking the entropy and inflection and maximum infection time instants and their reached values of the discussed multi-population structures to their counterparts of the maximum dissipation rate being associated to the formulation of a simpler model based on the log-normal distribution and one-dimensional infection dynamics. Some numerical tests are performed for comparisons of the entropies and its width of the basic model with two of the discussed examples in the previous sections which involve the presence of more than one integrated subpopulations. Finally, conclusions end the paper.

Notation

$${\mathit{R}}_{+}=\left\{\hspace{0.17em}r\in \mathit{R}\hspace{0.17em}:\hspace{0.17em}r>0\right\}; {\mathit{R}}_{0+}=\left\{\hspace{0.17em}r\in \mathit{R}\hspace{0.17em}:\hspace{0.17em}r\ge 0\right\}={\mathit{R}}_{+}\cup \left\{0\right\}$$

$${\mathit{Z}}_{+}=\left\{\hspace{0.17em}r\in \mathit{Z}\hspace{0.17em}:\hspace{0.17em}r>0\right\}; {\mathit{Z}}_{0+}=\left\{\hspace{0.17em}r\in \mathit{Z}\hspace{0.17em}:\hspace{0.17em}r\ge 0\right\}={\mathit{Z}}_{+}\cup \left\{0\right\}$$

$$\overline{n}=\left\{0,1,\cdots ,n\right\}$$

2. The Basic Model Description and Some Related Technical Results

Since disease propagation can be interpreted as a thermodynamic system, it can be assumed that the rate of increase or decrease is proportional to the infection at the previous day following the approach of modelling the rate of chemical reactions, [11]. Thus, assume that the infection evolution obeys the following time-varying differential equation:

$$\dot{I}\left(t\right)=\alpha \left(t\right)I\left(t\right); I\left(0\right)=I_0>0$$

(1)

where $\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is continuous and time differentiable on $\left(0,\hspace{0.17em}+\infty \right)$. The particular structure of the varying coefficient $\alpha \left(t\right)$ depends on the balances between the spreading mechanism and the exerted controls during the public intervention. Such a coefficient contains the available information related to the incorporation of all the control mechanisms and the coupling dynamics between the infectious populations and the remaining interacting ones such as the susceptible, immune or vaccinated ones. By taking time-derivatives with respect to time in (1), one gets:

$$\begin{array}{c}\ddot{I}\left(t\right)=\dot{\alpha }\left(t\right)I\left(t\right)+\alpha \left(t\right)\dot{I}\left(t\right)\\ =\left(\dot{\alpha }\left(t\right)/\alpha \left(t\right)+\alpha \left(t\right)\right)\dot{I}\left(t\right)\\ =\left(\dot{\alpha }\left(t\right)+{\alpha }^2\left(t\right)\right)I\left(t\right); \dot{I}\left(0\right)=\hspace{0.17em}{\dot{I}}_0=\alpha \left(0\right)I_0\end{array}$$

(2)

It is proposed in [11] to consider two relevant time instants in the disease evolution, namely:

(1)

The inflection time instant of I(t) which is the date in the infection evolution at which the controlling actions take effect on the evolution. Typically, this time instant is the undulation point date in the evolution of $I\left(t\right)$, that is the zero of $\ddot{I}\left(t\right)$, provided that the first non-zero derivative of $I^{\left(n\right)}\left(t\right)=\frac{d^{n}I\left(t\right)}{d{t}^n}$, $n\ge 3$ occurs for some even $n$ since this last condition ensures that the undulation time instant is the inflection time instant.

(2)

The critical time instant at which the spread rate turns from initial growing to decrease which can be empirically attributed to the global influence of the control interventions. This time instant is a relative maximum of I(t) and it satisfies the constraints $\dot{I}\left(t\right)=0$ and $\ddot{I}\left(t\right)<0$ under the reasonable assumption that ${\dot{I}}_0>0$.

It turns out that, along the whole disease evolution, several successive inflection points and relative maxima can happen. The subsequent result which is concerned with the non-negativity, boundedness and asymptotic vanishing property of the infection as time tends to infinity and its two first- time derivatives is immediate from the above expressions (1) and (2):

Theorem 1.

The following properties hold:

(i)

The infection population and its two first-time derivatives obey the following time evolution equations:

$$I\left(t\right)=e^{{\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }}{I}_0;\dot{I}\left(t\right)=\alpha \left(t\right)e^{{\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }}{I}_0;\ddot{I}\left(t\right)=\left(\dot{\alpha }\left(t\right)+{\alpha }^2\left(t\right)\right)e^{{\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }}{I}_0;\forall t\in {\mathit{R}}_{0+}$$

(3)

(ii)

$I\left(t\right)>0$;$\forall t\in {\mathit{R}}_{0+}$if and only if$I_0\ge 0$; and$I\left(t\right)=0$;$\forall t\in {\mathit{R}}_{0+}$if and only if$I_0=0$.

(iii)

If$+\infty >I_0\ge 0$then$I\left(t\right)\le K{I}_0<+\infty$;$\forall t\in {\mathit{R}}_{0+}$for some$K\in {\mathit{R}}_{+}$if and only if$\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$is such that${\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }\le K<+\infty$;$\forall t\in {\mathit{R}}_{0+}$.

(iv)

$I\left(t\right)\to 0$as$t\to +\infty$for any given finite$I_0$if and only if$\underset{t\to +\infty }{lim}{\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }=-\infty$.

(v)

If$+\infty >I_0\ge 0$and${\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }\le K<+\infty$;$\forall t\in {\mathit{R}}_{0+}$for some$K\in {\mathit{R}}_{0+}$then$\left|\dot{I}\left(t\right)\right|<+\infty$;$\forall t\in {\mathit{R}}_{0+}$if and only if, for some$K_1\in {\mathit{R}}_{+}$,$\left|\alpha \left(t\right)\right|\le K_1<+\infty$;$\forall t\in {\mathit{R}}_{0+}$. If$+\infty >I_0\ge 0$then$\left|\dot{I}\left(t\right)\right|<+\infty$;$\forall t\in {\mathit{R}}_{0+}$if and only if$\left|\alpha \left(t\right)\right|e^{{\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }}\le K_2<+\infty$;$\forall t\in {\mathit{R}}_{0+}$, for some$K_2\in {\mathit{R}}_{+}$provided that$\alpha \left(0\right)<+\infty$.

(vi)

$\dot{I}\left(t\right)\to 0$as$t\to +\infty$for any given finite$I_0$if and only if$\underset{t\to +\infty }{lim}\left(\alpha \left(t\right)e^{{\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }}\right)=0$.

If$\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$is bounded and$I\left(t\right)\to 0$as$t\to +\infty$then$\dot{I}\left(t\right)\to 0$as$t\to +\infty$.

(vii)

If$+\infty >I_0\ge 0$then$\left|\ddot{I}\left(t\right)\right|<+\infty$;$\forall t\in {\mathit{R}}_{+}$if and only if$\left|\dot{\alpha }\left(t\right)+{\alpha }^2\left(t\right)\right|e^{{\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }}\le K_3<+\infty$;$\forall t\in {\mathit{R}}_{0+}$, for some$K_3\in {\mathit{R}}_{+}$.$\ddot{I}\left(t\right)\to 0$as$t\to +\infty$for any given finite$I_0$if and only if$\underset{t\to +\infty }{lim}\left(\left(\dot{\alpha }\left(t\right)+{\alpha }^2\left(t\right)\right)e^{{\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }}\right)=0$.

If$\left(\alpha +\dot{\alpha }\right):{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$is bounded and$I\left(t\right)\to 0$as$t\to +\infty$then$\ddot{I}\left(t\right)\to 0$as$t\to +\infty$.

Note that $\alpha \left(t\right)$ (respectively, $\alpha \left(t\right)+\dot{\alpha }\left(t\right)$) is infinity at $t=0$ while it is bounded for $t>0$, as it happens for instance with the $\alpha$—function proposed in [11], then $\dot{I}\left(t\right)$ (respectively, $\ddot{I}\left(t\right)$) is still bounded under the conditions of Theorem 1 (v) (respectively, Theorem 1 (vii)) on ${\mathit{R}}_{+}$.

Note also that the vanishing infection condition of Theorem 1 typically occurs under convergence of the solution to the disease-free equilibrium point if the disease reproduction number is less than one [19,22,23,24,27,29,30,36]. However, it can happen that the infection oscillates around some stable equilibrium or that it converges to a nonzero positive constant defining the corresponding component of the endemic equilibrium steady-state as it is discussed in the next result.

Corollary 1.

The following properties hold:

(i)

Assume that there exists some$C\in {\mathit{R}}_{+}$such that${\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }\to C$as$t\to +\infty$and that$\alpha \left(t\right),\hspace{0.17em}\dot{\alpha }\left(t\right)\to 0$as$t\to +\infty$. Then,$I\left(t\right)\to e^C{I}_0$,$\dot{I}\left(t\right)\to 0$and$\ddot{I}\left(t\right)\to 0$as$t\to +\infty$.

(ii)

Assume that${\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }\to C$as$t\to +\infty$and that$\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$is uniformly continuous. Then,$\alpha \left(t\right)\hspace{0.17em}\to 0$,$I\left(t\right)\to e^C{I}_0$and$\dot{I}\left(t\right)\to 0$as$t\to +\infty$. Assume, in addition, that$\dot{\alpha }:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$is uniformly continuous. Then$\ddot{I}\left(t\right)\to 0$as$t\to +\infty$.

Proof of Property (i).

Follows directly from (1)–(3). On the other hand, since $\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is uniformly continuous and the limit $\underset{t\to +\infty }{lim}{\displaystyle {\int }_0^t\alpha \left(\tau \right)d\tau }=C$ exists and it is finite then $\alpha \left(t\right)\to 0$ as $t\to +\infty$ (Barbalat´s Lemma) and $I\left(t\right)\to e^C{I}_0$ as $t\to +\infty$ from (3), $\dot{I}:{\mathit{R}}_{+}\to {\mathit{R}}_{0+}$ is bounded, since being continuous, it cannot diverge in finite time, and $\dot{I}\left(t\right)\to 0$ as $t\to +\infty$ from (1). If, furthermore, $\dot{\alpha }:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is uniformly continuous and, since $\underset{t\to +\infty }{lim}{\displaystyle {\int }_0^t\dot{\alpha }\left(\tau \right)d\tau }=\underset{t\to +\infty }{lim}\alpha \left(t\right)-\alpha \left(0\right)=-\alpha \left(0\right)$ then $\hspace{0.17em}\dot{\alpha }\left(t\right)\to 0$ as $t\to +\infty$ (again from Barbalat´s Lemma). Since $\alpha \left(t\right),\hspace{0.17em}\dot{\alpha }\left(t\right)\to 0$ as $t\to +\infty$ then $\ddot{I}\left(t\right)\to 0$ as $t\to +\infty$ from (2). □

Let us introduce the following definitions and lemma of usefulness for the proof of the subsequent theorem [36]:

Definition 1.

Let$f:\mathit{R}\to \mathit{R}$be everywhere continuous and twice differentiable at$t_0\in \mathit{R}$. Then,$t_0$is an undulation point (or pre-inflection point) of$f$if$\ddot{f}\left(t_0\right)=0$.

Inflection points of the continuous and twice-differentiable$f:\mathit{R}\to \mathit{R}$are the undulation points of the function where the curvature changes its sign, that is, points of change of local convexity to local concavity or vice-versa. They are also the isolated extrema of$\dot{f}:\mathit{R}\to \mathit{R}$. A well-known technical definition and a related result on inflection points follow:

Definition 2.

Let$f:\mathit{R}\to \mathit{R}$be everywhere continuous and twice differentiable at$t_0\in \mathit{R}$which is an isolated extremum of$f$(that is, a local maximum or minimum, and also an undulation point of,$f$as a result).

Lemma 1.

The following properties hold:

(i) 

Let$f:\mathit{R}\to \mathit{R}$be everywhere continuous and twice differentiable at$t_0\in \mathit{R}$. Then,$t_0$is an inflection point of$f$if$\ddot{f}\left(t+\epsilon \right)\ddot{f}\left(t-\epsilon \right)<0$for some sufficiently small$\epsilon \in {\mathit{R}}_{+}$.

(ii) 

Let$f:\mathit{R}\to \mathit{R}$be everywhere continuous and an odd number$k\left(\ge 3\right)$-times differentiable, within a neighborhood of$t_0\in \mathit{R}$which is an undulation point of$f$satisfying$f^{\left(j\right)}\left(t_0\right)=0$for$j=2,3,\dots k-1$and$f^{\left(k\right)}\left(t_0\right)\ne 0$. Then,$t_0$is an inflection point of$f$.

The subsequent result has a very technical proof leading to the basic result that the zeros at finite time instants of $\dot{I}\left(t\right)$ and $\ddot{I}\left(t\right)$ alternate if I(t) is sufficiently smooth and $\alpha \left(t\right)$ is sufficiently smooth. In order to simplify the result proof, it is assumed, with no loss in generality, that the disease dynamics (1)–(2) has no equilibrium points such that the zeros under study are isolated.

Theorem 2.

Assume that the function$\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$defined by$\alpha \left(t\right)=-\frac{cln\left(g\left(t\right)/E\right)}{h\left(t\right)}$, where$c\hspace{0.17em},\hspace{0.17em}E\in {\mathit{R}}_{+}$and$g,h:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$are everywhere continuous and time-differentiable such that$g\left(0\right)=0$with$\underset{t\to 0}{lim}\frac{ln\left(g\left(t\right)/E\right)}{h\left(t\right)}\le -\epsilon$for some$\epsilon \in {\mathit{R}}_{0+}$, and furthermore,$\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$fulfills the constraints:

$$\alpha \left(D_i\right)=0; \dot{\alpha }\left(L_i\right)=-{\alpha }^2\left(L_i\right)$$

(4)

$$\frac{h\left(L_i\right)\dot{g}\left(L_i\right)-ln\left(g\left(L_i\right)/E\right)\dot{h}\left(L_i\right)g\left(L_i\right)}{g\left(L_i\right){ln}^2\left(g\left(L_i\right)/E\right)}=K>0; \forall L_i\in L_S\cap \left[0,\hspace{0.17em}\overline{L}\right]$$

(5)

for any given positive real number$\overline{L}$, with$D_i\in D_S$and$L_i\in L_S$, where$D_S=\left\{D\in {\mathit{R}}_{+}\hspace{0.17em}:\hspace{0.17em}\alpha \left(D\right)=g\left(D\right)=E\right\}\subset {\mathit{R}}_{0+}$and$L_S=\left\{L\in {\mathit{R}}_{+}\hspace{0.17em}:\hspace{0.17em}\dot{\alpha }\left(L\right)+{\alpha }^2\left(L\right)=0\right\}\subset {\mathit{R}}_{0+}$are assumed to be nonempty and of zero Lebesgue measure.

Then, the following properties hold:

(i)

$g\left(L_i\right)=E\iff h\left(L_i\right)\dot{g}\left(L_i\right)>0$, equivalently,$D_S\cap L_S=\varnothing$.

(ii)

(a)$card{L}_S=card{D}_S+\vartheta$with$\vartheta =\left\{0,\hspace{0.17em}1\right\}$, and (a) if$card\left(D_S\right)=card\left(L_S\right)\le {\aleph }_0$(with${\aleph }_0$denoting the infinite cardinality of denumerable sets) then$L_i<D_i<L_{i+1}$;$\forall i\in {\mathit{Z}}_{0+}$for any pairs$D_i\hspace{0.17em},\hspace{0.17em}{D}_{i+1}\in D_S$and$L_i\hspace{0.17em},\hspace{0.17em}{L}_{i+1}\in L_S$fulfilling$\left(D_i,\hspace{0.17em}{D}_{i+1}\right)\cap D_S=\varnothing$and$\left(L_i,\hspace{0.17em}{L}_{i+1}\right)\cap L_S=\varnothing$,

(b) if$1\le card\left(L_S\right)=card\left(D_S\right)-1=\ell <\infty$then$L_i<D_i<L_{i+1}$for$i\in \overline{\ell -1}$.$\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$is subject to the constraint$c=K$,$\dot{I}\left(D_i\right)=\ddot{I}\left(L_i\right)=0$;$\forall L_i\in L_S\cap \left[0,\hspace{0.17em}L\right]$and$D_i>L_i$.

(iii)

$\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$is subject to$\dot{I}\left(D_i\right)=\ddot{I}\left(L_i\right)=0$;$\forall D_i\in D_S\cap \left[0,\hspace{0.17em}L\right]$,$\forall L_i\in L_S\cap \left[0,\hspace{0.17em}L\right]$and$D_1>L_1$for any$I_0>0$.

Proof. 

First, note that $\dot{I}\left(D\right)=\ddot{I}\left(L\right)=0$; $\forall D\in D_S$, $\forall L\in L_S$ since $\alpha \left(D\right)=0$ even if $I\left(D\right)\ne 0$. On the other hand, $L_S$ is the set of undulation points of $I:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ and it is clear that $D_S$ is contained in the set of relative maximum and minimum points of $I\left(t\right)$. The properties (i)–(iii) are now proved:

Proof of Property (i).

It is now proved that $D_S$ is the set of extreme points of $I\left(t\right)$ which is disjoint to its set of undulation points $L_S$. Assume, on the contrary, that there is some $D\notin D_S$ such that $\dot{I}\left(D\right)=0$. Then, $I\left(D\right)=0$ since $\alpha \left(D\right)\ne 0$, and then the disease-free equilibrium point is reached in finite time contradicting the fact that $\alpha \left(t\right)=-\frac{cln\left(g\left(t\right)/E\right)}{h\left(t\right)}$ is only zero at finite time for a discrete set of time instants satisfying $g\left(t\right)=E$ so that $\dot{I}\left(D\right)=0$ if and only if $D\in D_S$. Then, $I\left(D\right)=\dot{I}\left(D\right)=\ddot{I}\left(D\right)=0$ is a disease-free equilibrium point which is reached in finite time which contradicts the given hypothesis. So, it is easy to see that $L_S$ and $D_S$ are discrete sets of non-negative real time instants which can be strictly ordered. Note also from (1)–(2) that:

$$\alpha \left(D_i\right)=-\frac{cln\left(g\left(D_i\right)/E\right)}{h\left(D_i\right)}=0; \forall D_i\in D_S$$

(6)

$$\dot{\alpha }\left(L_i\right)=-{\alpha }^2\left(L_i\right); \forall L_i\in L_S$$

(7)

If $D_i\le D<+\infty$ then $g\left(D_i\right)=E$ since $c\ne 0$. Also, $\alpha \left(t\right)=-\frac{cln\left(g\left(t\right)/E\right)}{h\left(t\right)}$ and, if $L_i\in L_S$ and since $h\left(L_i\right)>0$, one has:

$$\begin{array}{c}{\alpha }^2\left(L_i\right)=-\dot{\alpha }\left(L_i\right)=\frac{d}{dt}{\left[\frac{cln\left(g\left(t\right)/E\right)}{h\left(t\right)}\right]}_{t=L_i}\\ =c\frac{h\left(L_i\right)\dot{g}\left(L_i\right)/g\left(L_i\right)-ln\left(g\left(L_i\right)/E\right)\dot{h}\left(L_i\right)}{h^2\left(L_i\right)}=c^2\frac{{ln}^2\left(g\left(L_i\right)/E\right)}{h^2\left(L_i\right)}\end{array}$$

(8)

$$c=\frac{h\left(L_i\right)\dot{g}\left(L_i\right)-ln\left(g\left(L_i\right)/E\right)\dot{h}\left(L_i\right)g\left(L_i\right)}{g\left(L_i\right){ln}^2\left(g\left(L_i\right)/E\right)}>0$$

(9)

Now, if there is some $L_i\in L_S\cap D_S$, equivalently, $L_S\cap D_S\ne \varnothing$, then $g\left(L_i\right)=E\iff h\left(L_i\right)\dot{g}\left(L_i\right)\ne 0$ from (9) since $c\ne 0$ and $ln\left(g\left(L_i\right)/E\right)=0$ and, furthermore, one gets from (8) that $\dot{\alpha }\left(L_i\right)\ne 0$ since $g\left(L_i\right)=E$. But one also has that $\dot{\alpha }\left(L_i\right)=\alpha \left(L_i\right)=0$, since $\dot{\alpha }\left(L_i\right)=-{\alpha }^2\left(L_i\right)$; $\forall L_i\in L_S$ from the first identity of (8). Then, $0\ne \dot{\alpha }\left(L_i\right)=0$ is a contradiction so that $L_i\notin L_S\cap D_S$. Equivalently, $D_S\cap L_S=\varnothing$. Property (i) has been proved. □

Proof of Property (ii).

Since $\dot{\alpha }\left(L_i\right)=-{\alpha }^2\left(L_i\right)$ then $\ddot{\alpha }\left(L_i\right)=-2\alpha \left(L_i\right)\dot{\alpha }\left(L_i\right)$ so that:

$$\ddot{I}\left(L_i\right)=\left(\dot{\alpha }\left(L_i\right)+{\alpha }^2\left(L_i\right)\right)I\left(L_i\right)=0$$

$$\stackrel{⃛}{I}\left(L_i\right)=\left(\ddot{\alpha }\left(L_i\right)+2\alpha \left(L_i\right)\dot{\alpha }\left(L_i\right)\right)I\left(L_i\right)+\left(\dot{\alpha }\left(L_i\right)+{\alpha }^2\left(L_i\right)\right)\dot{I}\left(L_i\right)$$

$$=\left(\ddot{\alpha }\left(L_i\right)-2{\alpha }^3\left(L_i\right)\right)I\left(L_i\right)$$

Since the zeros of $\alpha \left(t\right)$ and those of its first time- derivative do not coincide since $D_S\cap L_S=\varnothing$ (from Property (i)), it turns out that the two sets of respective zeros alternate if there are not two zeros of $\alpha \left(t\right)$ within any open time interval of two consecutive zeros of $\dot{\alpha }\left(t\right)$ or vice-versa. One proceeds by contradiction arguments by assuming two cases which are both rebutted.

Case 1: Assume that there are two consecutive zeros of $\dot{I}\left(t\right)$ between two consecutive zeros of $\ddot{I}\left(t\right)$, then satisfying the constraint $0\le L_i<D_i<D_{i+1}<L_{i+1}$ for some two consecutive time instants $D_i\hspace{0.17em},D_{i+1}$ in $D_S$ and two consecutive time instants $L_i\hspace{0.17em},L_{i+1}$ in $L_S$ so that $\alpha \left(D_i\right)=\alpha \left(D_{i+1}\right)=\ddot{I}\left(L_i\right)=\ddot{I}\left(L_{i+1}\right)=0$. Assume that $I\left(t\right)=0$ for some $t\in \left(D_i,D_{i+1}\right)$ then $\dot{I}\left(t\right)=\alpha \left(t\right)I\left(t\right)=0$ so that $t\in D_S$ and then $D_i\hspace{0.17em},D_{i+1}$ are not consecutive time instants in $D_S$ and this case has to be excluded from further reasoning. Now, assume that $I\left(t\right)\ne 0$ for all $t\in \left(D_i,D_{i+1}\right)$ and $\dot{\alpha }\left(t\right)\ne 0$, otherwise, if $\dot{\alpha }\left(t\right)=0$ then $t\in D_s$ and $D_i\hspace{0.17em},D_{i+1}$ are not consecutive time instants in $D_S$. Thus, $\alpha \left(t\right)=\alpha \left(D_i\right)+{\displaystyle {\int }_{D_i}^t\dot{\alpha }\left(\tau \right)}d\tau ={\displaystyle {\int }_{D_i}^t\dot{\alpha }\left(\tau \right)}d\tau$ for all $t\in \left(D_i,D_{i+1}\right)$. Since $\dot{\alpha }\left(t\right)\ne 0$ for all $t\in \left(D_i,D_{i+1}\right)$, it has no sign change in $\left(D_i,D_{i+1}\right)$ so that $\underset{t\to D_{i+1}^{-}}{lim}\alpha \left(t\right)\ne 0$ and since $\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is continuous then $\alpha \left(D_{i+1}\right)\ne 0$ which contradicts that $D_{i+1}\in D_S$. It has been proved that Case 1 is impossible $0\le L_i<D_i<D_{i+1}<L_{i+1}$ cannot happen.

Case 2: Assume now that there are two consecutive zeros of $\ddot{I}\left(t\right)$ between two consecutive zeros of $\dot{I}\left(t\right)$, that is $0\le D_i<L_i<L_{i+1}<D_{i+1}$ for some consecutive time instants $D_i\hspace{0.17em},D_{i+1}$ in $D_S$ and some two consecutive time instants $L_i\hspace{0.17em},L_{i+1}$ in $L_S$. Then, $\alpha \left(t\right)\ne 0$ for all $t\in \left(L_i,\hspace{0.17em}{L}_{i+1}\right)$ since, otherwise, there exists some $t\in \left(L_i,\hspace{0.17em}{L}_{i+1}\right)$ such that $t\in D_S$, and then the previously claimed constraint $0\le D_i<L_i<L_{i+1}<D_{i+1}$ does not hold, and also $\dot{\alpha }\left(t\right)\ne -{\alpha }^2\left(t\right)<0$ for all $t\in \left(L_i,\hspace{0.17em}{L}_{i+1}\right)$ since, otherwise, there exists some $t\in \left(L_i,\hspace{0.17em}{L}_{i+1}\right)$ such that $t\in L_S$ and then $L_i$ and $L_{i+1}$ are not two consecutive time instants in $L_S$ as claimed. Also, note that.

$\dot{\alpha }\left(L_i\right)+{\alpha }^2\left(L_i\right)=\dot{\alpha }\left(L_i\right)+{\alpha }^2\left(L_i\right)=0$ with $\alpha \left(L_i\right)\ne 0$ and $\alpha \left(L_{i+1}\right)\ne 0$ since $L_i\hspace{0.17em},L_{i+1}\notin D_S$. But then, by continuity arguments on $\dot{\alpha }\left(t\right)+{\alpha }^2\left(t\right)$, there is a change of sign point $t\in \left(L_i,\hspace{0.17em}{L}_{i+1}\right)$ which zeroes this function which contradicts $\dot{\alpha }\left(t\right)\ne -{\alpha }^2\left(t\right)<0$ for all $t\in \left(L_i,\hspace{0.17em}{L}_{i+1}\right)$. Then, Case 2 is impossible so that $0\le D_i<L_i<L_{i+1}<D_{i+1}$ cannot happen and Property (ii) has been proved. □

Proof of Property (iii).

Assume that, contrarily to the statement, $D_1\le L_1$. If $L_1=D_1$ then $\dot{I}\left(L_1\right)=\ddot{I}\left(L_1\right)$ and the equilibrium point is reached in finite time what is impossible, since $I_0>0$, for a non-trivial solution of a continuous-time first-order differential equation with continuous-time parameterization. Then, $L_1=D_1$ is impossible. Now, assume that $L_1>D_1$ and $0=\dot{I}\left(L_1\right)=\dot{I}\left(D_1\right)+{\displaystyle {\int }_{L_1}^{D_1}\ddot{I}\left(\tau \right)}d\tau ={\displaystyle {\int }_{L_1}^{D_1}\ddot{I}\left(\tau \right)}d\tau$ with $\ddot{I}\left(L_1\right)=0$ and then it exists some $L_2\in \left(L_1,D_1\right)$ such that $\ddot{I}\left(L_2\right)=0$ and $L_2\in L_S$. As a result, there is $D_1>L_2>L_1$ and then there are two consecutive undulation time instants what contradicts Property (ii). As a result, $D_1>L_1$ as claimed. □

Remark 1.

In Theorem 2, note that the sets$D_S$and$L_S$have the following properties:

They are nonempty so that there is at least one$D\in D_S$such that$\alpha \left(D\right)=0$implying that$\dot{I}\left(D\right)=0$and at least one$L\in L_S$such that$\dot{\alpha }\left(L\right)=-{\alpha }^2\left(L\right)$implying that$\ddot{I}\left(L\right)=0$. Otherwise, the infection could converge asymptotically to zero as time goes to infinity but it would not have finite zeros,

They are sets of zero Lebesgue measure so that they are denumerable discrete sets of strictly ordered isolated real points, for any real numbers,

They fulfill that$card{L}_S=card{D}_S+\vartheta$with$\vartheta =\left\{0,\hspace{0.17em}1\right\}$so that they are of either identical finite or infinite cardinal or the cardinal of$L_S$is finite and exceeds that of$D_S$by one,

If$\vartheta =0$then$card\left(D_S\right)=card\left(L_S\right)\le {\aleph }_0$, that is, if both sets have infinity cardinal or identical finite one then any ordered points of$L_S$and$D_S$alternate.

On the other hand, note that:

Equation (4) establishes that $D_S$ is the set of zeros of $\alpha \left(t\right)$. At those zeros, the first-time derivative of the infection function is zeroed from (1) without such a function being necessarily zero while on the other hand, Equation (5) is a nonzero real constant for any finite undulation time instant $L_i\le \overline{L}$ of $I:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ zeroing the second derivative of the infection function according to (2) which holds if $c=K$ from (5). The fact that (5) is constant follows easily under periodicity conditions of the same or integer multiple/submultiple periods of $g\left(t\right)$ and $h\left(t\right)$.

Since $\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ has no finite zero coincident with a zero of its first time-derivative, by hypothesis, then $g\left(L_i\right)=E\iff h\left(L_i\right)\dot{g}\left(L_i\right)\ne 0$ since $c\ne 0$ from inspection of (8)–(9). This is equivalent to $D_S\cap L_S=\varnothing$, that is, the finite zeros which make zero $\dot{I}\left(t\right)$ and which do not make zero $I\left(t\right)$ do not make zero either $\ddot{I}\left(t\right)$. However, $\ddot{I}\left(t\right)=0$ if $I\left(t\right)=\dot{I}\left(t\right)=0$ from (2), provided that $\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is twice everywhere continuously differentiable in $\left[0,\hspace{0.17em}+\infty \right)$ but this can only happen as time tends to infinity for certain structures of $g\left(t\right)$ and $h\left(t\right)$. Note that the constraint (5) also implies that the auxiliary functions $g,\hspace{0.17em}h:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ used to define the function $\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ in (1) fulfill the constraint $h\left(L_i\right)\dot{g}\left(L_i\right)\ne ln\left(g\left(L_i\right)/E\right)\dot{h}\left(L_i\right)g\left(L_i\right)$; $\forall L_i\in L_S$.

By examining Definitions 1 and 3 and Lemma 1, it turns out that the set $L_S$ of undulation points of $I\left(t\right)$ includes but, maybe non-properly, the set of its inflection points. However, it suffices to give some further weak conditions on $\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$, that is, on $g,\hspace{0.17em}h:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ to guarantee that every undulation point of $I\left(t\right)$ is also an inflection point. Some such conditions are discussed in the next corollary.

Corollary 2.

The following properties hold:

(i)

Assume that:

$$\underset{\epsilon \to 0^{+}}{lim\hspace{0.17em}sup}\left[\theta \left(L_i+\epsilon \right)\theta \left(L_i-\epsilon \right)\right]<0; \forall L_i\in L_S$$

where:

$$\theta \left(t\right)=h\left(t\right)\dot{g}\left(t\right)-g\left(t\right)\dot{h}\left(t\right)ln\left(g\left(t\right)/E\right); \forall t\in {\mathit{R}}_{0+}$$

Then, the set$L_S$of undulation points of$I\left(t\right)$is the set of its inflection points.

(ii)

Assume that$f,g:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$are twice continuously differentiable at each undulation point$L_i\in L_S$. Then, the sets of undulation points and that of the inflection points of$I\left(t\right)$coincide if

$$\frac{h^3\left(L_i\right)\left(g\left(L_i\right)\hspace{0.17em}\ddot{g}\left(L_i\right)-{\dot{g}}^2\left(L_i\right)\right)}{h^4\left(L_i\right)g^2\left(L_i\right)}\ne \hspace{0.17em}\frac{2}{h^3\left(L_i\right)}\left({ln}^3\frac{g\left(L_i\right)}{E}+\dot{h}\left(L_i\right)ln\frac{g\left(L_i\right)}{E}\right); \forall L_i\in L_S$$

Proof. 

Note that $\ddot{I}\left(t\right)=\left({\alpha }^2\left(t\right)+\dot{\alpha }\left(t\right)\right)I\left(t\right)$, $\forall t\in {\mathit{R}}_{0+}$ so that $\ddot{I}\left(L_i\pm \epsilon \right)=\left({\alpha }^2\left(L_i\pm \epsilon \right)+\dot{\alpha }\left(L_i\pm \epsilon \right)\right)I\left(L_i\pm \epsilon \right)$. Since $L_i>0$, $g\left(t\right)h\left(t\right)>0$ if $t>0$ and $\underset{\epsilon \to 0}{lim}I\left(L_i\pm \epsilon \right)=I\left(L_i\right)$, since $I\left(t\right)$ is continuous, one gets that $\underset{\epsilon \to 0^{+}}{lim\hspace{0.17em}sup}\left[\ddot{I}\left(L_i+\epsilon \right)\ddot{I}\left(L_i-\epsilon \right)\right]<0$ if and only if $\underset{\epsilon \to 0^{+}}{lim\hspace{0.17em}sup}\left[\theta \left(L_i+\epsilon \right)\theta \left(L_i-\epsilon \right)\right]<0$. Property (i) has been proved.

On the other hand, if $f,g:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ are twice continuously differentiable at each undulation point $L_i\in L_S$ of $I\left(t\right)$, then $\ddot{f},\ddot{g}$ exist in $L_S$. Then, defining $\widehat{\alpha }\left(t\right)=-c^{-1}\alpha \left(t\right)=\frac{ln\left(g\left(t\right)/E\right)}{h\left(t\right)}$; $\forall t\in {\mathit{R}}_{0+}$ yields:

$$\dot{\widehat{\alpha }}\left(t\right)=\frac{h\left(t\right)\dot{g}\left(t\right)-g\left(t\right)\dot{h}\left(t\right)ln\left(g\left(t\right)/E\right)}{h^2\left(t\right)g\left(t\right)}; \forall t\in L_S$$

$$\ddot{\widehat{\alpha }}\left(t\right)=\frac{h^3\left(t\right)\left(g\left(t\right)\ddot{g}\left(t\right)-{\dot{g}}^2\left(t\right)\right)/g^2\left(t\right)+\dot{h}\left(t\right)\dot{g}\left(t\right)/g\left(t\right)-2h\left(t\right)\dot{h}\left(t\right)ln\left(g\left(t\right)/E\right)}{h^4\left(t\right)}; \forall t\in L_S$$

$$\ddot{I}\left(t\right)=\left({\alpha }^2\left(t\right)+\dot{\alpha }\left(t\right)\right)I\left(t\right)\Rightarrow \stackrel{⃛}{I}\left(t\right)=0 \mathrm{with} {\alpha }^2\left(t\right)=-\dot{\alpha }\left(t\right) \mathrm{and} I\left(t\right)>0; \forall t\in L_S$$

$$\stackrel{⃛}{I}\left(t\right)=\left({\alpha }^2\left(t\right)+\dot{\alpha }\left(t\right)\right)\dot{I}\left(t\right)+\left(2\alpha \left(t\right)\dot{\alpha }\left(t\right)+\ddot{\alpha }\left(t\right)\right)I\left(t\right)\Rightarrow \stackrel{⃛}{I}\left(t\right)=\left(2\alpha \left(t\right)\dot{\alpha }\left(t\right)+\ddot{\alpha }\left(t\right)\right)I\left(t\right)=\left(\ddot{\alpha }\left(t\right)-2{\alpha }^3\left(t\right)\right)I\left(t\right); \forall t\in L_S$$

Since $I\left(t\right)>0$; $\forall t\in {\mathit{R}}_{0+}$ then $\stackrel{⃛}{I}\left(t\right)\ne 0$; $\forall t\in L_S$ if and only if $\ddot{\alpha }\left(t\right)\ne 2{\alpha }^3\left(t\right)$; $\forall t\in L_S$, equivalently, if and only if $\ddot{\widehat{\alpha }}\left(t\right)\ne 2{\widehat{\alpha }}^3\left(t\right)$; $\forall t\in L_S$ which is fully equivalent to the condition of Property (ii). The proof is complete.  □

Remark 2.

Note that Theorem 2 applies, in particular, to the case when there are equilibrium points with the initial conditions being distinct from such points. It can be also extended by including the above case by redefining finite discrete sets of the zeros of$\dot{I}\left(t\right)$and$\ddot{I}\left(t\right)$$D_S\to D_S\cap \left[0,\hspace{0.17em}L\right]$¸$L_S\to L_S\cap \left[0,\hspace{0.17em}L\right]$for any given$L\in \left[0,\hspace{0.17em}\infty \right)$in the sense that the eventual zeros at finite time of$\dot{I}\left(t\right)$and$\ddot{I}\left(t\right)$alternate although an equilibrium points has not still been reached provided that it exists.

Inspired in Theorem 2, some conditions are discussed in the next result which imply that the first undulation point of the infection evolution function (i.e., the first zero of its second-time derivative) precedes the first zero of its first time-derivative. It is not required that the infection has necessarily a disease-free equilibrium point or that it might be oscillatory leading to successive zeros of its time- derivative along time.

Theorem 3.

Assume that the function$\alpha \left(t\right)=-\frac{cln\left(g\left(t\right)/E\right)}{h\left(t\right)}$, where$c\hspace{0.17em},\hspace{0.17em}E\in {\mathit{R}}_{+}$and$g,h:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$are everywhere continuous and time-differentiable and satisfy the constraints:

(1) 

$g\left(t\right)<E$;$\forall t\in \left[0,\hspace{0.17em}D\right)$,$g\left(D\right)=E$

(2) 

$\dot{g}\left(0\right)<\hspace{0.17em}\left(\frac{c{ln}^2\left(g\left(0\right)/E\right)}{h\left(0\right)}-\frac{\left|ln\left(g\left(0\right)/E\right)\right|}{h\left(0\right)\dot{h}\left(0\right)}\right)g\left(0\right)$

(3) 

$g\left(t\right)>0$and$h\left(t\right)>0$if$t>0$

(4) 

$\frac{ln\left(g\left(0\right)/E\right)}{h\left(0\right)}\ne 0$

Assume also that$I_0\hspace{0.17em}>0$. Then,$min\left(I\left(t\right),\hspace{0.17em}\dot{I}\left(t\right)\right)>0$;$\forall t\in \left[0,\hspace{0.17em}D\right)$;$\dot{I}\left(D\right)=0$and there is some$L\in \left(0,\hspace{0.17em}D\right)$such that$\ddot{I}\left(t\right)\ne 0$;$\forall t\in \left[0,\hspace{0.17em}L\right)$and$\ddot{I}\left(L\right)=0$.

Proof. 

Note from the definition of $\alpha \left(t\right)$, (1), (2) and the given constraints 1 and 2 that

$$\dot{\alpha }\left(0\right)=-\frac{c}{h\left(0\right)}\left(\frac{1}{h\left(0\right)\dot{h}\left(0\right)}\left|ln\frac{g\left(0\right)}{E}\right|+\frac{\dot{g}\left(0\right)}{g\left(0\right)}\right)$$

$\alpha \left(0\right)>0$, since $0\le g\left(0\right)<E$, $\alpha \left(D\right)=0$, since $g\left(D\right)=E$), ${\alpha }^2\left(0\right)+\dot{\alpha }\left(0\right)>0$, from the condition 2 since $\alpha \left(0\right)>0$ and since $\alpha :{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is continuous and time-differentiable since $g,h:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ are everywhere continuous and time-differentiable. Note also that, from the given assumptions and constraints, $min\left(I_0,\hspace{0.17em}{\dot{I}}_0,{\ddot{I}}_0\hspace{0.17em}\right)>0$ since $I_0>0$ by hypothesis, ${\dot{I}}_0=\alpha \left(0\right)\hspace{0.17em}{I}_0>0$ and ${\ddot{I}}_0=\left({\alpha }^2\left(0\right)+\dot{\alpha }\left(0\right)\right)I_0>0$. Furthermore, $\dot{I}\left(D\right)=\alpha \left(D\right)I\left(D\right)=0$. From the constraint 3 and the continuity of $g,h:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$, one has that $\alpha ,\dot{\alpha },\ddot{\alpha }:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ are continuous and bounded on $\left(0,+\infty \right)$, $\dot{I}\left(t\right)>0$; $\forall t\in \left[0,\hspace{0.17em}D\right)$ and $\ddot{I}\left(t\right)>0$; $\forall t\in \left[0,\hspace{0.17em}{L}_0\right)$ and some $L_0\in R_{+}$. Furthermore since $c>0$ and $\frac{ln\left(g\left(0\right)/E\right)}{h\left(0\right)}\ne 0$, from the constraint 4, $g\left(t\right)<E$; $\forall t\in \left[0,\hspace{0.17em}D\right)$, from the constraint 1, and $g\left(t\right)>0$ and $h\left(t\right)>0$ if $t>0$, from the constraint 3. Then $\alpha \left(t\right)>0$; $\forall t\in \left[0,\hspace{0.17em}D\right)$. Since $g,h:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ are continuous and positive on any bounded interval $\left[0,\hspace{0.17em}T\right)$ then $\alpha \left(t\right)$ is positive and finite on $\left[0,\hspace{0.17em}D\right)$. It is now proved that $t=D$ is the first zero of $\dot{I}\left(t\right)$. Assume that this is not the case so that there is some $D_1<D$ such that $\dot{I}\left(D_1\right)=0$, with $\alpha \left(D_1\right)\ne 0$, and $\dot{I}\left(t\right)>0$; $\forall t\in \left[0,\hspace{0.17em}{D}_1\right)$. Then $I\left(D_1\right)=\dot{I}\left(D_1\right)=\dot{I}\left(D_1\right)=0$ from (2) and the infection extinguishes in a finite time $D_1<D$. This leads to a contradiction since $I\left(D_1\right)=I_0+{\displaystyle {\int }_0^{D_1}\dot{I}\left(\tau \right)}d\tau >0$ since $I_0>0$ and $\dot{I}\left(t\right)>0$; $\forall t\in \left[0,\hspace{0.17em}{D}_1\right)$. Therefore, if $D_1<D$ such that $\dot{I}\left(D_1\right)=0$ then $I\left(D_1\right)>0$. But then $\alpha \left(D_1\right)=\dot{I}\left(D_1\right)/I\left(D_1\right)=0$ from (1) which contradicts that $\alpha \left(t\right)\ne 0$; $\forall t\in \left[0,\hspace{0.17em}D\right)$. As a result, $t=D$ is the first zero of $\dot{I}\left(t\right)$ and there is no $D_1<D$ such that $I\left(D_1\right)=0$. Since $I,\dot{I}:{\mathit{R}}_{0+}\cap \left[0,\hspace{0.17em}D\right]\to {\mathit{R}}_{0+}$ are continuous with $\dot{I}\left(t\right)>0$; $\forall t\in \left[0,\hspace{0.17em}D\right)$ and $\dot{I}\left(D\right)=0$ and $\ddot{I}\left(t\right)$; $\forall t\in \left[0,\hspace{0.17em}{L}_0\right)$ and some $L_0\in {\mathit{R}}_{+}$ then there is some $L\in \left(0,\hspace{0.17em}D\right)$ such that $\ddot{I}\left(L\right)=0$. Assume that this is not the case. Then, $0=\dot{I}\left(D\right)={\dot{I}}_0+{\displaystyle {\int }_0^D\ddot{I}\left(\tau \right)}d\tau >0$. Hence, a contradiction arises. Thus, there is some $L\in \left(0,\hspace{0.17em}D\right)$ such that $\ddot{I}\left(L\right)=0$. □

Remark 3.

Note that, under all the conditions of Theorem 3,$\alpha \left(t\right)>0$;$\forall t\in \left[0,D\right)$and$\alpha \left(D\right)=0$. Furthermore, the first zero of$\dot{I}\left(t\right)=0$occurs at$t=D$, there is no$t<D$such that$I\left(t\right)=\dot{I}\left(t\right)=0$and there is some$L<D$such that$\ddot{I}\left(L\right)=0$.

The following example describes the basic model proposed in [11] under a first-order differential equation for the infection evolution without any entropy considerations at this stage:

Example 1.

The function$\alpha \left(t\right)=-c\hspace{0.17em}ln\left(t/D\right)/t$, for some$D>0$, proposed in [11] satisfies all the conditions of Theorem 3 with$h\left(t\right)=g\left(t\right)=t$and$E=D$. It satisfies, in addition, that$\alpha \left(0\right)=+\infty$. This function satisfies also the given further conditions of Theorem 2$g\left(0\right)=h\left(0\right)=0$with$\underset{t\to 0}{lim}\frac{ln\left(g\left(t\right)/E\right)}{h\left(t\right)}\le -\epsilon$.

Note that the condition$\alpha \left(0\right)>0$of Theorem 3 avoids that${\dot{I}}_0=0$if$I_0\ne 0$so that$t=0$is a zero of$\dot{I}\left(t\right)$.

It can be argued that the proposed basic model (1) is a very simple time-varying differential equation of first-order which describes the infective population time-evolution. Note that the use of appropriate particular structures in the definition of the time-varying coefficient $\alpha \left(t\right)$ can take care of the eventual incorporation of the necessary supplementary environment information to make such an equation well-posed to practically describe a concrete disease evolution through time. The incorporation which can be incorporated is the eventual couplings of the infectious subpopulation with another ones (such as the susceptible, recovered or vaccinated subpopulations and their associated dynamics) or the information about the feedback information controls in more elaborated models. The next section develops some work in this direction.

3. Further Examples of Linking the Basic Model to Some Existing Epidemic Models Incorporating Other Subpopulations

The infection description via (1) assumes implicitly that it has a first-order dynamics. It has been argued that $\alpha \left(t\right)$ in (1) contains the information about the controls and other coupled subpopulations influencing the disease evolution through time. It can be of interest to discuss its application to infection descriptions described by differential equations of orders higher than one which is a very common situation in disease transmission mathematical models.

It is now seen how a well-known epidemic model can be also discussed under the point of view of Theorem 3. In the subsequent example, the above characterization, based on the first zero of infection evolution time-derivative and on the undulation point of the infection evolution, is used for a model with three subpopulations via an appropriate choice of $g\left(t\right)$ and $h\left(t\right)$ in the definition of $\alpha \left(t\right)$.

Example 2.

Consider the following SIR model without demography [30]:

$$\dot{S}\left(t\right)=-\beta S\left(t\right)I\left(t\right); \dot{I}\left(t\right)=\left(\beta S\left(t\right)-\gamma \right)I\left(t\right); \dot{R}\left(t\right)=\gamma I\left(t\right); \forall t\in {\mathit{R}}_{0+}$$

(10)

where$S\left(t\right)$,$I\left(t\right)$and$R\left(t\right)$are, respectively, the susceptible, infectious and recovered (or immune) subpopulations, under nonzero initial conditions being subject to$min\left(S\left(0\right),\hspace{0.17em}I\left(0\right),\hspace{0.17em}R\left(0\right)\right)\ge 0$, where$\beta$is the coefficient transmission rate and$\gamma$is the removal or recovery rate (its inverse${\gamma }^{-1}$being the average infectious period). The mathematical study of this model and their variants is not easy as seen in [30,40]. First, note that the total population$N\left(t\right)=S\left(t\right)+R\left(t\right)+I\left(t\right)=S_0+R_0+I_0$;$\forall t\in {\mathit{R}}_{0+}$is constant for all time. The basic reproductive ratio (or reproduction number) is$R_{{}_{*}}=\beta /\gamma$and, if$S_0\le R_{{}_{*}}^{-1}$, then${\dot{I}}_0\le 0$while if$S_0>R_{{}_{*}}^{-1}$, it becomes endemic for all time since${\dot{I}}_0>0$. The solution of (10) becomes in closed form:

$$S\left(t\right)=e^{-\beta {\displaystyle {\int }_0^{t}I\left(\tau \right)d\tau }}{S}_0; I\left(t\right)=e^{{\displaystyle {\int }_0^t\left(\beta S\left(\tau \right)-\gamma \right)d\tau }}{I}_0; R\left(t\right)=S_0+R_0+I_0-S\left(t\right)-I\left(t\right); \forall t\in {\mathit{R}}_{0+}$$

(11)

Note that by combining the above equations that:

$$S\left(t\right)=e^{-\beta I_0{\displaystyle {\int }_0^t{e}^{{\displaystyle {\int }_0^{\tau }\left(\beta S\left(\sigma \right)-\gamma \right)d\sigma }}d\tau }}{S}_0; I\left(t\right)=e^{{\displaystyle {\int }_0^t\left(\beta e^{-\beta {\displaystyle {\int }_0^{\tau }I\left(\sigma \right)d\sigma }}{S}_0-\gamma \right)d\tau }}{I}_0$$

(12)

Note from (11) that $S:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is non-increasing so that there exists a susceptible equilibrium subpopulation $S_e=\underset{t\to \infty }{lim}S\left(t\right)\le S_0$ for any given non-negative initial conditions. Note also from (10) that $\dot{N}\left(t\right)=0$ and then $N\left(t\right)=N_0$; $\forall t\in {\mathit{R}}_{0+}$ Note that If $I_0=0$ then $I\left(t\right)=0$, $S\left(t\right)=S_0$ and $R\left(t\right)=R_0=N_0-S_0$; $\forall t\in {\mathit{R}}_{0+}$. We examine three cases for $I_0>0$:

Case (a) if $S_0<R_{{}_{*}}^{-1}$ then $S\left(t\right)\le S_0$ and $\beta S\left(t\right)-\gamma <0$; $\forall t\in {\mathit{R}}_{0+}$, then $I\left(t\right)\to 0,S\left(t\right)\to S_e$ and $R\left(t\right)\to R_e=N_0-S_e$ as $t\to \infty$. Since $S:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is non-increasing, $S_e\le S_0<R_{{}_{*}}^{-1}$. This implies that $\underset{t\to \infty }{lim}{\displaystyle {\int }_0^t\left(\beta S\left(\tau \right)-\gamma \right)d\tau }=-\infty$ and $\dot{I}\left(t\right)=-\lambda \left(t\right)I\left(t\right)\le -{\lambda }_{a}I\left(t\right)$, $I\left(t\right)\to 0$ at exponential rate as $t\to \infty$ for some ${\lambda }_a>0$ from (10) and (11) since $I_0-I\left(t\right)\ge {\lambda }_a{\displaystyle {\int }_0^{t}I\left(\tau \right)}d\tau$ so that ${\displaystyle {\int }_0^{\infty }I\left(\tau \right)}d\tau \le I_0/{\lambda }_a<+\infty$. Then, $I:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is integrable on $\left[0,\infty \right)$. Thus, $C=\beta {\displaystyle {\int }_0^{\infty }I\left(t\right)dt}<+\infty$ so that $S_e=e^{-\beta {\displaystyle {\int }_0^{\infty }I\left(t\right)dt}}{S}_0=e^{-C}{S}_0>0$ (then there is a nonzero susceptible equilibrium level) and $R_e=N_0-S_e<N_0$.

Case (b) if $S_0=R_{{}_{*}}^{-1}$ then $S\left(t\right)\to S_e\le S_0=\gamma /\beta$ as $t\to \infty$ since $S:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is non-increasing and then it converges to $S_e$ satisfying $0\le S_e\le S_0$. By inspection of the second equation of (11), it also follows that $I\left(t\right)\to I_e$ and $R\left(t\right)\to R_e$ as $t\to \infty$ satisfying $I_e\ge 0$ and $R_e\ge 0$. Assume that $I_e>0$ then $S_e=0$ from the first equation of (11). But if $S_e=0$ then $I_e=0$ since then $I:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is strictly decreasing on $\left[t_a,\hspace{0.17em}\infty \right)$ for some finite $t_a>0$ from the second equation of (11). Hence, a contradiction to $I_e>0$ follows implying that $I_e=0$ if $S_e=0$. Now, assume that $\gamma /\beta >S_e>0$. Then, from the second equation of (11), $I\left(t\right)\to I_e=0$ as $t\to \infty$. But then $S_e>0$, from the first equation of (12), since $\gamma /\beta >S_e$ if $I_0>0$ and then $R_e=N_0-S_e$. From the second equation of (12) and, under a similar reasoning as that of Case a, $I:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is integrable on $\left[0,\infty \right)$ and $S_e>0$. In summary, if $S_0=R_{{}_{*}}^{-1}=\gamma /\beta$ and $I_0>0$ then $I\left(t\right)\to 0,S\left(t\right)\to 0$ and $R\left(t\right)\to N_0=S_0+R_0+I_0$ as $t\to \infty$ in the same way as in Case a if $S_0\le R_{{}_{*}}^{-1}$.

Case (c) if $S_0>R_{{}_{*}}^{-1}$ then ${\dot{I}}_0>0$ from (10) and $S:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is increasing on some interval $\left[0,t_0\right]$. The fact that $I:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is strictly increasing on some initial time interval is of interest from the point of view of hospital management of availability of beds and other sanitary specific means in the event that the disease might have a relevant number of seriously infected individuals. Since $S:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is non-increasing then either $I\left(t\right)\to I_e=S_0+I_0+R_0=N_0$, $S\left(t\right)\to S_e=0$ and $R\left(t\right)\to R_e=0$ as $t\to \infty$ or $S\left(t\right)\to S_e\in {\left(0,\hspace{0.17em}{R}_{{}_{*}}^{-1}\hspace{0.17em}\right]}^{}$ as $t\to \infty$ from (11) since $S:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is non-increasing. The firs possibility $I\left(t\right)\to I_e=N_0$ is unfeasible since from the first equation of (11) $I\left(t\right)\to \infty$ as $t\to \infty$. Then, $S\left(t\right)\to S_e\in {\left(0,\hspace{0.17em}{R}_{{}_{*}}^{-1}\hspace{0.17em}\right]}^{}$ as $t\to \infty$. Now, first, assume that $S_e\in \left(\gamma /\beta ,\hspace{0.17em}{R}_{{}_{*}}^{-1}\right]$. Then, from the first equation of (12), $S\left(t\right)\to 0$ as $t\to \infty$. Then, $S_e=0$ which contradicts that $S_e>\gamma /\beta \hspace{0.17em}$, As a result, $0\le S_e\le \gamma /\beta \hspace{0.17em}$. Now, assume that $S_e=0$. Then, from (11), $I\left(t\right)\to I_e=0$ and $I:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ being square-integrable, and following a similar argument as that of Cases a–b, one again concludes that $S_e>0$ so that $S_e\in \left(0,\hspace{0.17em}\gamma /\beta \hspace{0.17em}\right]$ and $R_e=N_0-S_e$, as a result. But, since $S_e\le \gamma /\beta$ then $I_e=0$ from (11) since $I:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is strictly decreasing after some finite time instant $t_0$ and integrable on $\left[0,\infty \right)$ and a following again the reasoning of Cases a–b, one concludes that $S_e>0$. As a result, if $S_0>R_{{}_{*}}^{-1}$ and $I_0>0$, then $I_e=0$, $S_e>0$ and $R_e=N_0-S_e$. Thus, the relevant conclusions on the disease- free equilibrium point which is a disease- free one are similar for the three above cases.

On the other hand, since $S:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ it exists a finite $t=D>0$ such that $S\left(D\right)=R_{{}_{*}}^{-1}=\gamma /\beta$ and $\dot{I}\left(D\right)=\alpha \left(D\right)I\left(D\right)=\left(\beta S\left(D\right)-\gamma \right)I\left(D\right)=0$, $I\left(D\right)=e^{{\displaystyle {\int }_0^D\left(\beta S\left(\tau \right)-\gamma \right)d\tau }}{I}_0\ne 0$, if $I_0\ne 0$ and, furthermore,

$$\begin{array}{c}\ddot{I}\left(D\right)=\left(\beta \dot{S}\left(D\right)-\gamma \right)I\left(D\right)+\left(\beta S\left(D\right)-\gamma \right)\dot{I}\left(D\right)\\ =\left(\beta \dot{S}\left(D\right)-\gamma \right)I\left(D\right)=-{\beta }^{2}S\left(D\right)I^2\left(D\right)-\gamma I\left(D\right)\\ =-\gamma \left(\beta \hspace{0.17em}{e}^{{\displaystyle {\int }_0^D\left(\beta S\left(\tau \right)-\gamma \right)d\tau }}{I}_0+1\right)e^{{\displaystyle {\int }_0^D\left(\beta S\left(\tau \right)-\gamma \right)d\tau }}{I}_0<0\end{array}$$

(13)

and also:

$${\ddot{I}}_0=-{\beta }^2{S}_0{I}_0^2+\left(\beta S_0-\gamma \right){\dot{I}}_0=I_0{\left[{\left(\beta S_0-\gamma \right)}^2-{\beta }^2{S}_0{I}_0\right]}^{}$$

(14)

and ${\ddot{I}}_0>0$ under the reasonable assumption that $I_0$ is sufficiently small (the initial numbers of infectious is usually very small in practice) satisfying $I_0<\frac{{\left(\beta S_0-\gamma \right)}^2}{{\beta }^2{S}_0}$. As a result, there is some time instant $L\in \left(0,\hspace{0.17em}D\right)$ such that $\ddot{I}\left(L\right)=0$ so that it is an undulation point of $I:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$. As a result, we find that if the basic reproduction number exceeds unity then the infection curve corresponding to the endemic solution has a minimum at a larger time instant that the one defining its undulation point. That situation corresponds to the situation of small initial infection force with reproduction number greater than one. On the other hand, if ${\ddot{I}}_0\le 0$, then ${\dot{I}}_0>0$ does not hold.

Comparing the infectious subpopulation evolution to (1) and the structure of the function in Theorem 3 yields:

$$\alpha \left(t\right)=\beta \hspace{0.17em}S\left(t\right)-\gamma =-\frac{cln\left(g\left(t\right)/E\right)}{h\left(t\right)}$$

(15)

$$\dot{\alpha }\left(t\right)=\beta \hspace{0.17em}\dot{S}\left(t\right)=-c\frac{d}{dt}\left(\frac{ln\left(g\left(t\right)/E\right)}{h\left(t\right)}\right)$$

(16)

$$=-{\beta }^{2}S\left(t\right)I\left(t\right)=-\frac{c}{h\left(t\right)}\left(\frac{1}{h\left(t\right)\dot{h}\left(t\right)}\left|ln\frac{g\left(t\right)}{E}\right|+\frac{\dot{g}\left(t\right)}{g\left(t\right)}\right);$$

(17)

$\forall t\in {\mathit{R}}_{0+}$. If one defines $g\left(t\right)=t$; $\forall t\in {\mathit{R}}_{0+}$ and $h\left(t\right)=\frac{cln\left(t/E\right)}{\gamma -\beta \hspace{0.17em}S\left(t\right)}$; $\forall t\in {\mathit{R}}_{0+}$, then $h\left(t\right)=\frac{c\left|ln\left(t/E\right)\right|}{\beta \hspace{0.17em}S\left(t\right)-\gamma }$; $\forall t\in {\mathit{R}}_{0+}$. It is easy to verify that these functions satisfy the conditions of Theorem 3.

In the case when the reproduction number is less than unity and it is an upper-bound of the normalized susceptible population, each primary infection generates, in average, less than one secondary one so that the infection extinguishes asymptotically. According to this particular model, also the susceptible subpopulation extinguishes asymptotically. See Case a referred to (11). Thus, the disease-free equilibrium point is ${\left(S_{df}^{*}\hspace{0.17em},\hspace{0.17em}{I}_{df}^{*}\hspace{0.17em},\hspace{0.17em}{R}_{df}^{*}\right)}^T={\left(0\hspace{0.17em},\hspace{0.17em}0\hspace{0.17em},\hspace{0.17em}N\right)}^T$. In this case, $I\left(t\right),\dot{I}\left(t\right),\ddot{I}\left(t\right)\to 0$ as $t\to \infty$ but there are no finite time instants of minimum and undulation of the infectious curve to the light of Theorem 3.

However, we can have a practical visualization of the disease removal by defining a design quadruple $\left(k_1,k_2,k_3,\epsilon \right)\in {\mathit{R}}_{+}^4$ and the following cut associate time instants:

$$t_{Ii}\left(k_i,\hspace{0.17em}\epsilon \right)=min\left(\tau \in {\mathit{R}}_{0+}:\hspace{0.17em}\left|\frac{d{I}^{\left(i-1\right)}}{dt}\right|\le k_i\epsilon :\hspace{0.17em}t\in \left[\tau \hspace{0.17em},\hspace{0.17em}+\infty \right)\right); i=1,2,3$$

(18)

Note that $t_{I2}\left(k_2,\hspace{0.17em}\epsilon \right)$ and $t_{I3}\left(k_3,\hspace{0.17em}\epsilon \right)$ generalize the roles of the time instants $D$ and $L$, that is, the finite minimum infection and undulation time instants, respectively, within prescribed margins when those time instants do not exist.

Example 3.

Consider Case a of Example 2 so that$S\left(t\right)\le S_0<\gamma /\beta$leading to$I\left(t\right)\to 0,S\left(t\right)\to S_e>0$and$R\left(t\right)\to R_e=N_0-S_e$as$t\to \infty$and$I\left(t\right)>0$,$\dot{I}\left(t\right)<0$and$\ddot{I}\left(t\right)<0$are strictly decreasing on$\left[0,\hspace{0.17em}+\infty \right)$. Take prescribed constants$\epsilon \in \left(0,1\right)$$k_i\ge 1$for$i=1,2,3$. The solution trajectory converges to the disease-free equilibrium point at exponential rate. Then, one gets by combining (10)–(12) and (18) that:

$$\left|{\displaystyle {\int }_0^{t_{I1}}\left(\gamma -\beta S\left(\tau \right)\right)d\tau }\right|=\left|{\displaystyle {\int }_0^{t_{I1}}\left(\gamma -\beta e^{-\beta {\displaystyle {\int }_0^{\sigma }I\left(\sigma \right)d\sigma }}{S}_0\right)d\tau }\right|\le ln{I}_0-ln{k}_1+\left|ln\epsilon \right|; \forall t\in {\mathit{R}}_{0+}$$

(19)

$$\left(\gamma -\beta e^{-\beta {\displaystyle {\int }_0^{t_{I2}}I\left(\tau \right)d\tau }}{S}_0\right)e^{-{\displaystyle {\int }_0^{t_{I2}}\left(\gamma -\beta e^{-\beta {\displaystyle {\int }_0^{\tau }I\left(\sigma \right)d\sigma }}{S}_0\right)d\tau }}{I}_0\le k_2\epsilon ; \forall t\in {\mathit{R}}_{0+}$$

(20)

$${\left[{\beta }^2\hspace{0.17em}S\left(t\right)I\left(t\right)-{\left(\beta S\left(t\right)-\gamma \right)}^2\right]}^{}I\left(t\right)\le k_3\epsilon ¸ \forall t\in {\mathit{R}}_{0+}$$

(21)

implying that:

$$t_{I1}=min\left(t\in {\mathit{R}}_{0+}\hspace{0.17em}:\gamma t-\beta S_0{\displaystyle {\int }_0^t{e}^{-\beta {\displaystyle {\int }_0^{\sigma }I\left(\sigma \right)d\sigma }}d\tau }=ln{I}_0-ln{k}_1+\left|ln\epsilon \right|\right)\ge \frac{1}{\gamma }\left(ln{I}_0-ln{k}_1+\left|ln\epsilon \right|\right)$$

(22)

$$\begin{array}{c}2min\left(k_2\epsilon I_0,\hspace{0.17em}\beta S_0\right)e^{-\beta {\displaystyle {\int }_0^{t_{I2}}I\left(\tau \right)d\tau }}\\ \le \eta \left(t_{I2}\right)=k_2\epsilon e^{{\displaystyle {\int }_0^{t_{I2}}\left(\gamma -\beta S\left(\tau \right)\right)d\tau }}{I}_0+\beta e^{-\beta {\displaystyle {\int }_0^{t_{I2}}I\left(\tau \right)d\tau }}{S}_0\\ \le \left(k_2\epsilon I_0+\beta S_0\right)e^{{\displaystyle {\int }_0^{t_{I2}}\left(\gamma -\beta S\left(\tau \right)\right)d\tau }}{I}_0\end{array}$$

(23)

which leads to:

$$\begin{array}{c}{e}^{{\displaystyle {\int }_0^{t_{I2}}\left(\gamma -\beta S\left(\tau \right)+\beta I\left(\tau \right)\right)d\tau }}\ge \frac{2min\left(k_2\epsilon I_0,\hspace{0.17em}\beta S_0\right)}{\left(k_2\epsilon I_0+\beta S_0\right)I_0}\\ \Rightarrow t_{I2}\ge max\left(t>0:{\displaystyle {\int }_0^t\left(\gamma -\beta S\left(\tau \right)+\beta I\left(\tau \right)\right)d\tau }\right)=ln\left[\frac{2min\left(k_2\epsilon I_0,\hspace{0.17em}\beta S_0\right)}{\left(k_2\epsilon I_0+\beta S_0\right)I_0}\right]\end{array}$$

(24)

$$\begin{array}{c}{e}^{{\displaystyle {\int }_0^{t_{I2}}\left(\beta S\left(\tau \right)-\beta I\left(\tau \right)-\gamma \right)d\tau }}\le \frac{\left(k_2\epsilon I_0+\beta S_0\right)I_0}{2min\left(k_2\epsilon I_0,\hspace{0.17em}\beta S_0\right)}\hfill \\ \hfill \Rightarrow t_{I2}\le min\left(t>0:{\displaystyle {\int }_0^t\left(\beta S\left(\tau \right)-\beta I\left(\tau \right)-\gamma \right)d\tau }\right)=ln\left[\frac{\left(k_2\epsilon I_0+\beta S_0\right)I_0}{2min\left(k_2\epsilon I_0,\hspace{0.17em}\beta S_0\right)}\right]\end{array}$$

(25)

and:

$$-k_3\epsilon \le \ddot{I}\left(t\right)=\left(\beta S\left(t\right)-\gamma \right)\dot{I}\left(t\right)+\beta \dot{S}\left(t\right)I\left(t\right)={\left[{\left(\gamma -\beta S\left(t\right)\right)}^2-{\beta }^{2}S\left(t\right)I\left(t\right)\right]}^{}I\left(t\right)\le k_3\epsilon$$

what implies that $\left|\ddot{I}\left(t\right)\right|\le k_3\epsilon$; $\forall t\in \left[t_{I3},\hspace{0.17em}\infty \right)$ such that:

$$t_{I3}\ge max\left(t>0:{\left[{\left(\gamma -\beta S\left(t\right)\right)}^2-{\beta }^{2}S\left(t\right)I\left(t\right)\right]}^{}I\left(t\right)\right)\ge -k_3\epsilon ,$$

$$t_{I3}\le min\left(t>0:{\left[{\left(\gamma -\beta S\left(t\right)\right)}^2-{\beta }^{2}S\left(t\right)I\left(t\right)\right]}^{}I\left(t\right)\right)\le k_3\epsilon$$

Example 4.

Consider the following SIS model with vaccination and antiviral or antibiotic controls:

$$\dot{S}\left(t\right)=\gamma I\left(t\right)-\beta S\left(t\right)I\left(t\right)-k_{V}S\left(t\right); \dot{I}\left(t\right)=\left(\beta S\left(t\right)-\gamma -k_T\right)I\left(t\right); \forall t\in {\mathit{R}}_{0+}$$

(26)

subject to$S\left(0\right)=S_0$,$I\left(0\right)=I_0$with$min\left(S_0,\hspace{0.17em}{I}_0\hspace{0.17em}\right)\ge 0$where the vaccination and treatment feedback controls on the susceptible and infectious are, respectively,$V\left(t\right)=k_{V}S\left(t\right)$and$T\left(t\right)=k_{T}I\left(t\right)$with$min\left(k_V,\hspace{0.17em}{k}_T\right)\ge 0$. If it is assumed that the total population$N\left(t\right)=N_0=S_0+I_0$;$\forall t\in {\mathit{R}}_{0+}$is constant through time then there is a complementary recovered (or immune) subpopulation present which obeys the differential equation$\dot{R}\left(t\right)=k_{V}S\left(t\right)+k_{T}I\left(t\right)$with$R\left(0\right)=R_0=0$. The solution is:

$$\begin{array}{c}S\left(t\right)=e^{-{\displaystyle {\int }_0^t\left(\beta I\left(\tau \right)+k_V\right)d\tau }}{S}_0+\gamma {\displaystyle {\int }_0^t{e}^{-{\displaystyle {\int }_{\tau }^t\left(\beta I\left(\sigma \right)+k_V\right)d\sigma }}I\left(\tau \right)}d\tau \\ =e^{-k_{V}t}{S}_0-{\displaystyle {\int }_0^t{e}^{-k_V\left(t-\tau \right)}}\left(\beta S\left(\tau \right)-\gamma \right)I\left(\tau \right)d\tau \end{array}$$

(27)

$$I\left(t\right)=e^{\beta {\displaystyle {\int }_0^{t}S\left(\tau \right)d\tau }}{e}^{-\left(\gamma +\hspace{0.17em}{k}_T\right)t}{I}_0$$

(28)

$$R\left(t\right)={\displaystyle {\int }_0^t\left(k_{V}S\left(\tau \right)+k_{T}I\left(\tau \right)\right)}d\tau$$

(29)

The following result links the above SIS model with a complementary recovered subpopulation to the generic one (1) under a minimum number of initial susceptible and sufficiently large number of initial infectious with initial growing rate.

Theorem 4.

Assume that$S_0>\hspace{0.17em}\frac{\gamma +k_T}{\beta }$,$I_0<1+\frac{1}{\gamma }\left(k_T+k_V{S}_0\right)$and${\dot{I}}_0>\frac{\beta \left|{\dot{S}}_0\right|I_0}{\beta S_0-\gamma -k_T}$.

Then, the following properties hold:

(i)

${\dot{S}}_0<0$and${\ddot{I}}_0>0$,

(ii)

$S\left(t\right)$is strictly decreasing on$\left[0,\hspace{0.17em}{t}_{Smin}\right]$with$t_{Smin}=min\left(t\in {\mathit{R}}_{0+}\hspace{0.17em}:\hspace{0.17em}S\left(t\right)=\gamma /\beta \right)$,

(iii)

$I\left(t\right)$is strictly increasing on$\left[0,\hspace{0.17em}{t}_{Imax}\right]$, and

$I_{max}=I\left(t_{max}\right)=max\left(I\left(t\right)\hspace{0.17em}:\hspace{0.17em}\hspace{0.17em}t\in \left[0\hspace{0.17em},\hspace{0.17em}{t}_{Imax}\right],\hspace{0.17em}{t}_I{}_{max}=min\left(t\in R_{0+}:S\left(t\right)=\left(\gamma +k_T\right)/\beta \right)\right)$with$t_{Imax}\ge t_{Smin}$,

(iv)

There is$t_{und}<t_{Imax}$which is an undulation and, furthermore, strict inflection time instant of$I\left(t\right)$,

(v)

Assume, in addition, that$I_0$is large enough to satisfy$I_0>\frac{\left(\gamma +k_T\right)k_V}{\left(\gamma -\beta \left(\gamma +k_T\right)\right)}{e}^{-\beta {\displaystyle {\int }_0^{t_{Imax}}S\left(\tau \right)d\tau }}{e}^{\left(\gamma +\hspace{0.17em}{k}_T\right)t_{Imax}}$. Then, the epidemic model (26) can be written in the form (1) on$\left[0,\hspace{0.17em}{t}_{Imax}\right]$with the following function$\alpha :\left[0,\hspace{0.17em}{t}_{Imax}\right]\to {\mathit{R}}_{0+}$:

$$\alpha \left(t\right)=\beta \left(e^{-{\displaystyle {\int }_0^t\left(\beta I\left(\tau \right)+k_V\right)d\tau }}{S}_0+\gamma {\displaystyle {\int }_0^t{e}^{-{\displaystyle {\int }_{\tau }^t\left(\beta I\left(\sigma \right)+k_V\right)d\sigma }}I\left(\tau \right)}d\tau \right)-\gamma -k_T; t\in \left[0,\hspace{0.17em}{t}_{Imax}\right]$$

(30)

which is of the form$\alpha \left(t\right)=-\frac{cln\left(g\left(t\right)/E\right)}{h\left(t\right)}$with$g:\left[0,\hspace{0.17em}{t}_{Imax}\right]\to \left[0,\hspace{0.17em}E\right]$;$\forall t\in \left[0,\hspace{0.17em}{t}_{Imax}\right]$and any given$E\in {\mathit{R}}_{+}$and$h\left(t\right)=\frac{c\left|ln\left(g\left(t\right)/E\right)\right|}{\beta \left(e^{-{\displaystyle {\int }_0^t\left(\beta I\left(\tau \right)+k_V\right)d\tau }}{S}_0+\gamma {\displaystyle {\int }_0^t{e}^{-{\displaystyle {\int }_{\tau }^t\left(\beta I\left(\sigma \right)+k_V\right)d\sigma }}I\left(\tau \right)}d\tau \right)-\gamma -k_T}$;$\forall t\in \left[0,\hspace{0.17em}{t}_{Imax}\right]$.

(vi)

The equilibrium points are$S_1^{*}=I_1^{*}=0$,$R_1^{*}=N_0$if$k_V\ne 0$and$k_T\ge 0$, and$S_2^{*}=\frac{\gamma +k_T}{\beta }$,$I_2^{*}=0$and$R_2^{*}=N_0-\frac{\gamma +k_T}{\beta }$which is only reachable if$k_V=0$since, otherwise,$I_2^{*}=-\frac{k_V}{k_T}\frac{\gamma +k_T}{\beta }<0$.

Proof. 

Since $S_0>\hspace{0.17em}\frac{\gamma +k_T}{\beta }$ and $I_0<1+\frac{1}{\gamma }\left(k_T+k_V{S}_0\right)$ then ${\dot{I}}_0>0$ and ${\dot{S}}_0<0$. Also, ${\ddot{I}}_0=\beta {\dot{S}}_0{I}_0+\left(\beta S_0-\gamma -k_T\right){\dot{I}}_0=\left(\beta S_0-\gamma -k_T\right){\dot{I}}_0-\beta \left|{\dot{S}}_0\right|I_0>0$ if ${\dot{I}}_0>\frac{\beta \left|{\dot{S}}_0\right|I_0}{\beta S_0-\gamma -k_T}$. Property (i) has been proved. Furthermore, $S_0>\hspace{0.17em}\frac{\gamma +k_T}{\beta }\ge \frac{\gamma }{\beta }$ implies from (27) that $S\left(t\right)$ is strictly decreasing on $\left[0,\hspace{0.17em}t\prime \right]$ where $t\prime =min\left(t\in {\mathit{R}}_{0+}\hspace{0.17em}:\hspace{0.17em}S\left(t\right)=\gamma /\beta \right)$ what proves Property (ii) with $t_{Smin}=t\prime$. On the other hand and since $S:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is continuous, there exists some $t″\in \left[0,\hspace{0.17em}t\prime \right]$ such that $S\left(t″\right)=\hspace{0.17em}\frac{\gamma +k_T}{\beta }$ with $t″=t\prime$ if and only if $k_T=0$. From (26), $\dot{I}\left(t^{″}\right)=0$ and $\dot{I}\left(t\right)>0$ for $t\in \left[0,t″\right)$ since ${\dot{I}}_0>0$. On the other hand, one has from (26) and (28) that:

$$\ddot{I}\left(t″\right)=\left(\beta S\left(t″\right)-\gamma -k_T\right)\dot{I}\left(t″\right)+\beta \dot{S}\left(t″\right)I\left(t″\right)$$

$$=\beta \left[\beta \left(\gamma -\beta S\left(t″\right)\right)I\left(t″\right)-k_{V}S\left(t″\right)\right]I\left(t″\right)$$

$$=-{\left[{\beta }^2{k}_{T}I\left(t″\right)+k_V\left(\gamma +k_T\right)\right]}^{}I\left(t″\right)$$

$$=-{\left[{\beta }^2{k}_T{e}^{\beta {\displaystyle {\int }_0^{t}S\left(\tau \right)d\tau }}{e}^{-\left(\gamma +\hspace{0.17em}{k}_T\right)t}{I}_0+k_V\left(\gamma +k_T\right)\right]}^{}{e}^{\beta {\displaystyle {\int }_0^{t″}S\left(\tau \right)d\tau }}{e}^{-\left(\gamma +\hspace{0.17em}{k}_T\right)t″}{I}_0<0$$

and $I\left(t\right)$ has a relative maximum $I_{max}$ at $t=t″=t_{Imax}$ which is also the absolute maximum on $\left[0,t_{max}\right]$. Property (iii) has been proved. Note also that since $\ddot{I}\left(t\right)$ is continuous and ${\ddot{I}}_0>0$, there exists some $t_{und}<t″$ such that $t_{und}$ is an undulation point of $I\left(t\right)$. Note furthermore that

$$\ddot{I}\left(t_{und}\right)=\left(\beta S\left(t_{und}\right)-\gamma -k_T\right)\dot{I}\left(t_{und}\right)+\beta \dot{S}\left(t_{und}\right)I\left(t_{und}\right)=0$$

From Lemma 1(i), $\ddot{I}\left(t_{und}-\epsilon \right)\ddot{I}\left(t_{und}+\epsilon \right)<0$; $\forall \epsilon \in B\left(0,r\right)$ and some $r\in {\mathit{R}}_{+}$ implies that $t_{und}$ is also an inflection time instant of $I\left(t\right)$. The equivalent logic contrapositive proposition establishes that:

$$\left[\forall r\in {\mathit{R}}_{+},\hspace{0.17em}\exists \hspace{0.17em}\epsilon \in \left[0,\hspace{0.17em}r\right]:\ddot{I}\left(t_{und}-\epsilon \right)\ddot{I}\left(t_{und}+\epsilon \right)\ge 0\right]\Rightarrow [t_{und} \mathrm{is} \mathrm{not} \mathrm{an} \mathrm{inflection} \mathrm{time} \mathrm{instant} \mathrm{of} I\left(t\right)]$$

Then, if $\ddot{I}\left(t_{und}-\epsilon \right)\ddot{I}\left(t_{und}+\epsilon \right)<0$; $\forall \epsilon \in B\left(0,r\right)$ and some $r\in {\mathit{R}}_{+}$ then $t_{und}$ is in fact an inflection time instant of $I\left(t\right)$. Assume that there is some arbitrarily small $\epsilon \in {\mathit{R}}_{+}$ such that $\ddot{I}\left(t_{und}-\epsilon \right)\ddot{I}\left(t_{und}+\epsilon \right)\ge 0$

Then:

$\dot{I}\left(t_{und}+\epsilon \right)=\ddot{I}\left(t_{und}\right)+{\displaystyle {\int }_0^{\epsilon }\ddot{I}\left(t_{und}+\tau \right)}d\tau$; $\dot{I}\left(t_{und}-\epsilon \right)=\ddot{I}\left(t_{und}\right)+{\displaystyle {\int }_0^{-\epsilon }\ddot{I}\left(t_{und}+\tau \right)}d\tau$.

Since $\ddot{I}\left(t\right)$ is continuous on $\left[t_{und}-\epsilon ,\hspace{0.17em}{t}_{und}+\epsilon \right]$ and one gets that

$$\dot{I}\left(t_{und}+\epsilon \right)-\dot{I}\left(t_{und}-\epsilon \right)={\displaystyle {\int }_0^{\epsilon }\ddot{I}\left(t_{und}+\tau \right)}d\tau -{\displaystyle {\int }_0^{-\epsilon }\ddot{I}\left(t_{und}+\tau \right)}d\tau$$

It is known that $0<{\epsilon }_I\le \dot{I}\left(t_{und}\right)<{\dot{I}}_0$ so that, for some arbitrarily small $\epsilon \in {\mathit{R}}_{+}$ such that $\ddot{I}\left(t_{und}-\epsilon \right)\ddot{I}\left(t_{und}+\epsilon \right)\ge 0$, there are ${\epsilon }_1\in \left[0,\hspace{0.17em}\epsilon \right]$ and ${\epsilon }_2\in {\mathit{R}}_{+}$ with $-{\epsilon }_2\in \left[-\epsilon ,\hspace{0.17em}0\right]$ such that the following joint constraints hold:

(1)

$\dot{I}\left(t_{und}+\tau \right)>0$; $\forall \tau \in \left[-{\epsilon }_2,\hspace{0.17em}{\epsilon }_1\right]\subset \left[-\epsilon ,\hspace{0.17em}\epsilon \right]$ with $\dot{I}\left(t\right)$ being strictly increasing on $\left[-{\epsilon }_2,\hspace{0.17em}{\epsilon }_1\right]$.

(2)

${\displaystyle {\int }_0^{{\epsilon }_1}\ddot{I}\left(t_{und}+\tau \right)}d\tau ={\displaystyle {\int }_0^{-{\epsilon }_2}\ddot{I}\left(t_{und}+\tau \right)}d\tau$

Then, one gets from Condition 2 that:

$$\dot{I}\left(t_{und}+{\epsilon }_1\right)-\dot{I}\left(t_{und}-{\epsilon }_2\right)={\displaystyle {\int }_0^{{\epsilon }_1}\ddot{I}\left(t_{und}+\tau \right)}d\tau -{\displaystyle {\int }_0^{-{\epsilon }_2}\ddot{I}\left(t_{und}+\tau \right)}d\tau =0$$

so that $\dot{I}\left(t\right)$ is not strictly increasing on $\left[-{\epsilon }_2,\hspace{0.17em}{\epsilon }_1\right]$, hence a contradiction. As a result, the undulation time instant $t_{und}$ of $I\left(t\right)$ is also a strict inflection time instant of $I\left(t\right)$ since $\dot{I}\left(t_{und}\right)\ne 0$ since Lemma 1 (ii) holds and the first zero of $\dot{I}\left(t\right)$ occurs at $t_{Imax}>t_{und}$. Property (iv) has been proved. To prove Property (v), note that Equation (30) follows from (26)–(27). Now, we equalize (30) to (1) to get admissible functions $g,h:{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ leading to:

$$\alpha \left(t\right)=\beta \left(e^{-{\displaystyle {\int }_0^t\left(\beta I\left(\tau \right)+k_V\right)d\tau }}{S}_0+\gamma {\displaystyle {\int }_0^t{e}^{-{\displaystyle {\int }_{\tau }^t\left(\beta I\left(\sigma \right)+k_V\right)d\sigma }}I\left(\tau \right)}d\tau \right)-\gamma -k_T=-\frac{cln\left(g\left(t\right)/E\right)}{h\left(t\right)}$$

(31)

and note that $\alpha \left(0\right)=\beta S_0-\gamma -k_T>0$. Note also that $\alpha \left(0\right)=\frac{+\infty }{h\left(0\right)}$ from the use of (31) in (30) implies that $h\left(0\right)=0$ irrespective of $g\left(t\right)$ while $g\left(t\right)$ is chosen arbitrary and continuous time-differentiable subject to $g\left(0\right)=0$ and $\alpha \left(t_{Imax}\right)=0$, $g\left(t_{Imax}\right)=E$ (so that $ln\left(g\left(t_{Imax}\right)/E\right)=0$) with $h\left(t\right)=\frac{c/E}{\hspace{0.17em}\beta \gamma I\left(t\right)-\beta \left(\beta I\left(t\right)+k_V\right)S\left(t\right)}$ for $t\in \left[0,\hspace{0.17em}{t}_{Imax}\right]$.

Now, note that $h\left(t_{Imax}\right)$ is a primary $\left(0/0\right)$—type indetermination which is resolved through L´H $\hat{\mathrm{o}}$ pital rule leading to:

$$\begin{array}{c}h\left(t_{Imax}\right)=\frac{c/g\left(t_{Imax}\right)}{\beta \dot{S}\left(t_{Imax}\right)}=\frac{c/E}{\hspace{0.17em}\beta \gamma I\left(t_{Imax}\right)-\beta \left(\beta I\left(t_{Imax}\right)+k_V\right)S\left(t_{Imax}\right)}\\ =\frac{c/\left(\beta E\right)}{\hspace{0.17em}\gamma I\left(t_{Imax}\right)-\left(\beta I\left(t_{Imax}\right)+k_V\right)\left(\gamma +k_T\right)}\end{array}$$

Since $I\left(t_{Imax}\right)=e^{\beta {\displaystyle {\int }_0^{t_{Imax}}S\left(\tau \right)d\tau }}{e}^{-\left(\gamma +\hspace{0.17em}{k}_T\right)t}{I}_0$ then for sufficiently large $I_0$ such that

$$I_0>\frac{\left(\gamma +k_T\right)k_V}{\left(\gamma -\beta \left(\gamma +k_T\right)\right)}{e}^{-\beta {\displaystyle {\int }_0^{t_{Imax}}S\left(\tau \right)d\tau }}{e}^{\left(\gamma +\hspace{0.17em}{k}_T\right)t_{Imax}}$$

then:

$$h\left(t\right)=\frac{c\left|ln\left(g\left(t\right)/E\right)\right|}{\beta \left(e^{-{\displaystyle {\int }_0^t\left(\beta I\left(\tau \right)+k_V\right)d\tau }}{S}_0+\gamma {\displaystyle {\int }_0^t{e}^{-{\displaystyle {\int }_{\tau }^t\left(\beta I\left(\sigma \right)+k_V\right)d\sigma }}I\left(\tau \right)}d\tau \right)-\gamma -k_T}$$

$$=\frac{cln\left(g\left(t\right)/E\right)}{\gamma +k_T-\beta \left(e^{-k_{V}t}{S}_0-{\displaystyle {\int }_0^t{e}^{-k_V\left(t-\tau \right)}}\left(\beta S\left(\tau \right)-\gamma \right)I\left(\tau \right)d\tau \right)}$$

fulfilling, in particular:

$$h\left(t_{Imax}\right)=\frac{c/\left(\beta E\right)}{\hspace{0.17em}\left(\gamma -\beta \left(\gamma +k_T\right)\right)I{t}_{Imax}-\left(\gamma +k_T\right)k_V}$$

$$=\frac{c/\left(\beta E\right)}{\hspace{0.17em}\left(\gamma -\beta \left(\gamma +k_T\right)\right)e^{\beta {\displaystyle {\int }_0^{t_{Imax}}S\left(\tau \right)d\tau }}{e}^{-\left(\gamma +\hspace{0.17em}{k}_T\right)t_{Imax}}{I}_0-\left(\gamma +k_T\right)k_V}>0$$

Property (v) has been proved. Property (vi) is obvious by zeroing (26). □

Example 4 is tested numerically in the sequel with the following data β = 30, γ = 50 years⁻¹, implying that the average infectious period is T_γ = 365/50 = 7.3 days, k_V = 1 and k_T = 50. The time scale of the figures is in a scale of years accordingly with the above numerical values. In Figure 1, the solution trajectories of all the subpopulation are shown with the constraints of Theorem 4 being fulfilled by the initial conditions, in particular $S_0>\hspace{0.17em}\frac{\gamma +k_T}{\beta }$, $I_0=1-S_0$ and $R_0=0$ so that $N_0$ is normalized to unity. It is seen that the infectious subpopulation trajectory has a maximum at a finite time and that the state trajectory solution converges asymptotically to an endemic equilibrium point. In Figure 2, the state trajectory solution is shown with $N_0=1$ when $S_0=\left(\gamma +k_T\right)/\beta$ which violates the conditions of Theorem 4 with ${\dot{I}}_0=0$. In this case, there is no relative maximum of the infectious subpopulation at finite time. In both situations, it has been observed by extending the overall simulation time that the susceptible and the infectious subpopulations converge asymptotically to zero while the recovered subpopulation converges to unity as time tends to infinity. The controls are suppressed in Figure 3 with N₀ = 1. In this case, the recovered subpopulation may be deleted from the model since it is unnecessary while being identically zero. The infectious and susceptible subpopulations are in an endemic equilibrium point for all time so that the infection results to be permanent in the sense that it cannot be asymptotically removed. See Theorem 4(vi) for the case k_V= 0. Figure 4 exhibits a trajectory solution which agrees with Theorem 4 while there is no normalization of the initial conditions to unity. In this case, the maximum of the infectious subpopulation at a finite time becomes very apparent.

4. Links with Entropy and Maximum Dissipation Mechanism Issues

4.1. Comparison of the Epidemic Model and Reference Model Information Entropies

Since (1) is a scalar equation, a valid solution for the particular model-dependent time-varying coefficient $\alpha \left(t\right)=-cln\left(g\left(t\right)/E\right)/h\left(t\right)$ of Theorem 2 and Theorem 3 is, according to Theorem 1:

$$I\left(t\right)=e^{-c{\displaystyle {\int }_0^t{h}^{-1}\left(\tau \right)\hspace{0.17em}ln\left(g\left(\tau \right)/E\right)d\tau }}{I}_0; t\in {\mathit{R}}_{0+}$$

(32)

Under the particular constraints $E=D$, $c=\left(1-ln\left(L/D\right)\right)/{ln}^2\left(L/D\right)$ and $g\left(t\right)=h\left(t\right)=t$, it is got in [11] that $\alpha \left(t\right)={\left[\left(ln\left(L/D\right)-1\right)/{ln}^2\left(L/D\right)\right]}^{}{t}^{-1}ln\left(t/D\right)$ and (32), namely:

$$I_p\left(t\right)=e^{\left(ln\left(L/D\right)-1\right)/{ln}^2\left(L/D\right){\displaystyle {\int }_0^t{\tau }^{-1}ln\left(\tau /D\right)d\tau }}{I}_0; t\in {\mathit{R}}_{+}$$

(33)

approaches the log-normal distribution:

$$I_r\left(t\right)=\frac{k}{\sqrt{2\pi }\hspace{0.17em}{\sigma }_{r}t}{e}^{-\frac{{\left(ln\hspace{0.17em}t-{\mu }_r\right)}^2}{2{\sigma }_r^2}}; t\in {\mathit{R}}_{+}$$

(34)

for reference values $D=D_r$ and $L=L_r$ of the maximum and inflection reference time instants where ${\mu }_r=ln{D}_r+{\sigma }_r^2$ and ${\sigma }_r$ is given by the principle of extreme entropy production rate, typically ${\sigma }_r\approx 0.408$ gives the width of the distribution function for the maximum dissipation rate for the usual definition of the Shannon entropy. The main reason for the limitation of such a width is that the medical and social interventions are a dissipation mechanism which controls and limits the disease propagation. Comparing (33) and (34), one gets that $k=\sqrt{2\pi }\hspace{0.17em}{\sigma }_r^3{I}_0$ after solving the indetermination $0/0$ at $t=0$ via L´ $\mathrm{H}\hat{\mathrm{o}}\mathrm{pital}$ rule leading to the “infection reference evolution” $I_r\left(t\right)=I_p\left(t\right)$, that is by equalizing (23) and (24), under the above set of particular constraints, where:

$$I_r\left(t\right)=t^{-1}{\sigma }_r^2{e}^{-\frac{{\left(lnt-ln{D}_r-{\sigma }_r^2\right)}^2}{2{\sigma }_r^2}}{I}_0; t\in {\mathit{R}}_{+}$$

(35)

Now, equalize $I\left(t\right)=I_r\left(t\right)+\tilde{I}\left(t\right)$; $\forall t\in {\mathit{R}}_{+}$ for some perturbation function $\tilde{I}:{\mathit{R}}_{+}\to {\mathit{R}}_{0+}$ resulting to be from (32) and (35) for $I_0>0$:

$$\tilde{I}\left(t\right)=\left(e^{-c{\displaystyle {\int }_0^t{h}^{-1}\left(\tau \right)\hspace{0.17em}ln\left(g\left(\tau \right)/E\right)d\tau }}-t^{-1}{\sigma }_r^2{e}^{-\frac{{\left(lnt-ln{D}_r-{\sigma }_r^2\right)}^2}{2{\sigma }_r^2}}\right)I_0; t\in {\mathit{R}}_{+}$$

(36)

The Shannon entropy of the infection $S_I\left(\eta \right)$ results to be given by the following Riemann- Stieljes integral which quantifies the entropy error ${\tilde{S}}_I\left(\eta \right)$ of that associated with any given model related to the entropy of the “infection reference evolution” given by the log- normal function $S_{I_r}\left(\eta \right)=S_{I_r}\left(\eta ,{\sigma }_r\right)$ for the given reference width value ${\sigma }_r=\sqrt{1/2\eta }$:

$$\begin{array}{c}{S}_I\left(\eta \right)=-{\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}I\left(t\right)ln\left(t^{1-\eta }I\left(t\right)\right)d{t}^{\eta }\hfill \\ =-{\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}I\left(t\right)\left(\left(1-\eta \right)lnt+lnI\left(t\right)\right)d{t}^{\eta }\hfill \\ =-{\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}\left(I_r\left(t\right)+\tilde{I}\left(t\right)\right)\left(\left(1-\eta \right)lnt+ln\left(I_r\left(t\right)\left(1+I_r^{-1}\left(t\right)\tilde{I}\left(t\right)\right)\right)\right)d{t}^{\eta }\hfill \\ =-{\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}\left(I_r\left(t\right)+\tilde{I}\left(t\right)\right)\left(\left(1-\eta \right)lnt+ln{I}_r\left(t\right)+ln\left(1+I_r^{-1}\left(t\right)\tilde{I}\left(t\right)\right)\right)d{t}^{\eta }\hfill \\ =S_{I_r}\left(\eta \right)-{\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}{I}_r\left(t\right)ln\left(1+I_r^{-1}\left(t\right)\tilde{I}\left(t\right)\right)d{t}^{\eta }-{\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}\tilde{I}\left(t\right)lnI\left(t\right)d{t}^{\eta }-\left(1-\eta \right){\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}\tilde{I}\left(t\right)lntd{t}^{\eta }\hfill \\ =S_{I_r}\left(\eta \right)-{\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}{I}_r\left(t\right)ln\left(1+I_r^{-1}\left(t\right)\tilde{I}\left(t\right)\right)d{t}^{\eta }-\left(1-\eta \right){\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}\tilde{I}\left(t\right)lntd{t}^{\eta }\hfill \\ +{\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}{I}_r\left(t\right)\left(1-I_r^{-1}\left(t\right)I\left(t\right)\right)lnI\left(t\right)d{t}^{\eta }\hfill \\ =S_{I_r}\left(\eta \right)+{\tilde{S}}_I\left(\sigma \right); t\in {\mathit{R}}_{+}\hfill \end{array}$$

(37)

after using $I\left(t\right)=I_r\left(t\right)\hspace{0.17em}\left(1+I_r^{-1}\left(t\right)\tilde{I}\left(t\right)\right)$ and its equivalent expression $\tilde{I}\left(t\right)=-I_r\left(t\right)\left(1-I_r^{-1}\left(t\right)I\left(t\right)\right)$, where the reference entropy based on the identification of the log-normal function (34)with the solution of (1), that is, (33), yields for ${\sigma }_r=\sqrt{1/2\eta }$:

$$\begin{array}{c}{S}_{I_r}\left(\hspace{0.17em}\eta \right)=-{\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}{I}_r\left(t\right)ln\left(t^{1-\eta }{I}_r\left(t\right)\right)d{t}^{\eta }\\ =\eta \left(ln\left(\sqrt{\frac{\pi }{\eta }}\right)+\eta \left(ln{D}_r+\frac{1}{2\eta }\right)+\frac{1}{2}\right)\end{array}$$

(38)

after converting the Riemann-Stieljes integral (39) in a Riemann integral via differentiation of $d{t}^{\eta }$ by using (35). Note that it is assumed that both current and reference entropies are evaluated for the same parameter $\eta$ which is typically chosen as $\eta =3$. At the same time, it is assumed that the maximum dissipation rate proportional to the maximum rate of entropy production is governed by the width of the distribution function $\sigma$. So the current model can potentially have a value $\sigma \ne {\sigma }_r$. See [11] for the normalized case obtained for $I_0=1$, and, also one gets the following entropy error:

$$\begin{array}{c}{\tilde{S}}_I\left(\eta \right)=-{\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}{\left[ln\left({\left(I\left(t\right)/I_r\left(t\right)\right)}^{I_r\left(t\right)}\right)+ln\left(I{\left(t\right)}^{I\left(t\right)-I_r\left(t\right)}\right)\right]}^{}d{t}^{\eta }-\left(1-\eta \right){\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}\tilde{I}\left(t\right)lntd{t}^{\eta }\hfill \\ =-{\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}{\left[ln\left({\left(I\left(t\right)/I_r\left(t\right)\right)}^{I_r\left(t\right)}\left(I{\left(t\right)}^{I\left(t\right)-I_r\left(t\right)}\right)\right)\right]}^{}d{t}^{\eta }+\left(1-\eta \right){\displaystyle {\int }_0^{\infty }{t}^{1-\eta }}{I}_r\left(t\right)\left(1-I_r^{-1}\left(t\right)I\left(t\right)\right)lntd{t}^{\eta }\hfill \\ =-\eta {\displaystyle {\int }_0^{\infty }ln\left(\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}}\right)dt}+\eta \left(\eta -1\right){\displaystyle {\int }_0^{\infty }ln\left(t^{I_r\left(t\right)-I\left(t\right)}\right)dt}\hfill \\ =-\eta {\displaystyle {\int }_0^{\infty }ln\left(\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}\right)dt}; t\in {\mathit{R}}_{+}\hfill \end{array}$$

(39)

It turns out obvious that the integrand of (39) is identically zero if $\tilde{I}\left(t\right)\equiv 0$, so that $I\left(t\right)\equiv I_r\left(t\right)$, leading to ${\tilde{S}}_I\left(\eta \right)\equiv 0$. The expression (37), subject to (38)–(39), parameterizes the incremental entropy with the same parameter $\eta$ which parameterizes the reference entropy $S_{I_r}\left(\hspace{0.17em}{\eta }_r\right)$. Now, define the error:

$$\delta \left(t\right)=\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}-1; t\in {\mathit{R}}_{0+}$$

(40)

so that ${\tilde{S}}_I\left(\eta \right)\equiv 0$ if $\delta \left(t\right)\equiv 0$ and, expanding $ln\left(\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}}\right)$ via the Newton- Mercator series for the logarithm, leads to:

$$ln\left(\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}\right)=ln\left(1+\delta \left(t\right)\right)=\delta \left(t\right)+{\displaystyle {\sum }_{n=2}^{\infty }\frac{{\left(-1\right)}^{n+1}}{n}\hspace{0.17em}}{\delta }^n\left(t\right); t\in {\mathit{R}}_{0+}$$

(41)

and such a series converges to $ln\left(1+\delta \left(t\right)\right)$ for all $t\in {\mathit{R}}_{0+}$ provided that $\delta \left(t\right)\in \left(-1,\hspace{0.17em}1\right]$, equivalently, $\tilde{I}\left(t\right)\in \left(-I_r\left(t\right),\hspace{0.17em}{I}_r\left(t\right)\right]$; $\forall t\in {\mathit{R}}_{0+}$; $\forall t\in {\mathit{R}}_{0+}$. Thus, the following description in linear and higher-order additive terms of the entropy error follows from (40)–(41) into (39):

$${\tilde{S}}_I\left(\eta \right)={\tilde{S}}_{IL}\left(\eta \right)+{\tilde{\tilde{S}}}_I\left(\eta \right); t\in {\mathit{R}}_{0+}$$

(42)

where:

$${\tilde{S}}_{IL}\left(\eta \right)=-\eta {\displaystyle {\int }_0^{\infty }\left(\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}-1\right)dt}; t\in {\mathit{R}}_{0+}$$

(43)

$${\tilde{\tilde{S}}}_I\left(\eta \right)=-\eta \left({\displaystyle {\sum }_{n=2}^{\infty }\frac{{\left(-1\right)}^{n+1}}{n}\hspace{0.17em}}{\left(\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}-1\right)}^{n}dt\right); t\in {\mathit{R}}_{0+}$$

(44)

The subsequent results hold related to the case when the error between the infectious functions of the model and the reference one associated to the log-normal function converges asymptotically to zero as time tends to infinity. The first result, stated separately by convenience concerned its proof, discusses the simplest case for $\eta =1$.

Proposition 1.

Assume that$\eta =1$and$\underset{t\to +\infty \hspace{0.17em}}{lim}\left|{\displaystyle {\int }_0^{t}ln\left(\frac{I{\left(\tau \right)}^{I\left(\tau \right)}}{I_r{\left(\tau \right)}^{I_r\left(\tau \right)}}\right)d\tau }\right|<+\infty$.

then,${\tilde{S}}_I\left(\hspace{0.17em}1\right)<+\infty$for all$t\in {\mathit{R}}_{0+}$and$\underset{t\to +\infty }{lim}\left(I\left(t\right)-I_r\left(t\right)\right)=0$.

Proof. 

Note from (39) that ${\tilde{S}}_I\left(1\right)={\displaystyle {\int }_0^{\infty }ln\left(\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}}\right)dt}<+\infty$. Since the function $\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}}$ is uniformly continuous on ${\mathit{R}}_{0+}$ and $\underset{t\to +\infty \hspace{0.17em}}{lim}\left|{\displaystyle {\int }_0^{t}ln\left(\frac{I{\left(\tau \right)}^{I\left(\tau \right)}}{I_r{\left(\tau \right)}^{I_r\left(\tau \right)}}\right)d\tau }\right|<+\infty$ then $ln\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}}\to 0$ as $t\to +\infty$ from Barbalat´s lemma and then $\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}}\to 1$ as $t\to +\infty$. It is clear that a limit solution which satisfies this constraint is $\underset{t\to +\infty }{lim}\left(I\left(t\right)-I_r\left(t\right)\right)=0$. It is now proved that no alternative limiting constraint on the pair $\left(I\left(t\right),\hspace{0.17em}{I}_r\left(t\right)\right)$ as $t\to +\infty$ is compatible with $\underset{t\to +\infty }{lim}\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}}=1$. Assume that $\underset{t\to +\infty }{lim\hspace{0.17em}\inf }\left|\hspace{0.17em}I\left(t\right)-I_r\left(t\right)\right|>0$ It can happen that:

(a)

$\underset{t\to +\infty }{lim\hspace{0.17em}\inf }\left(I\left(t\right)-I_r\left(t\right)\right)>0$. Then, $\underset{t\to +\infty }{lim\hspace{0.17em}\inf }ln\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}}=\underset{t\to +\infty }{lim\hspace{0.17em}\inf }\left(I\left(t\right)lnI\left(t\right)-I_r\left(t\right)ln{I}_r\left(t\right)\right)$$>\underset{t\to +\infty }{lim\hspace{0.17em}\inf }\left(I_r\left(t\right)lnI\left(t\right)-I_r\left(t\right)ln{I}_r\left(t\right)\right)>\underset{t\to +\infty }{lim\hspace{0.17em}\inf }\left(I_r\left(t\right)ln{I}_r\left(t\right)-I_r\left(t\right)ln{I}_r\left(t\right)\right)=0$ so that $\underset{t\to +\infty }{lim\hspace{0.17em}\inf }ln\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}}>0$. Hence, a contradiction to Barbalat´s lemma; or

(b)

$\underset{t\to +\infty }{lim\hspace{0.17em}\inf }\left(I_r\left(t\right)-I\left(t\right)\right)>0$. Under a similar reasoning to that of a), one gets that $\underset{t\to +\infty }{lim\hspace{0.17em}\inf }ln\frac{I_r{\left(t\right)}^{I_r\left(t\right)}}{I{\left(t\right)}^{I\left(t\right)}}>0$. Again, a contradiction to Barbalat´s lemma.

The second result discusses the simplest case for $\eta \ne 1$. It is seen that the basic limit result $\underset{t\to +\infty }{lim}\left(I\left(t\right)-I_r\left(t\right)\right)=0$ of Proposition 1 is still kept under the reasonable assumption that the infection and reference infection functions are bounded. □

Proposition 2.

Assume that$\eta \ne 1$,$I,I_r:\hspace{0.17em}{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$are bounded and$\underset{t\to +\infty \hspace{0.17em}}{lim}\left|{\displaystyle {\int }_0^{t}ln\left(\frac{I{\left(\tau \right)}^{I\left(\tau \right)}}{I_r{\left(\tau \right)}^{I_r\tau }{t}^{\left(\eta -1\right)\left(I_r\left(\varsigma \right)-I\left(\tau \right)\right)}}\right)d\tau }\right|<+\infty$.

Then,${\tilde{S}}_I\left(\hspace{0.17em}\eta \right)<+\infty$for all$t\in {\mathit{R}}_{0+}$and$\underset{t\to +\infty }{lim}\left(I\left(t\right)-I_r\left(t\right)\right)=0$.

Proof. 

Note that ${\tilde{S}}_I\left(1\right)<+\infty$ and that, from the uniform continuity of $\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}$ everywhere in ${\mathit{R}}_{0+}$, the boundedness of its integral on $\left[0,\infty \right)$ and Barbalat´s lemma, it follows that $\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}\to 1$ as $t\to +\infty$ what implies that:

$$\underset{t\to +\infty }{lim}\left(I\left(t\right)lnI\left(t\right)-I_r\left(t\right)ln{I}_r\left(t\right)+\left(1-\eta \right)\left(I_r\left(t\right)-I\left(t\right)\right)\hspace{0.17em}ln\hspace{0.17em}t\right)=0$$

If $\eta >1$ and $lnt\to \infty$ as $t\to \infty$ then there exists some strictly increasing real sequence ${\left\{t_i\right\}}_{i=0}^{\infty }$, such that $\underset{k\to \infty }{lim}\left|\left(1-\eta \right)\left(I_r\left(\left(t_k\right)\right)-I\left(\left(t_k\right)\right)\right)\hspace{0.17em}ln\hspace{0.17em}{t}_k\right|=\infty$ with $t_k\in {\left\{t_i\right\}}_{i=0}^{\infty }$ if $\underset{t\to +\infty }{lim}\left(I\left(t\right)-I_r\left(t\right)\right)\ne 0$. But this can hold only if $\underset{k\to +\infty }{lim}\left|\hspace{0.17em}I\left(t_k\right)lnI\left(t_k\right)-I_r\left(t_k\right)ln{I}_r\left(t_k\right)\right|=+\infty$. But, since $I_r:\hspace{0.17em}{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is bounded for all time, this implies that $I\left(t_k\right)\to +\infty$ as $t_k\left(\in {\left\{t_i\right\}}_{i=0}^{\infty }\right)\to +\infty$ and $I:\hspace{0.17em}{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is unbounded. But then

$$\underset{k\to +\infty }{lim}\left(I\left(t_k\right)lnI\left(t_k\right)-I_r\left(t_k\right)ln{I}_r\left(t_k\right)+\left(\eta -1\right)\left(I\left(t_k\right)-I_r\left(t_k\right)\right)\hspace{0.17em}ln\hspace{0.17em}{t}_k\right)=\infty +\infty =\infty$$

and a contradiction follows to the above limit to be zero. As a result, $\underset{t\to +\infty }{lim}\left(I\left(t\right)-I_r\left(t\right)\right)=0$ if $\eta >1$.

Now, assume that $\eta <1$. Since $\underset{k\to \infty }{lim}\left|\left(1-\eta \right)\left(I_r\left(t_k\right)-I\left(t_k\right)\right)\hspace{0.17em}ln\hspace{0.17em}{t}_k\right|=\infty$ for $t_k\left(\in {\left\{t_i\right\}}_{i=0}^{\infty }\right)\to +\infty$ and some strictly increasing real sequence ${\left\{t_i\right\}}_{i=0}^{\infty }$, provided that $\underset{t\to +\infty }{lim}\left(I\left(t\right)-I_r\left(t\right)\right)\ne 0$, then $I\left(t_k\right)\to +\infty$ as $t_k\left(\in {\left\{t_i\right\}}_{i=0}^{\infty }\right)\to +\infty$ since $I_r:\hspace{0.17em}{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is bounded. Since $I:\hspace{0.17em}{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$ is unbounded, because it has a divergent subsequence ${\left\{I\left(t_k\right)\right\}}_{k=0}^{\infty }$ and it is a solution of a unstable time-invariant linear differential system, it is of positive exponential order ${\varsigma }_0>0$ and there exists a real constant $\varsigma <{\varsigma }_0$ such that $I\left(t_k\right)\ge e^{\varsigma t_k}$; $\forall t_k\in {\left\{t_i\right\}}_{i=0}^{\infty }$ and $I\left(t_k\right)/ln{t}_k\left(\ge e^{\varsigma t_k}/ln{t}_k\right)\to \infty$ as $t_k\left(\in {\left\{t_i\right\}}_{i=0}^{\infty }\right)\to \infty$ and, furthermore,

$$\underset{k\to +\infty }{lim}\left(I\left(t_k\right)lnI\left(t_k\right)-I_r\left(t_k\right)ln{I}_r\left(t_k\right)\right)=\left(1-\eta \right)\underset{k\to +\infty }{lim}\left(I\left(t_k\right)-I_r\left(t_k\right)\right)\hspace{0.17em}ln\hspace{0.17em}{t}_k=\infty$$

but the expression below is an infinity limit (and not a $\infty -\infty$ indetermination since $I\left(t_k\right)/ln{t}_k\to \infty$):

$$\underset{k\to +\infty }{lim}\left(I\left(t_k\right)lnI\left(t_k\right)-I_r\left(t_k\right)ln{I}_r\left(t_k\right)-\left(1-\eta \right)\left(I\left(t_k\right)-I_r\left(t_k\right)\right)\hspace{0.17em}ln\hspace{0.17em}{t}_k\right)=\infty$$

which contradicts:

$$\underset{t\to +\infty }{lim}\left(I\left(t\right)lnI\left(t\right)-I_r\left(t\right)ln{I}_r\left(t\right)+\left(1-\eta \right)\left(I_r\left(t\right)-I\left(t\right)\right)\hspace{0.17em}ln\hspace{0.17em}t\right)=0$$

As a result, $\underset{t\to +\infty }{lim}\left(I\left(t\right)-I_r\left(t\right)\right)=0$ if $\eta \ne 1$.

It is now briefly discussed the fact that the boundedness hypothesis of Proposition 2 is not very restrictive for some of the given examples, like for instance, Examples 2,3, where the infectious subpopulation converges asymptotically to zero. For such a purpose, note from $\left(35\right)$ that $I_r\left(t\right)\to 0$ exponentially fast as $t\to \infty$. In example 2, $I\left(t\right)\to 0$ exponentially as $t\to \infty$ so their difference function also converges to zero exponentially as $t\to \infty$. The integral boundedness invoked in the assumption of Proposition 2 is of the form $F=\left|\hspace{0.17em}{\displaystyle {\int }_0^{\infty }ln\hspace{0.17em}x\left(t\right)dt}\right|<+\infty$, where $x\left(t\right)=\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_{r}t-I\left(t\right)\right)}}$ is everywhere differentiable with respect to time. In order to convert the elevant Riemann-Stieljes integral into a standard Riemann one, take $dx=\dot{x}\left(t\right)dt$ and, later on, perform the change of variable $x\to u$ defined by $u=lnx$, $du=dx/x$ to yield:

$$F=\left|\hspace{0.17em}{\displaystyle {\int }_{x_0}^1\frac{ln\hspace{0.17em}x}{\dot{x}}dx}\right|=\left|\hspace{0.17em}{\displaystyle {\int }_{x_0}^1\frac{ln\hspace{0.17em}x}{x}\frac{x}{\dot{x}}dx}\right|\le \left(\underset{0\le t\le +\infty }{sup}\left|\frac{x\left(t\right)}{\dot{x}\left(t\right)}\right|\right)\left|\hspace{0.17em}{\displaystyle {\int }_{x_0}^1\frac{ln\hspace{0.17em}x}{x}dx}\right|$$

$$=M\left(\eta \right)\left|\hspace{0.17em}{\displaystyle {\int }_{ln{x}_0}^{0}udu}\right|=M\left(\eta \right)\frac{ln{x}_0^2}{2}$$

where $x\left(0\right)=x_0$ and $M\left(\eta \right)=\underset{0\le t\le +\infty }{sup}\left|\frac{x\left(t\right)}{\dot{x}\left(t\right)}\right|\le +\infty$ for the given constant $\eta$. Note that $M\left(\eta \right)<+\infty$ if and only if the set of zeros of $\dot{x}\left(t\right)$ at any finite time instant is empty, that is, if and only if $Z_{xdot}\left(\eta \right)=\varnothing$, where $Z_{xdot}\left(\eta \right)=\left\{t\ge 0\hspace{0.17em}: \dot{x}\left(t\right)=0\hspace{0.17em}\right\}=\varnothing$ (equivalently, $M\left(\eta \right)=+\infty$ if and only if $Z_{xdot}\left(\eta \right)\ne \varnothing$. Rewriting $x\left(t\right)=\frac{y\left(t\right)}{t^{\eta -1}z\left(t\right)}$ it follows that $\dot{x}\left(t\right)=0$ for any $t\ge 0$ if and only the following constraint holds $t=\frac{\left(\eta -1\right)z\left(t\right)y\left(t\right)}{z\left(t\right)\dot{y}\left(t\right)-y\left(t\right)\dot{z}\left(t\right)}$. Therefore, $Z_{xdot}\left(\eta \right)=\left\{t\ge 0\hspace{0.17em}: t=\frac{\left(\eta -1\right)z\left(t\right)y\left(t\right)}{z\left(t\right)\dot{y}\left(t\right)-y\left(t\right)\dot{z}\left(t\right)}\hspace{0.17em}\right\}\ne \varnothing$ is an event of zero probability. Thus, the boundedness hypothesis of Proposition 2 happens almost surely in the event that the infectious subpopulation converges asymptotically to zero as time tends to infinity. □

Propositions 1 and 2 yield the direct joint result independently of the value of $\eta$:

Proposition 3.

Assume that$\eta \in {\mathit{R}}_{0+}$,$I,I_r:\hspace{0.17em}{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$are bounded and$\underset{t\to +\infty \hspace{0.17em}}{lim}\left|{\displaystyle {\int }_0^{t}ln\left(\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}\right)dt}\right|<+\infty$.

Then,${\tilde{S}}_I\left(\hspace{0.17em}\eta \right)<+\infty$for all$t\in {\mathit{R}}_{0+}$and$\underset{t\to +\infty }{lim}\left(I\left(t\right)-I_r\left(t\right)\right)=0$.

Concerning Proposition 3, note that the boundedness of ${\tilde{S}}_I\left(\hspace{0.17em}\eta \right)$ does not guarantee that the linear part and the remaining part of higher- order terms in the decomposition of (42), subject to (43) and (44), are both finite. It could “a priori” happen that they both tend to infinity with opposite signs. But if any of them is bounded, the other one should be bounded as well according to Proposition 3. Fortunately, this does not happen under weak extra assumptions. In particular, the following result holds:

Proposition 4.

Assume that$\eta \in {\mathit{R}}_{0+}$,$I,I_r:\hspace{0.17em}{\mathit{R}}_{0+}\to {\mathit{R}}_{0+}$are bounded, and

$${\displaystyle {\int }_0^{\infty }ln\left(\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}-1\right)dt}<+\infty .$$

Then,${\tilde{S}}_{IL}\left(\eta \right)<+\infty$;$t\in {\mathit{R}}_{0+}$.

If, in addition,${\displaystyle {\int }_0^{\infty }ln\left(\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}dt\right)}<+\infty$then$\left|{\tilde{\tilde{S}}}_I\left(\eta \right)\right|<+\infty$and${\tilde{S}}_I\left(\eta \right)<+\infty$for all$t\in {\mathit{R}}_{0+}$and$\underset{t\to +\infty }{lim}\left(I\left(t\right)-I_r\left(t\right)\right)=0$.

Proof. 

It is direct to see that ${\tilde{S}}_{IL}\left(\eta \right)<+\infty$. Also, and again from Barbalat´s lemma, $\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}\to 1$ as $t\to +\infty$. Thus, from Proposition 3, $\underset{t\to +\infty }{lim}\left(I\left(t\right)-I_r\left(t\right)\right)=0$. If, furthermore, ${\displaystyle {\int }_0^{\infty }ln\left(\frac{I{\left(t\right)}^{I\left(t\right)}}{I_r{\left(t\right)}^{I_r\left(t\right)}{t}^{\left(\eta -1\right)\left(I_r\left(t\right)-I\left(t\right)\right)}}\right)dt}<+\infty$ then, again from Proposition 3, ${\tilde{S}}_I\left(\eta \right)<+\infty$ and $-\infty <-\left|\tilde{S}\left(\eta \right)-{\tilde{S}}_{IL}\left(\eta \right)\right|\le {\tilde{\tilde{S}}}_I\left(\eta \right)\le \left|\tilde{S}\left(\eta \right)-{\tilde{S}}_{IL}\left(\eta \right)\right|<+\infty$. □

Note that the above results agree with the asymptotic results of Examples 1–4, where $I\left(t\right)\to 0$ as $t\to \infty$, and with Theorem 1, since the reference $I_r\left(t\right)\to 0$, jointly implying $\left(I\left(t\right)-I_r\left(t\right)\right)\to 0$ as $t\to \infty$.

Remark 4.

The rationale behind the definition of a time-varying coefficient in (1) is to reduce the higher-order epidemic model with two or more states to a single-order differential equation based on the assumption that the log-normal distribution is a sufficiently accurate model for the infectious evolution. It is apparent that the profile of the log-normal distribution remembers the behavior of the strong infections in their blowing–up evolution phase along time. However, it is obvious that the epidemic models have the concourse of several coupled subpopulations so that it the model is reduced to a first-order dynamics the influence of the remaining dynamics should be accounted for through a time-varying parameterization and dynamics uncertainty in (1) since the model order is reduced to unity. The accuracy of the modeling procedure is evaluated by means of the entropy through (37). Hence if the actual infectious population curve is close to the reference one, then we have$S_I\left(\eta \right)=S_{I_r}\left(\eta \right)$which generates the dissipation rate of the model. On the other hand, if the current system differs from the reference model, then the entropy becomes corrected with the additional term${\tilde{S}}_I\left(\eta \right)$. Therefore, the contributing terms in (37) provide an estimation of the modeling uncertainty based on the assumed log- normal reference distribution. As a result, the best approximation of the current model to the reference one is that which minimizes the error entropy${\tilde{S}}_I\left(\eta \right)$, i.e., the one which reduces as much as possible the uncertainty introduced by the approximation.

Remark 5.

Note that the entropy of the infection$I\left(t\right)$for$\eta =1$is defined as$S_I\left(1\right)=-{\displaystyle {\int }_0^{\infty }I}\left(\tau \right)lnI\left(\tau \right)d\tau$The entropy of the truncated function$I_t\left(\tau \right)=I\left(\tau \right)$for$\tau \in \left[0,t\right]$and$I_t\left(\tau \right)=0$for$\tau \notin \left[0,t\right]$is$S_{I_t}\left(1\right)=-{{\displaystyle {\int }_0^{\infty }I}}_t\left(\tau \right)ln{I}_t\left(\tau \right)d\tau =-{\displaystyle {\int }_0^{t}I\left(\tau \right)}lnI\left(\tau \right)d\tau$. Note also that${\dot{S}}_{I_t}\left(1\right)=-I\left(t\right)lnI\left(t\right)$and${\ddot{S}}_{I_t}\left(1\right)=-\dot{I}\left(t\right)\left(1+lnI\left(t\right)\right)=0$if$t=D$. That is, the inflection point of the truncated entropy occurs at the relative extreme values of$I\left(t\right)$. In particular, if the infection is in its first expanding phase, this occurs at its maximum$t=D$.

4.2. Estimation of Errors of the Distribution Widths between the Log-Normal Reference and Current Model Information Entropies

One gets from (38) for the usual reference entropy definition based on the log-normal distribution of width ${\sigma }_r=\frac{1}{\sqrt{2\eta }}$, [11,33,37], that:

$$S_{I_r}\left({\sigma }_r,D_r,\eta \right)=\eta \left(ln\left(\sqrt{2\pi }{\sigma }_r\right)+\eta \left(ln{D}_r+{\sigma }_r^2\right)+\frac{1}{2}\right)$$

(45)

and the particular value:

$${\sigma }_r={\sigma }_r\left(\eta \right)=\arg \left(r\in {\mathit{R}}_{0+}\hspace{0.17em}:\hspace{0.17em}\frac{d^{2}S\left({\sigma }_r,D_r,\eta \right)}{d{r}^2}=0\right)$$

(46)

is the width distribution maximum value which makes the reference entropy to cease to increase while giving the maximum dissipation rate which leads to:

$$S_{I_r}\left(\sqrt{\frac{1}{2\eta }},\hspace{0.17em}{D}_r,\eta \right)=\eta \left(ln\sqrt{\frac{\pi }{\eta }}+\eta \left(ln{D}_r+{\sigma }_r^2\right)+\frac{1}{2}\right)=\arg \hspace{0.17em}\left(S_{I_r}\left({\sigma }_r,\hspace{0.17em}\eta \right)\hspace{0.17em}:\hspace{0.17em}\frac{d^{2}S\left({\sigma }_r,\eta \right)}{d{\sigma }_r^2}=0\right)$$

(47)

Note that the above reference description is easily associated to an epidemic model given by a first-order differential equation involving only the infection evolution. Note also, in particular, that the infection curve solution is of exponential order as it is the log-normal function. Such an order is negative if the disease-free equilibrium point is globally asymptotically stable (that is, the reproduction number is less than one) so that the infection converges exponentially to zero. In other words, the curves (43) and (44) can be reasonably identified with each other as it has been made in the above subsection by considering the influence of the initial conditions. In more sophisticated models involving the concourse of more subpopulations (say susceptible, immune, etc.), like those discussed in the above section, the differential equation is of higher-order than one so that the $\alpha \left(t\right)$ -function describing the time evolution of $I\left(t\right)$ depends on the remaining subpopulations. This translates into the following facts:

(1)

Fact 1: It is known that, for $\eta =3$, ${\sigma }_r=\sqrt{\frac{1}{6}}\approx 0.408$; $\frac{D_r}{L_r}=1.649$; $\frac{I\left(D_r\right)}{I\left(L_r\right)}=\frac{f\left(D_r\right)}{f\left(L_r\right)}=2.120$, [11].

(2)

Fact 2: A modification of the relevant time instants $D$ and $L$ of maximum infection and previous inflection point with respect to $D_r$ and $L_r$, and the corresponding entropies as it has been discussed analytically in Section 4.1. Those parameters depend on each particular model. This also will translate, as a result, into a change of the distribution width $\sigma$ related to the reference width ${\sigma }_r$ for the maximum dissipation concerns.

(3)

Fact 3: Although the above reference values ${\sigma }_r$ and ratio $D_r/L_r$ are independent on $I_0$, since they are got from the log-normal distribution function, the current ones are, in general, dependent on $I_0$. The entropy of the current given multi-subpopulation model is given explicitly by (37), subject to (38)–(39). The time instants $D$ and $L$ of the respective maximum and inflection infection time instants and their values $I\left(D\right)$ and $I\left(L\right)$ are calculated from the first zeros of the curves $\dot{I}\left(t\right)$ and $\ddot{I}\left(t\right)$, respectively which also lead directly to their corresponding rates.

(4)

Fact 4: The entropy of the current model might be interpreted in terms of the maximum dissipation rate by assuming a description via a log-normal distribution. However, it is easy to verify that the log-normal function is zeroed as its argument is either zero or $+\infty$, although its profile is close, but not identical, to the solution of a first-order differential equation describing a decaying exponential infection evolution towards a disease- free equilibrium point. For this reason, and having in mind the comparison of the solution of models with more than one subpopulation (with associated differential system of order larger than one) to the log-normal distribution $f\left(t\right)$ which is zero at zero and at infinity and which satisfies ${\displaystyle {\int }_0^{\infty }f\left(t\right)}d\left(t\right)=1$, we first normalize the infectious subpopulation of the current model in order to get a comparable entropy to the reference one associated with the log-normal function, that is, we define:

$$I_n\left(t\right)=\frac{I\left(t\right)}{{\displaystyle {\int }_0^{\infty }I\left(\tau \right)d\tau }}; S_{I_n}\left(\eta \right)=-{\displaystyle {\int }_0^{\infty }{I}_n\left(\tau \right)ln\hspace{0.17em}\left(I_n\left(\tau \right)/t^{\eta -1}\right)d\tau }$$

(48)

4.3. Some Numerical Tests on Reference and Current Model Entropies

Now Example 2 and Example 4 are compared to the infection study of [11], by introducing the appropriate tools of normalized infection entropy (48) associated with the maximum dissipation rate for the choice $\eta =1$. Recall the basic notation $D_r$, $L_r$, $D$ and $L$ being the first time instants such that $\dot{I}\left(D_r\right)=0$, $\ddot{I}\left(L_r\right)=0$, $\dot{I}\left(D\right)=0$, $\ddot{I}\left(L\right)=0$ (Examples 2 and 4). One gets from (45) for $\eta =1$ and, correspondingly, ${\sigma }_r=\sqrt{1/2}$ that the parameterized reference entropy is:

$$S_{I_r}\left(\sqrt{\frac{1}{2}},D_r,1\right)=\left(ln\sqrt{\pi }+ln{D}_r+1\right)$$

(49)

and one gets for Example 2 that its associated normalized entropy for $\eta =1$ being un-parameterized in $\left(D,\hspace{0.17em}\sigma \right)$ becomes from (48):

$$S_{I_n}\left(1\right)=-{\displaystyle {\int }_0^{\infty }{I}_n\left(t\right)ln}{I}_n\left(t\right)dt$$

(50)

Numerical experimentation with Example 2: Note that $D$ is the first time instant such that $\dot{I}\left(D\right)=0$ and $I\left(D\right)$ is a relative maximum, which in practice, gives the maximum expected infectious numbers. Also, $L$ is the first time instant such that $\ddot{I}\left(L\right)=0$. Note also that the basic model, of response being close to a log- normal function, has only an infectious subpopulation while the examples of Section 3 have more subpopulations integrated in the models. Therefore, the reasonable condition that the initial conditions of the infectious subpopulation are the one percent of the total population, we consider a total population of $N_0=S_0+I_0+R_0=1$ for Example 2 in order to get a feasible comparison.

Thus, we perform several alternative experiments as follows:

(a)

We get the values of the time instants $D$ and $L$ and the corresponding infection numbers $I\left(D\right)$ and $I\left(L\right)$, from the solution trajectory of Example 2 and its first two-time derivatives trajectories through time, as well as the normalized entropy $S_{I_n}\left(1\right)$ from (50). Later on, by equalizing (50) to (49), one then gets the value of $D_{rm}$ which specifies the time instant given a maximum infectious subpopulation with a maximum dissipation rate in a log normal distribution. This equalization yields:

$$D_{rm}=e^{S_{I_n}\left(1\right)-1}/\sqrt{\pi }$$

(51)

(b)

We equalize again (49) by fixing $D_r=D$ in (50). Then, we get the necessary value ${\sigma }_{rm}$ for such an equality to hold.

(c)

We define the variance with distribution function $I_n\left(t\right)$ and log- normal distribution resulting to be:

$$var\left(\omega \right)={\displaystyle {\int }_0^{\infty }\omega }\left(t\right){\left(t-{\displaystyle {\int }_0^{\infty }\omega \left(t\right)tdt}\right)}^{2}dt$$

(52)

where $\omega \left(t\right)=I_n\left(t\right)$ or $\omega \left(t\right)=x\left(t,D_r,{\sigma }_r\right)$, the log-normal distribution. Then, we obtain the necessary ${\sigma }_{rmv}={\sigma }_{rmv}\left(var\left(I_n\right),\hspace{0.17em}D\right)$ got from

$$var\left(I_n\right)=var\left(x\left(D_r\hspace{0.17em}=D,{\sigma }_{rmv}\right)\right)$$

(53)

One observes that, in general, ${\sigma }_{rmv}\ne {\sigma }_r=\frac{1}{\sqrt{2}}$ which ensures that the variance of log- normal distribution is equal to $var\left(I_n\right)$ for such a value of ${\sigma }_{rmv}$. Some numerical data on Example 2 are now compared with the log-normal distribution function. The model parameters are $\beta =$13,065 and γ = 50.1 year⁻¹ what means that the average infectious period is T_γ= 1/γ = 365/50.1 = 7.29 days. The initial infectious subpopulation is the one percent of the normalized total one N₀ = 1. For those initialization, the quotient $S_0/I_0$ (percentage of initial susceptible subpopulation versus recovered subpopulation) is used to plot Figure 5, Figure 6, Figure 7 and Figure 8 whose time scale are in years. Figure 5 displays the time instants of maximum infection and inflection point versus different values of $S_0/R_0$. The values of $D_{rm}$ from (51) is also plotted. The corresponding infectious subpopulations are displayed in Figure 6. Figure 7 gives the entropies of (50) and (49). On the other hand, Figure 8 displays ${\sigma }_{rm}$, ${\sigma }_{rmv}$ and the variance of the normalized infectious $I_n\left(t\right)$ of (52). It is basically concluded that for the model of example 2 which has three subpopulations, the results are distinct from to those obtained from the log-normal distribution which we can recall that behave closely to the solution of a first-order differential equation involving the infectious only for initial infection being close to zero and small susceptible amounts. The above discrepancy increases as the quotient $S_0/I_0$ increases. The reason of the approximation discrepancy is that the couplings of the infectious subpopulations with the remaining ones becomes increasingly relevant to the transient responses evolution as the proportion of susceptible to infectious increases.

Numerical experimentation with Example 4: The initial values satisfy a normalization constraint $N_0=S_0+R_0=1$ with subpopulations S₀ = 0.99, i_o = 0.01 (that, is the initial infectious subpopulation is $1\%$ of the total one) and $R_0=0$ since the recovered populations is compensatory in the model in order to take into account the effects of the intervention controls. The parameters $\beta$ and $\gamma$ are fixed as in Example 2. In particular, Figure 9 and Figure 10 show the maximum infection and its previous value at the inflection time instant and the corresponding time instants without vaccination and with a vaccination effort rate of k_T = 290 for different values of the vaccination control gain. It is basically seen that the maximum and inflection amounts decrease as the treatment control gain gives a skip from zero to an important effort as that, in parallel, the above values also decrease as the vaccination control gain increases. Figure 10 and Figure 11 describe parallel experiments where the roles of the vaccination and treatment control gains are reversed with respect to the data of Figure 9 and Figure 10. The obtained conclusions are similar. The time instants of maximum infection and the inflection value are reached without and with vaccination control as the treatment control effort increases for Example 4 are plotted in Figure 12. The corresponding entropies for those to experiments compared to the reference entropy are displayed in Figure 13 and Figure 14. Note that the entropies (48) and (50) reach negative values because of the normalization of the infection by the total infection integral contribution (48) used to evaluate the normalized entropy (50). Note that the vaccination control does not affect to the entropy as significantly as the treatment control gains since it influences less significantly to the model dynamics.

5. Conclusions

This paper has investigated the extensions of a first-order differential system describing the infection propagation through time to epidemic models integrating more than one subpopulation. The main involved tool has been the consideration of the coupling of inter-populations dynamics and the control intervention information through the structure of the time-varying coefficient which drives the basic differential equation model of first-order. The control of the infection along its transient to fight more efficiently against a potential initial exploding transmission from a high initial growth rate is considered relevant. Special attention has been paid throughout the manuscript to the discussion of the profiles of the transients of the infection curve in terms of the time instants of its first relative maximum towards its previous inflection time instant, so the study is mainly focused on the transient behavior characterization rather than on the steady-state equilibrium points. The time instants leading to the maximum infection and inflection numbers have been investigated via the Shannon´s information entropy for the maximum dissipation rate linked to a previous background study for a first-order differential equation describing the infection propagation. Since it is relevant to know the time instants of maximum infection and inflection as well as its numbers in order to monitor the availability of hospitalization resources, some examples related to existing epidemic models integrated by more than a subpopulation have been studied. The obtained results have been compared, both via theoretical work and also by numerical experimentation, to the background results obtained from a reference model, just involving a single infectious population, which is based on a description via a log-normal distribution which has a close profile to the solution response of a first-order differential equation. In those examples, special attention is paid to the comparisons of the maximum infection and inflection time dates for different values of initial conditions and to the entropy discrepancies related to the reference one. It can be concluded that the influence of the couplings of the dynamics of other subpopulations in the model to the infectious one is relevant to the infection evolution, especially, in the cases when the initial amounts of the susceptible are significantly large compared to the initial amounts of the infectious.