Decoding how higher-order network interactions shape complex contagion dynamics

TL;DR

This paper develops a unified mean-field framework to analyze how higher-order interactions in complex networks influence contagion dynamics, revealing conditions for various bifurcation behaviors and multistability.

Contribution

It introduces a generalized model that characterizes bifurcation regimes in higher-order interaction models, including novel insights into multistability with three- and four-body interactions.

Findings

Three- and four-body interactions lead to different bifurcation regimes.

Multistability can occur with three-body interactions in single and multi-population models.

Two-population models exhibit multistability through symmetry breaking.

Abstract

Complex contagion models that involve contagion along higher-order structures, such as simplicial complexes and hypergraphs, yield new classes of mean-field models. Interestingly, the differential equations arising from many such models often exhibit a similar form, resulting in qualitatively comparable global bifurcation patterns. Motivated by this observation, we investigate a generalized mean-field-type model that provides a unified framework for analysing a range of different models. In particular, we derive analytical conditions for the emergence of different bifurcation regimes exhibited by three models of increasing complexity -- ranging from three- and four-body interactions to two connected populations with both pairwise and three-body interactions. For the first two cases, we give a complete characterisation of all possible outcomes, along with the corresponding conditions on network and epidemic parameters. In the third case, we demonstrate that multistability is possible despite only three-body interactions. Our results reveal that single population models with three-body interactions can only exhibit simple transcritical transitions or bistability, whereas with four-body interactions multistability with two distinct endemic steady states is possible. Surprisingly, the two-population model exhibits multistability via symmetry breaking despite three-body interactions only. Our work sheds light on the relationship between equation structure and model behaviour and makes the first step towards elucidating mechanisms by which different system behaviours arise, and how network and dynamic properties facilitate or hinder outcomes.