# Modelling Alcoholism as a Contagious Disease: A Mathematical Model with Awareness Programs and Time Delay

Authors: Shuang-Hong Ma; Hai-Feng Huo; Xin-You Meng
Journal: Discrete Dynamics in Nature and Society
 (2015)

---

## Abstract

A dynamic alcohol consumption model with awareness programs and time delay is formulated and analyzed. The aim of this model is to capture the effects of awareness programs and time delay in controlling the alcohol problems. We introduce awareness programs by media in the model as a separate class with growth rate of the cumulative density of them being proportional to the number of mortalities induced by heavy drinking. Susceptible population will isolate themselves and avoid contact with the heavy drinkers or become aware of risk of heavy drinking and decline to drink due to such programs. In particular, we incorporate time delay because the nonconsumer population will take a period of time to become an alcohol consumer. We find that the model has two equilibria: one without alcohol problems and one where alcohol problems are endemic in population. The model analysis shows that though awareness programs cannot eradicate alcohol problems, they are effective measures in controlling the alcohol problems. Further, we conclude that the time delay in alcohol consumption habit which develops in susceptible population may result in Hopf bifurcation by increasing the value of time delay. Some numerical simulation results are also given to support our theoretical predictions.


---

## Body

## 1. Introduction

Alcoholism, also known as alcohol dependence, is a function of social epidemic, environmental contexts, individuals’ preferences, and family history. Alcohol consumption has been identified as a major contributor to the global burden of chronic disease, injury, and economic cost [1–3]. Alcohol abuse can also lead to a range of negative social effects such as violence and antisocial and criminal behavior. The World Health Organization reports that the harmful use of alcohol causes approximately 3.3 million deaths every year (or 5.9% of all the global deaths), and 5.1% of the global burden of disease is attributable to alcohol consumption [4]. Therefore, alcoholism is among the main targeted health risk behaviors due to the high relevance of negative health and social effects. The study of alcoholism has become an important aspect of social epidemic.The spread of health risk behavior within a community can be viewed as a diffusion process with its own incidence rate. In this situation, the social interaction is considered to be the key factor in spreading the behavior which can result in adverse health effects. For this reason, alcoholism can be viewed as a treatable contagious disease. Mathematical model is a predictive tool which can mimic the process of infectious diseases and provide useful measures to analyze the spread and control of infectious diseases [5, 6]. Furthermore, mathematical models for human behaviors related to addictions have been developed from epidemiological models for the spread of infectious diseases, including drinking, smoking, and drug use (see [7–10] and the references contained therein). In particular, numerous mathematical models for alcohol problems have been formulated and studied recently [11–21]. Benedict [11] modelled alcoholism as a contagious disease and studied how “infected” drinking buddies spread problem drinking. Manthey et al. [12] studied campus drinking and suggested that the reproductive numbers are not sufficient to predict whether drinking behavior will persist on campus and that the pattern of recruiting new members plays a significant role in the reduction of campus alcohol problems. Environmental and peer influences combined to create a culture of drinking were studied in [13–16]. In addition, the two-stage models (one stage where people admit having alcohol problem and another stage where people do not admit having alcohol problem) have been developed in [17, 18]. Bhunu and Mushayabasa [19] studied the cointeraction of alcoholism and smoking in a community. Walters et al. [20] also discussed alcohol problems, and their results showed that an increase in the recovery rate decreased the proportion of binge drinkers in the population. Wang et al. [21] presented a deterministic SATQ-type mathematical model (the model including susceptible, alcoholism, treating, and quitting drinking compartments) for the spread of alcoholism with two control strategies to gain insights into this increasingly concerned-about health and social phenomenon. The optimal control strategies are derived by proposing an objective functional and using Pontryagin’s Maximum Principle.A rational step is to make people aware of the alcohol problems through the media. Media (e.g., Radio, Newspapers, Billboards, TV, and Internet), being the prime source of information, can not only influence the individuals’ behavior but also increase the governmental health care involvement to control the spread of heavy drinking. These behavioral responses can change the transmission patterns and declination to drink. In view of this, there is a need to incorporate the effect of awareness programs through the media in the mathematical models. In recent years, many mathematical models have been used for studying the impact of awareness programs by media on epidemic outbreaks (see [22–32] and the references cited therein). These studies suggested that education and media have a huge impact on controlling the spread of infectious diseases. However, most of these studies assumed that the media reduced the contact rate of susceptible with infective individuals. In addition, these studies assumed that the cumulative density of awareness programs is constant in models. Hence, it is more reasonable to consider that the cumulative density of awareness programs varies. Misra et al. [33] have explicitly incorporated the cumulative density of awareness programs in the modelling process, which is a separate dynamic variable whose growth depends on the size of epidemic. The results showed that awareness programs by media are implemented in accordance with the size of infective population. Considering the effect of susceptible population influenced by media to be limited, Misra et al. [34] have considered interaction between susceptible population and awareness programs as Holling-II functional response function.Recently, Huo and Wang [35] have studied drinking dynamics and focused on awareness programs and treatment in the modelling process. They have extended the model in [33] via including a treatment class and established some sufficient conditions for the stability of the alcohol-free and alcohol-present equilibria.It is well known that delay differential equations exhibit much more complicated dynamics than ordinary differential equations since the time delay may affect the stability of the system and even lead to instability, oscillation, or bifurcation phenomena [34, 36–39]. In most of the above studies, with susceptible individual contacts with heavy drinkers, he/she will become a heavy drinker at once. However, the alcoholism epidemic should consider the effects of the time delay. Since a nondrinker or moderate drinker is infected by heavy drinkers, there is a time delay during which a susceptible individual becomes heavy drinker. Hence, it is more realistic to consider the time delay in the modelling alcoholism process.In this paper, motivated by the above works, we present a nonlinear alcoholism model with awareness programs and time delay. First, we assume that the growth rate of the cumulative density of awareness programs is proportional to the number of mortalities induced by heavy drinking. Hence, the awareness about drinking will alert the susceptible individuals so that they isolate themselves and decline to drink or drink moderately (small intake of alcohol may be beneficial to health). Second, we assume that the heavy drinkers can recover from heavy alcohol drinking due to counselling, health reasons, treatment, prohibition, tax hiked on alcohol beverages, and so forth. A fraction of recovered population will join the unaware susceptible population whereas the rest will join the aware population. In particular, we incorporate time delayτ to describe the time needed for a susceptible individual to become a heavy alcohol user. The goal is to analyze the impact of awareness programs and time delay on the alcohol drinking distribution. We also analyze the effect of time delay on system stability. Our results show that awareness programs are an effective measure in controlling heavy drinking. Our results also show that the system will change its stability with the increase of the value of time delay.The remainder of this paper is organized as follows. In Section2, we formulate the delayed mathematical model. The condition for the existence of the equilibria is shown in Section 3. In Section 4, we discuss the stability of equilibria and the existence of Hopf bifurcation. To support our theoretical predictions, some numerical simulations are included in Section 5. A brief conclusion is also given in Section 6.


## 2. Model Formulation

The total population in the model isP
(
t
) at time t, and it is divided into three types: S
(
t
), X
(
t
), and A
(
t
). S
(
t
) denotes susceptible individuals who do not drink or drink only moderately. A
(
t
) is referred to as heavy drinkers. The class X
(
t
) is aware population who are aware of risk and decline to drink or drink moderately. It is supposed that the cumulative density of awareness programs driven by media at time t is M
(
t
). The growth rate of M
(
t
) is proportional to the number of mortalities induced by heavy drinking. It is also assumed that heavy drinking only spreads due to the direct contact between susceptible individuals and heavy drinkers. Susceptible population avoid being in contact with the heavy drinkers because the implementation of awareness programs results in awareness about heavy drinking. In view of this, a separate class X
(
t
) is formed. It is assumed that when a nondrinker or moderate drinker is “infected” by heavy drinkers, there is a time delay τ during which the alcohol consumption habit develops in his/her body. It is only after such time delay τ that a susceptible individual becomes a heavy drinker.Considering homogenous mixing, the transfer diagram is shown in Figure1.Figure 1
The transfer diagram of system (1).The transfer diagram leads to the following dynamic alcohol consumption model with awareness programs and time delay:(1)
d
S
t
d
t
=
Λ
-
β
S
t
-
τ
A
t
-
τ
-
λ
S
t
M
t
k
0
+
M
t
+
σ
X
t
+
1
-
p
γ
A
t
-
d
S
t
,
d
A
t
d
t
=
β
S
t
-
τ
A
t
-
τ
-
γ
+
d
+
α
A
t
,
d
X
t
d
t
=
λ
S
t
M
t
k
0
+
M
t
+
p
γ
A
t
-
d
+
σ
X
t
,
d
M
t
d
t
=
k
α
A
t
-
μ
M
t
,where Λ is the recruitment rate of susceptible individuals who do not drink or drink only moderately, β is the transmission coefficient of the infection for the susceptible individuals from the heavy drinkers, λ is the dissemination rate of awareness programs among susceptible population, k
0 limits the effect of awareness programs on susceptible population, known as half saturation point for interactions of Holling-II functional response function, σ is the transmission rate of transferring from aware population to unaware susceptible population due to fading of memory or certain factors, and so forth, γ is the recovery rate of heavy drinkers, p represents fraction of recovered individuals joining aware class, 1
-
p represents fraction of recovered individuals joining unaware susceptible class, d is the natural death rate of the general population in the model, α is the death rate due to heavy alcohol consumption, k is the implementation rate of awareness programs, μ is the depletion rate of awareness programs due to ineffectiveness such as social and psychological factors in population, and τ is the time lag due to development of alcohol consumption habit.All the above parameters are assumed to be positive constants. We suppose that the initial condition for system (1) takes the form S
(
θ
)
,
A
(
θ
)
,
X
(
θ
)
≥
0, θ
∈
[
-
τ
,
0
], and S
(
0
)
,
A
(
0
)
,
X
(
0
)
>
0, where S
(
θ
)
,
A
(
θ
)
,
X
(
θ
)
∈
C
(
[
-
τ
,
0
]
,
R
+
), R
+
=
x
∣
x
≥
0, and M
(
0
)
=
M
0
≥
0.Using the fact thatP
(
t
)
=
S
(
t
)
+
A
(
t
)
+
X
(
t
), system (1) is rewritten as(2a)
d
A
t
d
t
=
β
P
t
-
τ
-
A
t
-
τ
-
X
t
-
τ
A
t
-
τ
-
γ
+
d
+
α
A
t
,
(2b)
d
X
t
d
t
=
λ
P
t
-
A
t
-
X
t
M
t
k
0
+
M
t
+
p
γ
A
t
-
d
+
σ
X
t
,
(2c)
d
P
t
d
t
=
Λ
-
d
P
t
-
α
A
t
,
(2d)
d
M
t
d
t
=
k
α
A
t
-
μ
M
t
.Since system (1) is equivalent to system (2a), (2b), (2c), and (2d), we study system (2a), (2b), (2c), and (2d) in this paper.From (2c), we can obtain (3)
l
i
m
t
→
∞
sup
⁡
P
t
≤
Λ
d
; then (4)
0
≤
A
,
X
≤
P
≤
Λ
d
.From (2d), we have (5)
l
i
m
t
→
∞
sup
⁡
M
t
≤
k
α
Λ
μ
d
.Hence, the feasible region of system (2a), (2b), (2c), and (2d) is given by the set (6)
Ω
=
A
,
X
,
P
,
M
∈
R
+
4
:
0
≤
A
,
X
≤
P
≤
Λ
d
,
0
≤
M
≤
k
α
Λ
μ
d
,which is positively invariant with respect to system (2a), (2b), (2c), and (2d). We will consider dynamic behavior of system (2a), (2b), (2c), and (2d) on the set Ω.


## 3. Existence of Equilibria

It is obvious that system (2a), (2b), (2c), and (2d) has one alcohol-free equilibrium E
0
(
0,0
,
Λ
/
d
,
0
) and one alcohol-present equilibrium E
1
(
A
∗
,
X
∗
,
P
∗
,
M
∗
). The existence of equilibrium E
0 is trivial. In the following, we will prove the existence of alcohol-present equilibrium E
1.The basic reproduction number of the alcohol consumption model can be easily obtained by the next generation matrix method formulated in [40], which is given by the following expression: (7)
R
0
=
β
Λ
d
γ
+
d
+
α
,and it acts as a threshold as is shown in the following result.Theorem 1.
System (2a), (2b), (2c), and (2d) always has the alcohol-free equilibrium E
0
(
0,0
,
Λ
/
d
,
0
). Furthermore, system (2a), (2b), (2c), and (2d) has a unique alcohol-present equilibrium E
1
(
A
∗
,
X
∗
,
P
∗
,
M
∗
) if R
0
>
1.Proof.
For equilibriumE
1, A
∗, X
∗, P
∗, and M
∗ satisfy the following algebraic equations:(8a)
β
P
∗
-
A
∗
-
X
∗
A
∗
-
γ
+
d
+
α
A
∗
=
0
,
(8b)
λ
P
∗
-
A
∗
-
X
∗
M
∗
k
0
+
M
∗
+
p
γ
A
∗
-
d
+
σ
X
∗
=
0
,
(8c)
Λ
-
d
P
∗
-
α
A
∗
=
0
,
(8d)
k
α
A
∗
-
μ
M
∗
=
0
. From (8c), we obtain(9)
P
∗
=
Λ
-
α
A
∗
d
.The fourth equation (8d) yields(10)
M
∗
=
k
α
μ
A
∗
.Substituting (9) into (8a), we get(11)
X
∗
=
Λ
-
α
A
∗
-
d
A
∗
d
-
γ
+
d
+
α
β
.Then, substituting (9)–(11) into (8b), we have a quadratic equation of A
∗ as follows:(12)
a
1
A
∗
2
+
a
2
A
∗
+
a
3
=
0
,where(13)
a
1
=
α
k
p
γ
+
α
+
d
+
α
σ
d
+
σ
,
a
2
=
α
k
γ
+
d
+
α
λ
+
d
+
σ
β
+
μ
k
0
p
γ
+
α
+
d
+
α
σ
d
+
σ
-
α
k
Λ
1
+
σ
d
,
a
3
=
μ
k
0
d
+
σ
γ
+
d
+
α
β
-
Λ
d
=
μ
k
0
d
+
σ
γ
+
d
+
α
β
1
-
R
0
.
It is apparent that the existence of equilibriumE
1 is equivalent to the existence of A
∗. Obviously, a
1
>
0, and a
3
<
0 holds naturally when R
0
>
1. Hence, (12) has a unique positive solution, say A
∗. Using the value of A
∗, we can get positive solutions P
∗
,
M
∗, and X
∗, respectively. So system (2a), (2b), (2c), and (2d) has a unique alcohol-present equilibrium E
1 if R
0
>
1. This completes the proof.


## 4. Stability Analysis

In this section, we will investigate the stability of the equilibria of system (2a), (2b), (2c), and (2d).


### 4.1. Stability Analysis withτ
=
0

First, we investigate the stability under the condition ofτ
=
0. The stability of two equilibria E
0 and E
1 of system (2a), (2b), (2c), and (2d) can be given in the following theorem.Theorem 2.
Forτ
=
0, one has the following results:(i)
The alcohol-free equilibriumE
0 is globally asymptotically stable if R
0
<
1.
(ii)
The alcohol-free equilibriumE
0 is unstable and the alcohol-present equilibrium E
1 exists if R
0
>
1.
Furthermore, equilibriumE
1 is locally asymptotically stable if the coefficients of the characteristic equation of system (2a), (2b), (2c), and (2d) at E
1 satisfy the Routh-Hurwitz criterion.Proof.
Whenτ
=
0, the Jacobian matrix J of system (2a), (2b), (2c), and (2d) is (14)
J
=
β
P
-
2
A
-
X
-
γ
+
d
+
α
-
β
A
β
A
0
p
γ
-
λ
M
k
0
+
M
-
λ
M
k
0
+
M
-
d
+
σ
λ
M
k
0
+
M
λ
k
0
P
-
A
-
X
k
0
+
M
2
-
α
0
-
d
0
α
k
0
0
-
μ
.
(i) For the alcohol-free equilibriumE
0
(
0,0
,
Λ
/
d
,
0
), the matrix J at E
0 can be written as (15)
J
0
=
J
E
0
=
β
Λ
d
-
γ
+
d
+
α
0
0
0
p
γ
-
d
+
σ
0
λ
Λ
d
+
k
0
-
α
0
-
d
0
α
k
0
0
-
μ
.It is obvious that every eigenvalue of J
0 has negative real part whenever R
0
<
1. So the alcohol-free equilibrium E
0
(
0,0
,
Λ
/
d
,
0
) of system (2a), (2b), (2c), and (2d) is globally asymptotically stable if R
0
<
1.
(ii) IfR
0
>
1, then the alcohol-free equilibrium E
0 is a saddle point. So the equilibrium E
0 of system (2a), (2b), (2c), and (2d) is unstable. In this case, the alcohol-present equilibrium E
1 exists. Further, the matrix J at E
1 can be written as (16)
J
1
=
J
E
1
=
β
P
∗
-
2
A
∗
-
X
∗
-
γ
+
d
+
α
-
β
A
∗
β
A
∗
0
p
γ
-
λ
M
∗
k
0
+
M
∗
-
λ
M
∗
k
0
+
M
∗
-
d
+
σ
λ
M
∗
k
0
+
M
∗
λ
k
0
P
∗
-
A
∗
-
X
∗
k
0
+
M
∗
2
-
α
0
-
d
0
α
k
0
0
-
μ
=
J
11
J
12
J
13
0
J
21
J
22
J
23
J
24
J
31
0
J
33
0
J
41
0
0
J
44
.
The characteristic equation of system (2a), (2b), (2c), and (2d) at E
1 is d
e
t
(
φ
I
-
J
1
)
=
0, where I is the identity matrix of size 4 and φ is the eigenvalue of matrix J
1. By simple computation, the characteristic equation can be expressed: (17)
φ
4
+
b
1
φ
3
+
b
2
φ
2
+
b
3
φ
+
b
4
=
0
,where(18)
b
1
=
-
J
11
+
J
22
+
J
33
+
J
44
,
b
2
=
J
11
J
22
+
J
11
J
33
+
J
11
J
44
+
J
22
J
33
+
J
22
J
44
+
J
33
J
44
-
J
12
J
21
+
J
13
J
31
,
b
3
=
J
12
J
21
J
33
+
J
12
J
21
J
44
+
J
13
J
31
J
22
+
J
13
J
31
J
44
-
J
11
J
22
J
33
+
J
11
J
22
J
44
+
J
11
J
33
J
44
+
J
22
J
33
J
44
+
J
12
J
24
J
41
+
J
12
J
31
J
23
,
b
4
=
J
11
J
22
J
33
J
44
+
J
12
J
23
J
31
J
44
+
J
12
J
24
J
33
J
41
-
J
12
J
21
J
33
J
44
-
J
13
J
22
J
31
J
44
.
Notice thatR
0
>
1; it can be calculated that the coefficients b
i (i
=
1,2
,
3,4
.) all are positive according to J
i
i (i
=
1,2
,
3,4
.) and J
12
,
J
31 all are negative and J
13, J
21, J
23, J
24, and J
41 all are positive. Furthermore, by using the Routh-Hurwitz criterion [41], it is easy to see that the eigenvalues of the above characteristic equation have negative real parts if the conditions (19)
b
1
b
2
-
b
3
>
0
,
b
1
b
2
b
3
-
b
3
2
-
b
1
2
b
4
>
0hold. This shows that the alcohol-present equilibrium E
1 is locally asymptotically stable, and the proof is complete.Moreover, we know thatR
0 is a threshold value which determines whether alcohol problems will die out or become endemic in the population. The alcohol-free equilibrium E
0 undergoes a transcritical bifurcation at R
0
=
1. We are now interested in knowing which model parameter has the greatest effect on R
0 and hence has the greatest effect in controlling the alcohol consumption.Using the approach in Chitnis et al. [42], we calculate the normalized forward sensitivity indices of R
0. Let (20)
Υ
q
R
0
=
∂
R
0
∂
q
×
q
R
0denote the sensitivity index of R
0 with respect to the parameter q. We get(21)
Υ
β
R
0
=
Υ
Λ
R
0
=
+
1
,
Υ
d
R
0
=
-
1
+
d
γ
+
d
+
α
,
Υ
γ
R
0
=
-
γ
γ
+
d
+
α
,
Υ
α
R
0
=
-
α
γ
+
d
+
α
.Considering that the natural death rate d is a relatively stable constant, we find that Υ
β
R
0 and Υ
Λ
R
0 are the largest sensitivity indexes. We focus on reducing the value of β or Λ whose reduction by 1% would reduce R
0 by 1%. In contrast, there is a need to increase γ in order to control alcohol problems. In particular, increasing α results in promoting the dissemination rate λ of awareness programs; hence, the value of R
0 decreases. Therefore, media coverage for the negative effects and the number of mortalities induced by heavy drinking is very important. At the same time, this sensitivity analysis tells us that the efforts through awareness programs to increase prevention, which cause the number of aware populations X
(
t
) to increase, are effective in controlling the spread of habitual alcohol consumption. To clarify this, the numerical simulation results will be also performed in Section 5 (see Figures 2 and 3).Figure 2
The effect of increasingβ and Λ on R
0.Figure 3
Trajectories of populations and awareness program changes for different values ofλ.
(a)
(b)
(c)
(d)


### 4.2. Stability Analysis withτ
>
0

In this subsection, we will be concerned with the effect of time delay on the stability of the alcohol-present equilibriumE
1. Let u
(
t
)
=
(
a
(
t
)
,
x
(
t
)
,
p
(
t
)
,
m
(
t
)
)
T; by linearizing system (2a), (2b), (2c), and (2d) around the alcohol-present equilibrium E
1, we have(22)
d
u
d
t
=
B
u
t
+
C
u
t
-
τ
,where (23)
B
=
-
γ
+
d
+
α
0
0
0
p
γ
-
λ
M
∗
k
0
+
M
∗
-
λ
M
∗
k
0
+
M
∗
-
d
+
σ
λ
M
∗
k
0
+
M
∗
λ
k
0
P
∗
-
A
∗
-
X
∗
k
0
+
M
∗
2
-
α
0
-
d
0
α
k
0
0
-
μ
=
b
11
0
0
0
b
21
b
22
b
23
b
24
b
31
0
b
33
0
b
41
0
0
b
44
,
C
=
γ
+
d
+
α
-
β
A
∗
-
β
A
∗
β
A
∗
0
0
0
0
0
0
0
0
0
0
0
0
0
=
c
11
c
12
c
13
0
0
0
0
0
0
0
0
0
0
0
0
0
.By simple computation, the characteristic equation of the linearized system (22) is rewritten as follows:(24)
ψ
4
+
c
1
ψ
3
+
c
2
ψ
2
+
c
3
ψ
+
c
4
+
d
1
ψ
3
+
d
2
ψ
2
+
d
3
ψ
+
d
4
e
-
ψ
τ
=
0
,where(25)
c
1
=
-
b
22
-
b
33
-
b
44
,
c
2
=
b
11
b
22
+
b
11
b
33
+
b
22
b
33
+
b
22
b
44
+
b
33
b
44
,
c
3
=
-
b
11
b
22
b
33
-
b
11
b
22
b
44
-
b
11
b
33
b
44
-
b
22
b
33
b
44
,
c
4
=
b
11
b
22
b
33
b
44
,
d
1
=
-
c
11
,
d
2
=
b
44
c
11
-
b
21
c
12
-
b
31
c
13
,
d
3
=
b
22
c
11
+
b
33
c
11
+
b
21
b
33
c
12
+
b
22
b
31
c
13
+
b
21
b
44
c
12
+
b
31
b
44
c
13
-
b
22
b
33
c
11
-
b
24
b
41
c
12
-
b
31
b
23
c
12
,
d
4
=
b
24
b
33
b
41
c
12
+
b
22
b
33
b
44
c
11
+
b
23
b
31
b
44
c
12
-
b
22
b
44
c
11
-
b
33
b
44
c
11
-
b
22
b
44
b
31
c
13
-
b
21
b
33
b
44
c
12
.Now, we will derive the conditions for the stability ofE
1 and the occurrence of Hopf bifurcation. Let ψ
=
i
ω (ω
>
0) be a root of (24). By substituting it into (24) and separating the real and imaginary parts, we can get the following equations:(26)
d
2
ω
2
-
d
4
cos
⁡
ω
τ
+
d
1
ω
3
-
d
3
ω
sin
⁡
ω
τ
=
ω
4
-
c
2
ω
2
+
c
4
,
d
1
ω
3
-
d
3
ω
cos
⁡
ω
τ
+
d
4
-
d
2
ω
2
sin
⁡
ω
τ
=
-
c
1
ω
3
+
c
3
ω
.Using the fact that sin
2
⁡
(
ω
τ
)
+
cos
2
⁡
(
ω
τ
)
=
1 and denoting η
=
ω
2, we obtain(27)
f
η
=
η
4
+
e
1
η
3
+
e
2
η
2
+
e
3
η
+
e
4
=
0
,where(28)
e
1
=
d
1
2
-
2
c
2
,
e
2
=
c
2
2
+
2
c
4
+
2
d
1
d
3
-
d
2
2
,
e
3
=
2
d
2
d
4
-
2
c
2
c
4
-
d
3
2
,
e
4
=
c
4
2
-
d
4
2
.If the coefficientse
i (i
=
1,2
,
3,4
.) of (27) satisfy the Routh-Hurwitz criterion [41], then we do not get any positive value of η. In this case, (24) does not have any purely imaginary roots. In summary, we have the following result.Theorem 3.
IfR
0
>
1 and the coefficients e
i (i
=
1,2
,
3,4
.) of (27) satisfy the Routh-Hurwitz criterion, then the alcohol-present equilibrium E
1 is locally asymptotically stable for all τ
>
0.Next, we consider the case that the coefficientse
i (i
=
1,2
,
3,4
.) of (27) do not satisfy the Routh-Hurwitz criterion. For example, since l
i
m
η
→
+
∞
f
(
η
)
=
+
∞, we assume that e
4
<
0; that is, f
(
0
)
=
e
4
<
0. Then, (27) has at least one positive real root, defined as η
∗. In this case, we get (24) has a pair of purely imaginary roots; that is, ±
i
ω
∗
(
ω
∗
=
η
∗
>
0
). Furthermore, via combining this ω
∗ with (26), we can get(29)
τ
j
=
1
ω
∗
arctan
⁡
d
1
ω
∗
7
+
c
1
d
2
-
c
2
d
1
-
d
3
ω
∗
5
+
c
2
d
3
+
c
4
d
1
-
c
1
d
4
-
c
3
d
2
ω
∗
3
+
c
3
d
4
-
c
4
d
3
ω
∗
d
2
-
c
1
d
1
ω
∗
6
+
c
1
d
3
+
c
3
d
1
-
c
2
d
2
-
d
4
ω
∗
4
+
c
2
d
4
+
c
4
d
2
-
c
3
d
3
ω
∗
2
-
c
4
d
4
+
j
π
,
j
=
0,1
,
2
,
…
;then ψ
=
±
i
ω
∗ are a pair of purely imaginary roots of (24) with τ
=
τ
j (j
=
0,1
,
2
,
…). Thus, we give the following assumption:(H1) Equation (27) has at least one positive real root.Suppose that (H1) holds. Without loss of generality, we assume that (27) has four positive real roots, which are defined by η
k, k
=
1,2
,
3,4. Then, we will have four positive real roots ω
k
=
η
k, k
=
1,2
,
3,4
. Thus, we denote (30)
τ
k
j
=
1
ω
k
arctan
⁡
d
1
ω
k
7
+
c
1
d
2
-
c
2
d
1
-
d
3
ω
k
5
+
c
2
d
3
+
c
4
d
1
-
c
1
d
4
-
c
3
d
2
ω
k
3
+
c
3
d
4
-
c
4
d
3
ω
k
d
2
-
c
1
d
1
ω
k
6
+
c
1
d
3
+
c
3
d
1
-
c
2
d
2
-
d
4
ω
k
4
+
c
2
d
4
+
c
4
d
2
-
c
3
d
3
ω
k
2
-
c
4
d
4
+
j
π
,
k
=
1,2
,
3,4
,
j
=
0,1
,
2
,
…
.Then, ψ
=
±
i
ω
k is a pair of purely imaginary roots of (24) with τ
=
τ
k
{
j
}, k
=
1,2
,
3,4, j
=
0,1
,
2
,
….Define(31)
τ
0
=
τ
k
0
0
=
m
i
n
k
∈
1,2
,
3,4
τ
k
0
.In order to investigate whether the Hopf bifurcation occurs or not as time delayτ increases through τ
0, we will check whether the following transversality condition is satisfied. It is necessary to make the following assumption:(H2)f
′
(
ω
0
2
)
≠
0, where ω
0
=
ω
k
0.Assume that (H2) holds; the following result is true.Lemma 4.
Assume that (H2) holds; then the transversality condition(32)
sgn
⁡
d
R
e
ψ
τ
d
τ
τ
=
τ
0
>
0is satisfied.Proof.
Differentiating the two sides of (24) with respect to τ, we get (33)
d
ψ
τ
d
τ
-
1
=
4
ψ
3
+
3
c
1
ψ
2
+
2
c
2
ψ
+
c
3
+
3
d
1
ψ
2
+
2
d
2
ψ
+
d
3
e
-
ψ
τ
-
d
1
ψ
3
+
d
2
ψ
2
+
d
3
ψ
+
d
4
τ
e
-
ψ
τ
d
1
ψ
3
+
d
2
ψ
2
+
d
3
ψ
+
d
4
ψ
e
-
ψ
τ
.Combining (26) with (H2), we compute that (34)
d
R
e
ψ
τ
d
τ
ψ
=
i
ω
0
-
1
=
ω
0
2
f
′
ω
0
2
ω
0
2
d
1
ω
0
3
-
d
3
ω
0
2
+
d
2
ω
0
2
-
d
4
2
=
f
′
ω
0
2
d
1
ω
0
3
-
d
3
ω
0
2
+
d
2
ω
0
2
-
d
4
2
≠
0
.Suppose (35)
sgn
⁡
d
R
e
ψ
τ
d
τ
ψ
=
i
ω
0
-
1
<
0
;we let τ
<
τ
0 and τ
→
τ
0, which implies that (24) has a root ψ
(
τ
)
=
α
(
τ
)
+
i
ω
(
τ
) such that α
(
τ
)
>
0, but this contradicts the definition of τ
0. Hence, we obtain (36)
sgn
⁡
d
R
e
ψ
τ
d
τ
τ
=
τ
0
=
sgn
⁡
d
R
e
ψ
τ
d
τ
ψ
=
i
ω
0
-
1
>
0
. This completes the proof of the lemma.From the above discussion and Lemma4, we have the following theorem.Theorem 5.
IfR
0
>
1 and (H1) and (H2) hold, then the alcohol-present equilibrium E
1 is locally asymptotically stable for τ
∈
[
0
,
τ
0
) and unstable for τ
>
τ
0. System (2a), (2b), (2c), and (2d) exhibits Hopf bifurcation at the equilibrium E
1 when τ
=
τ
0.


## 4.1. Stability Analysis withτ
=
0

First, we investigate the stability under the condition ofτ
=
0. The stability of two equilibria E
0 and E
1 of system (2a), (2b), (2c), and (2d) can be given in the following theorem.Theorem 2.
Forτ
=
0, one has the following results:(i)
The alcohol-free equilibriumE
0 is globally asymptotically stable if R
0
<
1.
(ii)
The alcohol-free equilibriumE
0 is unstable and the alcohol-present equilibrium E
1 exists if R
0
>
1.
Furthermore, equilibriumE
1 is locally asymptotically stable if the coefficients of the characteristic equation of system (2a), (2b), (2c), and (2d) at E
1 satisfy the Routh-Hurwitz criterion.Proof.
Whenτ
=
0, the Jacobian matrix J of system (2a), (2b), (2c), and (2d) is (14)
J
=
β
P
-
2
A
-
X
-
γ
+
d
+
α
-
β
A
β
A
0
p
γ
-
λ
M
k
0
+
M
-
λ
M
k
0
+
M
-
d
+
σ
λ
M
k
0
+
M
λ
k
0
P
-
A
-
X
k
0
+
M
2
-
α
0
-
d
0
α
k
0
0
-
μ
.
(i) For the alcohol-free equilibriumE
0
(
0,0
,
Λ
/
d
,
0
), the matrix J at E
0 can be written as (15)
J
0
=
J
E
0
=
β
Λ
d
-
γ
+
d
+
α
0
0
0
p
γ
-
d
+
σ
0
λ
Λ
d
+
k
0
-
α
0
-
d
0
α
k
0
0
-
μ
.It is obvious that every eigenvalue of J
0 has negative real part whenever R
0
<
1. So the alcohol-free equilibrium E
0
(
0,0
,
Λ
/
d
,
0
) of system (2a), (2b), (2c), and (2d) is globally asymptotically stable if R
0
<
1.
(ii) IfR
0
>
1, then the alcohol-free equilibrium E
0 is a saddle point. So the equilibrium E
0 of system (2a), (2b), (2c), and (2d) is unstable. In this case, the alcohol-present equilibrium E
1 exists. Further, the matrix J at E
1 can be written as (16)
J
1
=
J
E
1
=
β
P
∗
-
2
A
∗
-
X
∗
-
γ
+
d
+
α
-
β
A
∗
β
A
∗
0
p
γ
-
λ
M
∗
k
0
+
M
∗
-
λ
M
∗
k
0
+
M
∗
-
d
+
σ
λ
M
∗
k
0
+
M
∗
λ
k
0
P
∗
-
A
∗
-
X
∗
k
0
+
M
∗
2
-
α
0
-
d
0
α
k
0
0
-
μ
=
J
11
J
12
J
13
0
J
21
J
22
J
23
J
24
J
31
0
J
33
0
J
41
0
0
J
44
.
The characteristic equation of system (2a), (2b), (2c), and (2d) at E
1 is d
e
t
(
φ
I
-
J
1
)
=
0, where I is the identity matrix of size 4 and φ is the eigenvalue of matrix J
1. By simple computation, the characteristic equation can be expressed: (17)
φ
4
+
b
1
φ
3
+
b
2
φ
2
+
b
3
φ
+
b
4
=
0
,where(18)
b
1
=
-
J
11
+
J
22
+
J
33
+
J
44
,
b
2
=
J
11
J
22
+
J
11
J
33
+
J
11
J
44
+
J
22
J
33
+
J
22
J
44
+
J
33
J
44
-
J
12
J
21
+
J
13
J
31
,
b
3
=
J
12
J
21
J
33
+
J
12
J
21
J
44
+
J
13
J
31
J
22
+
J
13
J
31
J
44
-
J
11
J
22
J
33
+
J
11
J
22
J
44
+
J
11
J
33
J
44
+
J
22
J
33
J
44
+
J
12
J
24
J
41
+
J
12
J
31
J
23
,
b
4
=
J
11
J
22
J
33
J
44
+
J
12
J
23
J
31
J
44
+
J
12
J
24
J
33
J
41
-
J
12
J
21
J
33
J
44
-
J
13
J
22
J
31
J
44
.
Notice thatR
0
>
1; it can be calculated that the coefficients b
i (i
=
1,2
,
3,4
.) all are positive according to J
i
i (i
=
1,2
,
3,4
.) and J
12
,
J
31 all are negative and J
13, J
21, J
23, J
24, and J
41 all are positive. Furthermore, by using the Routh-Hurwitz criterion [41], it is easy to see that the eigenvalues of the above characteristic equation have negative real parts if the conditions (19)
b
1
b
2
-
b
3
>
0
,
b
1
b
2
b
3
-
b
3
2
-
b
1
2
b
4
>
0hold. This shows that the alcohol-present equilibrium E
1 is locally asymptotically stable, and the proof is complete.Moreover, we know thatR
0 is a threshold value which determines whether alcohol problems will die out or become endemic in the population. The alcohol-free equilibrium E
0 undergoes a transcritical bifurcation at R
0
=
1. We are now interested in knowing which model parameter has the greatest effect on R
0 and hence has the greatest effect in controlling the alcohol consumption.Using the approach in Chitnis et al. [42], we calculate the normalized forward sensitivity indices of R
0. Let (20)
Υ
q
R
0
=
∂
R
0
∂
q
×
q
R
0denote the sensitivity index of R
0 with respect to the parameter q. We get(21)
Υ
β
R
0
=
Υ
Λ
R
0
=
+
1
,
Υ
d
R
0
=
-
1
+
d
γ
+
d
+
α
,
Υ
γ
R
0
=
-
γ
γ
+
d
+
α
,
Υ
α
R
0
=
-
α
γ
+
d
+
α
.Considering that the natural death rate d is a relatively stable constant, we find that Υ
β
R
0 and Υ
Λ
R
0 are the largest sensitivity indexes. We focus on reducing the value of β or Λ whose reduction by 1% would reduce R
0 by 1%. In contrast, there is a need to increase γ in order to control alcohol problems. In particular, increasing α results in promoting the dissemination rate λ of awareness programs; hence, the value of R
0 decreases. Therefore, media coverage for the negative effects and the number of mortalities induced by heavy drinking is very important. At the same time, this sensitivity analysis tells us that the efforts through awareness programs to increase prevention, which cause the number of aware populations X
(
t
) to increase, are effective in controlling the spread of habitual alcohol consumption. To clarify this, the numerical simulation results will be also performed in Section 5 (see Figures 2 and 3).Figure 2
The effect of increasingβ and Λ on R
0.Figure 3
Trajectories of populations and awareness program changes for different values ofλ.
(a)
(b)
(c)
(d)


## 4.2. Stability Analysis withτ
>
0

In this subsection, we will be concerned with the effect of time delay on the stability of the alcohol-present equilibriumE
1. Let u
(
t
)
=
(
a
(
t
)
,
x
(
t
)
,
p
(
t
)
,
m
(
t
)
)
T; by linearizing system (2a), (2b), (2c), and (2d) around the alcohol-present equilibrium E
1, we have(22)
d
u
d
t
=
B
u
t
+
C
u
t
-
τ
,where (23)
B
=
-
γ
+
d
+
α
0
0
0
p
γ
-
λ
M
∗
k
0
+
M
∗
-
λ
M
∗
k
0
+
M
∗
-
d
+
σ
λ
M
∗
k
0
+
M
∗
λ
k
0
P
∗
-
A
∗
-
X
∗
k
0
+
M
∗
2
-
α
0
-
d
0
α
k
0
0
-
μ
=
b
11
0
0
0
b
21
b
22
b
23
b
24
b
31
0
b
33
0
b
41
0
0
b
44
,
C
=
γ
+
d
+
α
-
β
A
∗
-
β
A
∗
β
A
∗
0
0
0
0
0
0
0
0
0
0
0
0
0
=
c
11
c
12
c
13
0
0
0
0
0
0
0
0
0
0
0
0
0
.By simple computation, the characteristic equation of the linearized system (22) is rewritten as follows:(24)
ψ
4
+
c
1
ψ
3
+
c
2
ψ
2
+
c
3
ψ
+
c
4
+
d
1
ψ
3
+
d
2
ψ
2
+
d
3
ψ
+
d
4
e
-
ψ
τ
=
0
,where(25)
c
1
=
-
b
22
-
b
33
-
b
44
,
c
2
=
b
11
b
22
+
b
11
b
33
+
b
22
b
33
+
b
22
b
44
+
b
33
b
44
,
c
3
=
-
b
11
b
22
b
33
-
b
11
b
22
b
44
-
b
11
b
33
b
44
-
b
22
b
33
b
44
,
c
4
=
b
11
b
22
b
33
b
44
,
d
1
=
-
c
11
,
d
2
=
b
44
c
11
-
b
21
c
12
-
b
31
c
13
,
d
3
=
b
22
c
11
+
b
33
c
11
+
b
21
b
33
c
12
+
b
22
b
31
c
13
+
b
21
b
44
c
12
+
b
31
b
44
c
13
-
b
22
b
33
c
11
-
b
24
b
41
c
12
-
b
31
b
23
c
12
,
d
4
=
b
24
b
33
b
41
c
12
+
b
22
b
33
b
44
c
11
+
b
23
b
31
b
44
c
12
-
b
22
b
44
c
11
-
b
33
b
44
c
11
-
b
22
b
44
b
31
c
13
-
b
21
b
33
b
44
c
12
.Now, we will derive the conditions for the stability ofE
1 and the occurrence of Hopf bifurcation. Let ψ
=
i
ω (ω
>
0) be a root of (24). By substituting it into (24) and separating the real and imaginary parts, we can get the following equations:(26)
d
2
ω
2
-
d
4
cos
⁡
ω
τ
+
d
1
ω
3
-
d
3
ω
sin
⁡
ω
τ
=
ω
4
-
c
2
ω
2
+
c
4
,
d
1
ω
3
-
d
3
ω
cos
⁡
ω
τ
+
d
4
-
d
2
ω
2
sin
⁡
ω
τ
=
-
c
1
ω
3
+
c
3
ω
.Using the fact that sin
2
⁡
(
ω
τ
)
+
cos
2
⁡
(
ω
τ
)
=
1 and denoting η
=
ω
2, we obtain(27)
f
η
=
η
4
+
e
1
η
3
+
e
2
η
2
+
e
3
η
+
e
4
=
0
,where(28)
e
1
=
d
1
2
-
2
c
2
,
e
2
=
c
2
2
+
2
c
4
+
2
d
1
d
3
-
d
2
2
,
e
3
=
2
d
2
d
4
-
2
c
2
c
4
-
d
3
2
,
e
4
=
c
4
2
-
d
4
2
.If the coefficientse
i (i
=
1,2
,
3,4
.) of (27) satisfy the Routh-Hurwitz criterion [41], then we do not get any positive value of η. In this case, (24) does not have any purely imaginary roots. In summary, we have the following result.Theorem 3.
IfR
0
>
1 and the coefficients e
i (i
=
1,2
,
3,4
.) of (27) satisfy the Routh-Hurwitz criterion, then the alcohol-present equilibrium E
1 is locally asymptotically stable for all τ
>
0.Next, we consider the case that the coefficientse
i (i
=
1,2
,
3,4
.) of (27) do not satisfy the Routh-Hurwitz criterion. For example, since l
i
m
η
→
+
∞
f
(
η
)
=
+
∞, we assume that e
4
<
0; that is, f
(
0
)
=
e
4
<
0. Then, (27) has at least one positive real root, defined as η
∗. In this case, we get (24) has a pair of purely imaginary roots; that is, ±
i
ω
∗
(
ω
∗
=
η
∗
>
0
). Furthermore, via combining this ω
∗ with (26), we can get(29)
τ
j
=
1
ω
∗
arctan
⁡
d
1
ω
∗
7
+
c
1
d
2
-
c
2
d
1
-
d
3
ω
∗
5
+
c
2
d
3
+
c
4
d
1
-
c
1
d
4
-
c
3
d
2
ω
∗
3
+
c
3
d
4
-
c
4
d
3
ω
∗
d
2
-
c
1
d
1
ω
∗
6
+
c
1
d
3
+
c
3
d
1
-
c
2
d
2
-
d
4
ω
∗
4
+
c
2
d
4
+
c
4
d
2
-
c
3
d
3
ω
∗
2
-
c
4
d
4
+
j
π
,
j
=
0,1
,
2
,
…
;then ψ
=
±
i
ω
∗ are a pair of purely imaginary roots of (24) with τ
=
τ
j (j
=
0,1
,
2
,
…). Thus, we give the following assumption:(H1) Equation (27) has at least one positive real root.Suppose that (H1) holds. Without loss of generality, we assume that (27) has four positive real roots, which are defined by η
k, k
=
1,2
,
3,4. Then, we will have four positive real roots ω
k
=
η
k, k
=
1,2
,
3,4
. Thus, we denote (30)
τ
k
j
=
1
ω
k
arctan
⁡
d
1
ω
k
7
+
c
1
d
2
-
c
2
d
1
-
d
3
ω
k
5
+
c
2
d
3
+
c
4
d
1
-
c
1
d
4
-
c
3
d
2
ω
k
3
+
c
3
d
4
-
c
4
d
3
ω
k
d
2
-
c
1
d
1
ω
k
6
+
c
1
d
3
+
c
3
d
1
-
c
2
d
2
-
d
4
ω
k
4
+
c
2
d
4
+
c
4
d
2
-
c
3
d
3
ω
k
2
-
c
4
d
4
+
j
π
,
k
=
1,2
,
3,4
,
j
=
0,1
,
2
,
…
.Then, ψ
=
±
i
ω
k is a pair of purely imaginary roots of (24) with τ
=
τ
k
{
j
}, k
=
1,2
,
3,4, j
=
0,1
,
2
,
….Define(31)
τ
0
=
τ
k
0
0
=
m
i
n
k
∈
1,2
,
3,4
τ
k
0
.In order to investigate whether the Hopf bifurcation occurs or not as time delayτ increases through τ
0, we will check whether the following transversality condition is satisfied. It is necessary to make the following assumption:(H2)f
′
(
ω
0
2
)
≠
0, where ω
0
=
ω
k
0.Assume that (H2) holds; the following result is true.Lemma 4.
Assume that (H2) holds; then the transversality condition(32)
sgn
⁡
d
R
e
ψ
τ
d
τ
τ
=
τ
0
>
0is satisfied.Proof.
Differentiating the two sides of (24) with respect to τ, we get (33)
d
ψ
τ
d
τ
-
1
=
4
ψ
3
+
3
c
1
ψ
2
+
2
c
2
ψ
+
c
3
+
3
d
1
ψ
2
+
2
d
2
ψ
+
d
3
e
-
ψ
τ
-
d
1
ψ
3
+
d
2
ψ
2
+
d
3
ψ
+
d
4
τ
e
-
ψ
τ
d
1
ψ
3
+
d
2
ψ
2
+
d
3
ψ
+
d
4
ψ
e
-
ψ
τ
.Combining (26) with (H2), we compute that (34)
d
R
e
ψ
τ
d
τ
ψ
=
i
ω
0
-
1
=
ω
0
2
f
′
ω
0
2
ω
0
2
d
1
ω
0
3
-
d
3
ω
0
2
+
d
2
ω
0
2
-
d
4
2
=
f
′
ω
0
2
d
1
ω
0
3
-
d
3
ω
0
2
+
d
2
ω
0
2
-
d
4
2
≠
0
.Suppose (35)
sgn
⁡
d
R
e
ψ
τ
d
τ
ψ
=
i
ω
0
-
1
<
0
;we let τ
<
τ
0 and τ
→
τ
0, which implies that (24) has a root ψ
(
τ
)
=
α
(
τ
)
+
i
ω
(
τ
) such that α
(
τ
)
>
0, but this contradicts the definition of τ
0. Hence, we obtain (36)
sgn
⁡
d
R
e
ψ
τ
d
τ
τ
=
τ
0
=
sgn
⁡
d
R
e
ψ
τ
d
τ
ψ
=
i
ω
0
-
1
>
0
. This completes the proof of the lemma.From the above discussion and Lemma4, we have the following theorem.Theorem 5.
IfR
0
>
1 and (H1) and (H2) hold, then the alcohol-present equilibrium E
1 is locally asymptotically stable for τ
∈
[
0
,
τ
0
) and unstable for τ
>
τ
0. System (2a), (2b), (2c), and (2d) exhibits Hopf bifurcation at the equilibrium E
1 when τ
=
τ
0.


## 5. Numerical Simulations

With the help ofMATLAB, some numerical results of system (2a), (2b), (2c), and (2d) are provided to substantiate the analytic results obtained in this section. We choose a set of values of parameters: Λ
=
5; β
=
0.00002; λ
=
0.02; k
0
=
100; σ
=
0.001; p
=
1; γ
=
0.2; d
=
0.00004; α
=
0.002; k
=
2.5; μ
=
0.05. It is easy to check that the basic reproduction number(37)
R
0
=
12.3738
>
1is satisfied and one positive root A
∗
=
291.6366 can be obtained by solving (12). So, the alcohol-present equilibrium E
1 of system (2a), (2b), (2c), and (2d) is obtained as follows: (38)
A
∗
=
291.6366
,
X
∗
=
100075.5350
,
P
∗
=
110468.1715
,
M
∗
=
29.0637
.First, we can see that the value ofR
0 increases as the value of β or Λ increases, which means that β and Λ play a key role in ensuring that the alcohol problem occurs (see Figure 2).Second, the effect of awareness programs on the alcohol consumption behavior is demonstrated in Figure3. It is apparent that awareness programs are an effective measure in controlling heavy drinking from this figure. We note that the aware population X
(
t
) increase whereas heavy drinkers A
(
t
) decrease with an increase in the value of dissemination rate λ of awareness programs among susceptible population while the other parameters remain unchanged. Hence, alcohol problems should be controlled through the execution of awareness programs.Moreover, according to the above given parameter values, the critical value of time delayτ
0
=
1.5767 can be obtained by solving the corresponding expression. From Theorem 5, the corresponding waveforms are shown in Figures 4–7. That is, when τ
=
1.50
<
τ
0
=
1.5767, the alcohol-present equilibrium E
1 of system (2a), (2b), (2c), and (2d) is locally asymptotically stable (see Figures 4 and 5, resp.); and Hopf bifurcation occurs once τ
=
2.50
>
τ
0
=
1.5767 (see Figures 6 and 7, resp.). These figures indicate that the system is stable for lower values of τ but for its higher values periodic solutions may arise and system may lose its stability.Figure 4
The alcohol equilibriumE
1 of system (2a), (2b), (2c), and (2d) is locally asymptotically stable with τ
=
1.50
<
τ
0
=
1.5767.
(a)
(b)
(c)
(d)Figure 5
This corresponding phase plot of system (2a), (2b), (2c), and (2d) in A
(
t
)
-
X
(
t
) plane.Figure 6
Hopf bifurcations occur withτ
=
2.50
>
τ
0
=
1.5767.
(a)
(b)
(c)
(d)Figure 7
Periodic orbit of system (2a), (2b), (2c), and (2d) with τ
=
2.50 in A
(
t
)
-
X
(
t
) plane.


## 6. Conclusions and Discussions

The goal of this paper is to analyze the impacts of awareness programs and time delay on the alcohol consumption behavior. A vast body of literature (see [11–35]) has been used for studying alcohol problems and describing the effects of awareness programs by media on the infectious disease or alcohol consumption behavior. However, most of these models are ODE models and do not incorporate the effects of the time delay. Based on the previous works, we establish an alcohol consumption model with awareness programs and time delay in Section 2. Compared to previous models, a key novelty of our model is that we introduce the delay in our model to describe the time needed for a susceptible individual to become an infectious alcohol user. To be more realistic, we considered interaction between susceptible population and awareness programs with Holling-II functional response function. In addition, we assume that the growth rate of the cumulative density of awareness programs is proportional to the number of mortalities induced by heavy drinking. We obtain the conditions for the existence and stability of two equilibria of our model in Section 3.At the same time, the condition of transcritical bifurcation is also given whenτ
=
0. In order to study the effect of time delay τ
>
0 on system (2a), (2b), (2c), and (2d), we compute the linear form of system (2a), (2b), (2c), and (2d) around the alcohol-present equilibrium E
1 and analyze the local stability of E
1. In addition, we obtain the existence of Hopf bifurcation by taking time delay τ as bifurcation parameter. Our results suggest that the alcohol-present equilibrium E
1 of system (2a), (2b), (2c), and (2d) is locally asymptotically stable for τ
∈
[
0
,
τ
0
) and unstable for τ
>
τ
0. System (2a), (2b), (2c), and (2d) undergoes a Hopf bifurcation at the equilibrium E
1 when τ
=
τ
0. Furthermore, we study the impacts of awareness programs on the alcohol consumption behavior. Our findings show that awareness programs are an effective measure in controlling heavy drinking.Not only awareness programs but also time delay is introduced in our model. It can help us to reduce the economic burden of disease caused by alcohol abuse. It is also beneficial to the social problems such as traffic accident and violent crime. Whether in theory or in the practical sense, these issues are very interesting and need further studies.



---

*Source: 260195-2015-10-19.xml*