Abstract

A more realistic two-stage model for binge drinking problem is introduced, where the youths with alcohol problems are divided into those who admit the problem and those who do not admit it. We also consider the direct transfer from the class of susceptible individuals towards the class of admitting drinkers. Mathematical analyses establish that the global dynamics of the model are determined by the basic reproduction number, . The alcohol-free equilibrium is globally asymptotically stable, and the alcohol problems are eliminated from the population if . A unique alcohol-present equilibrium is globally asymptotically stable if . Numerical simulations are also conducted in the analytic results.

1. Introduction

Young people’s binge drinking problem is a major concern to public health. Recently, US surveys indicate that approximately 90% of college students have consumed alcohol at least once [1], and more than of college students have engaged in binge drinking [2, 3]. The binge drinking refers to youths in 17–30 age group who drink a large amount of alcohol and become so drunk; they are likely to exhibit antisocial behavior [4]. Although there have been many attempts to reduce the problem, alcohol abuse by college students has persisted and in some cases increased over the past several decades [5]. Prior studies have indicated that heavy alcohol drinkers are likely to engage in risky sexual behaviours and more likely to get sexually transmitted infection than social drinkers [6, 7]. There is a strong medical evidence that treatment of individuals with alcohol problems is a major issue [810].

Thus, it is very important to use a mathematical method to study the binge drinking problems in youths. A simple model for alcohol treatment is presented by Sanchez et al. [11]. Since then, there have been numerous studies investigating campus drinking and the associated consequences [1215]. Manthey et al. [12] focus on a college campus, divide the student population into three classes: nondrinkers, social drinkers and problem drinkers, and show that campus alcohol abuse may be reduced by minimizing the ability of problem drinkers to directly recruit nondrinkers. Cintron-Arias et al. [13] focus on situations where relapse rates are high and conclude that the systematic removal of individuals from high-risk environments, or the development of programs that limit access or reduce the residence times in such environments (or both approaches combined), may reduce the level of alcohol abuse. Mubayi et al. [14] show that if the relative residence times of moderate drinkers are distributed randomly between low- and high-risk environments, then the proportion of heavy drinkers is likely to be higher than expected. Mulone and Straughan [15] investigate a model for binge drinking taking into account admitting and nonadmitting drinkers. But the global stability of binge drinking model is not discussed in the literature.

Motivated by the binge drinking model in [15], we develop a more realistic model with two stages. Drinking is often encouraged by peer pressure. A susceptible individual acquires alcohol problems through the direct contact with the admitting drinker or the nonadmitting drinker. The new drinker can become either the admitting drinker or the nonadmitting drinker. So we consider the direct transfer from the class of susceptible individuals towards the class of admitting drinkers; furthermore, we study the global dynamics of the model. The reason to introduce this new direct transfer is that about one-third of the American population admit to drinking problems, of the population admit to the binge drinking problem [16]. Therefore, this fact cannot be neglected in the binge drinking model.

The organization of this paper is as follows. In the next section, the binge drinking model with two stages and some basic properties are derived. In Section 3, the existence and the global stability of equilibria are investigated. Some numerical simulations are given in Section 4. Some discussions are given in Section 5.

2. The Model

2.1. System Description

The total population is divided into four compartments, namely, the susceptible compartment of those who do not drink or drink only moderately, denoted by , those who drink heavily at least some of the time but do not admit having a problem, denoted by , those who drink heavily and admit having a problem, denoted by , and those people in treatment, denoted by . The total number of population at time is given by The model structure is shown in Figure 1. The transfer diagram leads to the following system of ordinary differential equations: where is the number of individuals entering into the system in a given time interval (say each year), so represents the rate of entry. Since we are dealing with youths, we assume that the death rate is negligible, and so the leaving rate is also . A susceptible individual can be turned to drink through direct contact with an admitting drinker or a nonadmitting drinker. is the transmission coefficient for the individuals who do not admit to have alcohol problems, is the transmission coefficient for the individuals who admit to have these problems. We assume that due to the fact that they are more likely to be unaware of their condition. is the rate at which represents those with the alcohol problems admitting to have these problems and then transferring from class to , is that fraction of who go into treatment, is that fraction of who relapse into admitting drinkers (we adopt a linear relapse term rather than peer pressure since we argue that the relapse is primarily due to the person), is the probability of a new drinker to admit having the problem. We may show that is constant, and then we introduce the fractions of , , , and : with . Then the system (2.2) becomes As system (2.4) is equivalent to system (2.2), we only need to study system (2.4).

2.2. Basic Properties
2.2.1. Invariant Region

Adding all the equations of system (2.4) gives

Thus, the total population is a constant. Since system (2.4) monitors human population, it is plausible to assume that all its state variables and parameters are nonnegative for all . Further, it can be shown that the region is positively invariant. Thus, each solution of system (2.4), with initial conditions in , remains there for . Therefore, the -limit sets of the solutions of system (2.4), in , are contained in . Furthermore, in , the usual existence, uniqueness, and continuation results hold for the system, so that the system (2.4) is well posed mathematically and epidemiologically [17]. So we consider the dynamics of system (2.4) on the set in this paper.

2.2.2. Positivity of Solutions

For system (2.4), it is necessary to prove that all the state variables are positive, so that the solutions of the system with positive initial conditions remain positive for all . We thus state the following lemma.

Lemma 2.1. If and , the solutions , , , and of system (2.4) are positive for all .

Proof. Under the given initial conditions, it is easy to prove that the solutions of system (2.4) are positive; if not, we assume a contradiction: that there exists a first time such that there exists a , there exists a , and there exists a , In the first case, we have which is a contradiction, meaning that , . In the second case, we have which is a contradiction, meaning that , . Similarly, it can be shown that and for all .

Thus, the solutions , , ?, and of system (2.4) remain positive for all .

3. Analysis of the Model

There are one alcohol-free equilibrium and one alcohol-present equilibrium for system (2.4).

3.1. Alcohol-Free Equilibrium and the Reproduction Number

The model has an alcohol-free equilibrium given by In the following, the basic reproduction number of system (2.4) will be obtained by the next-generation matrix method formulated in [18].

Let , then system (2.4) can be written as where The Jacobian matrices of and at the alcohol-free equilibrium are, respectively, where The model reproduction number, denoted by , is thus given by Following Theorem 2 of [18], we have the following result on the local stability of .

Theorem 3.1. The alcohol-free equilibrium is locally asymptotically stable for and unstable otherwise.

3.2. Global Stability of

Theorem 3.2. For system (2.4), the alcohol-free equilibrium is globally asymptotically stable if .

Proof. We use the comparison theorem to prove the global stability of the alcohol-free equilibrium. The rate of change of the variables of system (2.4) can be rewritten as where and are defined in (3.5). Since for all in , then Since the eigenvalues of the matrix all have negative real parts (this comes from the local stability results in Lemma 1 in [18]), then system (2.4) is stable whenever . So, as . By the comparison theorem [19], it follows that and as . The as . So, is globally asymptotically stable for .

3.3. Alcohol-Present Equilibrium
3.3.1. Existence of the Alcohol-Present Equilibrium

If , system (2.4) has a unique alcohol-present equilibrium , where

3.3.2. Global Stability of the Alcohol-Present Equilibrium

Theorem 3.3. If , the alcohol-present equilibrium is globally asymptotically stable.

Proof. To study the global stability of the alcohol-present equilibrium, motivated by [2022], we use a Lyapunov function as follows: Applying the identity , the derivative of is given by
The positive constants , ?, ?, and are chosen such that the coefficients of , ?, ?, ?, and are equal to zero, that is, From (3.12), we have We regroup terms in such that , where Using the values for , ?, ?, and in (3.13), relations in (3.12), and the equilibrium relations, We can rewrite as Similarly, we can rewrite as where Let , with each term representing the expression enclosed in a pair of big square brackets. Using the arithmetic-geometric mean inequality, we obtain From (3.19), we have From (3.16) and (3.20), we have . Therefore, . Furthermore, if and only if and . We can show that . Thus, is a negative definite with respect to . So the alcohol-present equilibrium is globally asymptotically stable.

4. Numerical Simulation

To illustrate the analytic results obtained above, we give some simulations using the parameter values in Table 1. Numerical results are displayed in the following figures. First, we choose ,?; numerical simulation gives , then Theorem 3.2 indicates that youths do not have alcohol problems. Figure 2 confirms this conclusion. Second, we choose , ?; numerical simulation gives ; then Theorem 3.3 shows that alcohol problems persist. Figure 3 further validates the conclusion. Finally, we choose ; numerical simulation gives the relation between and . Figure 4 confirms that shows a decline, while the number of admitting drinkers increases. From the figures above, we find that making more people admit having alcohol problems can reduce the alcohol problems. Then, the treatment of admitting drinkers is an effective measure in alcohol problems.

5. Discussion

We have formulated a binge drinking model with two stages and investigated their dynamical behaviors. Depending on the basic reproduction number , the steady state is either the alcohol free or the alcohol-present. By using the comparison theorem, we prove that all solutions converge to when the basic reproduction number is less than one, that is, the alcohol problems disappear eventually. By constructing the Lyapunov function, we prove that the unique alcohol-present equilibrium is globally stable, that is, the alcohol problems will persist in the population, and the number of binge drinking individuals tends to be a positive constant when the basic reproduction number exceeds one.

To better understand the binge drinking model with two stages, we consider the more realistic model and investigate the global stability of equilibria of the model for binge drinking problem. The global stability is not discussed in the literature [15].

The basic reproductive number in this paper is given by

If the model has only one transfer from the class of susceptible individuals towards the class of nonadmitting drinkers , then the basic reproductive number is It is similar to the basic reproductive number in [15]. Comparing with the two expressions, we see that is a function on . In Figure 4, we find that the basic reproductive number decreases when the number of admitting drinkers increases. In general, the nonadmitting drinkers are also unlikely to seek help. Making more people admit to have alcohol problems through the awareness programs, such as education and media programs is an effective measure to control the alcohol problems.

In deriving our model (2.2), we consider the direct transfer from the class of susceptible individuals towards the class of admitting drinkers and assume the relapse from the class of treatment individuals towards the class of admitting drinkers. However, if we hope to include another relapse from the class of treatment individuals towards the class of nonadmitting drinkers, we could modify (2.2) to the following model: where represents relapse from the class of treatment individuals towards the class of nonadmitting drinkers.

In another way, we might consider the class of quit drinkers and also consider the following model: where represents the population that has stopped drinking permanently at time , say quit drinkers. The term is the rate from the class of treatment individuals towards the class of quit drinkers. We leave these works for the future.

Acknowledgments

This work was partially supported by the NNSF of China (10961018), the NSF of Gansu Province of China (1107RJZA088), the NSF for Distinguished Young Scholars of Gansu Province of China (1111RJDA003), the Special Fund for the Basic Requirements in the Research of University of Gansu Province of China, and the Development Program for HongLiu Distinguished Young Scholars in Lanzhou University of ?Technology.