Emergent Hyperbolic Network Geometry.

A large variety of interacting complex systems are characterized by interactions occurring between more than two nodes. These systems are described by simplicial complexes. Simplicial complexes are formed by simplices (nodes, links, triangles, tetrahedra etc.) that have a natural geometric interpretation. As such simplicial complexes are widely used in quantum gravity approaches that involve a discretization of spacetime. Here, by extending our knowledge of growing complex networks to growing simplicial complexes we investigate the nature of the emergent geometry of complex networks and explore whether this geometry is hyperbolic. Specifically we show that an hyperbolic network geometry emerges spontaneously from models of growing simplicial complexes that are purely combinatorial. The statistical and geometrical properties of the growing simplicial complexes strongly depend on their dimensionality and display the major universal properties of real complex networks (scale-free degree distribution, small-world and communities) at the same time. Interestingly, when the network dynamics includes an heterogeneous fitness of the faces, the growing simplicial complex can undergo phase transitions that are reflected by relevant changes in the network geometry.

Simplicial complexes are the many-body generalization of networks123456 and they can encode interactions occurring between two or more nodes7891011121314. While networks are formed exclusively by nodes and links, simplicial complexes include higher dimensional simplices i.e. triangles, tetrahedra etc. As such they are fundamental to study a large variety of real complex interacting systems, including brain functional networks7, protein interaction networks15, collaboration networks12. Because simplices have a natural topological and geometrical interpretation, simplicial complexes are ideal to investigate the underlying geometry and topology of networks7161718 and for these reasons they are extensively used in quantum gravity1920212223.

One of the fundamental quests of quantum gravity is to describe the emergence of a continuous, finite dimensional space, using pre-geometric models, where space is an emergent property of a network or of a simplicial complex21222425. This fundamental mathematical problem has its relevance also in the field of network theory16 where one of the major aim of network geometry is to characterize the continuous hidden metric behind the inherently discrete structure of complex networks. In fact, it is believed that most complex networks have a continuous hidden network geometry2627282930 such that any two connected nodes are also close in the hidden metric29303132. In this context, there is increasing evidence that the hidden geometry of a large variety of networks including the Internet, airport networks, the brain functional networks, and metabolic networks27282930 is hyperbolic. Characterizing the hyperbolicity of networks is not only a fundamental theoretical question, but it can also have practical implications as it can be used to improve significantly the navigability on such networks2930. While the mathematical definition of the curvature of networks is a hot mathematical subject for which different definitions have been given163334353637, most of the results obtained so far are related to the embeddings of complex networks in hyperbolic spaces2728293032.

The underlying assumptions of several models2627 of complex hyperbolic networks is that nodes are sprinkled randomly in the hidden hyperbolic metric and links are established according to their hyperbolic distance. Interestingly this type of models can be related to causal sets38 in de Sitter space and they have been used to describe a “network cosmology”39. From the complexity point of view, if we want for example to use this type of models for describing the evolution of the World-Wide-Web, the sprinkling of the nodes in hyperbolic space reflects some distribution of interest of the webpage owners, and links between the webpages are established depending on the similarities between the interests of webpage owners. Although this is a very plausible mechanism for network evolution, it cannot be adopted if we aim at describing the emergence of the underlying hyperbolic geometry as the result of the endogenous dynamics of the network.

Here we will show in the framework of a very simple, stylized model, that the hyperbolic network geometry can be an emergent property of growing simplicial complexes that share the universal properties of complex network structures. Specifically, we will propose a model in which the hidden hyperbolic metric is not causing the network dynamics but it is instead the outcome of the network evolution.

Our model of emergent geometry is based on a growing simplicial complex. Metric spaces satisfy the triangular inequality, therefore a network with non-trivial geometry should include a high clustering coefficient and high density of triangles, ensured by building the network using simplicial complexes. Additionally growing networks have been extensively used as a non-equilibrium framework1340414243 for the emergence of complex statistical properties of networks such as the power-law degree distribution. By extending the well established framework of growing network models1340414243 to simplicial complexes we will provide significant new insights into emergent geometry. Importantly, the emergent hidden geometry of growing simplicial complexes is hyperbolic, i.e. the hyperbolic geometry emerges spontaneously from the evolution of the simplicial complexes. In this way we provide evidence that hyperbolic network geometry emerges from growing simplicial complexes whose temporal evolution is purely combinatorial, i.e. it does not take into account the hidden geometry.

Interestingly the properties of the network geometry change significantly with the dimension of the simplicial complex and the network geometry can be strongly affected by phase transitions occurring when a fitness parameter9104243 is associated to each face of the simplicial complexes describing intrinsic local heterogeneities.

Results

We consider simplicial complexes formed by gluing together d-dimensional simplices. A d-dimensional simplex (or d-simplex) is a topological object including the set of d + 1 nodes and all its subsets. The underlying network structure of a d-simplex is constituted by a fully connected network, or a clique, of d + 1 nodes, such as links (1-simplices), triangles (2-simplices), tetrahedra (3-simplices) etc. The δ-faces of a d-dimensional simplex are all the δ-dimensional simplices that can be built by a subset of (δ + 1) of its nodes. For example the faces of a triangle (2-simplex) are its three links (1-simplices) and its three nodes (0-simplices), the faces of a tetrahedron are its four triangular faces (2-simplices), its six links (1-simplices) and its four nodes (0-simplices), etc. As long as we are concerned exclusively with the network properties of simplicial complexes, the use of simplicial complexes is equivalent with the use of hypergraphs and hypernetworks that are recently attracting increasing attention1244.

The simplicial complexes that we are considering in this paper are constructed by gluing d-simplices along their (d − 1)-faces. To every (d − 1)-face α of the simplicial complex, (i.e. a link for d = 2, or a triangular face for d = 3) we associate an incidence number n given by the number of d-dimensional simplices incident to it minus one. The simplicial complex dynamics is dictated by the following algorithm and depends on a parameter s = −1, 0, 1 called flavor. We start from a single d-dimensional simplex, i.e a triangle for d = 2, a tetrahedron for d = 3. At each time we add a d-dimensional simplex to a (d − 1)-face α. The face α is chosen randomly with probability Π given by

The new d-dimensional simplex is induced by a new node and all the nodes of the chosen (d − 1)-face α. For this type of dynamics, the combinatorial condition to obtain a discrete manifold is that n can take exclusively the values n = 0, 1.

For s = −1 it is possible to attach a simplex only to faces with n = 0. In fact for n = 0 we have but for n = 1 we have Π = 0. As a consequence of this, the resulting network is a manifold, with each (d − 1)-face incident at most to two d-dimensional simplicial complexes, i.e. n = 0, 1. For s = 0 the d–dimensional simplices are attached with uniform probability to any (d − 1)-face, while for s = 1 the dynamics follows a generalized preferential attachment and the new simplex is attached to a (d − 1)-face α proportionally to the number of simplicies already attached to the face. i.e. 1 + n. Therefore for s = 0 as for s = 1 the incidence number n can take values n = 0, 1, 2, 3 ….

Simplices are topological objects that can be turned into geometrical entities when we attribute a given length to their links. Here, in order to describe the emergent geometry of our model of growing simplicial complexes, we assume that every simplicial complex is built by simplices that have links of equal length across the entire simplicial complex.

The resulting networks are small world for every flavor s and any dimension d except from the special case s = −1, d = 1 in which the resulting network is a chain. This implies that the number of nodes in the network N increases exponentially with its diameter D, i.e. . Therefore, if all the links have equal length, the hidden geometry of these networks cannot be the one of a Euclidean space of finite Hausdorff dimension d because this would imply a power-law scaling . As a consequence of this the small-world property suggests that the natural embedding of these networks is hyperbolic. Nevertheless the small-world property might not be sufficient to guarantee an embedding in the hyperbolic space. Here we show that for our class of growing simplicial complexes the hidden geometry, corresponding to the embedding where all the links have the same distance, are the hyperbolic spaces , and specifically the Poincaré ball model45.

The great advantage of the present class of models with respect to general small-world networks is that their dual is a tree. The dual network can be constructed by associating to every d-simplex a node of the dual network, and to every pair of d-simplicies sharing a (d − 1)-face a link of the dual network. Since in our simplicial complex evolution at each time we glue a new d-simplex to a (d − 1)-face, the resulting structure of the dual network is a tree. Taking advantage of this simple structure of the dual the present class of models admits several embeddings in the hyperbolic space model. Between the different possible embeddings, only one embedding can fill the entire space in the asymptotic limit t → ∞. Therefore this embedding defines the emergent geometry of our simplicial complexes.

Let us consider a Poincaré ball model of . The Poincaré ball model includes all the points of the unit ball , with |…| indicating the Euclidean norm. The Poicaré ball model is associated to the hyperbolic metric d assigning to each pair of points the distance

Here we identify every d-dimensional simplex of our simplicial complex with an ideal simplex of the Poincaré ball model. An ideal simplex has all its nodes at the boundary of the hyperbolic ball, so all the nodes i have a position satisfying |r| = 1. This allows to have all the nodes of each simplex at equal hyperbolic distance. Note that interestingly this distance is actually infinite but this is the “cost” required for having an embedding that asymptotically in time fills the entire hyperbolic space. In order to fully characterize the hidden geometry of the studied networks, we need also to determine further the position of the ideal nodes at the boundary of the ball. To this end we start from a d-dimensional simplex whose (d + 1) nodes have the same maximum (Euclidean) distance from each other. Therefore the positions r of the initial (d + 1) nodes i = 1, 2, … d + 1 satisfy

Each new node i of the network has a position at the boundary of the ball determined by the position of its “ancestors”, i.e. the nodes of the face α connected to the new node. In particular the new node i is placed at equal (Euclidean) distance from all the nodes j of the α face to which it is attached, i.e.

In this way the angular position of the new node is fully determined by the stochastic dynamics of the network (see Fig. 1 for details).

Figure 1

The first steps of the temporal evolution of growing simplicial complex with N nodes is shown here in the emergent hyperbolic space for d = 2 (panel A) and d = 3 (panel B). The flavor is s = −1.

The resulting networks have a rich geometrical structure, which is linked to the mathematics of Farey sequences in d = 246. Additionally, the simplicial complexes in dimension d = 3 are characterized by a boundary with notable geometrical features.The induced geometry on this boundary can be studied by placing the nodes of the network in a (d − 1)-dimensional space characterized by the angular coordinates of the nodes. The network resulting from the projection on the boundary of the ball B is a random Apollonian network4748 for d = 3.

Let us make three important observations related to the geometric nature of the proposed class of growing simplicial complexes.

First of all we note that the hyperbolic nature of the emergent geometry it is a consequence of the assumption that each link of the simplicial complex must have equal length. This assumption implies that distances of different links can be compared. Therefore strictly speaking here the network geometry is actually a consequence of a kind of “proto-geometry” that allows comparison of length of different links. If we allow, instead, to have links of different lengths the curvature of the hidden geometry is not determined and it is even possible to tune the length of the links such that the same simplicial complex can be embedded in a d-sphere. For instance this can be achieved simply by taking the embedding on the Poicaré ball model described above, and considering instead of the hyperbolic metric on the ball the Euclidean metric.

Secondly we note that the natural hyperbolic embedding of growing simplicial complexes that we discussed above, works particularly well for flavor s = −1 while some caution is required when using this embedding for flavors s = 0 and s = 1. In fact, for s = −1 the simplicial complexes are manifolds and as a consequence of this, links do not cross and each node of the simplicial complex has a distinct position in the hyperbolic space. However for flavor s = 0 and s = 1 the proposed embedding implies that some nodes of the simplicial complex (the nodes that are immediate “descendant” of the same face) acquire the same position in the hyperbolic space. As a result in a geometrical embedding links are effectively weighted. For this reason the growing simplicial complexes with flavor s = −1 play a very special role with respect to the other flavors s = 0 and s = 1.

Finally we note that for s = −1 and d = 3 the growing simplicial complex model presented here belongs to the class of stacked polytopes that are equivalent to Apollonian packings, whose discrete Lorentzian geometry is raising recent interest in the mathematical community48495051. In fact, these stacked polytopes in d = 3 have a symmetry group G that is a noncompact discrete subgroup of the Lorentz group SO(3, 1) = SL(2, )/. Therefore these results provide an additional important insight on the hyperbolic nature of the underlying geometry of the class of models proposed in this paper.

The networks resulting from the proposed model of growing simplicial complexes are strongly affected by their dimensionality. In fact they are scale-free for dimension

while for dimension d = 1 − s they have exponential degree distribution. In particular, the degree distribution of growing simplicial complexes of dimension d and flavor s is given for for d + s = 1 by (see Methods for details)

with k ≥ d while for d + s > 1 it is given by (see Methods for details)

with k ≥ d.

Therefore for d + s > 1 the degree distribution is scale-free and has a power-law scaling

for k ≫ 1, with power-law exponent

Finally for d + s = 0 the simplicial complexes reduce to chains (see Methods for details).

In Fig. 2 we show the perfect agreement between the predicted degree distribution and simulation results for d = 1, 2, 3 and flavor s = −1, 0, 1. For dimensions d > 1 these networks display a significant community structure (high modularity) and high average clustering coefficient (see Methods) as most complex networks. The values of the modularity and the clustering coefficient are modulated by the dimension d and the flavor s of the growing simplicial complex (see Table 1).

Figure 2

The effect of dimensionality and flavor in the degree distribution .

The symbols (blue dimonds) indicate the simulation result of a single realization of the growing simplicial complex with dimension d and flavor s with N = 105 nodes. The solid lines indicate the theoretical predictions.

Table 1

Modularity M and average clustering coefficient C of the growing simplicial complex with N = 104 nodes averaged over 20 realizations are reported here for the dimension d = 2, 3 and flavor s = −1, 0, 1.

M	s = −1	s = 0	s = 1	C	s = −1	s = 0	s = 1
d = 2	0.97	0.94	0.90	d = 2	0.65	0.74	0.79
d = 3	0.91	0.85	0.80	d = 3	0.77	0.81	0.84

The modularity M is obtained by running the generalized Louvain method5253.

The emergent geometry of growing simplicial complexes is strongly affected by phase transitions occurring in the network evolution. In order to show evidence for this statement we study a variation of the model including fitness of the faces of the simplicial complex. First we associate to each node an energy drawn from a distribution. To any δ-face α with 0 < δ < d we associate an energy

and a fitness

where β = 1/T > 0 is a global parameter called inverse temperature, such that for β = 0 all the nodes have the same fitness while for β ≫ 1 small differences in energy yield big differences in the fitness of the nodes. The model remains the same with the exception that the probability Π that a (d − 1)-face α is selected is no longer given by Eq. (1) but is given by

This model displays a structural phase transition at high value of β. In order to characterize this phase transition we make use of the hyperbolic embedding. While above the phase transition the simplicial complex is growing in all directions of the hyperbolic space and has statistical properties related to quantum statistics910, below the phase transition there is symmetry breaking and the network evolves asymmetrically in the hyperbolic space. In Fig. 3 we show a visualization of the model above and below the phase transition for dimension d = 2, 3 showing that also the geometry of the boundary of simplicial complexes in d = 3 is strongly affected by the geometrical phase transition occurring in the model.

Figure 3

The hyperbolic emergent network geometry is shown here for the growing simplicial complex with fitness of the nodes and flavor s = −1 above and below the phase transition in dimension d = 2 (panel a-b) and dimension d = 3 (panel c-d). Panel (e-f) display the projection of the network on the boundary of the d = 3 hyperbolic space. The energy distribution is uniform over discrete values of the energies of the nodes 0 ≤ ε < 10. The inverse temperature is β = 0.01 for panels a, c, e, β = 50 for panel b, β = 20 for panels d, f. The number of nodes is N = 200. The color of the links indicates the different values of their energies.

In order to numerically study the phase transition occurring in the growing simplicial complexes with fitnesses we define a vector R given by

where r is the (Euclidean) position vector of the node i in the Poincaré ball. We study the Euclidean norm R = |R| and the standard deviation and the maximum degree k as a function of the inverse temperature β. Here we focus on the result for d = 2, 3, flavor s = −1 (see Fig. 4). As β increases across the phase transition, σ develops relevant finite size effects and becomes vanishingly small in the large network limit. For the same values of the parameters R approaches one indicating that the simplicial complex grows in a well defined preferential direction. The phase transition in the network geometry for flavor s = −1 has these characteristics both for d = 2 and d = 3. However the behavior of the maximum degree k across the phase transition shows major differences between the case d = 2 and d = 3, displaying a clear maximum at the transition point for d = 3. This implies that the transition might affect the degree distribution in different ways depending on the dimension d. Interestingly similar transitions are observed for different flavors s = 0, 1.

Figure 4

The phase transition in the network geometry of the growing simplicial complexes with fitnesses is characterized here by showing R, σ and the maximum degree k as a function of β.

Here simulation results are reported for growing simplicial complexes with flavor s = −1, dimension d = 2 (left panels), and d = 3 (right panels) and network sizes N = 2500 (blue dashed line), 5000 (red dot dashed line), 10000 (green solid line). All the data are averaged over 500 realizations.

Discussion

In conclusion, this paper shows that the study of simplicial complexes allows great advances in our understanding of complex networks. In fact by extending the framework of growing complex networks to simplicial complexes our simple model produces networks that display most of the universal properties of complex networks including scale-free degree distribution, small-world properties and significant modular structure. These networks have statistical and geometrical properties that are a function of their dimension d and their flavor s that modulate the values of their modularity and their clustering coefficient. These non-equilibrium models of simplicial complexes are ideal frameworks to show the emergence of the hyperbolic network geometry. Specifically they can explain how real hyperbolic networks might result from purely combinatorial rules that do not take into consideration the hidden metric of the network.This network geometry has a very interesting structure linked to Farey sequences and Apollonian tilings. Additionally these network geometries can undergo relevant changes following phase transitions in the network evolution.

We believe that this paper, showing that growing simplicial complexes give rise to a complex emergent hyperbolic geometry, related to Apollonian packings, is opening new perspectives for understanding the origin of the emergent hyperbolic geometry of complex networks. On the same time, our hope is that further research in this direction could also indicate a path for establishing a cross-fertilization between network theory and quantum gravity.

Methods

Degree distribution of growing simplicial complexes with flavor s

In order to derive the degree distribution of the growing simplicial complexes with flavor s, we use the master equation approach34. It can be easily shown that the average number of nodes of degree k that at each time increase their degree by one is given by

as long as d + s ≠ 0, i.e. (d, s) ≠ (1, −1) for which the growing simplicial complex reduces to a chain. To derive the exact degree distribution of the simplicial complex, we consider the master equation for the average number of nodes that at time t have degree k in a growing d dimensional simplicial complex of flavor s. The master equation [3, 4] for reads

with k ≥ d.

Solving this equation we get both exponential and power-law degree distribution. In particular, the degree distribution of growing simplicial complexes of dimension d and flavor s is given for for d + s = 1 by

with k ≥ d while for d + s > 1 it is given by

with k ≥ d.

Therefore for d + s > 1 the degree distribution is scale-free and has a power-law scaling

for k ≫ 1, with power-law exponent

Community structure of growing simplicial complexes

The modularity M evaluates the significance of the community structure of a network. It is defined4 as

where a denotes the adjacency matrix of the network, L the total number of links, and {q}, where q = 1, 2 … Q, indicates to which community the node i belongs. Finding the network community structure that optimizes modularity is NP hard. One of the most popular greedy algorithms to find the community structure is the generalized Louvain method5253 that is able to determine a lower bound on the maximum modularity of the network. As shown in Table 1 of the main text the growing simplicial complexes with d ≥ 2 are characterized by a large modularity.

Average clustering coefficient of growing simplicial complexes

The local clustering coefficient C measures the density of triangles passing through node i, i.e.

where a denotes the adjacency matrix of the network, and where k = ∑a is the degree of node i.

As shown in Table 1 of the main text the growing simplicial complexes with d ≥ 2 are characterized by a large clustering coefficient.

Additional Information

How to cite this article: Bianconi, G. and Rahmede, C. Emergent Hyperbolic Network Geometry. Sci. Rep. 7, 41974; doi: 10.1038/srep41974 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.