Effect of density dependence on coinfection dynamics

In this paper we develop a compartmental model of SIR type (the abbreviation refers to the number of Susceptible, Infected and Recovered people) that models the population dynamics of two diseases that can coinfect. We discuss how the underlying dynamics depends on the carrying capacity K: from a simple dynamics to a more complex. This can also help in understanding the appearance of more complicated dynamics, for example, chaos and periodic oscillations, for large values of K. It is also presented that pathogens can invade in population and their invasion depends on the carrying capacity K which shows that the progression of disease in population depends on carrying capacity. More specifically, we establish all possible scenarios (the so-called transition diagrams) describing an evolution of an (always unique) locally stable equilibrium state (with only non-negative compartments) for fixed fundamental parameters (density independent transmission and vital rates) as a function of the carrying capacity K. An important implication of our results is the following important observation. Note that one can regard the value of K as the natural ‘size’ (the capacity) of a habitat. From this point of view, an isolation of individuals (the strategy which showed its efficiency for COVID-19 in various countries) into smaller resp. larger groups can be modelled by smaller resp. bigger values of K. Then we conclude that the infection dynamics becomes more complex for larger groups, as it fairly maybe expected for values of the reproduction number R0≈1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0\approx 1$$\end{document}. We show even more, that for the values R0>1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0>1$$\end{document} there are several (in fact four different) distinguished scenarios where the infection complexity (the number of nonzero infected classes) arises with growing K. Our approach is based on a bifurcation analysis which allows to generalize considerably the previous Lotka-Volterra model considered previously in Ghersheen et al. (Math Meth Appl Sci 42(8), 2019).