A two-sex demographic model with single-dependent divorce rate

Divorce appears to be one of the least studied demographic processes, both empirically and in two-sex demographic models. In this paper, we study mathematical as well as biological implications of the assumption that the divorce rate is positively affected by the amount of single (i.e., unmarried/unpaired) individuals in the population. We do that by modifying the classical exponential two-sex model accounting for pair formation and separation. We model the divorce rate as an increasing function of the single population size and show that the single population pressure on the established couples alters the exponential behavior of the classical model in which the divorce rate is assumed constant. In particular, the total population size becomes bounded and a unique positive equilibrium exists. In addition, a Hopf bifurcation analysis around the positive equilibrium shows that the modified model may exhibit sustained oscillations.