Global asymptotic stability, extinction and ergodic stationary distribution in a stochastic model for dual variants of SARS-CoV-2.

Several mathematical models have been developed to investigate the dynamics SARS-CoV-2 and its different variants. Most of the multi-strain SARS-CoV-2 models do not capture an important and more realistic feature of such models known as randomness. As the dynamical behavior of most epidemics, especially SARS-CoV-2, is unarguably influenced by several random factors, it is appropriate to consider a stochastic vaccination co-infection model for two strains of SARS-CoV-2. In this work, a new stochastic model for two variants of SARS-CoV-2 is presented. The conditions of existence and the uniqueness of a unique global solution of the stochastic model are derived. Constructing an appropriate Lyapunov function, the conditions for the stochastic system to fluctuate around endemic equilibrium of the deterministic system are derived. Stationary distribution and ergodicity for the new co-infection model are also studied. Numerical simulations are carried out to validate theoretical results. It is observed that when the white noise intensities are larger than certain thresholds and the associated stochastic reproduction numbers are less than unity, both strains die out and go into extinction with unit probability. More-over, it is observed that, for weak white noise intensities, the solution of the stochastic system fluctuates around the endemic equilibrium (EE) of the deterministic model. Frequency distributions are also studied to show random fluctuations due to stochastic white noise intensities. The results presented herein also reveal the impact of vaccination in reducing the co-circulation of SARS-CoV-2 variants within a given population.