Estimating the dynamics of the drinking epidemic model with control interventions: A sensitivity analysis

: This article presents a non-linear mathematical model that captures the dynamics of drinking prevalence within a population. The model is analyzed under an optimal control framework, dividing the total population into four compartments: susceptible, heavy drinker, drinker in treatment, and recovered classes. The model ’ s validity is a ﬃ rmed through considerations of positivity, boundedness, reproduction number, stability, and sensitivity analysis. Stability theory is employed to explore both local and global stabilities. Sensitivity analysis identi ﬁ es parameters with a signi ﬁ cant impact on the reproduction number ( R 0 ), with maximum sensitivity observed in parameters related to drinking transmission and transitions from heavy drinking to treatment stages. These parameters exhibit sensitivity indices of ( ) 0.538, 1 , indicating that a 10% increase in these parameters would result in a ( ) 5.38, 1 increase in the threshold quantity. The study introduces an optimal control strategy that involves awareness campaigns and treatment as control variables. These controls aim to minimize the number of heavy drinkers while maximizing the number of recovered individuals. Pontryagin ’ s maximum principle is used to solve optimal control problems. Additionally, the research explores various parametric settings for each compartment, enriching the study environment. The e ﬀ ective-ness of the proposed control scheme is evaluated through rigorous numerical simulations, highlighting its competitive edge. The results, validated using MATLAB simulations, are detailed throughout the article.


Introduction
The ever-evolving research literature nominates the association of several physical, psychological, and socio-economic issues with the access use of alcohol.For example, Lees et al. [1] provided an enlightened account documenting the linkages of excessive alcohol and brain functioning ranging from compromised decision-making to memory loss.Similarly, World Health Organization [2,3] endorsed an intriguing association between chronic alcohol use and the development of medical complications, including a higher risk of liver disease, increased odds of cardiovascular issues, heart problems, high blood pressure, weakened immune system, and an increased risk of certain types of cancer.Furthermore, a thought-provoking investigation of Ferrari et al. [4] found vivid links connecting alcohol use with the streaming of negative emotions, such as depression, self-harm, suicidal attitude, anxiety, and alcohol dependence or addiction.Heavy drinking can also impair cognitive function, leading to difficulties with memory, attention, and decision-making.The main factors of alcohol drinking include income levels, education, employment status, and social norms.Research has found that individuals with lower incomes and lower levels of education are more likely to engage in heavy drinking [5].Unemployment and job instability have also been associated with increased alcohol consumption [6].Social norms and peer pressure can also influence drinking behaviors, as alcohol consumption is often seen as a way to fit in or cope with stress.The negative impulses of alcohol consumption are known to have drastic effects on more fragile groups of society, i.e., children.Excess drinking in mothers during pregnancy is found to be related to fetal alcohol spectrum disorders in children.Moreover, in a wider perspective, WHO Global Status Report [7] reported an astonishing estimate of 30% of absenteeism associated with alcohol, in Costa Rica.However, in the UK, the share of alcohol usage in fatal accidents is found to be 60%.The degree of burden of alcohol consumption on the health care systems can be witnessed from the documentation of Johnson et al. [8] reporting a 47% increase in the rate of all alcohol-related emergency visits from 2006 to 2014.From a global perspective, the unchecked use of alcohol is found to be related to the instigation of more than 200 medical complexities, including liver cirrhosis, malignancies, and cardiovascular disease [9,10].The resulting number of deaths is counted as higher as 3.3 million fatalities per annum, contributing to a notable share of 6% of all deaths, globally [3].More alarmingly, excessive alcohol consumption is becoming an emerging trend in adolescents aged 15-19 years [11].Almost 46% of the world's adolescents aged 15-19 years reported having overused alcohol.In 2016, the National Survey on Drug Use and Health reported that 85.6% of people aged 18 and older drank alcohol, in the USA.
Dynamic modeling, employing mathematical frameworks such as differential equations, unveils the temporal evolution of systems, articulating variable interactions for precise predictions and analysis across diverse fields.Mathematical representations unravel patterns, identify parameters, and enhance our understanding of complex systems.The impact of dynamic modeling spans scientific progress and technological innovations, underscoring its significance across disciplines [12].Chen has made significant contributions to dynamic modeling, particularly in integrated energy systems and natural gas dynamics [13].His work involves developing robust state estimators and Kalman filter-based approaches to enhance the accuracy and reliability of state predictions [14].Focusing on complex energy networks and pipeline systems, Chen's research addresses critical challenges and offers insights for advancing dynamic modeling methodologies in the energy sector.
Mathematical studies serve as invaluable tools for evaluating, testing, and implementing strategies on both short-and long-term scales, particularly in addressing chronic relapsing diseases such as alcoholism.There has been limited research on mathematical modeling methods for social problems such as alcohol and drug use, despite these issues often being referred to as epidemics, in contrast to research on other epidemic problems.However, in recent years, an increasing number of mathematicians have dedicated their efforts to developing mathematical models aimed at understanding and addressing alcoholrelated concerns.Khan and Khan [15] proposed a mathematical model for drinking using numerical analysis of fractional order.Chinnadurai et al. [16] focused on employing mathematical modeling to analyze strategies for controlling alcohol consumption, particularly examining their impact on socioeconomically disadvantaged populations.Through the optimal control approach, Imken and Fatmi [17] introduced a novel mathematical model to investigate alcohol consumption within the diabetic population, especially under antidiabetic drug treatments.ur Rahman et al. [18] proposed a fractional mathematical model for drinking under Atangana-Baleanu Caputo derivatives.Bunonyo et al. [19] research focused on using the eigenvalue method to develop a mathematical model for the time-dependent concentration of alcohol in the human bloodstream, contributing to a deeper understanding of alcohol metabolism dynamics and potentially informing strategies for alcoholrelated risk assessment and management.Huo and Wang's nonlinear mathematical model [20] illustrated the impact of awareness campaigns on binge drinking, demonstrating their effectiveness in mitigating alcohol-related issues.Ma et al. [21] tested a mathematical model of alcoholism as a communicable disease, incorporating awareness campaigns and a time delay, using optimal control techniques.Wang et al. [22] proposed and investigated a nonlinear model of alcoholism with optimal control to prevent interactions between susceptible and infected individuals.Sharma et al. [23] established a mathematical model of alcohol consumption, examining the stability and existence of an endemic equilibrium without drinking, and analyzing the sensitivity of R 0 [24,25].Manthey et al. [26] proposed a mathematical model to study the dynamics of campus drinking as an epidemiological model.Huo et al. [27] proposed a new social epidemic model of alcoholism with media coverage aimed at encouraging people to abstain from drinking.Huo and Song [28] distinguished between heavy drinkers who admit to drinking and those who do not, developing a two-stage model for problematic drinking that accounts for the shift in drinkers' status from susceptible individuals to admitting drinkers.This article discusses optimal control analysis of disease transmission within a community, as presented in the study by Anjam et al. [29].
In recent times, there has been a growing emphasis on employing optimal control strategies to improve viral transmission management.These applications extend beyond medical contexts and encompass various areas such as policy development, engineering risk assessment, emergency preparedness, and control program refinement [30,31].Optimal control strategies are designed to identify effective protocols for addressing complex problems, utilizing mathematical techniques such as dynamic programming, Pontryagin's principle, or model predictive control [32].Utilizing optimal control techniques in mathematical models offers benefits by enabling the identification of efficient control strategies to achieve desired outcomes, including improved system performance, maximized benefits, or reduced costs [33].Verma et al. proposed a mathematical model to investigate coronavirus dynamics, considering the impact of lockdown as an epidemiological measure [34].They also studied a nonlinear smoking model to analyze interactions between smokers and smoking quitters [35].Moreover, optimal control strategies have been introduced for coronavirus disease 2019, human immunodeficiency virus, and Dengue through mathematical models that capture relationships among these infectious diseases [36].For further insights, references such as Omame et al. [37] sightsee backward bifurcation and optimal control in a co-infection model for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and zika virus, while Omame et al. [38] delve into SARS-CoV-2 and HBV co-dynamics modeling.
Motivated by the importance of the aforementioned issue, this study delves into the application of optimal control theory to devise a new mathematical framework capable of illustrating drinking prevalence trends in a community.The goal of the optimal control scheme is to establish feasible control laws that optimize specific performance metrics.This article introduces a proposed mathematical model for a drinking epidemic, categorizing the entire population into four classes: susceptible drinkers, individuals in treatment, heavy drinkers, and recovered drinkers.Initially, we validate the proposed model by examining fundamental properties such as boundedness, positivity, reproduction number, equilibrium points, stability, and conducting sensitivity analyses.We conduct comprehensive stability analyses at equilibrium points to assess both local and global stabilities, outlining the legitimacy of the model through detailed sensitivity analysis and documenting local and global features.Subsequently, we develop the drinking epidemic model using control variables derived from qualitative optimal control principles, aiming to maximize susceptible and recovered drinkers while minimizing heavy drinkers and optimizing the aforementioned compartments.The optimal control involves two control variables: a social awareness campaign ( ) u t 1 and treatment assistance for reducing the drinking population ( ) u t

2
. Employing Pontryagin's maximum principle (Pontryagin et al. [39]), we establish a proof for the existence of the developed optimal control problem and analyze it to identify the necessary optimality conditions.The parametric setting covers a range of scenarios.
This article is organized into the following sections.Section 2 concentrates on model analysis, while Section 3 delves into sensitivity analysis.Section 4 discusses the stability aspects of the model under consideration.Section 5 elaborates on the integration of optimal control, emphasizing the implementation of more effective health interventions.Finally, Section 6 summarizes the key findings of this study.

Model formulation
Understanding the transmission patterns of diseases and creating effective control strategies heavily relies on mathematical models.Therefore, it is essential to focus on delineating the epidemiological facets of the disease and pinpointing crucial, adjustable parameters to bolster disease control measures.
This study delves into exploring a mathematical model aimed at understanding the drinking epidemic within the human population [40].It highlights several key limitations in the model analysis, including simplification of real-world complexities, potential oversights in parameter selection, reliance on assumed relationships, and the absence of dynamic external factors.These constraints may limit the model's ability to fully capture the complexities of human behavior and societal dynamics, potentially impacting the generalizability of the findings.Therefore, caution is advised in interpreting the predictions, as real-world applicability may be constrained by the inherent simplifications in the modeling approach.Additionally, certain assumptions were made in developing the proposed drinking epidemic model.The epidemic occurs within a closed environment, where factors such as sex, race, and social status do not influence the likelihood of developing heavy drinking habits.Members interact homogeneously, meaning they have an equal degree of interaction with one another.Heavy drinking behavior can be transmitted to non-drinkers when they come into contact with individuals who regularly consume large amounts of alcohol.Individuals undergoing treatment for their drinking habits may transition through stages of recovery and susceptibility, potentially leading to a return to heavy drinking.Conversely, those who have successfully stopped drinking enter the recovery compartment.
At any time ≥ t 0, the mathematical model divides the entire population ( ) t , signifying individuals who successfully recover from alcoholism through treatment and subsequently maintain sustained abstinence from alcohol.This structured approach enables us to account for the varying levels of susceptibility, severity of drinking habits, treatment interventions, and recovery rates within the population.It allows us to model the transitions between these compartments over time, reflecting the real-world complexities of drinking behavior and recovery processes.By incorporating specific parameters such as recruitment rates, transmission rates, death rates, and recovery rates into the model, we can simulate and analyze the interactions and flow patterns between these compartments more accurately.Moreover, the model formulation incorporates specific parameters such as the recruitment rate of the susceptible class b, the transmission rate from the susceptible class to the heavy drinker class α, the transmission rate from the recovered class to the susceptible class η, the natural death rate μ, the drinking- induced death rate in the class of heavy drinkers δ 1 , the drinking-induced death rate in the treatment class δ 2 , the proportion of individuals entering the drinker in treatment class ϕ, and the recovery rate γ.Considering the aforemen- tioned factors, the schematic diagram in Figure 1 depicts the flow pattern of interactions between different states related to drinking behavior.Hence, the total population at any given time t is expressed as: Therefore, the fundamental model for the drinking epidemic is governed by a system of nonlinear differential equations, which include: subject to the initial conditions, The specifics of each term that impact the system, as outlined in Model (1), are identified in Table 1.
2 Qualitative analysis of the proposed model In the forthcoming section, our analysis will focus on the mathematical formulation encapsulated in Model (1).Our examination will delve into the aspects of boundedness and positivity within the model.Additionally, we will determine the equilibrium points and compute the basic reproduction number.

Positivity
Positivity constraints are imposed on both the initial conditions and parameters in the devised model, ensuring that their values remain non-negative or greater than zero throughout the modeling process.This imposition enhances the model's reliability by minimizing the occurrence of unrealistic outcomes, thereby improving its applicability and making it more representative of real-world scenarios and behaviors.Consequently, we rearrange Model (1) in accordance with these constraints.
where ( ) , ( ) ( ( ) In a stress-free scenario, Consistent with the seminal discovery by Nagumo [41], the solution of Model (1) with an initial condition

Boundedness
Boundedness in the drinking model refers to inherent constraints within the study, such as the exclusive focus on specific demographics or contexts.These limitations risk compromising the model's reliability by narrowing its perspective and potentially overlooking broader behavioral patterns.As a result, reduced generalizability restricts the model's applicability to specific circumstances and populations, which may compromise its ability to capture the nuanced complexity of real-world drinking behaviors.
To understand the aforementioned model, we need the bounds on the dependent variables involved.To achieve this, we must determine the feasible region of attraction as outlined in the following theorem.
Theorem 2.1.There is a positive ψ for nonzero such that all solutions meet for a long time t.
Proof.The solutions of System (1) are greater than zero, now in the first compartment of System (1) as Thus, for undue time t, let us say > t t 0 for an excessive time t.Define ( ) ( ) . Now, differentiating of R 1 with the respect to the solution of Model (1) yields where ( ( 2 .Remember ( ) ≤ t 1 + b μ for any > t t 0 .ψ 1 exists, reliant merely on the para- meters of the considered model, such that ( ) ≤ R t ψ 1 1 for eventually > t t 0 and ( ) t and ( ) t are bounded above.Subsequently, the second and third compartments of Model (1), are eventually bounded above, and let ψ be the max- imum.This displays that Model (1) is destructive.

Equilibrium points
Equilibrium points denote the stable-state prevalence of individuals engaged in drinking behavior, wherein the influx of new drinkers is counterbalanced by the rates of recoveries and mortality.There exist two equilibrium points: the drinking-free equilibrium point, and the drinkingpresent equilibrium point.For the computation of the drinking model, the process involves setting the left-hand side of Model (1) equations equal to zero. (2) In the modeling of drinking behaviors, a drinking-free equilibrium point (E 0 ) represents a stable state with sustained absence of alcohol consumption, offering insights into effective prevention or mitigation measures.Understanding the dynamics around this equilibrium enhances comprehension of influential factors governing the spread of drinking behaviors.Analyzing these points is crucial for devising precise strategies in public health and intervention efforts.Therefore, the drinking-free equilibrium point (E 0 ) is derived by setting = 0, = 0, and = 0 in Model (1) as follows: An endemic equilibrium point (E ≠ 0. Upon solving the equations from Model (1), we arrive at the following outcomes:

Reproduction number
The basic reproduction number R 0 is a significant metric for assessing the transmission potential of heavy drinkers.It quantifies the total number of secondary heavy drinkers generated from the introduction of a single drinker individual into a susceptible population.In biological terms, it indicates the infectiousness and transmissibility of a pathogen.We utilize the next-generation matrix method [42,43] to calculate the reproduction number.The matrices F and V are employed to represent the recruitment of new drin- kers and the internal and external transmission terms associated with compartments of individuals engaged in drinking behavior, one obtains, The Jacobian of F and V matrix are given as The basic reproduction number R 0 is calculated as the spectral radius of the next-generation matrix − FV 1 , defined as follows: Hence, the basic reproduction number R 0 is obtained through the following derivation: , it indicates that the disease is unlikely to propagate extensively in the population and might eventually diminish.Conversely, if > R 1 0 , the disease can spread, potentially causing an outbreak.

Sensitivity analysis
Sensitivity analysis plays a crucial role in assessing the impact of input parameters on the dynamics of the epidemic drinking model.The parameters selected for sensitivity analysis in this study align closely with the model's objectives, which aim to understand the dynamics of drinking behavior and its transmission within a population.These carefully chosen parameters represent key factors influencing the spread of drinking behaviors, including the transmission rate of heavy drinkers, natural death rate, drinking-induced death rates, and the proportion of individuals entering treatment.The sensitivity indices used in the drinking model to calculate the reproduction number (R 0 ) are derived using the methodology proposed by Chitnis et al. [44].We specifically compute the normalized forward sensitivity indices of the parameters m in relation to R 0 .Let us delve into this analysis further The sensitivity estimation for the reproduction number concerning various parameters is presented as follows: Sensitivity analysis is commonly used to identify parameters significantly impacting the reproduction number R 0 .A positive value indicates a direct proportional relationship with the reproduction number, while a negative value indicates an inversely proportional relationship.In Table 2, the parameter (α) has a positive sign, indicating a directly proportional relationship with the reproduction number.Thus, an escalation in the reproduction number corresponds to an increase in parameters like the transmission rate of heavy drinkers (α).Reducing the parameter (α) can effectively control drinking in the population.Conversely, the reproduction number is negatively associated with parameters such as the natural death rate (μ), the proportion of drinkers entering treatment at the heavy drinker stage (ϕ), and the drinking-induced death rate of heavy drinker individuals (δ 1 ).Increasing either the natural death rate of the population (μ), the drinking-induced death rate of heavy drinker individuals (δ 1 ), or the proportion of drinkers entering treatment (ϕ) is neither ethical nor practical.Considering the sensitivity of R 0 to the drinking transmission rate of heavy drinkers (α), it is sensible to prioritize efforts on reducing this rate due to its direct impact on reproduction.In essence, this sensitivity analysis highlights the efficacy of prevention over cure.Enhancing preventive measures proves more effective in curbing the proliferation of habitual drinking compared to expanding access to treatment for individuals.
A notable revelation from this sensitivity analysis is the substantial influence of the transmission rate of heavy drinkers on the reproduction number (R 0 ).This discovery underscored the pivotal role that heavy drinkers have in propagating drinking behaviors within the population.Moreover, the sensitivity analysis underscored the significance of treatment interventions and their ability to lower heavy drinking rates, highlighting the importance of focused intervention strategies in addressing alcohol consumption.In simpler terms, this analysis suggests that prevention is more effective than trying to fix the issue later.Prioritizing prevention efforts is better for controlling the spread of drinking than solely focusing on increasing treatment numbers.The initial conditions for the model compartments, such as ( ) = 0 0.50, ( ) = 0 0.15, ( ) = 0 0.25, and ( ) = 0 0.1, as well as the values of compartmental parameters, are provided in Table 3 for numerical simulations and graphical representation of the model's dynamic behavior.
The effect of different significant model parameters on the value of R 0 is shown graphically in Figure 2. Figure 2(a) describes the reproduction number R 0 sensitivity versus drinking transmission rate (α) and drinking including death rate (δ 1 ).The influence of the drinking transmission rate (α) on the reproduction number is positive, whereas the impact of the drinking-related death rate (δ 1 ) is negative.Although the effect of the drinking-related death rate is almost negligible, the influence of the drinking transmission rate (α) is quite significant, so it increased the value of R 0 by almost 0.8.A similar interpretation can be seen in Figure 2(f) for ( ) α and ( ) μ on R 0 .Figure 2(b) depicts the sensitivity of the reproduction number (R 0 ) concerning the drinking-related death rate (δ 1 ) and the natural death rate (μ).Both parameters are decreased since they negatively impact the reproduction number (R 0 ).Therefore, both ( ) δ 1 and μ have a similar effect on R 0 , reaching a value close to 0.25.Figure 2(c) illustrates the impact of the parameters ( ) α and ( ) ϕ on R 0 .It is observed that R 0 increases significantly with these parameters, reaching a maximum value close to 0.8 or less than 1.This implies that the disease is unlikely to occur.Figure 2(d) illustrates how variations in the drinking-related death rate (δ 1 ) and the proportion of heavy drinkers in the treatment group (ϕ) influence the reproduction number (R 0 ).Both values decrease because they have a negative impact on R 0 .So, both ( ) δ 1 and ( ) ϕ have a similar effect on R 0 , making it approach a value near 0.5.Similarly, interpretation can be seen in Figure 2(e) for the proportional of heavy drinkers in the treatment compartment (α) and natural death rate (μ) on (R 0 ).

Analysis of stability
In this segment, we explore the local stability and global stability of the developed model at equilibrium points.The local stability is assessed by analyzing the eigenvalues of the Jacobian matrix at the equilibrium point.An equilibrium point is deemed locally stable if all eigenvalues have negative real parts.Global stability entails studying the system's behavior across its whole range, frequently necessitating progressive mathematical techniques such as Lyapunov analysis.

Local stability
A locally stable system tends to revert to its equilibrium state after experiencing minor disruptions, showcasing resilience and a tendency to maintain its original configuration near the equilibrium point.Conducting a local stability analysis within a drinking model aims to assess how minor disturbances can affect the persistence or reduction of drinking behavior.In the context of drinking models, the concept of local stability at equilibrium points becomes a crucial framework for understanding the complex dynamics of alcohol consumption and its spread within a population.In this discussion, we explore the locally stable drinking-free equilibrium point and the endemic equilibrium point in the devised model, utilizing the following theorem [29].for the purposed model is given as Evaluating the Jacobian matrix at the drinking-free equilibrium yields The resulting eigenvalues are  .λ λ , , ; otherwise, it is unstable.
Proof.The Jacobin matrix is computed as where Substituting C ij into equation ( 4), we obtain (5) The auxiliary equation of ( 5) can be calculated as follows: ) The aforementioned equation can be written as ) Applying the Routh-Hurwitz criterion to fourth-order polynomials [46], we ensure that . The non-positive eigenvalue of the auxiliary equation confirms the local asymptotic stability of E 1 as per the Routh-Hurwitz criterion.□

Global stability
Global stability in a drinking model investigates whether, under various conditions, the system converges to a stable state regarding drinking behavior across its entire range.Global stability analysis in the drinking model is essential for understanding the sustained, long-term dynamics of population drinking behavior.By examining the system's equilibrium conditions and considering the influence of broader societal factors, this analysis provides a comprehensive understanding of how the population's drinking patterns evolve over an extended period.Such insights are crucial for developing enduring public health strategies and policies tailored to address persistent trends and dynamics in alcohol consumption within the population.
Global stability analyses offer a holistic framework for effective long-term planning and management of drinking dynamics.In this discussion, we investigate the globally stable equilibrium points for both drinking-free and endemic states in the proposed model, employing the methodology outlined in the study by Anjam et al. [29].Proof.The global stability of the devised model at point E 0 is instigated through the construct of the Lyapunov function such as Upon computing the time derivative of equation ( 8) and applying it to the system of equations in (1), the following result is obtained: It is verifiable that if = 0 and = = 0.The Lyapunov LaSalle's invariant principle [48,49] indicates that E 0 achieves global asymptotic stabi- lity.□ Upon calculating the time derivatives of the aforementioned equation and utilizing Model (1), the resulting expression is as follows: As a consequence, ≤ 0

Results from numerical stability analysis
Runge-Kutta (RK) method is highly favored in epidemiological modeling because of its accuracy and robustness in solving the differential equations that govern disease dynamics.These models typically involve a system of coupled differential equations that describe interactions among different population compartments.The fourthorder RK method is particularly noteworthy for its high accuracy in approximating solutions to these equations, guaranteeing reliable and precise numerical results.Its robust convergence properties enhance both stability and efficiency, making it highly effective in solving diverse mathematical problems, especially those involving intricate and dynamic infectious disease dynamics.In this context, we employ the RK method to solve our proposed deterministic model, leveraging its ability to compute solutions through a series of approximations based on weighted function evaluations.Notably, the fourth-order RK method delivers accurate results with minimal model evaluations, making it a preferred choice for obtaining numerical solutions in Model (1).Hence, we obtain Drinking epidemic model with control interventions  11 To numerically solve the proposed model, we initially focus on the equation for the susceptible compartment in Model (12).This process results in Similarly, the remaining compartments of the posed model are solved, respectively, and give, 1 1 1

Algorithm
Step 1: Step 3: For the graphical depiction of the stability analysis, MATLAB software is utilized to run the aforementioned findings.The graphs presented in the figures below are generated using compartmental initial conditions for the devised model: ( ) = 0 0.50, ( ) = 0 0.15, ( ) = 0 0.25, and ( ) = 0 0.1, along with the parametric values specified in Table 3.These graphs serve to visually characterize the dynamics and stability characteristics of the model under varying initial conditions and parameter settings.Moreover, the time is taken between 0 and 30 units in terms of days with the initial population of susceptible individuals ( ) t , drinkers in treatment individuals ( ) t , heavy drinkers individuals ( ) t , and recovered individuals ( ) t .Figures 3-6 depict the stability dynamics of all population classes under various initial conditions.The susceptible class varies in magnitude with values of 0.5, 0.35, 0.2, and 0.05, showing a direct relationship among the size of the susceptible population and the time required to achieve stability (Figure 3).Similarly, Figure 4 illustrates stability achievements across different sizes of the heavy drinker class, with values of 0.25, 0.20, 0.15, and 0.10.The heavy drinker population decreases over time, stabilizing after 10 days.The population in the drinkers in treatment class is shown for various prevalence levels, represented by values such as 0.15, 0.13, 0.11, and 0.09 in Figure 5.The number of individuals in treatment for drinking initially rises over time, followed by a rapid decline as time progresses.Figure 6 illustrates the correlation between stability and the varying extent of the recovered population, considering values such as 0.1, 0.09, 0.08, and 0.07.The population experiences an initial rapid decrease, stabilizing with constant numbers after 15 days.

Optimal control theory
This section delves into the application of the drinking model with regard to optimal control strategies.Optimal control is utilized to intercede in the transmission of drinking within a specific community.The integration of an optimal control strategy, focusing on awareness campaigns ( ( ) u t

1
) and treatment ( ( ) u t 2), enhances the model's effectiveness in reducing heavy drinking rates by orchestrating a coordinated and targeted approach.Awareness campaigns play a pivotal role in educating the population about the risks and consequences of heavy drinking, promoting healthier behaviors, and reducing the social acceptability of excessive alcohol consumption.This can lead to a shift in social norms and attitudes towards drinking, thereby influencing individual behaviors and reducing heavy drinking rates over time.On the other hand, treatment interventions target individuals who are already engaged in heavy drinking behaviors.By providing access to effective treatment options, such as counseling, therapy, or medication, the model can help individuals reduce their alcohol consumption, manage cravings, and address underlying issues to their drinking habits.This approach not only aids in reducing heavy drinking rates but also promotes overall health and wellbeing among affected individuals.Furthermore, the combination of awareness campaigns and treatment interventions in the optimal control strategy allows for a comprehensive approach to addressing heavy drinking rates.Awareness campaigns create a supportive environment for individuals seeking help, while treatment interventions offer practical solutions for those struggling with heavy drinking.Together, these strategies work synergistically to reduce heavy drinking rates and promote healthier behaviors within the population.Moreover, the control variable ( ) u t 1 pertains to the implementation of a public awareness program targeting susceptible and treated individuals.This program aims to safeguard them against excessive alcohol consumption and facilitate successful avoidance.On the other hand, the ( ) u t 2 control variable represents the treatment applied to heavy drinker individuals.Therefore, heavy drinkers represent the proportion of individuals transitioning from regular alcohol consumption.Upon applying the ( ) u t 2 variable, they transition to the group of those permanently abstaining from alcohol.This control variable determines   the optimal values for decision variables, representing the adjustable action variable to achieve the desired result.By selecting the optimal value for the control variable, the best course of action can be determined.For emerging a control strategy, we bring into usage the optimal control theory [50][51][52].The adaptation of the control scheme interning susceptible drinkers class ( ) t , drinkers in treatment class ( ) t , heavy drinkers class ( ) t , and recovered class ( ) t is dis- cussed in next lines.The differential formulation of the model's equation (1), incorporating the control variable, is now expressed as follows: whereas the initial conditions are, Therefore, we consider an optimal control problem, wherein the objective function is specified as Here, is a control set with Lebesgue measurable onto [ ] 0, 1 and permissible range such as 1, 2 .Furthermore, Pontryagin's maximum principle [53,54] is used to find the optimal control strategy u that mini- mizes (or maximizes) the objective functional J by consid- ering the behavior of the Hamiltonian function H at each point in the state variables.Proof.A variety of techniques are considered to quantify the effectiveness of the control scheme [55,56], and consequently, nonnegative values are witnessed for control variables along with state variables.The convexity and the closeness of control variables describe the degree of solidity necessary to validate the control system.Moreover, it is noteworthy that the integrand exhibits a convex nature, representing an objective function in the form of ( ) Moreover, the presence of non-trivial vector functions like becomes apparent, particularly when ( ) y u * , * is regarded as the optimal procedure for the implemented control problem: The outcomes are achieved by the employment of the required condition to the Hamiltonian function.□

Existence of the optimal control problem
In this section, we establish the existence of optimal control by ensuring that all conditions of the control system are met at the initial time = t 0. To achieve this, we main- tain a bounded Lebesgue measurable control [56,57] with initial conditions and an upward bounded solution of the system.The problem of optimal control is tackled through an exploration of Lagrangian and Hamiltonian methods.The optimal control problem in Lagrangian formulation is represented by the following equation: over the region U is given by Proof.The application of optimal conditions provided outcomes such as ( ) ( ) The control variable is also solved such as One may note that the condition of control space remains expressible as under In simplification, the control variables are documented as follows: Drinking epidemic model with control interventions  15 The optimal system is now given as: The state variable and control variable in the optimal formation are achieved by the employment of an adjoint variable, the optimal system along with the initial conditions.The positivity of the second-order derivative of the objective function with respect to control variables is obvious.Therefore, the resultant Hamiltonian is documented as

Computational simulation and discussion
In this section, we aim to validate the proposed control strategy through numerical evaluation.We enrich the strategic framework by integrating a wide range of parametric settings.The numerical analysis commences with initializing the compartmental initial conditions from the study by Adu et al.The steeper rise of the susceptible population with respect to control is obvious in Figure 7. Similarly, the class of heavy drinkers is noted to be reduced more sharply with control as compared to without control in Figure 8.Moreover, the effectiveness of employed control can be seen for the compartment of drinkers in treatment from Figure 9. Finally, Figure 10 depicts the utility of the control scheme

Conclusion
This study is related to the development of a novel deterministic mathematical model elaborating the dynamics of the drinking epidemic model under the launch of an optimal control scheme.The investigation focuses on stratifying the population into four compartments: the susceptible class, heavy drinker class, drinkers in treatment class, and recovered drinker class.To ensure generalizability, a wide range of parametric settings are employed, encompassing the recruitment rate of the susceptible class, transmission rates from the susceptible class to the heavy drinker class, and from the recovered class to the susceptible class, natural death rate, drinking-induced death rates in the heavy drinker and treatment classes, the proportion of individuals entering the treatment class, and recovery rate.Throughout the investigation, a comprehensive analysis is provided, covering aspects such as positivity, boundedness, equilibrium points, basic reproduction number, stability, and sensitivity.The model ensures boundedness    and positivity by incorporating constraints and selecting appropriate functional forms.This precautionary approach prevents unrealistic growth of variables, aligning the model with real-world constraints and enhancing its reliability.The careful management of these aspects improves the model's applicability, offering policymakers an accurate representation that supports well-informed decision-making within practical limitations.The analysis of stability in the model indicates that when < R 1 0 , the system is both locally and globally asymptotically stable at the drinking-free equilibrium E 0 .Conversely, if > R 1 0 , the system exhibits an endemic equilibrium E 1 , at which point it is locally and globally asymptoti- cally stable.Furthermore, the sensitivity analysis to see the impact of different compartmental parameters on the reproduction number is discussed.The sensitivity analysis showed that the drinking transmission rates of heavy drinkers (α) are highly sensitive to R 0 , indicating their significant influence on the spread of drinking behaviors.This suggests that these parameters could be considered appropriate for inclusion as control variables in optimal control analysis.Furthermore, this finding underscores the critical role of heavy drinkers in driving the spread of drinking behaviors within the population.Additionally, the sensitivity analysis underscores the importance of treatment interventions and their potential to reduce heavy drinking rates, emphasizing the value of targeted intervention strategies in curbing alcohol consumption.
We initiate the optimal control process by integrating two time-dependent variables: a social awareness campaign and treatment, both acting as control variables.An optimal control problem is formulated by defining an objective function aimed at determining the optimal values for these control variables to minimize overall costs.By applying Pontryagin's maximum principle, we derive significant conditions for the optimal solution.Our research investigates the viability of two different optimal control methodologies.Numerical simulations conducted using MATLAB validate the effectiveness of the proposed control strategies.Graphical results reveal that the simultaneous use of both control variables is more effective in reducing the transmission flow of drinking behavior.These findings are extensively discussed, considering scenarios with and without optimal control.Through rigorous analytical procedures, we establish the utility of the control strategy in achieving optimal conditions across all four compartments.The susceptible population is estimated to show a gradual increase with control implementation, while a notable decrease in the proportions of heavy drinkers and individuals in treatment is observed with optimal control.Additionally, the recovery rate is significantly faster with the use of optimal control.Furthermore, our research advocates that simulations incorporating time-dependent controls are supplementary cost-effective equated to those utilizing time-independent controls.
In future endeavors, extending the model to include fractional-order dynamics alongside optimal control remains an intriguing avenue.This approach is expected to provide deeper insights into the dynamics of drinking behavior in society by leveraging available information more effectively.Additionally, there will be a strong emphasis on assessing the proposed scheme's applicability across diverse geographical regions, socio-economic strata, and varying healthcare access levels.Understanding how these strategies interact or counteract each other will help formulate more effective plans for disease control on a global scale.Moreover, future research may also focus on integrating real-time data and advanced machine-learning techniques to enhance the accuracy and predictive capabilities of the model.Author contributions: All authors have accepted responsibility for the entire content of this manuscript and have approved its submission.YNA and KI designed the concept, while KI and SAC carried out the methodology.MF and SM developed the software model code and conducted the simulations.YNA, SM, and NS prepared the manuscript with contributions from all coauthors.The authors applied the SDC approach to determine the sequence of authors.

Conflict of interest:
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this article.

Figure 1 :
Figure 1: The schematic diagram elucidates the dynamics of transmission in drinking behavior.

Theorem 4 . 1 .
The drinking-free equilibrium point E 0 is locally asymptotically stable with the condition <

Theorem 4 . 3 .
The drinking-free equilibrium point E is globally asymptomatically stable with the condition <

Theorem 4 . 4 . 1 0
The endemic equilibrium point E 1 of Model (1) exhibits global asymptotic stability when > R ; otherwise, it is unstable.Proof.The global stability is established by formulating the Lyapunov function at the endemic equilibrium points ( ) = *.Consequently, the endemic equilibrium E 1 is the only positively invariant set contained in [( As a result, the positive E 1 exhibits global asymptotic stability.□

Figure 3 :
Figure 3: Simulation results illustrate how changes in the susceptible class ( ) t manifest their impacts.

Figure 4 :
Figure 4: Simulation results depict how variations in the heavy drinker class ( ) t manifest their effects.

Figure 5 :
Figure 5: Simulation results depict how variations in the drinker in treatment class ( ) t manifest their effects.

Figure 6 :
Figure 6: Simulation results illustrate how changes in the recovered class ( ) t manifest their impacts.

1 . 6 .
[40] for our proposed model, comprising ( ) Furthermore, we ensure consistency by adopting literature-based parameter values as outlined in Table3.We address the optimization problem by employing the RK fourth-order technique and applying the transversality criterion within the time interval of[0,25].Additionally, we assume the weight constants to be Figures7-10illustrate the dynamic effectiveness of our model under both controlled and uncontrolled scenarios across all four compartments over time in days.The advantages of the optimal control strategy are evident across all population compartments.

Figure 7 :
Figure 7: Dynamics behavior of susceptible individual ( ) t with and without control are depicted in the graph.

Figure 8 : 1 and treatment for reducing the drinking population ( ) u t 2 ,
Figure 8: Dynamics behavior of heavy drinker individual ( ) t with and without control are depicted in the graph.

Figure 10 :
Figure 10: Dynamics behavior of recovered individual ( ) t with and without control are depicted in the graph.

Figure 11 :
Figure 11: Dynamics of control variable u 1 are depicted in the graphic.

Figure 12 :
Figure 12: Dynamics of control variable u 2 are depicted in the graphic.

Figure 9 :
Figure 9: Dynamics behavior of drinker in treatment individual ( ) t with and without control are depicted in the graph.

Funding information :
This research was supported by the Deanship of Research Development, Prince Mohammad Bin Fahd, University, Alkhubar 31952, Saudi Arabia.

Table 1 :
Explanation for each term in the equations of Model (1) is presented Terms of each equation Explanation b Recruitment rate of the susceptible people α The rate of transmission from the susceptible class to individuals with heavy drinking habits η Transmission rate of recovered people μThe rate of natural mortality among individuals in the susceptible population.

Table 2 :
Indices of sensitivity and parameters affecting the reproduction number R 0