Arxiv:2109.01099V1 [Physics.Soc-Ph] 2 Sep 2021
Total Page:16
File Type:pdf, Size:1020Kb
Hypergraph dynamics: assortativity and the expansion eigenvalue Nicholas W. Landry∗ and Juan G. Restrepoy Department of Applied Mathematics, University of Colorado at Boulder, Boulder, Colorado 80309, USA (Dated: September 3, 2021) The largest eigenvalue of the matrix describing a network's contact structure is often important in predicting the behavior of dynamical processes. We extend this notion to hypergraphs and motivate the importance of an analogous eigenvalue, the expansion eigenvalue, for hypergraph dynamical processes. Using a mean-field approach, we derive an approximation to the expansion eigenvalue and its associated eigenvector in terms of the degree sequence for uncorrelated hypergraphs. We introduce a generative model for hypergraphs that includes degree assortativity, and use a per- turbation approach to derive an approximation to the expansion eigenvalue and its corresponding eigenvector for assortative hypergraphs. We validate our results with both synthetic and empirical datasets. We define the dynamical assortativity, a dynamically sensible definition of assortativity for uniform hypergraphs, and describe how reducing the dynamical assortativity of hypergraphs through preferential rewiring can extinguish epidemics. Complex social systems often exhibit assortative mix- hypergraph assortative mixing patterns to other dynam- ing [1, 2], where individuals with similar characteristics ical processes. connect with each other more often than it would be ex- Our approach is as follows: first, we define and moti- pected if they were connected at random. Assortativity vate the importance of the expansion eigenvalue on dy- has been extensively studied in network science [3] and namical processes; second, we derive a mean-field approx- found to have significant effects on synchronization [4], imation of this eigenvalue and its corresponding eigen- epidemic dynamics [5, 6], stability [7], game theory [8], vector for hypergraphs without assortativity; third, we and general diffusion processes [9]. present a generative model for assortative hypergraphs; Recently, there has been much work on using hyper- fourth, we employ a perturbation approach to derive the graphs to describe systems with interactions involving effect of degree-degree mixing on these quantities and de- multiple agents [10, 11]. Hypergraphs are useful to de- fine the dynamical assortativity; and lastly, we show how scribe multi-way interactions in biology [12], social conta- our results can be used to modify hypergraph dynamics gion [13{15], synchronization [16], opinion models, infec- through preferential rewiring of hyperedges. tious disease spread [17], and real data [18]. It has been We start by defining terminology. A hypergraph is shown that contagion on a hypergraph cannot be pro- a mathematical object that describes group interactions jected onto a pairwise network if the higher-order conta- among a set of nodes. We represent it as H = (V; E), gion function is non-linear [19]. Recently the pairwise no- where V is the set of nodes and E is the set of hyper- tion of assortativity has been extended to hypergraphs for edges, which are subsets of V and represent unordered categorical node labels [20{22] and continuous attributes interactions of arbitrary size. We call a hyperedge with [18]. Assortativity on hypergraphs can provide different cardinality m an m-hyperedge and a hypergraph with insights on the structure of the interactions than assorta- only m-hyperedges an m-uniform hypergraph. It is useful tivity on the pairwise network projection [18] and, as we to consider weighted hypergraphs, where each edge e has will show, affect the outcome of hypergraph dynamical an associated positive weight βe. We define the hyperde- processes. gree sequence as in Ref. [14], where the mth order hyper- degree of node i, k(m), is the number of m-hyperedges to A fundamental problem when studying dynamics on i which it belongs. networks is to determine how structural characteristics We now define the expansion eigenvalue and discuss of the network affect the dynamical behavior. Many dy- its relevance to dynamical processes on hypergraphs. For namical properties such as the onset of epidemic spread- a weighted hypergraph, the expansion eigenvalue λ and arXiv:2109.01099v1 [physics.soc-ph] 2 Sep 2021 ing [23], synchronization [24], and percolation [25] are de- associated eigenvector u are defined by the eigenvalue termined by the largest eigenvalue of the network's ad- equation jacency matrix. In turn, this eigenvalue is affected by the network's degree distribution and assortative mixing X λui = βe(ui + ··· + ui ); (1) properties [26] as well as other structural characteristics. 1 m−1 e=fi;i1;:::;im−1g2E In this paper we show how the expansion eigenvalue, a suitably generalized eigenvalue for hypergraphs, is simi- where λ and u are the Perron-Frobenius eigenvalue and larly modified by the assortative properties of the hyper- eigenvector of the associated nonnegative matrix. graph. This eigenvalue has been shown to determine the In the unweighted case (i.e., βe = 1 for every hyper- extinction threshold for the SIS model on hypergraphs edge e), the eigenvector u corresponds to the Clique motif [27], and we believe it will also prove useful in relating Eigenvector Centrality [28]. Just as the largest eigenvalue 2 (0) of a network's adjacency matrix is determinant for net- By symmetry of the function fm , we get work dynamics, the expansion eigenvalue plays an impor- " # tant role in dynamical processes on hypergraphs. For ex- X P (k1) k1 uk λu = (m − 1) 1 k; (3) ample, consider an SIS process on a hypergraph, where a k Nhki healthy node can get infected via a hyperedge e to which k1 it belongs at rate βe if at least one other node in e is and multiplying both sides by kP (k)=(Nhki) and sum- infected (the case referred to as individual contagion in ming over k, we obtain for the uncorrelated case Ref. [14]) and heals spontaneously at rate γ. Then one can show (see Ref. [27] and Supplemental Material) that hk2i λ ≈ (m − 1) ; (4) a sufficient condition for epidemic extinction is given by hki γ > λ. In addition, we can determine the relative impor- tance of a node i with respect to this contagion model and uk / k from Eq. (3). For non-uniform hypergraphs, (in terms of its probability of infection at the onset of assuming that βe = βjej, the expansion eigenvalue is the the epidemic) from the ith entry of the associated eigen- largest eigenvalue of K, the degree-size correlation ma- (i) (j) (i) vector. More generally, the importance of the expansion trix, with entries Kij = βi(i − 1)hk k i=hk i (see the eigenvalue in spreading processes can be understood from Supplemental Material for the derivation). the fact that in the unweighted case the number of nodes In contrast with the uncorrelated case, we now assume reachable via hyperedges from a given starting node in ` that nodes are connected with an arbitrary function fm steps grows asymptotically as λ` [28]. determining the connection probability. We define (0) Having motivated the importance of the expansion fm(k1; : : : ; km) = fm (k1; : : : ; km) [1 + gm(k1; : : : ; km)] ; eigenvalue λ, we now focus on how it is affected by (5) assortative mixing by degrees in the hypergraph. In this paper we will consider random hypergraphs that where is a parameter and gm an assortativity function are constructed from a prescribed hyperdegree sequence for m-uniform hypergraphs. The assortativity function gm(k1; : : : ; km) determines how likely it is that nodes fk1;:::; kN g, where N is the number of nodes, ki = (2) with degrees k ; : : : ; k are joined by a m-hyperedge; [k ; : : : ; k(M)] is the target hyperdegree of node i, 1 m i if g > 0 (g < 0) it is more (less) likely than it and M is the maximum hyperedge size. The hy- m m would be expected if they were connected at random. In pergraph is then constructed by creating a hyperedge order to preserve the degree sequence, g must satisfy fi ; : : : ; i g with probability f (k ;:::; k ), where the m 1 m m i1 im P (0) fm (k1; : : : ; km)gm(k1; : : : ; km) = 0. functions fm specify the assortative mixing properties k1;:::;km We present two examples of the assortativity function of the hypergraph. The function fm(ki1 ;:::; kim ) = (0) (m) (m) (m) (m) m−1 gm: the aligned degrees and large degrees assortativity fm (k ; : : : ; k ) = (m − 1)!k : : : k =(Nhki) , i1 im i1 im functions. The aligned degrees function measures the PN where we define hxi = i=1 xi=N, corresponds to the average pairwise correlation on a hyperedge. It is defined case where nodes are connected with hyperedges com- as pletely at random if the hyperdegree of node i is ki (as 1 X (ki − hki) (kj − hki) in the hypergraph configuration model [14, 18, 29, 30]). g (k ; : : : ; k ) = We start by deriving a mean-field approximation for the m 1 m m hki hki 2 i<j expansion eigenvalue λ in this case, which we call the hk2i − hki2 2 uncorrelated case, before considering hypergraphs with − : (6) degree assortativity. For simplicity, from now on we will hki2 consider an unweighted m-uniform hypergraph, and will This function is similar to the expected value of the uni- denote k(m) by k and refer to it as the degree of node i i form function described in Ref. [18]. Likewise, the large i. Now we assume that all nodes with the same degree degrees function is defined as are statistically equivalent and that the eigenvector en- m m try of node i depends only on its degree, i.e., ui ! uki . 2 Y ki hk i Defining P (k) to be the number of nodes with degree g (k ; : : : ; k ) = − ; (7) m 1 m hki hki2 k, the equation defining the expansion eigenvalue can be i=1 written as and measures whether the product of degrees on a hy- peredge is larger than expected for a configuration null model.