Exact solution of generalized cooperative susceptible-infected-removed (SIR) dynamics

In this paper, we introduce a general framework for coinfection as cooperative susceptible-infected-removed (SIR) dynamics. We first solve the SIR model analytically for two symmetric cooperative contagions [L. Chen et al., Europhys. Lett. 104, 50001 (2013)10.1209/0295-5075/104/50001] and then generalize and solve the model exactly in the symmetric scenarios for three and more cooperative contagions. We calculate the transition points and order parameters, i.e., the total number of infected hosts. We show that the behavior of the system does not change qualitatively with the inclusion of more diseases. We also show analytically that there is a saddle-node-like bifurcation for two cooperative SIR dynamics and that the transition is hybrid. Moreover, we investigate where the symmetric solution is stable for initial fluctuations. We finally explore sets of parameters which give rise to asymmetric cases, namely, the asymmetric cases of primary and secondary infection rates of one pathogen with respect to another. This setting can lead to fewer infected hosts, a higher epidemic threshold, and also continuous transitions. These results open the road to a better understanding of disease ecology.