Abstract
It has been reported that COVID-19 patients had an increased neutrophil count and a decreased lymphocyte count in the severe phase and neutrophils may contribute to organ damage and mortality. In this paper, we present the bifurcation analysis of a dynamical model for the initial innate system response to pulmonary infection. The model describes the interaction between a pathogen and neutrophilis (also known as polymorphonuclear leukocytes). It is shown that the system undergoes a sequence of bifurcations including subcritical and supercritical Bogdanov–Takens bifurcations, Hopf bifurcation, and degenerate Hopf bifurcation as the parameters vary, and the model exhibits rich dynamics such as the existence of multiple coexistent periodic oscillations, homoclinic orbits, bistability and tristability, etc. Numerical simulations are presented to explain the theoretical results.
