Our objective in this section is to extend the initial model to include two intervention methods, called controls, represented as functions of time and assigned reasonable upper and lower bounds, each representing a possible method of influenza intervention. Using optimal control theory and numerical simulations, we determine the benefit of vaccination and media coverage when the latter has positive or negative effect on the former.

We will integrate the essential components into one SIVR-type model to accommodate the dynamics of an influenza outbreak determined by population-specific parameters such as the effect of contact reduction when infectious and vaccinated individuals are reported in the media.

Let

*u*_{v} and

*u*_{m} be the control variables for vaccination and media coverage respectively. Thus, model (1)-(4) now reads

A balance of multiple intervention methods can differ between populations. A successful mitigation scheme is one which reduces influenza-related deaths with minimal cost. A control scheme is assumed to be optimal if it maximizes the objective functional

The first two terms represent the benefit of the susceptible and vaccinated populations. The parameters

*B*_{1} and

*B*_{2} represent the weight constraints for the infected population and the control, respectively. They can also represent balancing coefficients transforming the integral into dollars expended over a finite time period of

*T* days [

34]. The goal is to maximize the populations of susceptible and vaccinated individuals, minimize the population of infectives, maximize the benefits of media coverage and vaccination, while minimizing the systemic costs of both media coverage and vaccination. The value

*u*_{v}(

*t*) =

*u*_{m}(

*t*) = 1 represents the maximal control due to vaccination and media coverage, respectively. The terms

and

represent the maximal cost of education, implementation and campaigns on both vaccination and media coverage.

*S*(

*t*) and

*V*(

*t*) account for the fitness of the susceptible and the vaccinated groups as a result of a reduction in the rate at which the vaccine wanes, and vaccination and treatment efforts are implemented [

35]. We thus seek optimal controls

and

such that

where

*U* = {(

*u*_{v},

*u*_{m})|

*u*_{v},

*u*_{m} measurable, 0 ≤

*a*_{11} ≤

*u*_{v} ≤

*b*_{11} ≤ 1, 0 ≤

*a*_{22} ≤

*u*_{m} ≤

*b*_{22} ≤ 1,

*t* [0,

*t*_{f}]} is the control set, with

*t* [

*t*_{0},

*t*_{f}]. The basic framework of this problem is to characterize the optimal control.

Existence of an optimal control

The existence of an optimal control can be obtained by using a result by Joshi [

36] and Fister

*et al.*[

37].

**Theorem 3***Consider the control problem with the system of Equations (4.1)-(4.4). There exists an optimal control**such that max***Proof.** To prove this theorem on the existence of an optimal control, we use a result from Fleming and Rishel [

38] (Theorem 4.1 pp. 68-69), where the following properties must be satisfied.

1. The set of controls and corresponding state variables is nonempty.

2. The control set *U* is closed and convex.

3. The right-hand side of the state system is bounded above by a linear function in the state and control.

4. The integrand of the functional is concave on *U* and is bounded above by *c*_{2} – *c*_{1}(*|u*_{v}|^{k} + *|u*_{m}|^{k}), where *c*_{1}, *c*_{2} > 0 and *k* > 1.

An existence result in Lukes [

39] (Theorem 9.2.1) for the system of equations (6)-(9) for bounded coefficients is used to give the first condition. The control set is closed and convex by definition. The right-hand side of the state system (Equations (4.1)-(4.4)) satisfies Condition 3 since the state solutions are a priori bounded. The integrand in the objective functional,

, is concave on

*U*. Furthermore,

*c*_{1},

*c*_{2} > 0 and

*k* > 1, so

Therefore, the optimal control exists, since the left-hand side of (11) is bounded; consequently, the states are bounded.

Since there exists an optimal control for maximizing the functional (10) subject to equations (6)-(9), we use Pontryagin’s Maximum Principle to derive the necessary conditions for this optimal control. Pontryagin’s Maximum Principle introduces adjoint functions that allow us to attach our state system (of differential equations), to our objective functional. After first showing existence of optimal controls, this principle can be used to obtain the differential equations for the adjoint variables, corresponding boundary conditions and the characterization of an optimal control

and

. This characterization gives a representation of an optimal control in terms of the state and adjoint functions. Also, this principle converts the problem of minimizing the objective functional subject to the state system into minimizing either the Lagrangian or the Hamiltonian with respect to the controls (bounded measurable functions) at each time

*t*[

40]. The Lagrangian is defined as

where

*w*_{11}(

*t*) ≥ 0,

*w*_{12}(

*t*) ≥ 0 are penalty multipliers satisfying

*w*_{11}(

*t*)(

*a*_{11} –

*u*_{v}(

*t*)) +

*w*_{12}(

*t*)(

*u*_{v}(

*t*) –

*b*_{11}) at the optimal

, and

*w*_{21}(

*t*) ≥ 0,

*w*_{22}(

*t*) ≥ 0 are penalty multipliers satisfying

*w*_{21}(

*t*)(

*a*_{22} –

*u*_{m}(

*t*)) +

*w*_{22}(

*t*)(

*u*_{m}(

*t*) –

*b*_{22}) at the optimal

.

Given optimal controls

and

, and solutions of the corresponding state system (6)-(9),
there exist adjoint variables λ

_{i}, for i = 1, 2, 3, 4 satisfying the following equations

with transversality conditions

*λ*_{i}[

*t*_{f}] = 0, for

*i* = 1

*,* 2

*,* 3

*,* 4. To determine the interior maximum of our Lagrangian, we take the partial derivatives of L with respect to

*u*_{v} and

*u*_{m}, respectively, and set it to zero. Thus,

To determine an explicit expression for our controls

,

(without

*w*_{11}*, w*_{12},

*w*_{21},

*w*_{22}), a standard optimality technique is utilized. The following cases are considered to determine a specific characterization of the optimal control.

**Case 1:** Optimality of

1. On the set

. Hence, the optimal control is

2. On the set

.We have

since *w*_{12} ≥ 0.

3. On the set

. Hence

Combining all the three sub-cases in a compact form gives

**Case 2:** Optimality of

1. On the set

. We have

2. On the set

. We have

since *w*_{22} ≥ 0.

3. On the set

. Hence

Combining all the three sub-cases in a compact form gives

The optimal system

The optimality system consists of the state system coupled with the adjoint system, with the initial conditions, the transversality conditions and the characterization of the optimal control:

where

and

are given by expressions (12) and (13), respectively, with

*S*(0) =

*S*_{0},

*I*(0) =

*I*_{0},

*V*(0) =

*V*_{0},

*R*(0) =

*R*_{0} and

*λ*_{i}[

*t*_{f}] = 0 for

*i* = 1

*,··· ,*4. Due to the a priori boundedness of the state and adjoint functions and the resulting Lipschitz structure of the ODEs, we obtain the uniqueness of the optimal control for small [

*t*_{f}] [

36]. The uniqueness of the optimal control follows from the uniqueness of the optimality system.

The state system of differential equations and the adjoint system of differential equations together with the control characterization above form the optimality system solved numerically and depicted in Figures , , , .

The model with pulse vaccination

The general model with pulse vaccination is given as

for *t* ≠ *t*_{k}, where *t*_{k} is the time of the *k*th vaccination. We may have *t*_{k+1} – *t*_{k} either constant or not, as we choose. The impulsive effect is given by

Δ*S* = –*θS*

Δ*V* = *θS*

when

*t* =

*t*_{k}. Here,

is the change in state at the impulse time.

In this model, vaccination occurs at fixed times, not continuously. This is closer to reality, since vaccination centres are only open at certain times, when people may get vaccinated in waves. Similarly, media stories tend to clump together, so that a big news story occurs on one day, which may trigger a short period of intense vaccination. We shall use a simplified version of this model to illustrate the possibility that media may have an adverse effect.

Adverse effects

Consider the following scenario. At the onset of the outbreak, the media - and hence the general population - is unaware of the disease and thus nobody gets the vaccine, allowing the disease to spread in its initial stages. At some point, there is a critical number of infected individuals, whereupon people are sufficiently aware of the infection to change their behaviour. We suppose that, initially, new infected people arrive at fixed times.

We further assume that vaccinated people mix more than susceptibles. In this case, people who are vaccinated feel confident enough to mix with the infected, even though they may still have the possibility to contract the virus. This might be the case for health-care workers, for instance, who get vaccinated and then have to tend to the sick.

Mathematically, we have a threshold for the critical number of infectives, *I*_{crit}*.*

For

*I* <

*I*_{crit}, this model would look like

For

*I* >

*I*_{crit}, the model becomes

with *β*_{4} – *β*_{6} ≥ 0.

However, to illustrate the adverse affect, we shall simplify the model even further. For a short timescale, we can assume recovery is permanent, so *σ* = 0. Thus, we can ignore the *R* equation.

For

*I* <

*I*_{crit}, we assume that there is no mixing, but rather that new infectives arrive impulsively into the system at fixed times

*t*_{k} and in numbers

*I*^{i}, where

*I*^{i} *I*_{crit}. (If the new infectives arrive at irregular times, then the broad results will be unchanged.)

For *I* >*I*_{crit}, fear of the disease keeps susceptibles from mixing with the infected, but the vaccinated will.

Thus

*β*_{4} =

*β*_{6} = 0. Since

*I*^{i} *I*_{crit}, we can assume that, for

*I* >

*I*_{crit}, the effects of new infectives are negligible.

for *I* >*I*_{crit}.

Thus, the effects of the media are to trigger a vaccinating panic whenever the number of infectives is large enough. We kept the model with impulse vaccination as simple as possible since even this simplified version shows that media reports could have an adverse effect.

Suppose new infectives appear regularly, so that

*t*_{k}_{+1} –

*t*_{k} =

*τ*. (If not, the analysis generalizes quite easily.) For

*t*_{k} <

*t* <

*t*_{k}_{+1}, we have

where

is the value immediately after the

*k*th impulse. Then, since the period is constant, we have

This is a recursion relation with solution

Thus, if

*m*^{+} >

*I*_{crit}, then eventually the system will switch from model (14)-(17) to model (18)-(20). The endemic equilibrium in model (18)-(20) satisfies

At the endemic equilibrium,

. Thus, we have

The characteristic equation is

It follows that the endemic equilibrium is stable if *Î* >*I*_{crit}. Thus, even in an extremely simplified version of the model, the media may make things significantly worse than if no media effect were included. We kept this model deliberately simple, partly for mathematical tractability and partly to show that the media effects apply even in this idealised scenario.

Note that, in reality, the fluctuations would apply in the upper region as well, making the actual value even

larger. In the lower region, we ignored interaction between susceptibles and infectives (ie we assume *β*_{4} = *β*_{6} = 0). The effect of including these terms would be to slow the exponential decay between impulses (or possibly cause it to increase). This would only increase the effect seen here.

In summary, a small series of outbreaks that would equilibrate at some maximum level *m*^{+} >*I*_{crit} will, as a result of the media, instead equilibrate at a much larger value *I* >*m*^{+} >*I*_{crit}. The driving factor here is if an imperfect vaccine causes overconfidence, so that people who have been vaccinated mix significantly more with infectives than susceptibles do. If this happens (as would be quite likely; most people who have been vaccinated feel invulnerable, even if the vaccine is not perfect, largely thanks to media oversimplifications), then the media effect is likely to be adverse. A simplified version of the model with pulse vaccination shows that the media can make things worse, if the vaccine is imperfect because the vaccinated mix over-confidently with the infectives.