Long-Term Bifurcation and Stochastic Optimal Control of a Triple-Delayed Ebola Virus Model with Vaccination and Quarantine Strategies

: Despite its high mortality rate of approximately 90%, the Ebola virus disease (EVD) has not received enough attention in terms of in-depth research. This illness has been responsible for over 40 years of epidemics throughout Central Africa. However, during 2014–2015, the Ebola-driven epidemic in West Africa became, and remains, the deadliest to date. Thus, Ebola has been declared one of the major public health issues. This paper aims at exploring the effects of external ﬂuctuations on the prevalence of the Ebola virus. We begin by proposing a sophisticated biological system that takes into account vaccination and quarantine strategies as well as the effect of time lags. Due to some external perturbations, we extend our model to the probabilistic formulation with white noises. The perturbed model takes the form of a system of stochastic differential equations. Based on some non-standard analytical techniques, we demonstrate two main approach properties: intensity and elimination of Ebola virus. To better understand the impact of applied strategies, we deal with the stochastic control optimization approach by using some advanced theories. All of this theoretical arsenal has been numerically conﬁrmed by employing some real statistical data of Ebola virus. Finally, we mention that this work could be a rich basis for further investigations aimed at understanding the complexity of Ebola virus propagation at pathophysiological and mathematics levels.


Introduction
In 2013, EVD was declared an epidemic. Guinea, Liberia and Sierra Leone were the most affected regions of 2013-2016, where the virus caused huge losses and disrupted the socioeconomic balance [1]. On 8 August 2014, WHO declared the EVD pandemic in West Africa a Public Health Emergency of International Concern (PHEIC), which is designated only for events with a risk of potential international spread or that require a coordinated international response. Over the duration of the epidemic, EVD spread to seven more countries: Italy, Mali, Nigeria, Senegal, Spain, the United Kingdom and the United States [2,3]. Ebola viruses belong to the Filoviridae family like Marburg and Lloviu viruses [4,5]. Currently, five species of Ebola virus have been described: Zaire, Reston, Bundibugyo, Sudan and Taï Forest [6][7][8]. They are highly pathogenic viruses, classified as level 4. They are responsible, for the most part, for viral hemorrhagic fevers (VHF) most often fatal in human and non-human primates, during transmission by contact with biological fluids [9,10]. This virus, known since 1976, had already caused several episodes in Africa but none of such rapid evolution and such serious impact on public health [11]. Indeed, before the year 2014, there were 1580 deaths with all epidemics combined. The month of December On 8 August 2014, faced with the scale of the epidemic, the Director General of the WHO (World Health Organization) formally announced that this surge of the Ebola virus is an emergency at the international level, which is related to public health [17]. An international mobilization was then formed. Means of prevention have been put in place to prevent the spread of the virus from person to person [18]. Common practices such as funeral rites have been stopped. No treatment or vaccine was known at the time [19]. Faced with the virulence of the Ebola virus and in view of the increasing number of deaths, researchers have tried to find solutions to end the epidemic. Search procedures have been accelerated. Experimental treatments and vaccines have gradually emerged [20]. The census of cases in Europe has reinforced the concern of Westerners about the epidemic. European countries, and in particular France, have prepared for the possibility of a national epidemic [21].
The foundations of epidemiology based on compartmental models were laid by physician Sir Ronald Ross who, in 1911, wrote the first compartmental model of malaria using [22]. Mathematical models can therefore help to determine the appropriate response and the effort required to control an infectious disease. More specifically, they help the scientific community to easily understand the transmission processes of infectious diseases and to evaluate the effectiveness of the different control strategies implemented to contain an epidemic [23]. Several mathematical models have been proposed to explain the dynamics of infectious diseases, some of which have included control strategies such as quarantine, isolation, contact tracing, and vaccination [24][25][26][27][28]. Berge et al. [29] suggested a deterministic model for explaining the transmission mechanism of Ebola virus that included a commonality between the epidemics and endemic phases. The model included both the direct and indirect modes of interactions within the three distinct populations: animals, humans and fruit bats. Actually, the paper [29] is considered to be an extension of the study [23], which presented a model that includes both the indirect and direct spread of Ebola virus and a source that continuously provides the Ebola virus. Further, models have been proposed that attempt to comprehend different intervention policies aimed at limiting the transmission of Ebola virus [30]. Through mathematical models, refs. [31][32][33] studied the impact of vaccines on the dynamics of the epidemics, and [34] looked at the issue of quarantining. In another work, Berge et al. [35] proposed a system of equations describing the effectiveness of the tracing of interactions as a control measure in case of Ebola virus. Researchers have recently shown a strong interest in optimal control theory of epidemic problems particularly in controlling the Ebola virus. Numerous strategies and prevention programs have been used to investigate the optimal control analysis of problems that reflects the dynamics of Ebola virus [36,37]. Particularly, Area et al. [38] introduced an Ebola virus model to explain how vaccination of the vulnerable population would affect the spread of diseases. Subsequently, they studied two different control problems for understanding the transmission mechanism of Ebola virus and its control through vaccination.
In fact, when dealing with biological systems, more realistic properties should be included such as the time delay. This property is used to describe the evolution of certain phenomena, which does not only depend on current states but also on the past [39][40][41]. Specifically, in the Ebola situation, time lag is widely employed to understand the viral interactions. In this paper, we focus on the effect of vaccination plan, quarantine strategy and time delay on the long-run behavior of Ebola virus. From mathematical modeling, the Ebola dynamics can be modeled in the general form by a system of four ordinary differential equations: dI e = hI s I e − (u + α 1 + q + w)I e dt, dI q = qI e − (u + α 2 + z)I q dt, where I s (t) denotes the number of susceptible class, I e (t) is the number of infected people, I q (t) stands for the number of quarantined population and I r (t) describes the size of recovered class. The parameter A is the constant recruitment into I s class; h is the rate at which the Ebola virus is spreading in the population; v is the vaccination rate of the susceptible individuals and u represents the natural death rate, which is constant for all compartments. The parameter w denotes the recovery rate of individuals that remain infectious for a period of 1/w; the notations α 1 and α 2 are the respective disease-related death rates in the infected and quarantined populations. The rate at which the infected population becomes quarantined is denoted by q, whereas z is the rate at which the quarantine population is achieving recovery. For biological purposes, we shall assume that all of the parameters used in the model are constant and positive. In system (1), the positive quantity ς 1 stands for the length of immunity period of the vaccinated, and the term e −uς 1 represents the probability that the susceptible individuals are vaccinated at time t − ς 1 and still alive at time t. The positive quantities ς 2 and ς 3 reflect, respectively, the length of immunity period of the recovered from infected and quarantined individuals. e −uς 2 and e −uς 2 denote the probability that the infected and quarantined individuals have been immunized at times t − ς 2 and t − ς 3 , and the probability they have lost immunity but are still alive at time t.
In the complex and real world, environmental fluctuations influence the spread of infectious diseases and make it more complicated to foresee their behavior [42][43][44][45]. In such cases, deterministic systems, while able to make very informative forecasts and previsions, are not appropriate enough [46][47][48][49]. So, there is a pressing need for a developed mathematical model that can take into account the randomness effect, especially in the context of a harmful infectious disease such as the Ebola virus. In this context, many scholars have suggested and evolved a large number of stochastic models that describe the Ebola virus propagation dynamics from different angles and viewpoints [50]. In this research, we will use the proportional white noise approach to better reflect the reality of Ebola spreading. In an uncertain and constantly changing environment, the dynamic system can be disturbed by the white noises. As a result, the following system of SDEs (being a counterpart of model (1)) may be presented: dI e = hI s I e − (u + α 1 + q + w)I e dt + e I e (t)dP e (t), where s , e , q and r are the positive linear intensities associated, respectively, to the mutually independent Brownian motions P s (t), P e (t), P q (t) and P r (t). These latter, and all the random variables that will be met in this paper, are defined on a filtered probability space (Ω, E , {E t } t≥0 , P Ω ) endowed with a filtration {E t } t≥0 that verifies the usual hypotheses [51]. By employing the identical techniques presented in (ref. [46], Theorem 2), we can facilely demonstrate that for every positive started value, there is a single positive solution to the probabilistic system (2). This indicates that the model (2) is well-constructed scientifically.
Analogous to the non-probabilistic framework, the primary goal of examining the dynamics of the probabilistic epidemic systems is to establish the conditions that ensure disappearance and continuation of illness [43,[52][53][54][55]. In this research, we try to cope with the disappearance of the Ebola virus and the existence of a single stationary distribution of system (2). Of course, ergodic steady distribution implies the permanence of the Ebola infection. Technically, one of the usual methods to ensure ergodicity is the Lyapunov candidate function, which gives just sufficient conditions in the majority of cases [56][57][58]. Therefore, the first problematic of this research can be formulated as follows: is it feasible to supply the acute threshold (sufficient and necessary condition) for the stationarity of the system Ebola and the disappearance of the Ebola virus? Specifically, the present study suggests a novel method for dealing with biological systems perturbed by white noises. We offer the sufficient and almost required criterion for the ergodic property of our model and the disappearance of Ebola virus. Using a probabilistic auxiliary equation, we found the acute threshold S • . In other words, if S • > 1, then model (2) has a single steady ergodic distribution, and if S • < 1, then the density of infected class will quickly converge to zero.
On the other hand, optimization theory of dynamical systems, especially as it relates to epidemic illnesses, has received a lot of interest from researchers in recent years. Several control methods have been suggested and successfully implemented to various epidemic diseases. The primary goal of the theory of optimal control, particularly in epidemiological problems, is to present a preventive measure that restricts the spread of the disease and to portray the transmission mechanism of the infections [59][60][61]. This theory has filled many loopholes in different fields, notably those that depends on dynamic system(s), such as physics and business [62,63]. Interested readers are advised to consult [64] and the bibliography therein for the details of step-by-step optimality criteria. In general, when dealing with controlling of contagious diseases, the ideal approaches are to mitigate infectious disease to a minimum standard while spending the least amount of money on the overall control program. Regarding real life scenarios, a much lower number of articles are available in the literature, and usually the investigated data are retrieved from the laboratories and clinics. One must cite many works and make many attempts in various mathematical formulas to develop a model that suits the data points well. After developing a very basic system, it is possible to modify it by incorporating various features of the disease that may fit the clinical data as well as experimentally. Liu and Meng [65] worked on the dynamics of an optimal stochastic harvesting problem with delay and derived a few insightful results, including the maximum limit of sustainability. As Ebola virus contains both stochasticity and delay, we intend to study this disease while considering these facts in modeling, which has not been done yet to the best of our knowledge.
The remaining parts of the work are as follows. In Section 2, we treat the dynamical bifurcation of the stochastic system (2) by proving that S • is the sharp threshold between stationarity and extinction of the Ebola virus. In Section 3, we deal with the Ebola system (2) under stochastic control by using the well-known methods of delayed stochastic control theory. In Section 4, we present the parameter estimation approach for Ebola virus with real statistical data of Western Guinea. In Section 5, we support our findings with some computer simulations, and we derive the main conclusions of the article in Section 6.

Stochastic Long-Run Dynamics of Ebola Model
This section exhibits a new approach to belay the dynamical bifurcation of our perturbed system. Since the dynamic of the recovered compartment has no influence on the infection transmission behavior, we can neglect the last equation of (2). We consider the following reduced model: dI e = hI s I e − (u + α 1 + q + w)I e dt + e I e (t)dP e (t), with initial historical data I s (t), I e (t), To probe the long-run dynamics of our proposed Ebola model (3), we will use an additional perturbed equation associated with the differential equation of the total class T (t) = I s (t) + I e (t) + I q (t) + I r (t). This new equation illustrates the biological situation when Ebola does not exist I e (t) = 0. Let J (t) be a Markov process that verifies the following probabilistic formulation: where the probabilistic part is given by and T (0) = I s (0) + I e (0) + I q (0) + I r (0). From Lemma 2.11 of [66], we infer that Equation (4) is scientifically well-posed and the time average of its unique solution is estimated as follows: In line with the probabilistic comparison result [67], we show that T (t) ≤ J (t) a.s. So, asymptotically, we can use inequality (5) to obtain information about the mean long-time behavior of T (t). Concerning the latter, we also present the following result. Lemma 1. Let I s , I e , I q , I r be the unique positive solution of (2). Then, for any ϕ ∈ 0.5, 0.5 + u m , where m = max{ 2 s , 2 e , 2 q , 2 r }, we have the following estimate: For a detailed demonstration of the above lemma, we refer the reader to the proof of Lemma 2.2 in [66].
Before stating the main theorem of this section, we present in Table 1 some notations used later and define the following quantity:

Subset Definition
To exhibit possible future scenarios for the spread of Ebola, we offer the next principal result.
Theorem 1. The Ebola propagation process described by system (3) has two possible scenarios: 1.

2.
The eradication case (S • < 1)-that is, the Ebola epidemic-will disappear with full probability.
Proof. First case: when S • > 1. We apply Itô's rule to the function ln I e (t), then Once again, Itô's rule leads to Now, we merge the above two equations as follows: In line with the positivity of the solution, we obtain For simplicity of notation, we set An integration from 0 to t on both sides of (8) gives By employing the strong law of large numbers for martingale [51] and (5), we have lim inf By using the definitions of the subsets W a , W b and W c appearing in Table 1, we easily obtain lim inf where E {Ω,P} denotes the mathematical expectation. In accordance with (10), we have lim inf We choose then, (11) becomes lim inf Let ϕ ∈ 0.5, 0.5 + u m and ψ is given by ψ −1 + ϕ −1 = 1. By the use of Young inequality, we obtain lim inf where κ is a positive constant such that According to Lemma 1, we infer that lim inf From Table 1 and Markov's inequality, we can conclude that To go further, we suppose that x verifies that Then, we obtain lim sup From (12), we show that lim inf Identical to the demonstration of (Lemma 3.2, ref. [68]) and the mutually limited possibilities of lemma [69], we establish the existence, uniqueness and ergodicity of a single invariant distribution for the perturbed model (3).

So, lim
t→∞ I e (t) = 0 a.s. Alternatively, we can say that infection, whose dynamics is governed by model (3), could be eradicated with a rate of deterioration (u + α 1 + q + w) S • − 1 (this is the minimum benchmark value). This ends the demonstration.

Stochastic Optimal Control Strategies
In this section, we intend to obtain the optimal strategies for the proposed stochastic system based on comparing different preventive measures. To accomplish this, first of all in the following, we will identify the time-dependent control variables: This strategy aims to reduce the number of infectives by shifting them into the quarantine class.
Taking into account the abovementioned two control variables (i.e., u 1 and u 2 ) in the stochastic system (2), the desired control problem can be written in the form of dI e (t) = hI s (t)I e (t)(1 − u 1 (t)) − (u + α 1 + q + w + u 2 (t))I e (t) dt + e I e (t)dP e (t), dI q (t) = qI e (t) − (u + α 2 + z + u 2 (t))I q (t) dt + q I q (t)dP q (t), dI r (t) = vI s (t) + (w + u 2 (t))I e (t) along with the subsidiary conditions Let us proceed by introducing the vectors listed below and The vectors, which consist of the data at t = 0, are as follows: Θ 0 = I s (0), I e (0), I q (0), I r (0) .
Here, ψ and φ are vectors that comprise the following components: The desired strategy, given the control problem (14), is to minimize the cost (objective) function defined by here, A l , B m , K n are positive constants for l = 1, 2, 3; m = 1, 2; and n = 1, · · · , 4. Physically, the foremost three integrand' terms stand for the decreasing size of vulnerable, latent and infectious individuals; whereas, the related constants A 1 , A 2 and A 3 are used to balance the size of the underlying populations I s (t), I e (t) and I q (t), respectively. As usual, we assumed the quadratic terms 0.5B i u 2 , = 1, 2 in the cost functional for the preventive measures. To balance the control parameters, we utilized the B in the cost functional.
To conclude the above-stated discussion, our aim in the current research is to find the minimum value of the objective functional (19) by employing the strategies of reducing the number of vulnerable, latent and Ebola-infected people; maximizing the size of quarantine compartment; and doing all this while utilizing the least resources. In other words, the control vector (u 1 (t), u 2 (t)) ∈ U aims to decrease the interactions of the vulnerable and infected people and quarantining a portion of the latent and EVD-infected people in a unit time. Thus, we are searching for control variables, which has the property where U is the feasible set of controls, which is defined by here, the notions u ,max for = 1, 2 are constant and positive. To use the well-known stochastic version of maximum principle, first of all, it is necessary to express the Hamiltonian H m (x, u, y, z) as where ·, · represent the usual Euclidean inner-product while y = [y 1 , y 2 , y 3 , y 4 ] and z = [z 1 , z 2 , z 3 , z 4 ] vectors' components denote the associated adjoint variables. We may reproduce the stochastic minimum/maximum principle in the following form dΘ (t) = ∂H(Θ , u , y, z)(t + τ) ∂y dt + g(Θ (t))dP (t), where θ (t) describes the optimal trajectory of the state variable. The terminal and initial condition related to Equations (23) and (24) are respectively. Since Equation (25) indicates that control u (t) (optimal) is a function of y(t), z(t) and x (t), we can write where Φ can be obtained from Equation (25). Hence, Equations (23) and (24) may be written in another way as As a result, the corresponding Hamiltonian is of the form According to the stochastic maximum principle, By taking the respective derivative of H with respect to I s , I e , I q and I r , we will obtain y 1 , y 2 , y 3 , and y 4 in the form of where the notation Due to the free terminal states (i.e., I s (T ), I e (T ), I q (T ), I r (T ) ∈ R 4 + ), the transversality conditions in association with the auxiliary conditions, we have and h I s , I e , I q , where y 1 (T ) = −K 1 I s , y 2 (T ) = −K 2 I e , y 3 (T ) = −K 3 I q , y 4 (T ) = −K 4 I r . The next step is to consider the derivative of the function H with respect to u 1 ,u 2 , yielding the following expressions that denote the characterization of the control measures (optimal) u 1 and u 2 as To obtain the required control schemes, we developed an optimum objective functional, which contains both the control variables as well as the underlying states of the system. The goal of this study is to minimize the value of cost functional subject to the proposed stochastic differential equation system (2) and the included time-delay. It is observed that the whole control problem relies on the objective functional; thus, care should be taken in its selection. Weights should be assigned to the factors that are more important when balancing the coefficient, especially in the case of two or more factors in the cost functional. Before applying the well-known Pontryagin' principle [70], it is required to check the theory of existence for such control measures. To demonstrate such existence, one must demonstrate that within the feasible region, the control measures are bounded, by considering a sequence of control measures and their possible dependence on the states to demonstrate its convergence. By following [70], we may define the Hamiltonian associated to the problem in the following form H = integrand in the objective functional + adjoint variables × RHS of the differential equations.
To obtain the optimality criterion, we will differentiate the function H with respect to the control variables at the optimum values of the controls u * . Further, by considering the derivative of the function H with respect to the states, one can easily obtain the adjoin equations and similarly the transversality conditions as given by relation (34).

Parameters Estimation and Curve Fitting
Parameter estimation is perhaps the most important step in any research related to epidemics. It is considered an inverse process-that is, given a set of state variables and a model, we identify the values of parameters that produce a good fit of the data with our model. This step in any research is difficult as there is no single numerical or analytical method that can be used. If a method gives a unique solution of the problem, still we need a good approximation of the initial data for computing a better estimates of the parameters. This is due to solving the minimization problem by matching the actual data with the model's solution and for a good approximation, optimization solvers need initial data very close to the real data. There is a limited knowledge of the system, obtaining good parameter values for the type of biological system studied in this research can be difficult. Keeping in view the proximity and feasibility of the proposed problem, we will utilize the MATLAB routine lsqcurvefit and will perform the estimation of parameters as shown in Table 2. This technique is proffered over other available methods as this tool can handled very easily both small and large scale problems and a detailed explanations, the readers are advised to see [71]. To do so, in this research work, we considered the real Ebola cases reported in Western Guinea from the first forty weeks of the 2015 Ebola epidemic [18]. Let us express the proposed model (1) in the following functional differential equation: Here, F denotes the functional relationship between the independent variable, state variables and parameter θ. The strategy will be to find the value of θ such that it minimizes the difference between the reported data and the data predicted by the model. For n number of data points, actually, we want to minimize the following function Thus, the complete optimization problem is given by minθ, subject to Equation (37).
By assuming the estimated values of the parameters, we simulated both the stochastic and deterministic models against the reported Ebola virus data from the first forty weeks of 2015 in Figure 2. By an inspection, one can see that the stochastic system explains the data with good agreement. The expression 1 12 ∑ 12 k=1 x real k −x approximate k x real k ≈ 1.5685e − 01 was used to measure the average relative error of fitness, which actually describes the goodness of fit value. To further strengthen our findings, we calculated the relative error in small scale by using (1.5685e − 01), which suggests that this fitness is very accurate compared with any other fitting tools.  Table 2.

Stochastic Long-Run Behavior of Ebola Virus Model
This subsection is devoted to the numerical confirmation of the theoretical results of this article. We seek to verify that S • is the acute threshold of the model (3). Additionally, we probe the complex impact of independent white noises on the long-term behavior of Ebola virus. We mention that the parameters of the deterministic model are based on the real data provided in Table 2. For stochastic intensities, we select it theoretically in order to precisely audit the results obtained in the cases of stationarity and extinction of Ebola virus and to show that a dynamic bifurcation occurs at certain noise values. The solution of our model (3) is simulated in our case with the initial data: I s (0) = 700, I e (0) = 100, I q (0) = 300 and I r (0) = 20. Henceforth, the units adopted for time and number of individuals are, respectively, one day and one thousand population.
Firstly, we choose the parameter values presented in Table 2, and also select s = 0.05, e = 0.09, q = 0.04 and r = 0.03. By performing simple calculations, we obtain S • = 3.8451 > 1. From Theorem 1, we infer the existence of a steady distribution. By way of explanation, stationarity indicates the continuation of Ebola virus in the population over time, and the ergodic property established by Theorem 1 implies the persistence in the mean of the virus provided by the Ebola infection-that is, the time average of all classes of model (3) reach a strictly positive value, which can be seen in the Figure 3. Now, we increase the noise amplitude as follows: s = 0.8, e = 2.19, q = 0.44 and r = 0.3. In this case, we obtain S • = 0.9422 < 1. From Theorem 1, we deduce that Ebola infection will not be present in the population. Numerically, this is shown in Figure 4. Biologically, we can establish that a large amount of outside noise leads to the eradication of Ebola infection. This means that a long-term bifurcation occurs at the same critical noise values.    Table 2 and the random intensities as follows: s = 0.05, e = 0.09, q = 0.04 and r = 0.03. Now, we shall explore the effect of triple delays on the long-term behavior of Ebola. From Figures 5-7, we show that increasing the length of the time delay leads to a significant decrease in the number of affected classes. Biologically speaking, a long delay has the great-est effect on reducing transmission. Control and tracking strategies, maximizing tracking delays-for example, using media-based technology-and increasing the effectiveness of contact tracing, have the potential to prevent contact transmissions. All this makes the control of Ebola virus possible among the population.

Stochastic Optimal Control Strategies
We will demonstrate the numerical (approximate) solution to the assumed stochastic control problem; we will show how the assumed control measures work by presenting the numerical example. A stochastic numerical method will be used to simulate the control system. Parallel to that, we will solve numerically the adjoint equations by using the same technique and utilizing the conditions at final time. We will code the proposed model (1) in Step 1 without using control variables; then, we will find the numerical solution of the control system. For obtaining the approximate solution of the control system, we need to simulate the adjoint system (14) with the help of a backward procedure and by considering the conditions at final time T. Since the adjoint and state variables depend upon the control parameters, the convex relationships and characterization of the control measures will be utilized in the simulation. Whereas, during the entire simulation, we assumed the parameter values from Table 2. Based on these procedures and values of the parameters, we illustrate graphical results in Figure 8, which show the usefulness of the controls in reducing the spread of Ebola. Figure 8 graphically shows the curves after and before the execution of the control strategies by employing the deterministic model (1) of the system. As the current strategy focuses on the minimization of the individuals in infected and quarantined classes, the plot shows the outcomes of the present study. Besides, we emphasized the increase in size of the recovered compartment, which is well-supported in the same figure. Further, we show that without control, the population of the quarantined class declines and will reach zero within a finite time; however, the same population tends to increase at a faster rate as we implement the control policy. Finally, Figure 9 explains the dynamics of the control variable during the entire control program both in the presence and absence of the environmental noises.

Conclusions
The Ebola virus has received much attention from worldwide scholars. It is caused by the Ebola virus, which takes the specific lymphocytes in the bodily immune system as the prime base and destroys its mechanism. This epidemic has a short incubation period with a high fatality rate for patients, and its treatment increases the predicted time for its transmission. This treatment strategy guarantees a more reasonable probability of affecting others.
The current paper analyzed a generalized stochastic Ebola epidemic model that includes Brownian motions as well as stochastic time-varying delays. More explicitly, the treated system represents a four compartmental system that takes the form of an Itô stochastic differential equations system. After having presented the deterministic system, a mathematical analysis is carried out to offer an insight into Ebola propagation-notably, its long-run behavior. The main mathematical and biological results of this paper are arranged as follows:

1.
We have given the global threshold of the Ebola model (2) based on some dynamical properties of a one-block boundary system (4) perturbed by white noises.

2.
In Theorem 1, we have proved the main result related to the stationarity and extinction of the Ebola virus. It is worthy to mention that the analysis of these longtime properties is very significant for the underlying perturbed systems. Especially, in the case of epidemiological models, the ergodicity offers a general idea of the infection permanence.

3.
For controlling the rapid spread of the disease, we assumed two control parameters, and these were incorporated into the stochastic model. The stochastic model was analyzed with the help of Pontryagin principle, and the required optimality conditions were derived therein. To check the validity of the theoretical results and effectiveness of the control parameters, we plotted the models by simulating the same in MATLAB. Through simulations, the obtained results for persistence/extinction of the Ebola infection are verified. It was concluded that by using both the control variables, one can easily eliminate the Ebola from the community with minimum cost.
To control emerging tropical infectious diseases, we need to strengthen our intervention strategies. The overall public health goal is to reduce the burden of disease by reducing transmission or mitigating its severity. There are at least two basic public health guidelines for managing the spread of infectious diseases, such as Ebola, for which there are no vaccines or effective treatment. These is (i) effective isolation of symptomatic persons and (ii) asymptomatic contact tracing of groups of exposed persons and their quarantine for monitoring. In this regard, our paper highlights these control strategies in the context of stochasticity. In general, we pointed out that this study enhances many previous works regarding white noise perturbation. Furthermore, it offers new insights into understanding Ebola spread with complex real-world assumptions. In other words, the approach described in this article opens up many possibilities for future research.