OalibOpen Access Library Journal2333-97212333-9705Scientific Research Publishing10.4236/oalib.1115488Oalib-153106ArticleBiomedicalLife SciencesBusinessEconomicsChemistryMaterials ScienceComputer ScienceCommunicationsEarthEnvironmental SciencesEngineeringMedicineHealthcarePhysicsMathematicsSocial SciencesHumanitiesCatastrophe on Epidemic and Endemic Transitions in SEIQRV Dynamics0000-0002-7617-722XKakiMohammed Nokas Murad10009-0000-7154-6084AlshkaneHemn Omer Ahmed2 Department of Civil Engineering, College of Engineering, Al-Qalam University, Kirkuk, Iraq College of Pharmacy, Al-Qalam University, Kirkuk, Iraq
Methodology: In this paper, we aim to investigate endemic and epidemic transitions in the SEIQRV model using catastrophe theory. By incorporating nonlinear interactions among susceptible, exposed, infected, quarantined, recovered, and vaccinated populations, the model reveals that disease dynamics may undergo abrupt qualitative changes under small variations in epidemiological or control parameters. Using catastrophe-theoretic analysis, we show that epidemic onset, persistence, and elimination correspond to structural changes in the equilibrium landscape, giving rise to fold and cusp catastrophes. Results:These structures explain the coexistence of disease-free and endemic equilibria, hysteresis effects, and delayed epidemic collapse under gradual parameter adjustment. The results demonstrate that traditional threshold-based analysis may be insufficient to capture critical transitions in SEIQRV systems. Discussion:Catastrophe geometry provides a unified analytical and visual framework for identifying epidemic tipping points and for designing robust intervention strategies. This approach highlights the nonlinear nature of epidemic dynamics and the potential limitations of conventional modeling methods in predicting sudden outbreaks or disease persistence. Conclusion:The methodology enhances both theoretical understanding and practical prediction of complex epidemic behavior, offering valuable insights for public health planning and intervention strategies.
Infectious diseases continue to pose significant challenges to public health worldwide, necessitating rigorous mathematical frameworks to understand, predict, and control their spread [1]. Compartmental epidemic models, formulated as systems of nonlinear differential equations, provide a structured approach for describing disease progression in populations. The SEIQRV model extends the classical SEIR framework by incorporating quarantine Q and vaccination V compartments, enabling more realistic modeling of contemporary public-health interventions, including isolation, mass immunization, and behavioral modifications [2]. The SEIQRV model captures the interplay between biological disease processes and intervention strategies [3]. Susceptible individuals S may acquire infection through contact with infectious individuals I or exposed carriers E , while quarantined individuals Q are isolated to reduce transmission. Recovered individuals R may acquire temporary or permanent immunity, and vaccinated individuals V are directly moved from the susceptible pool according to vaccination policies. Advanced formulations demonstrate how vaccination and quarantine control influence epidemic outcomes and system equilibria [4]. Nonlinear incidence rates, quarantine efficacy, and vaccination coverage collectively determine the dynamical behavior of the system, producing complex outcomes such as epidemic outbreaks, endemic persistence, or abrupt transitions between disease states.
Matrix-based methods and eigenvalue analysis provide systematic tools for investigating these dynamics. The disease-free equilibrium DFE serves as a critical reference point: its stability determines whether small introductions of infection fade or escalate into epidemics. Linearization of the SEIQRV system at the DFE yields the Jacobian matrix, whose eigenvalues reveal local stability properties. The basic reproduction number R 0 , obtained as the spectral radius of the next-generation matrix, acts as a sharp epidemic threshold: R 0 < 1 ensures local stability of the DFE, whereas R 0 > 1 indicates potential epidemic invasion.
When R 0 > 1 , the system admits a positive endemic equilibrium, where infection persists at nonzero steady-state levels [5][6]. Eigenvalue analysis at this equilibrium provides insights into the persistence and robustness of endemic states under small perturbations. The presence of vaccination and quarantine introduces additional control parameters that modulate both epidemic thresholds and long-term dynamics, potentially leading to multistability, hysteresis, and abrupt regime shifts [7][8].
The SEIQRV framework also forms a natural foundation for nonlinear and catastrophe-theoretic analyses [9][10]. Fold, cusp, swallowtail, and butterfly catastrophes [11]-[14] can emerge when multiple interventions interact with inherent nonlinearities, explaining sudden epidemic outbreaks, delayed suppression, and extreme sensitivity to small changes in parameters. By combining matrix modeling, eigenvalue derivation, and nonlinear analysis, the SEIQRV model provides a comprehensive platform for understanding epidemic transitions, evaluating control strategies, and predicting public health outcomes under complex intervention scenarios [15][16].
2. Methodology
The analysis is conducted within a deterministic SEIQRV framework, where disease transmission and control processes are modeled by a system of nonlinear differential equations. The model incorporates key biological and public-health mechanisms, including exposure, infection, quarantine, recovery, and vaccination. Equilibrium states are derived analytically, and their local stability is examined using Jacobian eigenvalue analysis. Control parameters such as transmission rate, quarantine efficiency, and vaccination rate are treated as bifurcation parameters [17]-[20].
Catastrophe theory is applied by reducing the system near critical equilibria to low-dimensional normal forms. Fold and cusp catastrophes are identified by analyzing degeneracy conditions in which equilibrium branches merge or disappear. Numerical continuation and parameter sweeps are employed to construct bifurcation diagrams and catastrophe surfaces, revealing regions of multistability and hysteresis. These geometric structures are then interpreted in epidemiological terms, linking mathematical singularities to epidemic outbreaks, endemic persistence, and sudden regime shifts. The combined analytical and numerical approach provides insight into critical transitions that govern epidemic and endemic dynamics in SEIQRV systems [21]-[24].
3. Hierarchical Catastrophe Progression in SEIQRV Dynamics
The nonlinear structure of the SEIQRV epidemic model supports a natural hierarchy of catastrophic transitions as model complexity and control interactions increase. This hierarchy—comprising fold, cusp, swallowtail, and butterfly catastrophes—provides a unified framework for understanding epidemic onset, persistence, and abrupt regime shifts under multiple intervention strategies [25]-[27].
3.1. Fold Catastrophe: Epidemic Threshold
The fold catastrophe represents the simplest transition in SEIQRV dynamics and corresponds to the classical epidemic threshold. As a single control parameter, typically the effective transmission rate, crosses a critical value, a stable disease-free equilibrium collides with an unstable equilibrium and disappears. This mechanism explains sudden epidemic outbreaks or abrupt disease elimination and aligns with basic reproduction number analysis. However, the fold catastrophe captures only single-threshold behavior and cannot account for coexistence or memory effects [28]-[31].
3.2. Fold Catastrophe in the SEIQRV Epidemic Model
Consider a deterministic SEIQRV model of the form:
d S d t = Λ − β S ( I + η E ) − ν S + ω R − μ S , d E d t = β S ( I + η E ) − ( σ + μ ) E , d I d t = σ E − ( γ + δ + μ ) I , d Q d t = δ I − ( γ q + μ ) Q , d R d t = γ I + γ q Q − ( ω + μ ) R , d V d t = ν S − μ V , where all parameters have their standard epidemiological meanings.
3.3. Disease-Free Equilibrium and Threshold Parameter
The disease-free equilibrium (DFE) is given by E 0 = ( S 0 , 0 , 0 , 0 , 0 , V 0 ) , S 0 = Λ μ + ν , V 0 = ν μ S 0 .
Linearization at the DFE yields a Jacobian matrix whose dominant eigenvalue determines local stability. Using the next-generation matrix approach, the basic reproduction number is
R0=βS0σ(σ+μ)(γ+δ+μ).3.4. Fold Catastrophe and Endemic Equilibria
When R 0 < 1 , the DFE is locally asymptotically stable, and the infection dies out. At the critical threshold R 0 = 1 the Jacobian has a simple zero eigenvalue, indicating loss of normal hyperbolicity. For R 0 > 1 , the system admits a positive endemic equilibrium E * = ( S * , E * , I * , Q * , R * , V * ) , I * > 0 .
Near the critical threshold, the dynamics on the center manifold reduce to the fold normal form
dudt=μ−u2
where u is a reduced measure of the prevalence of infection.
3.5. Eigenvalue Structure at the Fold Point
At the fold point R 0 = 1 , the Jacobian matrix J ( E 0 ) satisfies λ 1 = 0 , Re ( λ i ) < 0 ( i = 2 , ⋯ , 6 ) .
The zero eigenvalue corresponds to the infection subspace, while all remaining eigenvalues are strictly negative, ensuring a one-dimensional center manifold. As R 0 increases past unity, the eigenvalue becomes positive, destabilizing the DFE and forcing the system onto the endemic branch.
3.6. Epidemiological Interpretation
In the SEIQRV framework, the fold catastrophe formalizes the classical epidemic threshold while embedding it in a geometric context. The stable disease-free equilibrium and the emerging endemic equilibrium collide at R 0 = 1 , producing an abrupt transition. This explains why gradual changes in transmission, quarantine efficiency, or vaccination coverage may yield sudden outbreaks or rapid disease elimination. However, because the fold catastrophe involves only a single control parameter, it cannot account for bistability or hysteresis, which arise when vaccination and quarantine interact to produce higher-codimension catastrophes.
Note:In the SEIQRV model, the epidemic threshold R 0 = 1 corresponds to a fold catastrophe in which a simple zero eigenvalue of the Jacobian induces an abrupt transition between disease-free and endemic equilibria.
The Canonical Cusp Catastrophe:Consider the potential function V ( x ) = x 4 4 + a 2 x 2 + b x , where: x is the state variable, a and b are control parameters. The equilibrium points satisfy d V d x = x 3 + a x + b = 0 .
Thus, equilibria are the roots of the cubic equation x 3 + a x + b = 0 [9]-[14].
Number of Equilibria:The discriminant is Δ = − 4 a 3 − 27 b 2 .
If Δ < 0 , there is one real equilibrium.If Δ = 0 , two equilibria merge at a fold bifurcation.If Δ > 0 , there are three real equilibria.
The cusp boundary is therefore 4 a 3 + 27 b 2 = 0 .
Stability of Equilibria:For the gradient system x ˙ = − ( x 3 + a x + b ) , an equilibrium x * is locally stable when d d x ( x 3 + a x + b ) | x = x * = 3 ( x * ) 2 + a > 0 .
Inside the cusp region ( Δ > 0 ): the lower branch is stable, the middle branch is unstable, and the upper branch is stable. Hence, the system exhibits bistability.
What Causes a Catastrophe?
As a control parameter changes continuously, a stable equilibrium branch may reach a fold point where it disappears.
At that instant:
1) The current stable state no longer exists.
2) The system cannot remain nearby.
3) It jumps suddenly to another stable branch.
This abrupt transition is called a catastrophe.
Hysteresis:
Suppose a < 0 is fixed and b is slowly increased.
The system follows one stable branch.At the upper fold, that branch vanishes.The state jumps to the other stable branch.
If b is then decreased, the jump back occurs at a different fold point.
Thus, the forward and backward paths differ, producing hysteresis.
When two control parameters interact—such as transmission intensity and vaccination coverage—the SEIQRV system exhibits a cusp catastrophe. In this regime, disease-free and endemic equilibria may coexist, producing bistability and hysteresis. Epidemic emergence and collapse follow different paths depending on initial conditions and parameter history. This explains delayed epidemic suppression and resurgence despite gradual control improvements [31][32].
The swallowtail catastrophe arises from higher-order nonlinearities associated with multiple feedback mechanisms, including asymptomatic transmission and quarantine saturation. It introduces layered equilibrium branches and sequential bifurcations, allowing the SEIQRV system to transition through multiple epidemic phases. These dynamics account for prolonged outbreaks, temporary stabilization, and sudden re-escalation under modest parameter variations [8][11][21].
3.9. Butterfly Catastrophe: Extreme Sensitivity and Cascading Regime Shifts
At the highest level of complexity, the butterfly catastrophe governs SEIQRV dynamics when several control mechanisms—vaccination, quarantine efficiency, transmission heterogeneity, and behavioral response—act simultaneously. Governed by quintic nonlinearities, this catastrophe generates highly intricate equilibrium landscapes with multiple coexisting endemic states and nested bifurcation structures. Small perturbations can trigger cascading transitions, leading to explosive outbreaks or abrupt epidemic collapse. Strong hysteresis and path dependence dominate, emphasizing the difficulty of reversing high-prevalence endemic states once established [21][24][25].
3.10. Unified Interpretation
Together, these catastrophes form a progressive structure in which increasing model realism leads to richer and more realistic epidemic behavior. The fold captures basic epidemic thresholds, the cusp explains bistability and hysteresis, the swallowtail describes multi-stage epidemic evolution, and the butterfly reveals extreme sensitivity and cascading transitions. This hierarchy generalizes Thom’s catastrophe classification to SEIQRV epidemiological systems and provides a powerful analytical framework for identifying tipping points (Tipping points are critical thresholds in a dynamical system at which a small, gradual change in parameters or external conditions causes a sudden, qualitative, and often irreversible shift in the system’s behavior or state), designing early-warning indicators, and developing robust public-health interventions [33]-[36]. In epidemiological systems, tipping points mark transitions between disease-free, epidemic, and endemic states. For example, a slight increase in the effective transmission rate or a small reduction in vaccination coverage may push the system past a tipping point, leading to a rapid outbreak. Conversely, crossing a control-related tipping point may result in sudden disease elimination. These transitions are often accompanied by hysteresis and multistability, meaning that reversing the parameter change does not immediately restore the previous state.
4. SEIQR Epidemic Model
Let
S ( t ) : susceptible individuals. E ( t ) : exposed (infected but not yet infectious). I ( t ) : infectious individuals. Q ( t ) : quarantined infectious individuals. R ( t ) : recovered individuals.
The SEIQR dynamics are governed by the following nonlinear system:
Λ : recruitment rate into the susceptible population. β : effective transmission rate. σ : progression rate from exposed to infectious. δ : quarantine rate of infectious individuals. γ , γ Q : recovery rates from I and Q , respectively. μ : natural removal rate (birth/death, not disease-induced).
Key Properties
Nonlinearity arises from the bilinear incidence term β S I .The model reduces to SEIR when δ = 0 .The quarantine class Q modifies infectious duration and epidemic thresholds.The system preserves positivity and admits a biologically feasible region.
Disease-Free Equilibrium (DFE) and Epidemic Eigenvalues
Consider a generic compartmental epidemic model with compartments (susceptible S , exposed E , infected I , and other disease states that give rise to new infections). Using the next-generation matrix method, we can systematically derive the epidemic eigenvalues at the disease-free equilibrium (DFE). Essentially, one partitions the model into infected compartments and the rest. Let x = ( E , I , ⋯ ) represent the vector of infected states. Then the system of differential equations can be written in the form d x d t = F ( x ) − V ( x ) , where F ( x ) represents new infection terms and V ( x ) collects other transitions (progression, recovery, removal) without new infections. The next-generation matrix is constructed by linearizing these terms around the DFE.
Define the Jacobian matrices at the DFE:
F=[∂Fi∂xj]DFE,V=[∂Vi∂xj]DFE.
The next-generation matrix is then K = F V − 1 .
The epidemic eigenvalues are the eigenvalues of K . The basic reproduction number R 0 is the spectral radius (dominant eigenvalue) of K . If R 0 < 1 , all epidemic eigenvalues have negative real parts and the DFE is locally asymptotically stable, meaning the disease dies out. If R 0 > 1 , at least one eigenvalue is positive, indicating a potential epidemic outbreak.
For a simple SEIR model, for example, the new infection rate into the exposed class is β S I / N and other transitions include progression σ and recovery. Linearizing around ( S * , E * , I * ) = ( N , 0 , 0 ) yields a next-generation matrix with entries reflecting the expected number of secondary infections produced by one infectious individual during its infectious period. The resulting eigenvalues determine whether perturbations in the infected classes grow or decay—in other words, whether the system is in an epidemic regime.
5. Standard Derivation for the SEIQR Model
Below, I give a complete, standard derivation for the SEIQR model, including:
i) Disease-Free Equilibrium (DFE)
ii) Jacobian matrix at the DFE
iii) DFE (epidemic) eigenvalues and R 0
iv) Endemic equilibrium structure
v) Endemic eigenvalue interpretation
We will now start with the foregoing subsection titles:
1)Disease-Free Equilibrium (DFE)
At the disease-free state, E = I = Q = R = 0 . From d S d t = 0 , we obtain S 0 = Λ μ .
This matrix is block-triangular, so eigenvalues can be read systematically.
3)DFE (Epidemic) Eigenvalues and R 0
Trivial Eigenvalues
The following eigenvalues are always negative: λ 1 = − μ , λ 2 = − ( γ Q + μ ) , λ 3 = − μ .
These correspond to S , Q , R and do not affect epidemic onset.
Epidemic Subsystem ( E , I )
The epidemic dynamics are governed by the 2 × 2 block:
(−(σ+μ)βS0σ−(δ+γ+μ)).
The characteristic equation is λ 2 + a 1 λ + a 0 = 0 , where
a1=(σ+μ)+(δ+γ+μ),a0=(σ+μ)(δ+γ+μ)−βS0σ.
Basic Reproduction Number
Define
R0=βS0σ(σ+μ)(δ+γ+μ)=βΛσμ(σ+μ)(δ+γ+μ).
If R 0 < 1 : all eigenvalues have negative real parts → DFE is locally asymptotically stable.If R 0 > 1 : one eigenvalue becomes positive → epidemic outbreak occurs.
This is the epidemic eigenvalue condition.
4)Endemic Equilibrium
When R 0 > 1 , a unique endemic equilibrium E * = ( S * , E * , I * , Q * , R * ) exists with all components positive.
The infected level I * is proportional to R 0 − 1 , confirming a forward (transcritical) bifurcation at R 0 = 1 .
5)Endemic Eigenvalue Structure
The Jacobian evaluated at E * yields:
One eigenvalue associated with population balance (negative).Remaining eigenvalues determined by the infected-quarantine subsystem.
Key Properties
For R 0 > 1 , all eigenvalues have negative real parts → endemic equilibrium is locally asymptotically stable.No oscillatory instability arises unless additional delays or nonlinear feedback are introduced.The transition from DFE to the endemic state occurs via a transcritical bifurcation.
6)Interpretation
DFE eigenvalues determine epidemic invasion.Endemic eigenvalues determine persistence and robustness of infection.Quarantine reduces R 0 through δ , delaying or preventing outbreaks.This eigenvalue structure aligns naturally with matrix-based epidemic modeling, as emphasized in [37].
6. Endemic Dynamics and Catastrophe Analysis in SEIQRV Systems
In this section, we examine the long-term behavior of the SEIQRV system by analyzing the endemic equilibrium states and their associated eigenvalue structures, which provide insight into the local stability of the disease within the population. The analysis reveals how small perturbations around these equilibria can evolve, highlighting conditions under which the disease persists or dies out. Furthermore, by applying catastrophe theory, we explore the nonlinear transitions between epidemic and endemic states, illustrating how gradual changes in key parameters—such as transmission rate, recovery rate, and vaccination coverage—can lead to abrupt shifts in disease prevalence. These insights are crucial for understanding the thresholds at which control interventions may fail or succeed, providing a mathematical framework to guide effective public health strategies in mitigating sudden outbreaks [37]-[40].
6.1. Endemic Equilibrium and Endemic Eigenvalue Structure
An endemic equilibriumexists when the infected compartments attain positive steady-state values. Its stability can be examined by evaluating the full Jacobian matrix of the system at that equilibrium. In contrast to theDFE analysis, at the endemic state, one must account for all compartments contributing to infection, recovery, and other demographic processes. After substituting equilibrium values into the Jacobian and computing its eigenvalues, one typically finds that if R 0 > 1 , the endemic equilibrium is locally asymptotically stable, meaning perturbations from this steady state decay over time, and the disease persists [41].
In practice, the combination of matrix methods (like next-generation matrices) and classical stability analysis offers a robust framework for characterizing the conditions under which an epidemic either fades out or settles into an endemic state. This approach aligns with the matrix modeling [37] perspective that advocates making eigenvalue analysis central to both understanding and controlling infectious disease dynamics.
6.2. Catastrophe on Epidemic and Endemic Transitions in SEIQRV Dynamics
The SEIQRV epidemic model provides a comprehensive framework for describing infectious disease transmission by incorporating susceptible S , exposed E , infected I , quarantined Q , recovered R , and vaccinated V populations. Beyond classical threshold analysis, the dynamics of SEIQRV systems are inherently nonlinear and may exhibit catastrophictransitions, where small changes in epidemiological or control parameters produce sudden and qualitative shifts in disease behavior. These transitions cannot be fully explained by linear stability or smooth bifurcation theory alone [3][42][43].
Within the catastrophe-theoretic perspective, epidemic onset and elimination correspond to structural changes in the equilibrium landscape of the SEIQRV system. Variations in transmission rate, quarantine effectiveness, vaccination coverage, or immunity loss can move the system across critical bifurcation manifolds, triggering abrupt transitions between disease-free, epidemic, and endemic states. In particular, fold and cusp catastrophesnaturally arise in SEIQRV dynamics, capturing the emergence or disappearance of equilibria and the coexistence of multiple stable states. Such structures explain why epidemics may persist even when parameters suggest control, or why sudden outbreaks may occur despite gradual parameter variation [44]-[48]. Catastrophe geometry further reveals the presence of hysteresis and multistability in epidemic-endemic transitions. Once an epidemic state is established, reversing control parameters does not necessarily restore the disease-free equilibrium along the same path, leading to delayed epidemic collapse or rebound outbreaks. This behavior is especially pronounced when quarantine compliance or vaccination uptake changes slowly, while the underlying nonlinear interactions amplify small perturbations into large epidemiological consequences.
From a public health perspective, catastrophe analysis provides a powerful interpretive and predictive tool. By identifying tipping points and bifurcation surfaces in SEIQRV dynamics, it becomes possible to anticipate critical transitions and design intervention strategies that avoid catastrophic regime shifts. Thus, catastrophe theory not only enriches the mathematical understanding of epidemic and endemic dynamics but also offers actionable insight for managing infectious diseases under uncertainty and complex control interactions [37][38].
7. Key Findings
1)Hierarchical Catastrophe Dynamics
The SEIQRV model exhibits a hierarchy of catastrophic transitions—fold, cusp, swallowtail, and butterfly—corresponding to increasing nonlinear complexity and control interactions.
2)Epidemic Thresholds and R 0
The basic reproduction number determines whether the disease-free equilibrium is stable or unstable. DFE stability ( R 0 < 1 ) ensures disease extinction, while ( R 0 > 1 ) leads to epidemic invasion and the potential emergence of endemic states.
3)Endemic Multistability
The system admits multiple stable equilibria under certain parameter combinations, explaining hysteresis effects, delayed epidemic suppression, and potential resurgence after initial decline.
4)Impact of Vaccination and Quarantine
Control measures influence both epidemic eigenvalues and endemic equilibrium structure. Even small variations in vaccination coverage or quarantine efficiency can trigger abrupt transitions, particularly near critical bifurcation surfaces.
5)Nonlinear and Catastrophe-Theoretic Insights
Higher-order nonlinearities generate complex equilibrium landscapes. Swallowtail and butterfly catastrophes illustrate cascading bifurcations and extreme sensitivity to control parameters, highlighting the importance of coordinated interventions.
Practical Public-Health Implications
Understanding epidemic tipping points and potential abrupt transitions supports the design of robust intervention strategies, early-warning indicators, and policies to prevent catastrophic outbreaks.
8. Conclusions
The SEIQRV model provides a comprehensive framework for analyzing the spread and control of infectious diseases in populations where both quarantine and vaccination interventions are implemented. By incorporating susceptible, exposed, infected, quarantined, recovered, and vaccinated compartments, the model captures essential biological processes and public-health interventions that shape epidemic trajectories. Matrix-based analysis and eigenvalue derivation at the disease-free equilibrium (DFE) allow for the explicit computation of the basic reproduction number ( R 0 ) and establish rigorous thresholds for epidemic invasion. When ( R 0 < 1 ), the DFE is locally stable, ensuring disease extinction, whereas ( R 0 > 1 ) leads to potential outbreaks and the emergence of a positive endemic equilibrium.
The analysis of endemic equilibria reveals the conditions under which infection persists and highlights the influence of quarantine and vaccination on long-term disease dynamics. Eigenvalue structures at endemic states demonstrate local stability and the system’s sensitivity to parameter variations, emphasizing that even small changes in control measures can substantially alter epidemic outcomes. Moreover, the nonlinear interactions inherent in the SEIQRV framework may give rise to multistability, hysteresis, and abrupt transitions, suggesting that gradual policy adjustments can produce disproportionate epidemiological consequences.
Overall, the SEIQRV model unifies classical epidemic threshold theory with modern nonlinear and catastrophe perspectives, offering both qualitative and quantitative insights into epidemic and endemic transitions. This framework provides a robust analytical basis for designing intervention strategies, anticipating tipping points, and understanding the complex interplay between infection dynamics, quarantine, and vaccination. By combining matrix modeling, eigenvalue analysis, and nonlinear considerations, the SEIQRV approach enhances predictive capability and supports evidence-based public health planning in the control of infectious diseases.
The SEIQRV model provides a unified framework for understanding both epidemic onset and endemic persistence in populations subject to quarantine and vaccination interventions. Matrix-based analysis and eigenvalue derivation at the disease-free equilibrium yield explicit thresholds for epidemic invasion, while endemic equilibrium analysis demonstrates long-term infection persistence and stability. The inclusion of multiple control parameters and nonlinear interactions reveals complex phenomena, including multistability, hysteresis, and abrupt regime shifts, which are naturally interpreted through catastrophe theory.
By integrating fold, cusp, swallowtail, and butterfly catastrophes into the analysis, this framework captures a hierarchy of epidemic behaviors from simple outbreak thresholds to extreme sensitivity and cascading transitions. The results emphasize that marginal changes in vaccination or quarantine strategies can lead to disproportionately large epidemiological consequences. Overall, SEIQRV modeling, enhanced by matrix and catastrophe-theoretic methods, provides both qualitative insight and quantitative predictive capability, supporting robust public-health decision-making in controlling infectious diseases.
9. Future Work
Future work on SEIQRV dynamics could explore several directions to enhance the understanding and applicability of the model. One promising avenue is the inclusion of age-structured or spatially heterogeneous populations, which would allow the assessment of localized outbreaks and targeted interventions. Incorporating stochastic effects and real-time data assimilation could improve predictions under uncertainty and capture random epidemic fluctuations. Finally, integrating behavioral, social, and environmental factors into the model could provide a more comprehensive framework for designing adaptive and robust public health strategies.
ReferencesXu, W., Liu, H. and Qin, C. (2025) Dynamics of Nonlinear Stochastic SEIR Infectious Disease Model with Isolation and Latency Period. Symmetry, 17, Article 155. https://doi.org/10.3390/sym17020155 10.3390/sym17020155https://doi.org/10.3390/sym17020155Xu, W.Liu, H.Qin, C.2025Dynamics of Nonlinear Stochastic SEIR Infectious Disease Model with Isolation and Latency PeriodSymmetry1715510.3390/sym17020155Kiseleva, O., Yakovlev, S., Prytomanova, O. and Kuzenkov, O. (2025) Mathematical Modeling of Regional Infectious Disease Dynamics Based on Extended Compartmental Models. Computation, 13, Article 187. https://doi.org/10.3390/computation13080187 10.3390/computation13080187https://doi.org/10.3390/computation13080187Kiseleva, O.Yakovlev, S.Prytomanova, O.Kuzenkov, O.2025Mathematical Modeling of Regional Infectious Disease Dynamics Based on Extended Compartmental ModelsComputation1318710.3390/computation13080187Li, W. and Wu, K. (2025) Dynamical Behavior of a Stochastic SEIQRV Infectious Model with an Ornstein-Uhlenbeck Process and General Incidence. Frontiers in Applied Mathematics and Statistics, 11, Article 1687991. https://doi.org/10.3389/fams.2025.1687991 10.3389/fams.2025.1687991https://doi.org/10.3389/fams.2025.1687991Li, W.Wu, K.2025Dynamical Behavior of a Stochastic SEIQRV Infectious Model with an Ornstein-Uhlenbeck Process and General IncidenceFrontiers in Applied Mathematics and Statistics11168799110.3389/fams.2025.1687991Gupta, A. and Katarya, R. (2023) A Deep-SIQRV Epidemic Model for COVID-19 to Access the Impact of Prevention and Control Measures. Computational Biology and Chemistry, 107, Article ID: 107941. https://doi.org/10.1016/j.compbiolchem.2023.107941 10.1016/j.compbiolchem.2023.10794137625364https://doi.org/10.1016/j.compbiolchem.2023.107941Gupta, A.Katarya, R.2023A Deep-SIQRV Epidemic Model for COVID-19 to Access the Impact of Prevention and Control MeasuresComputational Biology and Chemistry107107941ID10.1016/j.compbiolchem.2023.10794137625364Bavithra, S.A.R., Chatzarakis, G.E. and Padmasekaran, S. (2025) Stability Analysis of Delay SEIQRV Mathematical Model for COVID-19. International Journal of Ap- pliedMathematics, 55, 2025-2032. https://www.iaeng.org/IJAM/issues_v55/issue_7/IJAM_55_7_14.pdfBavithra, S.A.R.Chatzarakis, G.E.Padmasekaran, S.2025Stability Analysis of Delay SEIQRV Mathematical Model for COVID-19International Journal of Ap-plied Mathematics55Ahmed, I., Tariboon, J., Muhammad, M. and Ibrahim, M.J. (2024) A Mathematical and Sensitivity Analysis of an HIV/AIDS Infection Model. International Journal of Mathematics and Computer in Engineering, 3, 35-46. https://doi.org/10.2478/ijmce-2025-0004 10.2478/ijmce-2025-0004https://doi.org/10.2478/ijmce-2025-0004Ahmed, I.Tariboon, J.Muhammad, M.Ibrahim, M.J.2024A Mathematical and Sensitivity Analysis of an HIV/AIDS Infection ModelInternational Journal of Mathematics and Computer in Engineering310.2478/ijmce-2025-0004Pathak, S. and Kota, V.R. (2025) Mathematical Analysis of an Infectious Disease Model: Stability, Persistence, and Equilibrium Analysis in a Non-Linear Epidemiological System. Journal of Nonlinear Mathematical Physics, 32, Article No. 64. https://doi.org/10.1007/s44198-025-00333-1 10.1007/s44198-025-00333-1https://doi.org/10.1007/s44198-025-00333-1Pathak, S.Kota, V.R.Stability, P2025Mathematical Analysis of an Infectious Disease Model: Stability, Persistence, and Equilibrium Analysis in a Non-Linear Epidemiological SystemJournal of Nonlinear Mathematical Physics32No10.1007/s44198-025-00333-1Bliman, P. and Fang, M. (2025) Reinfection Induced Multistability in an Epidemic Model. Nonlinear Science, 5, Article ID: 100057. https://doi.org/10.1016/j.nls.2025.100057 10.1016/j.nls.2025.100057https://doi.org/10.1016/j.nls.2025.100057Bliman, P.Fang, M.2025Reinfection Induced Multistability in an Epidemic ModelNonlinear Science5100057ID10.1016/j.nls.2025.100057Murad, M.N. (2011) On the Cusp Catastrophe Model and Stability. General Mathematics Notes, 2, 73-82.Murad, M.N.2011On the Cusp Catastrophe Model and StabilityGeneral Mathematics Notes2Kaki, M.N.M. (2012) Mathematical Catastrophe with Applications. General Mathematics Notes, 11, 35-46.Kaki, M.N.M.2012Mathematical Catastrophe with ApplicationsGeneral Mathematics Notes11Kaki, M.N.M. (2012) On the Catastrophic Model and Stability. International Journal of Basic & Applied Sciences, 12, 64-67.Kaki, M.N.M.2012On the Catastrophic Model and StabilityInternational Journal of Basic & Applied Sciences12Kaki, M.N.M. (2014) Catastrophic Types Depending on Degree of Non-Linearity. Pure and Applied Mathematics Journal, 3, 24-27.Kaki, M.N.M.2014Catastrophic Types Depending on Degree of Non-LinearityPure and Applied Mathematics Journal3Kaki, M.N.M. and Aziz, S.A. (2013) Stability and Existence of Periodic Solutions in Non-Linear Differential Equations. International Journal of Emerging Technology and Advanced Engineering, 3, 574-577.Kaki, M.N.M.Aziz, S.A.2013Stability and Existence of Periodic Solutions in Non-Linear Differential EquationsInternational Journal of Emerging Technology and Advanced Engineering3Kaki, M.N.M. (1985) Treatment of Phenomena of Instability by Methods of Catastrophe Theory. Master’s Thesis, University of Baghdad.Kaki, M.N.M.Thesis, U1985Treatment of Phenomena of Instability by Methods of Catastrophe TheoryMaster’s ThesisGuastello, S.J. (2025) Catastrophe Modelling for Time Series of Reported Cases of COVID-19: Workload Effects in the Health Care System. Nonlinear Dynamics, Psychology, and Life Sciences, 29, 453-472. https://pubmed.ncbi.nlm.nih.gov/41117151/Guastello, S.J.Dynamics, P2025Catastrophe Modelling for Time Series of Reported Cases of COVID-19: Workload Effects in the Health Care SystemNonlinear Dynamics29Montes-Olivas, S., Kucharski, A.J., Gravenor, M.B. and Frost, S.D.W. (2026) Exploring Epidemic Control Policies Using Nonlinear Programming and Mathematical Models. PLOS Computational Biology, 22, e1014238. https://doi.org/10.1371/journal.pcbi.1014238 10.1371/journal.pcbi.1014238https://doi.org/10.1371/journal.pcbi.1014238Montes-Olivas, S.Kucharski, A.J.Gravenor, M.B.Frost, S.D.W.2026Exploring Epidemic Control Policies Using Nonlinear Programming and Mathematical ModelsPLOS Computational Biology2210.1371/journal.pcbi.1014238Izadi, M. and Waezizadeh, T. (2024) Stability Analysis and Numerical Evaluations of a COVID-19 Model with Vaccination. BMC Medical Research Methodology, 24, Article No. 97. https://doi.org/10.1186/s12874-024-02209-2 10.1186/s12874-024-02209-238678207https://doi.org/10.1186/s12874-024-02209-2Izadi, M.Waezizadeh, T.2024Stability Analysis and Numerical Evaluations of a COVID-19 Model with VaccinationBMC Medical Research Methodology24No10.1186/s12874-024-02209-238678207Chen, Y. and Zhao, W. (2020) Asymptotic Behavior and Threshold of a Stochastic SIQS Epidemic Model with Vertical Transmission and Beddington-DeAngelis Incidence. Advances in Difference Equations, 2020, Article 353. https://doi.org/10.1186/s13662-020-02815-6 10.1186/s13662-020-02815-6https://doi.org/10.1186/s13662-020-02815-6Chen, Y.Zhao, W.2020Asymptotic Behavior and Threshold of a Stochastic SIQS Epidemic Model with Vertical Transmission and Beddington-DeAngelis IncidenceAdvances in Difference Equations202035310.1186/s13662-020-02815-6Zinihi, A., Ehrhardt, M. and Sidi Ammi, M.R. (2026) Spatiotemporal SEIQR Epidemic Modeling with Optimal Control for Vaccination, Treatment, and Social Measures. Journal of Mathematical Biology, 93, Article No. 10. https://doi.org/10.1007/s00285-026-02430-1 10.1007/s00285-026-02430-1https://doi.org/10.1007/s00285-026-02430-1Zinihi, A.Ehrhardt, M.Ammi, M.R.Vaccination, T2026Spatiotemporal SEIQR Epidemic Modeling with Optimal Control for Vaccination, Treatment, and Social MeasuresJournal of Mathematical Biology93No10.1007/s00285-026-02430-1Zhou, W., Zhao, T., Wang, A. and Tang, S. (2024) Bifurcations and Dynamics of a Filippov Epidemic Model with Nonlinear Threshold Control Policy and Medical-Resource Constraints. Chaos, Solitons & Fractals, 184, Article ID: 114992. https://doi.org/10.1016/j.chaos.2024.114992 10.1016/j.chaos.2024.114992https://doi.org/10.1016/j.chaos.2024.114992Zhou, W.Zhao, T.Wang, A.Tang, S.Chaos, S2024Bifurcations and Dynamics of a Filippov Epidemic Model with Nonlinear Threshold Control Policy and Medical-Resource ConstraintsChaos184114992ID10.1016/j.chaos.2024.114992Wang, X.L. and Wu, K.L. (2025) Bifurcation Analysis for a SIRS Model with a Nonlinear Incidence Rate. arXiv: 2510.27292.Wang, X.L.Wu, K.L.2025Bifurcation Analysis for a SIRS Model with a Nonlinear Incidence Rate2510Scarabel, F., Coldwell, H. and Cassidy, T. (2025) Genuine and Spurious Bistability in a Simple Epidemic Model with Waning Immunity. arXiv: 2511.09196.Scarabel, F.Coldwell, H.Cassidy, T.2025Genuine and Spurious Bistability in a Simple Epidemic Model with Waning Immunity2511Öztürk, Z., Bilgil, H. and Sorgun, S. (2023) Application of Fractional SIQRV Model for SARS-CoV-2 and Stability Analysis. Symmetry, 15, Article 1048. https://doi.org/10.3390/sym15051048 10.3390/sym15051048https://doi.org/10.3390/sym15051048Bilgil, H.Sorgun, S.2023Application of Fractional SIQRV Model for SARS-CoV-2 and Stability AnalysisSymmetry15104810.3390/sym15051048Scarabel, F., Polner, M., Wylde, D., Barbarossa, M.V. and Röst, G. (2025) Bistability and Complex Bifurcation Diagrams Generated by Waning and Boosting of Immunity. Journal of Mathematical Biology, 91, Article No. 30. https://doi.org/10.1007/s00285-025-02264-3 10.1007/s00285-025-02264-340844643https://doi.org/10.1007/s00285-025-02264-3Scarabel, F.Polner, M.Wylde, D.Barbarossa, M.V.2025Bistability and Complex Bifurcation Diagrams Generated by Waning and Boosting of ImmunityJournal of Mathematical Biology91No10.1007/s00285-025-02264-340844643Drake, J.M., Brett, T.S., Chen, S., Epureanu, B.I., Ferrari, M.J., Marty, É., et al. (2019) The Statistics of Epidemic Transitions. PLOS Computational Biology, 15, e1006917. https://doi.org/10.1371/journal.pcbi.1006917 10.1371/journal.pcbi.100691731067217https://doi.org/10.1371/journal.pcbi.1006917Drake, J.M.Brett, T.S.Chen, S.Epureanu, B.I.Ferrari, M.J.2019The Statistics of Epidemic TransitionsPLOS Computational Biology1510.1371/journal.pcbi.100691731067217Chakraborty, A.K., Gao, S., Miry, R., Ramazi, P., Greiner, R., Lewis, M.A., et al. (2024) An Early Warning Indicator Trained on Stochastic Disease-Spreading Models with Different Noises. Journal of the Royal Society Interface, 21, Article ID: 20240199. https://doi.org/10.1098/rsif.2024.0199 10.1098/rsif.2024.019939118548https://doi.org/10.1098/rsif.2024.0199Chakraborty, A.K.Gao, S.Miry, R.Ramazi, P.Greiner, R.Lewis, M.A.2024An Early Warning Indicator Trained on Stochastic Disease-Spreading Models with Different NoisesJournal of the Royal Society Interface21202401ID10.1098/rsif.2024.019939118548Wang, X. and Wu, K. (2025) Bifurcation Analysis for a SIRS Model with a Nonlinear Incidence Rate. arXiv:2510.27292.Wang, X.Wu, K.2025Bifurcation Analysis for a SIRS Model with a Nonlinear Incidence Rate2510Öztürk, Z., Bilgil, H. and Sorgun, S. (2023) Fractional SIQRV Model for COVID-19 and Numerical Solutions. Konuralp Journal of Mathematics, 11, 131-140.Bilgil, H.Sorgun, S.2023Fractional SIQRV Model for COVID-19 and Numerical SolutionsKonuralp Journal of Mathematics11Li, W., Huang, L. and Wang, J. (2020) Global Dynamics of Filippov-Type Plant Disease Models with an Interaction Ratio Threshold. Mathematical Methods in the Applied Sciences, 43, 6995-7008. https://doi.org/10.1002/mma.6450 10.1002/mma.6450https://doi.org/10.1002/mma.6450Li, W.Huang, L.Wang, J.2020Global Dynamics of Filippov-Type Plant Disease Models with an Interaction Ratio ThresholdMathematical Methods in the Applied Sciences4310.1002/mma.6450Sabbar, Y. (2024) Exploring Threshold Dynamics of a Behavioral Epidemic Model Featuring Two Susceptible Classes and Second-Order Jump-Diffusion. Chaos, Solitons & Fractals, 186, Article ID: 115216. https://doi.org/10.1016/j.chaos.2024.115216 10.1016/j.chaos.2024.115216https://doi.org/10.1016/j.chaos.2024.115216Sabbar, Y.Chaos, S2024Exploring Threshold Dynamics of a Behavioral Epidemic Model Featuring Two Susceptible Classes and Second-Order Jump-DiffusionChaos186115216ID10.1016/j.chaos.2024.115216Opoku-Sarkodie, R., Bartha, F.A., Polner, M. and Röst, G. (2024) Bifurcation Analysis of Waning-Boosting Epidemiological Models with Repeat Infections and Varying Immunity Periods. Mathematics and Computers in Simulation, 218, 624-643. https://doi.org/10.1016/j.matcom.2023.12.006 10.1016/j.matcom.2023.12.006https://doi.org/10.1016/j.matcom.2023.12.006Opoku-Sarkodie, R.Bartha, F.A.Polner, M.2024Bifurcation Analysis of Waning-Boosting Epidemiological Models with Repeat Infections and Varying Immunity PeriodsMathematics and Computers in Simulation21810.1016/j.matcom.2023.12.006Al-Jebouri, M.M., Kaki, M.M. and Mutlek, T. (2026) Quantitative Assessment of Wound Healing of Mice Treated with Combinations of Laser Irradiations and Antibiotics via Experimental Mathematical Modelling. World Journal of Pharmacy and Pharmaceutical Sciences, 15, 1165-1177.Al-Jebouri, M.M.Kaki, M.M.Mutlek, T.2026Quantitative Assessment of Wound Healing of Mice Treated with Combinations of Laser Irradiations and Antibiotics via Experimental Mathematical ModellingWorld Journal of Pharmacy and Pharmaceutical Sciences15Bar, T., Banerjee, A., Casals, B., Catalan, G. and Rodríguez-Viejo, J. (2025) Disorder-Aided Early Warning Signals: Predicting Catastrophic Shifts in Athermal Systems. Physical Review Research, 7, Article ID: 043228. https://doi.org/10.1103/t84q-m6zh 10.1103/t84q-m6zhhttps://doi.org/10.1103/t84q-m6zhBar, T.Banerjee, A.Casals, B.Catalan, G.Viejo, J.2025Disorder-Aided Early Warning Signals: Predicting Catastrophic Shifts in Athermal SystemsPhysical Review Research7043228ID10.1103/t84q-m6zhMiyauchi, Y., Ikeda, M. and Kawahara, Y. (2026) Generalized Stochastic Resilience for Early Warning Signals Based on Koopman Operator. Nonlinear Dynamics, 114, Article No. 246. https://doi.org/10.1007/s11071-025-12072-5 10.1007/s11071-025-12072-5https://doi.org/10.1007/s11071-025-12072-5Miyauchi, Y.Ikeda, M.Kawahara, Y.2026Generalized Stochastic Resilience for Early Warning Signals Based on Koopman OperatorNonlinear Dynamics114No10.1007/s11071-025-12072-5Sanders, P.A. and Bastiaansen, R. (2026) Early Warning of Critical Transitions: Distinguishing Tipping Points from Turing Destabilizations. Physica D: Nonlinear Phenomena, 495, Article ID: 135300. https://doi.org/10.1016/j.physd.2026.135300 10.1016/j.physd.2026.135300https://doi.org/10.1016/j.physd.2026.135300Sanders, P.A.Bastiaansen, R.2026Early Warning of Critical Transitions: Distinguishing Tipping Points from Turing DestabilizationsPhysica D: Nonlinear Phenomena495135300ID10.1016/j.physd.2026.135300Scheffer, M., Carpenter, S.R., Foley, J.A., Folke, C. and Walker, B. (2009) Critical Transitions in Nature and Society. Princeton University Press.Scheffer, M.Carpenter, S.R.Foley, J.A.Folke, C.Walker, B.2009Critical Transitions in Nature and SocietyAl-Jebouri, M.M. and Kaki, M.N.M. (2025) Application of Matrices Modelling for Infectious Diseases of Humans. Open Journal of Applied Sciences, 15, 2733-2758. https://doi.org/10.4236/ojapps.2025.159184 10.4236/ojapps.2025.159184https://doi.org/10.4236/ojapps.2025.159184Al-Jebouri, M.M.Kaki, M.N.M.2025Application of Matrices Modelling for Infectious Diseases of HumansOpen Journal of Applied Sciences1510.4236/ojapps.2025.159184Al-Jebouri, M.M. and Kaki, M.N.M. (2025) Dynamics and Phase Portraits of the SEIQR Model. Annals of Public Health & Epidemiology, 3, 1-8.Al-Jebouri, M.M.Kaki, M.N.M.2025Dynamics and Phase Portraits of the SEIQR ModelAnnals of Public Health & Epidemiology3Lee, H. and Li, V.C. (2025) Phase Boundaries and Critical Transitions in Coupled Epidemic-Behavioral Systems. Bulletin of Mathematical Biology, 87, Article No. 179. https://doi.org/10.1007/s11538-025-01546-7 10.1007/s11538-025-01546-7https://doi.org/10.1007/s11538-025-01546-7Lee, H.Li, V.C.2025Phase Boundaries and Critical Transitions in Coupled Epidemic-Behavioral SystemsBulletin of Mathematical Biology87No10.1007/s11538-025-01546-7Gul, N., Shah, I., Ahmad, S., Ullah, I. and De la Sen, M. (2025) Stability and Optimal Control Analysis of an SEIQR Epidemic Model with Saturated Incidence Rate. European Journal of Pure and Applied Mathematics, 18, Article No. 6596. https://doi.org/10.29020/nybg.ejpam.v18i4.6596 10.29020/nybg.ejpam.v18i4.6596https://doi.org/10.29020/nybg.ejpam.v18i4.6596Gul, N.Shah, I.Ahmad, S.Ullah, I.Sen, M.2025Stability and Optimal Control Analysis of an SEIQR Epidemic Model with Saturated Incidence RateEuropean Journal of Pure and Applied Mathematics18No10.29020/nybg.ejpam.v18i4.6596Earn, D.J.D. and McCluskey, C.C. (2025) Global Stability of Epidemic Models with Uniform Susceptibility. Proceedings of the National Academy of Sciences of the United States of America, 122, e2510156122. https://doi.org/10.1073/pnas.2510156122 10.1073/pnas.2510156122https://doi.org/10.1073/pnas.2510156122Earn, D.J.D.McCluskey, C.C.2025Global Stability of Epidemic Models with Uniform SusceptibilityProceedings of the National Academy of Sciences of the United States of America12210.1073/pnas.2510156122Xing, Y. and Li, H. (2025) Global Dynamics of a Reaction-Diffusion SEIVQR Epidemic Model in Almost Periodic Environments. Journal of Applied Analysis & Computation, 15, 762-785. https://doi.org/10.11948/20240121 10.11948/20240121https://doi.org/10.11948/20240121Xing, Y.Li, H.2025Global Dynamics of a Reaction-Diffusion SEIVQR Epidemic Model in Almost Periodic EnvironmentsJournal of Applied Analysis & Computation1510.11948/20240121Jeelani, M.B., Alnahdi, A.S., Ud Din, R., et al. (2025) Mathematical Model to Investigate Transmission Dynamics of COVID-19 with Vaccinated Class. AIMS Mathematics, 8, 29932-29955.Jeelani, M.B.Alnahdi, A.S.Din, R.2025Mathematical Model to Investigate Transmission Dynamics of COVID-19 with Vaccinated ClassAIMS Mathematics8Jing, S., Xue, L. and Yang, J. (2024) Backward Bifurcation Arising from Decline of Immunity against Emerging Infectious Diseases. Applied Mathematics Letters, 158, Article ID: 109241. https://doi.org/10.1016/j.aml.2024.109241 10.1016/j.aml.2024.109241https://doi.org/10.1016/j.aml.2024.109241Jing, S.Xue, L.Yang, J.2024Backward Bifurcation Arising from Decline of Immunity against Emerging Infectious DiseasesApplied Mathematics Letters158109241ID10.1016/j.aml.2024.109241Carlsson, R.-M., Childs, L.M., Feng, Z., Glasser, J.W., Heffernan, J.M., Li, J. and Röst, G. (2020) Modeling the Waning and Boosting of Immunity from Infection or Vaccination. Journal of Theoretical Biology, 497, Article 110265. https://doi.org/10.1016/j.jtbi.2020.110265 10.1016/j.jtbi.2020.11026532272134https://doi.org/10.1016/j.jtbi.2020.110265Carlsson, R.Childs, L.M.Feng, Z.Glasser, J.W.Heffernan, J.M.Li, J.2020Modeling the Waning and Boosting of Immunity from Infection or VaccinationJournal of Theoretical Biology49711026510.1016/j.jtbi.2020.11026532272134Faeq, I.R. and Abdalrahman, S.O. (2023) The Stability and Catastrophic Behavior of Finite Periodic Solutions in Non-Linear Differential Equations. Tikrit Journal of Pure Science, 28, 146-152. https://doi.org/10.25130/tjps.v28i6.1382 10.25130/tjps.v28i6.1382https://doi.org/10.25130/tjps.v28i6.1382Faeq, I.R.Abdalrahman, S.O.2023The Stability and Catastrophic Behavior of Finite Periodic Solutions in Non-Linear Differential EquationsTikrit Journal of Pure Science2810.25130/tjps.v28i6.1382Kiseleva, O., Yakovlev, S., Chumachenko, D. and Kuzenkov, O. (2024) Exploring Bifurcation in the Compartmental Mathematical Model of COVID-19 Transmission. Computation, 12, Article 186. https://doi.org/10.3390/computation12090186 10.3390/computation12090186https://doi.org/10.3390/computation12090186Kiseleva, O.Yakovlev, S.Chumachenko, D.Kuzenkov, O.2024Exploring Bifurcation in the Compartmental Mathematical Model of COVID-19 TransmissionComputation1218610.3390/computation12090186Al-Shbeil, I., Djenina, N., Jaradat, A., Al-Husban, A., Ouannas, A. and Grassi, G. (2023) A New COVID-19 Pandemic Model Including the Compartment of Vaccinated Individuals: Global Stability of the Disease-Free Fixed Point. Mathematics, 11, Article 576. https://doi.org/10.3390/math11030576 10.3390/math11030576https://doi.org/10.3390/math11030576Al-Shbeil, I.Djenina, N.Jaradat, A.Al-Husban, A.Ouannas, A.Grassi, G.2023A New COVID-19 Pandemic Model Including the Compartment of Vaccinated Individuals: Global Stability of the Disease-Free Fixed PointMathematics1157610.3390/math11030576