Noname manuscript No.
(will be inserted by the editor)
Guillermo Garcı́a-Pérez · M. Ángeles Serrano · Marián Boguñá
arXiv:1707.09610v1 [physics.soc-ph] 30 Jul 2017
Soft communities in similarity space
the date of receipt and acceptance should be inserted later
Abstract The S1 model has been a central geometric
model in the development of the field of network geometry. It has been mainly studied in its homogeneous
regime, in which angular coordinates are independently
and uniformly scattered on the circle. We now investigate
if the model can generate networks with targeted topological features and soft communities, that is, heterogeneous angular distributions. Under these circumstances,
hidden degrees must depend on angular coordinates and
we propose a method to estimate them. We conclude that
the model can be topologically invariant with respect to
the soft-community structure. Our results might have
important implications, both in expanding the scope of
the model beyond the independent hidden variables limit
and in the embedding of real-world networks.
1 Introduction
Complex networks have been widely studied in the last
twenty years in many different contexts, from biological
to social and technological [1, 2]. There seems to be some
Guillermo Garcı́a-Pérez
Departament de Fı́sica de la Matèria Condensada, Universitat de Barcelona, Martı́ i Franquès 1, 08028 Barcelona, Spain
Universitat de Barcelona Institute of Complex Systems
(UBICS), Universitat de Barcelona, Barcelona, Spain
E-mail: guille.garcia@ub.edu
M. Ángeles Serrano
Departament de Fı́sica de la Matèria Condensada, Universitat de Barcelona, Martı́ i Franquès 1, 08028 Barcelona, Spain
Universitat de Barcelona Institute of Complex Systems
(UBICS), Universitat de Barcelona, Barcelona, Spain
ICREA, Pg. Lluı́s Companys 23, E-08010 Barcelona, Spain
E-mail: marian.serrano@ub.edu
Marián Boguñá
Departament de Fı́sica de la Matèria Condensada, Universitat de Barcelona, Martı́ i Franquès 1, 08028 Barcelona, Spain
Universitat de Barcelona Institute of Complex Systems
(UBICS), Universitat de Barcelona, Barcelona, Spain
E-mail: marian.boguna@ub.edu
universal features common to the topology of many networks. For instance, in most cases, they are scale-free,
meaning that their degrees are power-law distributed, a
phenomenon that was explained in early times of network theory by the preferential attachment mechanism:
as the network grows, new nodes connect to highly connected —or popular— nodes with higher probability [3].
However, preferential attachment alone cannot explain the high level of clustering coefficient —the fraction of existing triangles— observed in real systems. To
explain clustering, the concept of similarity was introduced [4]. The basic idea is that nodes connect not only
because they are popular, but also because they are similar in some sense. Thus, if node A connects to nodes B
and C because they are similar to A, B and C should
also be similar and therefore have a high probability of
being connected. This transitivity of similarity suggests
encoding similarities between nodes as distances in metric spaces, since the triangle inequality is one of their
defining properties: if the distance dBC in the underlying metric space measures the dissimilarity between B
and C, it must be bounded by dBC ≤ dAB + dAC , therefore inducing the observed transitive connections.
The S1 model was proposed based on these ideas [4].
In this model, N nodes are randomly scattered into a
circle of radius R = N/2π. Every node i is also assigned
a hidden degree κi from any distribution (for instance, a
power-law P (κ) ∼ κ−γ ), and every pair of nodes i and j
is connected with probability
1
pij =
1+
dij
µκi κj
β ,
(1)
where dij is the distance along the circle; µ and β are two
global parameters controlling the average degree and the
clustering coefficient, respectively. Notice that this connection probability takes the form of a gravity law, as it
increases with the product of hidden degrees (popularities) and decreases with the distance between them (dissimilarity). Despite its apparent simplicity, this model
2
generates networks that resemble very much real networks; they are scale-free, small-world and have high levels of clustering. In fact, the degrees ki are proportional
to the hidden degrees κi , so the model is versatile enough
to generate networks with different degree distributions1 .
The possibilities of the S1 model go beyond generating realistic networks. The similarity space coordinates
of the nodes of a real network can be inferred by finding the coordinates that maximise the likelihood for the
real network to be generated by the model [6, 7, 8]. This
embedding process yields a map of the network strikingly meaningful. For instance, it allows to navigate the
network efficiently by mapping the coordinates to hyperbolic space. Moreover, being able to access to the similarity space coordinates of nodes opens the path to a
completely new way of analysing complex networks. For
example, in Refs. [6, 7, 8] it was shown that the angular
coordinates of nodes in real-world networks are not uniformly distributed. Instead, they are distributed in a heterogeneous manner, with angular regions more densely
populated than others. These dense regions reveal the
community structure of the network [9, 10, 11]. Indeed, by
partitioning the network using the largest gaps between
consecutive nodes along the circle as community boundaries, the partitions obtained have a modularity comparable to that of other community detection methods
currently available in the literature. Furthermore, this
geometric method seems to have higher resolution [8].
In Ref. [12], the authors introduced the Geometric
Preferential Attachment (GPA) model, a generalised version of the growing geometric model Popularity vs. Similarity Optimization [13] in which soft communities, as
they named these denser angular regions, emerge from
the growth dynamics of the network without altering
topological properties like the degree distribution or the
clustering spectrum of the resulting network. In this paper, we address the question of whether the S1 model
can generate networks with given target topological features and soft communities, that is, heterogeneous angular distributions. To that end, any angular distribution
could be considered in principle. For instance, one could
impose some non-uniform distribution function a priori.
However, the angular distribution from Ref. [12] is not
an imposition, but it rather emerges from a preferential
attachment process in similarity space that seems to be
a plausible explanation for the nature of communities in
real systems. We focus on that particular angular distribution for our study.
1
The original model defined in [4] is in fact more general,
allowing for any connection probability pij as long as it depends on the argument dij /(κi κj )1/D , where now the space
is the D-dimensional sphere and dij the geodesic distance on
the sphere. The particular functional form in Eq. (1) allows us
to interpret the network as a set of non-interacting fermions
(the links) embedded in the hyperbolic plane, with the hyperbolic length of a link playing the role of its energy and β
playing the role of the inverse of the bath temperature [5].
2 Results
In the bare S1 model, hidden degrees and similarity coordinates are typically assumed to be uncorrelated, so every node’s hidden variables are withdrawn independently
from some joint distribution ρ(κ, θ) that factorises [4, 14].
Nevertheless, the GPA is a growing model, so the angular coordinates and hidden degrees of different nodes are
correlated. We are therefore forced to drop such simplifying assumptions.
In the GPA growth process, the degree of a node is
determined by its age —the older the node is, the higher
its degree. Moreover, when a new node t is added to
the system, the probability for it to be placed at polar
coordinate θt depends on the number of nodes s < t
γ−2
1
at angular distance ∆θst < 2/(s γ−1 t γ−1 ), where γ is
the exponent of the power-law degree distribution. This
implies a very particular dependence between similarity coordinates and degrees: the angular coordinate of
a node must depend on the angular coordinates of all
nodes with higher degree. Hence, we must include the
implicit ordering in the sequence of nodes induced by
the degree sequence in our heterogeneous version of the
S1 model. That could be done by first assigning a hidden
degree κi from a power-law distribution P (κ) ∼ κ−γ to
every node, ordering the nodes according to their hidden
degrees and reproducing the angular preferential attachment from the GPA model with that particular ordering.
At the end of the process, we would obtain a set of N
nodes with hidden degrees power-law distributed with
exponent γ and the same angular distribution as the
GPA model for that value of γ. However, if we then connected every pair of nodes with the probabilities given
by Eq. (1), degrees and hidden degrees would not be proportional; the reason for such deviation from the usual
behaviour of the model is that a homogeneous angular
distribution is required for the proportionality between
hidden and observed degrees [14], which is not fulfilled
here by construction; hidden degrees must depend on
the spatial distribution of nodes as well. In the following
subsection, we address this issue. We explore the heterogeneous regime of the S1 model and show that it is
capable of generating networks with power-law degree
distributions, high clustering and soft communities.
2.1 Geometric Preferential Attachment in the S1 model
From the previous discussion, we see that hidden degrees
and angles are considerably entangled in the modelling of
geometric networks with soft communities. In this context, the S1 model requires the following steps:
1. Assigning angular coordinates: Angular coordinates are assigned according to the Geometric Preferential Attachment, which requires an ordering. Therefore, assign a label i = 1, . . . , N to every node. Then,
for every value of i from 1 to N :
3
γ = 2.1
Λ = 10−1
γ = 2.5
Λ = 10−1
γ = 3.0
Λ = 10−1
γ = 2.1
Λ = 100
γ = 2.5
Λ = 100
γ = 3.0
Λ = 100
γ = 2.1
Λ = 101
γ = 2.5
Λ = 101
γ = 3.0
Λ = 101
Fig. 1 Geometric layout of the networks generated by the S1 model with Geometric Preferential Attachment. In all cases,
N = 1000 and β = 2.5. Every column corresponds to a value of γ and every row to a value of Λ. As in Ref. [12], soft
communities emerge for low values of the initial attractiveness Λ. In order to clarify the figure, every node’s target degree is
tar
tar
tar 2
represented as a radial coordinate ri = R−2 ln kitar /kN
, where kN
is the smallest target degree and R = 2 ln N/(πµ(kN
) ) .
1
When using the hidden degrees instead of the target degrees, this mapping constitutes the isomorphism between the S model
and the H2 model in hyperbolic space [14, 6, 7, 8].
a. Sample i candidate angular positions φj , j = 1, . . . , i
for node i from U (0, 2π).
b. For every candidate position, define the attractiveness A(φj ) of candidate j as the number of
nodes s with an already defined angular position,
that is, with s < i, at angular distance ∆θis <
γ−2
1
2/(s γ−1 i γ−1 ), where γ is the exponent of the target power-law degree distribution.
c. Assign to node i the angular coordinate of candidate j, i. e. set θi = φj , with probability
Π(φj ) = Pi
A(φj ) + Λ
n=1
(A(φn ) + Λ)
.
(2)
The initial attractiveness Λ ≥ 0 is a parameter
that sets the strength of the geometric preferential
4
Pc(k)
100
γ = 2.1
10−1
10
−1
10−2 0
10
100
γ = 2.5
10−1
10−2
10−3 0
10
100
c̄(k)
100
10−2
10
1
10
2
10−3 0
10
10
100
3
Λ = 10−2
Λ = 10−1
Λ = 100
Λ = 101
Λ = 102
10−2
10
1
10
2
10−3 0
10
10
100
3
10−1
101
102
103
10−2 0
10
γ = 3.0
10−1
101
102
101
k
102
10−1
101
k
102
103
10−2 0
10
k
Fig. 2 Topological properties of the networks generated by the S1 model with Geometric Preferential Attachment with
N = 1000 and β = 2.5. Every color corresponds to a different
value of the initial attractiveness Λ. The top row shows the
R∞
complementary cumulative degree distribution Pc (k) = k P (k)dk, which behaves as Pc (k) ∼ k1−γ for P (k) ∼ k−γ . Black
dashed lines indicate such scaling. In the bottom row, the mean local clustering spectra c̄(k) are drawn. To avoid fluctuations
in the target degrees, all networks with the same exponent γ have been generated with the same target degree sequence
{kitar }. Despite their angular distributions being completely different (Fig. 1), their topologies are extremely similar.
attachment. For very high values of Λ, all candidate angles become equally likely, so the resulting
angular distribution is homogeneous (see Fig. 1).
This process generates a distribution of nodes in the
circle analogous to the angular distribution of the
GPA model. However, notice that, in the GPA model,
connections are established at the same time as positions are decided, whereas in the former steps, no
connections have yet been made.
2. Assigning hidden degrees: Once every node has
a defined angular position, we need to determine its
hidden degree such that the resulting observed degrees, that is, after the connections have been actually established, are power-law distributed with exponent γ. As mentioned earlier in this paper, we must
take into account that the spatial distribution is heterogeneous (especially for low values of Λ). We propose the following method:
tar
a. Generate a set of N target degrees k
from a
power-law distribution with exponent γ. Order
the target degrees such that k1tar > k2tar > · · · >
tar
kN
.
b. Assign to every node i a hidden degree κi , initially
set to κi = kitar .
c. Repeat N times:
i. Choose some node i randomly.
ii. Compute the expected degree k̄i of node i as
k̄i =
X
j6=i
1
1+
dij
µκi κj
β .
(3)
iii. Correct the value of κi so that the expected
tar
degree k̄i matches the target degree
ki . We
propose to reset |κi + kitar − k̄i δ| → κi ,
where δ is a random variable withdrawn from
the uniform distribution U (0, 0.1). Other numerical methods could be used with the same
end.
d. Compute all relative deviations
ǫi =
|kitar − k̄i |
.
kitar
(4)
If max {ǫi }i < η, where η is a tolerance which we
set to η = 10−2 , continue to step 3. Otherwise, go
back to step 2c.
3. Generating the network with the S1 model: In
this last step, we simply connect every pair of nodes
with the probabilities in Eq. (1). Since step 2 assigns
a hidden degree to every node such that its expected
degree matches its target degree, the resulting observed degrees in the network must be similar to the
target degrees as well.
5
Figure 1 shows the networks generated by the model for
different values of γ and Λ. As in Ref. [12], the angular distribution has an evident soft-community structure
for low values of Λ, whereas for high values of the initial attractiveness, the angular density resembles that of
the homogeneous S1 model. Despite the considerable differences in the similarity space distances between nodes
for different values of Λ, the displayed networks are extremely similar from a topological perspective (see Fig. 2),
with almost undistinguishable degree distributions and
clustering spectra. Notice that step 2 is not specifically
designed for the GPA angular distribution2 . In principle,
it should be valid for other distributions as well.
angular distribution requires correcting hidden degrees
in such a way that they depend on the hidden variables
of all other nodes. This is a very important result, since
it suggests that inferring hidden degrees via likelihood
maximisation as well might noticeably improve the quality of embeddings of real-world networks with community structure.
Acknowledgements We acknowledge support from a James
S. McDonnell Foundation Scholar Award in Complex Systems; the ICREA Academia prize, funded by the Generalitat
de Catalunya; Ministerio de Economı́a y Competitividad of
Spain projects no. FIS2013-47282-C2-1-P and no. FIS201676830-C2-2-P (AEI/FEDER, UE); the Generalitat de Catalunya
grant no. 2014SGR608.
3 Discussion
References
There is abundant evidence of the geometric origin of
many properties of complex networks, not only regarding
their topology [6, 7, 15, 8, 16], but also their weighted organisation [17]. The field of network geometry has therefore attracted much attention recently, and the S1 model
is one of its cornerstones. On the one hand, it provides
an intuitive and plausible explanation for clustering in
real networks by introducing the concept of similarity
space. On the other hand, it allows to build geometric maps of real networks by embedding them. These
maps are remarkably meaningful, to the extent of predicting symmetries in real systems [16]. In addition, they
are very useful; they can be used to navigate the network efficiently [6], to detect communities [6, 8] or even
to construct smaller-scale replicas of real networks for
efficiently testing dynamics on real networks [16].
So far, the S1 model has only been studied under several simplifying premises, like power-law degree distributions or independent hidden variables. Yet, it has been
able to explain many observed phenomena in complex
networks. However, it can be exploited beyond these assumptions, since the correlation between hidden degrees
and angles might clarify many more topological features
of real-world networks. This work opens the path towards
such line of study by showing that the model does not
require those simplifying assumptions, as it is capable
of generating topologically similar networks with highly
correlated hidden variables.
Moreover, the results presented in this paper might
also have an impact on the embedding of real networks.
Typically, the likelihood maximisation procedure only
seeks the best angular coordinates, whereas hidden degrees are considered to be a function of degree only and
known from the start [6, 7]. This hypothesis is a direct
consequence of the aforementioned simplifying assumptions usually contemplated in the S1 model. Nevertheless, as we have shown in this work, a heterogeneous
2
The ordering of the target degrees might not be necessary in a more general situation where, for instance, hidden
degrees are not correlated with angles.
1. M.E.J. Newman, Networks: An Introduction (Oxford
University Press, Oxford, 2010)
2. S.N. Dorogovtsev, J.F.F. Mendes, Handbook of Graphs
and Networks (Wiley-VCH, Berlin, 2003), chap. Accelerated Growth of Networks
3. A.L. Barabási, R. Albert, Science 286, 509 (1999)
4. M.Á. Serrano, D. Krioukov, M. Boguñá, Phys Rev Lett
100, 078701 (2008)
5. D. Krioukov, F. Papadopoulos, A. Vahdat, M. Boguñá,
Phys. Rev. E 80, 035101 (2009).
DOI 10.1103/
PhysRevE.80.035101.
6. M. Boguñá, F. Papadopoulos, D. Krioukov, Nature communications 1, 62 (2010). DOI 10.1038/ncomms1063.
7. M.A. Serrano, M. Boguna, F. Sagues, Mol. BioSyst. 8,
843 (2012). DOI 10.1039/C2MB05306C
8. G. Garcı́a-Pérez, M. Boguñá, A. Allard, M.A. Serrano, 6
(2016). DOI 10.1038/srep33441.
9. M.E.J. Newman, M. Girvan, Phys. Rev. E 69, 026113
(2004). DOI 10.1103/PhysRevE.69.026113
10. F. Radicchi, C. Castellano, F. Cecconi, V. Loreto,
D. Parisi, Proc. Natl. Acad. Sci. USA 101, 2658 (2004)
11. A. Arenas, A. Fernández, S. Gómez, New Journal of
Physics 10(5), 053039 (2008).
12. K. Zuev, M. Boguñá, G. Bianconi, D. Krioukov, Scientific
Reports 5, 9421 EP (2015).
13. F. Papadopoulos, M. Kitsak, M.A. Serrano, M. Boguna,
D. Krioukov, Nature 489(7417), 537 (2012).
14. D. Krioukov, F. Papadopoulos, M. Kitsak, A. Vahdat,
M. Boguñá, Phys Rev E 82, 036106 (2010)
15. A. Gulyás, J.J. Bı́ró, A. Kőrösi, G. Rétvári, D. Krioukov,
Nature communications 6, 7651 (2015)
16. G. Garcı́a-Pérez, M. Boguñá, M.A. Serrano, Multiscale
unfolding of real networks by geometric renormalization.
arXiv:1706.00394
17. A. Allard, M.Á. Serrano, G. Garcı́a-Pérez, M. Boguñá,
Nature Communications 8 (2017)