Next Article in Journal / Special Issue
Past and Present Trends in the Development of the Pattern-Formation Theory: Domain Walls and Quasicrystals
Previous Article in Journal
Recent Progress in Gamow Shell Model Calculations of Drip Line Nuclei
Previous Article in Special Issue
Gain-Assisted Optical Pulling Force on Plasmonic Graded Nano-Shell with Equivalent Medium Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Polygon-Based Hierarchical Planar Networks Based on Generalized Apollonian Construction

by
Mikhail V. Tamm
1,2,*,
Dmitry G. Koval
1 and
Vladimir I. Stadnichuk
1
1
Faculty of Physics, Moscow State University, 119992 Moscow, Russia
2
CUDAN Open Lab, Tallinn University, 10120 Tallinn, Estonia
*
Author to whom correspondence should be addressed.
Physics 2021, 3(4), 998-1014; https://doi.org/10.3390/physics3040063
Submission received: 1 July 2021 / Revised: 3 October 2021 / Accepted: 11 October 2021 / Published: 8 November 2021
(This article belongs to the Special Issue Dedication to Professor Michael Tribelsky: 50 Years in Physics)

Abstract

:
Experimentally observed complex networks are often scale-free, small-world and have an unexpectedly large number of small cycles. An Apollonian network is one notable example of a model network simultaneously having all three of these properties. This network is constructed by a deterministic procedure of consequentially splitting a triangle into smaller and smaller triangles. In this paper, a similar construction based on the consequential splitting of tetragons and other polygons with an even number of edges is presented. The suggested procedure is stochastic and results in the ensemble of planar scale-free graphs. In the limit of a large number of splittings, the degree distribution of the graph converges to a true power law with an exponent, which is smaller than three in the case of tetragons and larger than three for polygons with a larger number of edges. It is shown that it is possible to stochastically mix tetragon-based and hexagon-based constructions to obtain an ensemble of graphs with a tunable exponent of degree distribution. Other possible planar generalizations of the Apollonian procedure are also briefly discussed.

1. Introduction

It is often convenient to present big volumes of data as a graph, i.e., as a set of objects and binary relations (bonds) between them. This approach naturally arises in numerous contexts ranging from physics of disordered systems [1] and biology [2] to sociology [3] and linguistics (see, e.g., [4,5,6]). The rapid growth in information technology ensures that larger and larger datasets of this type are becoming available. This naturally stimulates interest in the tools to analyze these datasets and simple (or not so simple) reference mathematical models, which can be used to probe their properties. Thus, a rapid development in the last 20 years of a new interdisciplinary field on the boundary of the random graph theory, the data analysis and the statistical physics, known as complex network theory [7,8,9], has occurred.
Among the structural characteristics typical for many experimentally observed networks, there are three especially common and striking (see, e.g., [7]): (i) the small-world property (a very small average node-to-node distance measured along the network), (ii) extremely wide, approximately a power-law distribution of the node degrees (the networks with this property are often called ‘scale-free’), and (iii) large, as compared to referent randomized networks, the concentration of the short circles (e.g., triangles). It is reasonably easy to construct a model network that has one or two of these characteristics, e.g., Erdos–Rényi graphs [10,11] are small-world, random geometrical graphs [12] that have a large clustering coefficient. The Barabási–Albert model [13] generates small-world scale-free networks. The Watts–Strogatts model [14] generates small-world graphs with a large clustering coefficient, etc. Generating all three properties simultaneously is much harder. Random geometric networks in a hyperbolic space [15,16,17] constitute one example of networks with these properties. Another one is the Apollonian network.
The Apollonian network [18,19] is a planar graph that arises naturally as a network representation of the Apollonian gasket, a remarkable object, which is, apparently, the first known fractal (interestingly, its exact fractal dimension is still unknown) [20,21]. The construction of this network can be explained recursively as follows; see Figure 1. Take a triangle, pick a point inside it and connect it to the three corners of the triangle. As a result, one obtains a set of 3 adjacent triangles that form the first-generation Apollonian network. Now, pick a point inside each of the three triangles, and connect it to its corners, this gives a second-generation Apollonian network, then repeat ad infinitum. The resulting network has been studied extensively in recent years, and it has been shown to have many beautiful properties. For example, the degree distribution and the clustering coefficient have been calculated [18], as well as the average path length [22]. Notably, there is an interesting non-planar interpretation of the Apollonian network. Namely, it can be thought of as a simplicial complex in the following way [23]. A first-generation Apollonian network is a tetrahedron (3-simplex). A second-generation Apollonian network consists of four tetrahedra: the original one and another three, each having a common two-face with the original one. A third-generation Apollonian network consists of 13 tetrahedra: one produced in the first generation, three produced in the second generation and nine new ones attached to each free face of the three second-generation tetrahedra, etc. Thus, one can think of an Apollonian network as a regular rooted tree of tetrahedra touching each other by common faces. This construction is easy to visualize in a 3-dimensional (3D) space (see Figure 1b,c), and it makes the Apollonian network a natural discretization of the 3D hyperbolic space in the same way as a regular tree is a natural discretization of the hyperbolic plane. Many properties of the Apollonian network can be calculated exactly, which makes it a nice toy model for the study of various properties of real scale-free networks. As a result, there have been a significant number of papers in recent years studying percolation [23], spin models [24], signal spreading [25], synchronization [26], traffic [27], random walks [28,29], etc., on the Apollonian network.
Despite being such a beautiful and well-studied object, the Apollonian network has certain drawbacks as a model of real networks. Most importantly, it is a single deterministic object with certain fixed properties, e.g., a fixed degree distribution with a fixed power law exponent γ = ln 3 / ln 2 . Importantly, that degree distribution is not a true power law but rather a log-periodic distribution consisting of a sequence of atoms at points 3 × 2 n and a power-law envelope. This means that the network is scale-invariant only with respect to certain discrete renormalizations and thus do not have the full set of properties of a true power law distribution; see [30] for a recent discussion. One natural generalization is a random Apollonian network [31,32,33], which is constructed, instead of a regular generation-by-generation process, by sequential partitioning of arbitrarily chosen triangles. The average degree distribution in such network is a true power law with exponent γ R = 3 [32]. Notably, to the best of our knowledge, random Apollonian graphs remain the only scale-free planar graph model with a continuously growing size for which the exact degree distribution exponent is known. Another way of generalizing the network is to consider the k-simpliceswith k > 3 as building blocks of the network construction procedure. This gives rise to multidimensional Apollonian networks [34,35].
In this paper, the authors suggest another way of generalizing the Apollonian network construction. As a result, a novel famility of small-world, scale-free planar networks is obtained. The main idea is to construct an Apollonian-style iteration procedure based on polygons with different numbers of edges. The paper is organized as follows. In Section 2, a tetragon-based Apollonian-style network is constructed and the corresponding degree distribution is explicitely calculated. Then, the suggested procedure is generalized to polygons with an arbitrary even number of edges. In Section 4, it is shown that it is possible to construct a continuous one-parametric family of models interpolating between the tetragon- and hexagon-based models and demonstrate that the models in this family have a power low degree distribution with an exponent depending on the parameter, so it is possible to adjust it to fit the desired degree distribution (note that the adjustable exponent of the degree distribution can be obtained by different means in the so-called Evolving Apollonian networks [26,33]). Last section, summarizes the results of the paper and discusses further open questions and possible generalizations.

2. Tetragon-Based Network

2.1. Definition

Among several possible ways of generalizing the procedure described above to the case of polygons, consider the following procedure defined here for tetragons but easy to generalize for any polygon with an even number of edges. Note that given that we make such a generalization in further sections, we prefer to use the term ‘tetragon’ rather than ‘quadrilateral’ for a polygon with 4 sides in order to make the terminology more uniform.
Take a tetragon and pick a point inside it; then choose (at random) a pair of non-adjacent vertices of the original tetragon and connect them with a polyline with one new intermediate point. One now has two adjacent tetragons, for which one can repeat this construction, as shown in Figure 2. Importantly, contrary to the standard Apollonian network, which is a deterministic object, the network resulting from this procedure is stochastic. Indeed, already in the second generation, there are three topologically different realizations of the network, see Figure 2B. Notably, at any generation, this network has no triangles and is, in fact, bipartite.

2.2. Degree Distribution

The natural question to ask about this newly introduced class of planar Apollonian-like networks is what is the degree distribution of the nodes G n ( k ) in the n-th generation of the network and what distribution G ( k ) it converges to for n (here and in what follows the term “degree distribution” is used to mean the probability density function, i.e., the probability for the node to have a degree equal to k, as opposed to the cumulative distribution function, i.e., the probability for the node to have a degree larger or equal to k). By analogy with the Apollonian networks, one expects G ( k ) to be scale-free, i.e.,
G ( k ) C k α , k 1 ,
with some yet unknown constants C and α .
To calculate the degree distribution G n ( k ) , note that the degree of any given node is a random variable, whose distribution F n m ( k ) for all nodes except the four initial ones depends only on the number of generations between the generation m at which it was created and current generation n. Indeed, each node with degree k has exactly k adjacent tetragons ( k 1 for the four initial nodes), and at every step of the recurrent procedure, each of these tetragons is split in two, which results in the creation of a new edge adjacent to the node with the probability 1/2 (in the other half of the cases, the splitting path does not go through the given node). These splitting events happen independently for all tetragons. The overall degree distribution is therefore calculated by averaging over degree distributions of different generations:
G n all ( k ) = 4 F n ( 0 ) ( k ) + m = 1 n Q m F n m ( k ) 4 + m = 1 n Q m = 4 2 n + 3 F n ( 0 ) ( k ) + 2 n 1 2 n + 3 G n ( k ) , G n ( k ) = m = 1 n Q m F n m ( k ) m = 1 n Q m ,
where Q m = 2 m 1 is the number of nodes created in the m-th generation, F n ( 0 ) ( k ) is the degree distribution of the four initial nodes, and it is convenient to introduce G n ( k ) , the degree distribution of all nodes except four initial ones.

2.3. Recurrence Relation for F n ( K )

To construct the recurrence relation for F n ( k ) proceed as follows. Let l be the degree of a node in the ( n 1 ) -th generation. This means that this node has l tetragons adjacent to it, and when constructing the n-th generation of the network, each of them will be split in half, and with probability 1/2, the splitting path will go through the node under consideration. Every such path increases the degree of the node by one. Thus, the overall degree may increase by l , 0 l l , with the probability 2 l l l , leading to
F n ( k ) = l = ( k + 1 ) / 2 k 2 l l k l F n 1 ( l ) for n 1 , F 0 ( k ) = δ k , 2 ,
where the fact that all nodes are created with degree 2 is taken into account and the notation x is introduced for the integer part of x (i.e., greatest integer less or equal to x). This equation can be written down in a simpler form in terms of a generating function,
f n ( λ ) = k = 2 λ k F n ( k ) .
Indeed, after substituting Equation (3), one obtains f 0 ( λ ) = λ 2 and
f n ( λ ) = k = 2 l = ( k + 1 ) / 2 k λ k 2 l l k l F n 1 ( l ) = l = 2 m = 0 l ( λ / 2 ) l l m λ m F n 1 ( l ) = l = 2 λ ( 1 + λ ) 2 l F n 1 ( l ) = f n 1 λ ( 1 + λ ) 2 ,
where the order of summation is changed, m = k l is introduced, and the binomial formula,
( 1 + λ ) l = m = 0 l l m λ m ,
is used. The recurrence relation for the four initial nodes is a bit different because in their case, the node of degree l has only l 1 adjacent tetragons:
F n ( 0 ) ( k ) = l = k / 2 + 1 k 2 l + 1 l 1 k l + 1 F n 1 ( 0 ) ( l ) for n 1 , F 0 ( k ) = δ k , 2 ,
which leads to the following equation for the generating function
f n ( 0 ) ( λ ) = k = 2 λ k F n ( 0 ) ( k ) = 2 1 + λ f n 1 ( 0 ) λ ( 1 + λ ) 2 ,
In the n limit, both f n ( λ ) and f n ( 0 ) ( λ ) converge to zero for all | λ | < 1 . Indeed, the probability to have any finite degree many generations after the creation of a node is exponentially small.

2.4. Generating Function of the Degree Distribution

Combining Equation (2) for G n ( k ) , G n all ( k ) and the equations for the generating functions (5) and (8), one gets the recurrence relation for the full degree distributions in terms of generating functions:
g n ( λ ) = k G n ( k ) λ k , g n all ( λ ) = k G n all ( k ) λ k .
For g n ( λ ) , one gets:
( 2 n + 1 1 ) g n + 1 ( λ ) = 2 n λ 2 + ( 2 n 1 ) g n λ ( 1 + λ ) 2 for n 0 ; g 0 ( λ ) = λ 2 ,
which in the limit of large n reduces to
g n + 1 ( λ ) = 1 2 λ 2 + 1 2 g n λ ( 1 + λ ) 2 .
Contrary to Equation (5) and (8), Equation (10) has a non-trivial limiting solution for n . Indeed, if G n ( k ) converges to a limiting form G ( k ) , then
g ¯ ( λ ) = G ( k ) λ k = lim n G n ( k ) λ k = lim n G n ( k ) λ k = lim n g n ( λ ) .
where the summation and the limit are transposed, as one can do for convergent positive series. Thus, g ¯ ( λ ) is a solution of the functional equation,
2 g ¯ ( λ ) = λ 2 + g ¯ λ ( 1 + λ ) 2 , g ¯ ( λ ) = k G ( k ) λ k .
It seems impossible to solve this equation for all λ ; however, it is possible to extract most important information about G ( k ) directly from the equation. Indeed, the behavior of the distribution for the small and large k is controlled by the behavior of the generating function in the vicinity of λ = 0 and λ = 1 , respectively. For small λ , substituting
g ¯ ( λ ) = p 2 λ 2 + p 3 λ 2 + p 4 λ 4 +
into Equation (13), one obtains:
p 2 = 4 7 , p 3 = 16 105 , p 4 = 16 155 , p 5 = 64 1519 , etc .
for the limiting probabilities of having a node degree equal to 2, 3, 4, 5, ⋯.
In turn, from the behavior of g ¯ ( λ ) in the vicinity of λ = 1 , one can extract both the value of α defined in Equation (1) and the values of existing moments of G ( k ) . Indeed, define
G reg ( k ) = G ( k ) C k α ,
and, accordingly,
g reg ( λ ) = g ¯ ( λ ) C Li α ( λ ) .
Equation (1) implies that G reg ( k ) converges to zero with growing k faster than k α , which guarantees that g reg ( λ ) is less singular then the first one in the vicinity of λ = 1 . Thus, at λ 1 , the function g ¯ ( λ ) has a singularity of the type ( 1 λ ) α 1 and has smooth derivatives up to the order α 1 . Thus, in the lowest orders in ϵ = 1 λ ,
g ¯ ( λ ) = i = 0 α 1 a i ϵ i + C Γ ( 1 α ) ϵ α 1 + o ( ϵ α 1 ) ,
where values of a i depend on the small-k behavior of G ( k ) and contain information about the momenta of the distribution:
a 0 = k G ( k ) ; a 1 = k k G ( k ) = k , etc .
Now, substituting λ ( λ + 1 ) / 2 = 1 3 ϵ / 2 + ϵ 2 / 2 into Equation (13) and equating coefficients in front of different powers of ϵ , one obtains:
2 a 0 = 1 + a 0 , a 0 = 1 , 2 a 1 = 2 + 3 a 1 / 2 , a 1 = 4 , 2 C Γ ( 1 α ) = C Γ ( 1 α ) ( 3 / 2 ) α 1 , α = 1 + ln 2 ln ( 3 / 2 ) = ln 3 ln 3 ln 2 2.70951
Thus, α 1 = 1 , and only zeroth and first moments of the distribution converge:
k G ( k ) = a 0 = 1 ; k = a 1 = 4 ,
while all the higher moments, starting from k 2 , diverge with growing n.
It is instructive to calculate the exact values of moments k n , k 2 n for all finite n. To do this, note that
k n = d g n ( λ ) d λ λ = 1 ; k 2 n = k n + d 2 g n ( λ ) d λ 2 λ = 1 .
Equation (10) implies
( 2 n + 1 1 ) g n + 1 ( λ ) = 2 n + 1 λ + ( 2 n 1 ) 2 λ + 1 2 g n λ ( 1 + λ ) 2 ,
which, for λ = 1 , leads to
1 2 n 1 k n + 1 = 1 + 3 4 1 2 n k n .
Substituting
b n = 4 1 2 n k n ,
and allowing for the initial condition k 1 = 2 , b 1 = 3 , one obtains:
b n = 3 4 b n 1 = 4 3 4 n
and, thus,
k n = 4 1 ( 3 / 4 ) n 1 ( 1 / 2 ) n .
This is the average degree of all nodes except the original four at the n-th step of the network-generation process. Given Equation (2), one obtains the average degree of all nodes:
k n all = k G n all ( k ) = ( 2 n 1 ) k n + 4 ( 1 + ( 3 / 2 ) n ) 2 n + 3 = 4 ( 1 ( 3 / 4 ) n ) 2 n + 1 + ( 3 / 2 ) n 2 n + 3 = 4 2 n + 1 2 n + 3 .

2.5. Second Moment of the Finite-Generation Distribution

To calculate the second moment, take the second derivative of Equation (10):
( 2 n + 1 1 ) g n + 1 ( λ ) = 2 n + 1 + ( 2 n 1 ) g n λ ( 1 + λ ) 2 + ( 2 n 1 ) 2 λ + 1 2 2 g n λ ( 1 + λ ) 2
and take into account Equation (55). Substituting λ = 1 and allowing for the fact that
g n ( 1 ) = k n ; g n ( 1 ) = k 2 n k n
leads to
( 2 n + 1 1 ) ( k 2 n + 1 k n + 1 ) = 2 n + 1 + ( 2 n 1 ) k n + 9 4 ( 2 n 1 ) ( k 2 n k n )
or
( 2 n + 1 1 ) k 2 n + 1 9 4 ( 2 n 1 ) k 2 n = 2 n + 1 + ( 2 n + 1 1 ) k n + 1 5 4 ( 2 n 1 ) k n ,
which, after substituting Equation (27), simplifies to
( 2 n + 1 1 ) k 2 n + 1 9 4 ( 2 n 1 ) k 2 n = 2 n 5 3 4 n .
Now define the sequence
a n = k 2 n 2 n 1 2 n
and its generating function F ( s ) = a n s n . The recurrency for a n reads:
a n + 1 = 9 8 a n + 5 2 1 2 3 4 n
for n 1 , and a 1 = 2 . Then:
F ( s ) 1 9 8 s = 2 s + 5 2 s 2 1 s 3 8 s 2 1 3 s / 4 .
In order to proceed further, note that
2 s 1 9 s / 8 = 16 9 + 16 9 1 1 9 s / 8 , s 2 ( 1 s ) ( 1 9 s / 8 ) = 8 9 8 1 1 s + 64 9 1 1 9 s / 8 , s 2 ( 1 3 s / 4 ) ( 1 9 s / 8 ) = 32 27 32 9 1 1 3 s / 4 + 64 27 1 1 9 s / 8 .
Thus,
F ( s ) = 56 3 1 1 9 s / 8 20 1 1 s + 4 3 1 1 3 s / 4 ; a n = 56 3 9 8 n 20 + 3 4 n 1 ,
and
k 2 = ( 1 2 n ) 1 56 3 9 8 n 20 + 3 4 n 1 .
Proceeding in the same way, one obtains:
k 2 ( 0 ) = 4 3 9 4 n + 5 3 3 2 n + 1 .
Thus, the total average degree is:
k 2 all = 2 n 1 2 n + 3 k 2 + 4 2 n + 3 k 2 ( 0 ) = 1 + 3 2 n 1 24 9 8 n 20 + 8 3 4 n + 4 1 2 n 24 9 8 n 20 ,
where the approximal equality holds for n 1 . Comparing Equations (39) and (41) shows that, interestingly, the four initial nodes contribute a finite fraction to the overall value of k 2 all , which converges to 2/9 for large n.

2.6. Scaling Form of the Degree Distribution

For large n, the degree distribution G n ( k ) converges to G ( k ) . Typically (see, e.g., [36]), one expects the ratio of these functions Φ n ( k ) = G n ( k ) / G ( k ) to attain a universal shape for large n. More precisely, it means that there exists a scaling function ϕ ( x ) and a sequence K n for which
Φ n ( k ) ϕ k / K n 1 for n .
Here, K n is the characteristic scale of the distribution of the n-th generation graph, and it diverges as n . One must demand ϕ ( 0 ) = 1 in order for G n ( k ) to converge to G ( k ) and ϕ ( ) = 0 in order for k to be bound from above for any finite n. There are various ways of obtaining the scaling factor K n , e.g., one can use the large-n behavior of k 2 n :
k 2 n = k 2 G n ( k ) K n ( 3 α ) ,
where it is taken into account that (contrary to the lower moments) k 2 n is controlled by the tail of the distribution. Substituting Equation (41), one obtains:
K n 3 α = ( 9 / 8 ) n ; K n = ( 3 / 2 ) n ,
which is to be expected since the typical maximal degree of the network increases by a factor 3 / 2 on each step.
In order to check the predictions of the model, 2 × 10 5 realizations of the networks of up to the 14th generation were generated. Figure 3a shows the resulting degree distributions for sequential generations of the network. It can be seen from Figure 3b that after renormalization of the abscisse and ordinate axes by the factors ( 3 / 2 ) n and k α , respectively, the data collapse perfectly on a single scaling curve ϕ ( x ) .

3. Polygon-Based Networks for Polygons with Any Even Number of Edges

The procedure suggested in Section 2 can be easily generalized for any even number of edges 2 m ( m 2 ) in the generating polygon; see generalization for m = 3 in Figure 4. This procedure results in a sequence of planar scale-free network models with degree distributions converging to
G m ( k ) C m k α m , k 1 ,
with m-dependent exponents α m . At each generation, each polygon is split by a path connecting directly opposite nodes. There are m different ways of such a splitting, so each node of a polygon participates in the splitting with probability 1 / m . This allows generalizing the recurrence relations (3) and (7) for the degree distribution of a node n generations after its creation in the following way:
F n , m ( 0 ) ( k ) = l = k / 2 + 1 k l 1 k l + 1 1 m k l + 1 m 1 m 2 l k 2 F n 1 , m ( 0 ) ( l ) for n 1 F 0 ( 0 ) ( k ) = δ k , 2 ,
for the original 2 m nodes, and
F n , m ( k ) = l = ( k + 1 ) / 2 k l k l 1 m k l m 1 m 2 l k F n 1 , m ( l ) for n 1 ; F 0 , m ( k ) = δ k , 2 ;
for all the rest. The number of nodes created at n-th generation ( n 1 ) is ( m 1 ) 2 n . Proceeding in exactly the same way as before, gives
f n , m ( 0 ) ( λ ) = k = 2 λ k F n ( 0 ) ( k ) = m m 1 + λ f n 1 ( 0 ) λ ( m 1 + λ ) m ,
g n , m ( λ ) = 1 ( m 1 ) ( 2 n 1 ) l = 1 n k = 2 ( m 1 ) 2 l F n l , m ( k ) λ k = 2 n 1 2 n 1 λ 2 + 2 n 1 1 2 n 1 g n 1 , m λ ( m 1 + λ ) m for n 1 ; g 1 , m ( λ ) = λ 2 ,
and
g n , m all ( k ) = 2 m ( m 1 ) 2 n + m + 1 f n , m ( 0 ) ( λ ) + ( m 1 ) ( 2 n 1 ) ( m 1 ) 2 n + m + 1 g n , m ( λ )
for the generating functions of the degree distributions of the original, newly created and all nodes of the network, f n , m ( 0 ) ( λ ) , g n , m ( λ ) and g n , m all ( λ ) , respectively. The limiting function,
g ¯ m ( λ ) = lim n g n , m ( λ ) ,
satisfies
g ¯ m ( λ ) = λ 2 2 + 1 2 g ¯ m λ ( m 1 + λ ) m ,
and its behavior is easy to analyze both in the vicinity of λ = 0 and λ = 1 . Expanding Equation (52) for small λ , one obtains:
p 2 m = m 2 K 2 , m , p 3 m = 2 m 3 ( m 1 ) K 2 , m K 3 , m , p 4 m = m 3 ( 2 m 3 + 5 ( m 1 ) 3 ) l = 2 4 K l , m , p 5 m = m 5 ( m 1 ) 2 ( 12 m 4 + 8 m 3 ( m 1 ) + 14 ( m 1 ) 4 ) l = 2 5 K l , m , etc . ,
where a short-hand notation, K l , m = 2 m l ( m 1 ) l , is introduced.
In turn, in the vicinity of λ = 1 g ¯ m ( λ ) , it takes the form of Equation (18). Substituting the ansatz (18) into Equation (52), one obtains:
2 a 0 = 1 + a 0 , a 0 = 1 , 2 a 1 = 2 + ( m + 1 ) a 1 / m , a 1 = 2 m / ( m 1 ) , 2 a 2 = 1 a 1 / m + ( m + 1 ) 2 a 2 / m 2 , a 2 = ( m + 1 ) m 2 ( m 1 ) ( m 2 2 m 1 ) , 2 C Γ ( 1 α m ) = C Γ ( 1 α m ) ( m + 1 / m ) α m 1 , α m = 1 + ln 2 ln ( m + 1 ) ln m .
Thus, for any m 3 α m 1 2 , the second moment of G m ( k ) converges. The moments are controlled by the coefficients a i :
k G ( k ) = a 0 = 1 , k k G ( k ) = a 1 = 2 m m 1 , k k 2 G ( k ) = 2 a 2 a 1 = 2 m ( 2 m 2 m 1 ) ( m 1 ) ( m 2 2 m 1 ) .
Since the second moment of G m ( k ) is now controlled by the values of distribution at small k, the initial 2 m nodes do not contribute to the second moment.
Once again, in order to check the predictions 2 × 10 5 realizations of the networks of up to the 12th generation were generated. The results are shown in Figure 5a, and in Figure 5b the plot in the renormalized coordinates is shown. The collapse of the data onto a single master curve is apparent, although it is somewhat worse than in Figure 3b. Presumably, this happens because the typical degrees in the hexagon-based network are much smaller than in the tetragon-network of the same size, and finite size effects are therefore more important.

4. Polygon-Based Networks with Smoothly Changing Exponent of the Degree Distribution

As a result of Section 3, one now has a sequence of Apollonian-like models that generate planar scale-free networks with a discrete sequence of degree-distribution exponents α m = 1 + ln 2 / ( ln ( m + 1 ) ln m ) , m = 2 , 3 , Is it possible to further generalize the model to make α change continuously and take any intermediate values, including, for example, α = 3 , corresponding to the point where the second moment of the degree distributions diverges for the first time?
It turns out that this is indeed possible. One way to do that is as follows. Assume that when introducing a new shortcut dividing a polygon in two, one makes the resulting partition to be a pair of tetragons with probability p and a pair of hexagons with probability 1 p . That is to say, if the original polygon is a tetragon, then with probability p introduce a 2-step path connecting opposite vertices, and with probability q = 1 p , a 4-step path; if the original polygon is a hexagon, the new path connecting two opposite vertices is a 1-step path with probability p and a 3-step path with probability q.
Here, we restrict ourselves to this simplest construction, although it is possible to create more complicated rules. For example, one can introduce correlations between generations in a Markovian way so that there is a matrix p i j of probabilities for a tetragon to give birth to a couple of tetragons, a hexagon to give birth to a couple of tetragons, etc. As a result, it might be possible to construct a network that is, for example, tetragon-dominated at large scales (early generations) and hexagon-dominated at small scales (later generations).
Once again, consider a node, which is created at generation n 0 , and let us study its degree distribution at generation n 0 + n . This degree distribution depends only on n and on the number of edges of the initial two faces adjacent to the node where tetragons or hexagons.
Let the average fraction of tetragons at a given generation be p and the fraction of hexagons be q = 1 p . Then, for each face adjacent to a given node, the probability that this face is a tetragon is
π ( p ) = 4 p 4 p + 6 q = 2 p 3 p .
Assume now that different faces adjacent to a node are tetragons (hexagons) independently from each other. Generally speaking, that is not true: when a new edge is created, the two faces on the sides of it have a similar number of edges. However, one might expect that as the degree of the node grows the correlations become less and less relevant. In this approximation, the probability for a node of degree k to have exactly l adjacent tetragons is
k l π l ( 1 π ) k l .
when the next generation is created, a new edge adjacent to the node under consideration is created with probability 1 / 2 for each tetragon face and with probability 1 / 3 for each hexagon face. Therefore, one can write down the following approximate equation for the probability P ( k + r | k , p ) of the node that has degree k + r at the next generation given that it had degree k in the previous one:
P ( k + r | k , p ) = l = 0 k s = 0 r k l π l ( 1 π ) k l l s 1 2 l k l r s 1 3 r s 2 3 k l r + s ,
where the binomial coefficients m n are assumed to be zeros if n > m or n < 0 . Now, introduce the probability F n ( k | p ) for a node to have degree k n generations after its creation, and the corresponding generation function,
f n ( λ | p ) = k = 2 F n ( k ) λ k .
Then, f 0 ( λ | p ) = λ 2 ,
F n ( k ) = k = k + 1 2 k P ( k | k , p ) F n 1 ( k ) ,
and
f n ( λ | p ) = k = 2 k = k + 1 2 k λ k P ( k | k , p ) F n 1 ( k ) = k = 2 F n 1 ( k ) λ k k = k 2 k λ k k P ( k | k , p ) .
Using Equation (58), it is easy to calculate the second sum on the right-hand side:
r = 0 k λ r P ( k + r | k , p ) = r = 0 k l = 0 k s = 0 r k ! s ! ( l s ) ! ( r s ) ! ( k l r + s ) ! π l ( 1 π ) k l 1 2 l 1 3 r s 2 3 k l r + s λ r = r = 0 k l = 0 k s = 0 r k ! s ! ( r s ) ! ( l s ) ! ( k l r + s ) ! π λ 2 s π 2 l s ( 1 π ) λ 3 r s 2 ( 1 π ) 3 k l r + s = π λ 2 + π 2 + ( 1 π ) λ 3 + 2 ( 1 π ) 3 k ,
which leads to the following equation for the generating function:
f n ( λ ) = f n 1 λ 4 π 6 + 2 + π 6 λ = f n 1 2 p + λ 3 p λ ,
where Equation (56) is taken into account to obtain to the last expression. Proceeding as before, one gets the equation for the generating function of the full limiting degree distribution g ( λ ) (except for the initial set of nodes):
g ( λ ) = λ 2 2 + 1 2 g 2 p + λ 3 p λ .
Expanding g ( λ ) for λ = 1 ϵ , ϵ 1 in the form (compare Equation (18)):
g ¯ ( λ ) = i = 0 α ( p ) 1 a i ϵ i + C Γ ( 1 α ( p ) ) ϵ α ( p ) 1 + o ( ϵ α ( p ) 1 ) ,
and equating the coefficients term by term exactly in the same way as in Section 2, one gets the following equation for the degree distribution exponent α ( p ) :
2 = 4 p 3 p α ( p ) 1 , α ( p ) = 1 + ln 2 ln ( 4 p ) ln ( 3 p ) .
Thus, for example, the interesting case α ( p ) = 3 when the second moment of the degree distribution diverges for the first time corresponds to
p | α = 3 = 2 2 0.58579
In Figure 6, the numerical data for the degree distribution of the mixed networks are shown. One can clearly see that the slope of the distribution gradually changes with changing p. Moreover, after rescaling the degree distribution with the power law prescribed by Equation (66), all curves are approximately flat for small k, validating the approximation of independent phases.

5. Concluding Remarks and Open Questions

This paper presents one possible class of planar random graphs constructed from polygons with an even number of edges using a procedure similar to the construction of Apollonian graphs [18]. The 2 m -polygon-based graphs have a limiting power law degree distribution with the exponent,
α m = 1 + ln 2 ln ( m 1 ) ln m ,
and the moments of the degree distribution are given by Equations (21) and (55). The second moment of the degree distribution diverges as ( 9 / 8 ) n with the number of generations n in the case of tetragon-based graphs (see Equation (41)) and converges to a finite value in Equation (55) for the polygons with a larger number of edges. Moreover, as described in Section 4, it is possible to construct a mixed model based on two different polygons (tetragons and hexagons in our example) so that on all stages of construction, tetragons are formed with probability p and hexagons with probability 1 p . By varying p, one can adjust the slope of the degree distribution in order to achieve a desired value in a way reminiscent of evolving Apollonian networks [26].
Clearly, all graph classes presented here are small world. Indeed, the diameter of the graphs grows at most linearly with the number of generations:
d n + 1 d n + 2 m / 2 ,
where d n is the diameter of the n-generation graph. In turn, the total number of nodes grows exponentially with the number of generations; thus, the diameter is, at most, proportional to the logarithm of the number of nodes.
The shortest cycles in the graphs presented here are 2 m , and, in particular, there are no triangles in them, so, generally speaking, the clustering coefficient is zero. However, this should not obscure the fact that there is actually a huge number of short cycles in these graphs. Indeed, consider the following auxiliary construction: let the polygon-based construction be exactly as presented above up to n-th generation, but then connect all the nodes belonging to the same face on the last generation of the procedure, so the smallest faces (i.e., faces constructed on the last step) are considered to be complete graphs K 2 m ( ( 2 m 1 ) -simplices). The large-scale structure of the resulting graph (including, e.g., the slope of the degree distribution) will be the same as in the original polygon-based procedure, but a finite fraction of nodes (those created in the n-th generation of the construction) will have clustering coefficient 1, guaranteeing that the average clustering coefficient of the whole graph remains finite as n . In order to use polygon-based graphs as a toy model for experimental systems, it might be reasonable to add a random fraction of them in order to fit the observed clustering coefficient instead of adding all possible links connecting the vertices in the smallest faces.
Interestingly, polygon-based graphs with even m are bipartite, see Figure 7a. In that sense, the tetragon-based graph seems to be a natural generalization of the Apollonian construction for the case of bipartite graphs. We expect that there might be a connection between bipartite polygon-based graphs and space-filling bearings, which allow only cycles of even lengths [37] in a way similar to the connection between original Apollonian networks and space-filling systems of embedded disks. Exploring this question goes, however, beyond the scope of this paper.
This paper restricts itself to just one class of possible generalizations of the Apollonian construction based on polygons of arbitrary sizes. It is quite easy to suggest various other generalizations. The most obvious example is, probably, the random polygon construction where new graphs are constructed not generation-by-generation by splitting all the polygons of the previous generation at once, but rather by randomly choosing on each step a face to split. Figure 7b presents an example realization of such a tetragon-based random graph. In the standard Apollonian case, it is known that the exponent of the degree distribution is different for the regular and random constructions. Calculating the degree distribution of random Apollonian-like polygon-based graphs remains an interesting open question.
Another, this time a completely deterministic generalization, is as follows. Consider a polygon with an odd number of edges 2 m + 1 , m 1 . Put a point inside the polygon and connect it with all vertices of the polygon by chains of m edges and ( m 1 ) nodes. This splits a polygon into 2 m + 1 faces, each having exactly 2 m + 1 edges. On the next step, repeat this procedure for each of the faces and proceed ad infinitum. Figure 7c shows the second-generation pentagon-based graph obtained via such procedure. Clearly, this construction is an even more direct generalization of the Apollonian graph construction (indeed, m = 1 case is just the Apollonian graph itself). However, it means that it has standard drawbacks of the Apollonian graph in a sense that it is a single deterministic object rather than a stochastic ensemble of graphs and that its limiting degree distribution is not a power law but rather a log-periodic function with a power law envelope.
We think that the classes of graphs presented here are a useful addition to the toolkit of toy models to model scale-free graphs. Indeed, while having the main advantages of the Apollonian networks, they have additional flexibility in a sense that one might regulate the slope of the degree distribution and the average clustering coefficient in the way described above. In particular, such graphs might be, in our opinion, useful in the applications where graph planarity is essential [38], for example, in quantitative geography, such as the study of the formation of the systems of interconnected cities. On the other side, studying percolation, spectral properties, diffusion, synchronization, epidemic spreading, etc., on these generalized graphs might allow to systematically study the influence of varying degree exponents on these phenomena, which, to the best of our knowledge, have not yet been done for the scale-free planar networks.

Author Contributions

Conceptualization, M.V.T.; Investigation, M.V.T., D.G.K. and V.I.S.; Methodology, M.V.T.; Software, V.I.S.; Supervision, M.V.T.; Visualization, M.V.T.; Writing—original draft, M.V.T.; Writing—revised version, M.V.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by RSF, grant number 21-11-00215.

Acknowledgments

The authors are grateful to H. Hermann for encouraging comments on the idea of this work and to G. Bianconi, S. Nechaev and P. Krapivsky for many interesting discussions. The research presented here was undertaken in the Laboratory of Nonlinear, Nonequilibrium and Complex Systems at Moscow State University, which, for the last 10 years, has been headed by Professor M.I. Tribelsky. We are very grateful for his valuable scientific and personal advice, kindness and support over all these years and would like to dedicate this work to him on the occasion of his 70th birthday.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tikhonov, K.S.; Mirlin, A.D. From Anderson localization on random regular graphs to many-body localization. Ann. Phys. 2021, in press. [Google Scholar] [CrossRef]
  2. Sneppen, K.; Zocchi, G. Physics in Molecular Biology; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar] [CrossRef]
  3. Jackson, M.O. Social and Economic Networks; Princeton University Press: Princeton, NJ, USA, 2010. [Google Scholar]
  4. Kenett, Y.N.; Anaki, D.; Faust, M. Investigating the structure of semantic networks in low and high creative persons. Front. Hum. Neurosci. 2014, 8, 407. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Stella, M.; Beckage, N.M.; Brede, M.; Domenico, M.D. Multiplex model of mental lexicon reveals explosive learning in humans. Sci. Rep. 2018, 8, 2259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Valba, O.V.; Gorsky, A.S.; Nechaev, S.K.; Tamm, M.V. Analysis of English free association network reveals mechanisms of efficient solution of Remote Association Tests. PLoS ONE 2021, 16, e0248986. [Google Scholar]
  7. Newman, M. Networks; Oxford University Press: Oxford, UK, 2018. [Google Scholar] [CrossRef]
  8. Barabasi, A.L. Network Science; Cambridge University Press: Cambridge, UK, 2016. [Google Scholar]
  9. Dorogovtsev, S. Complex Networks; Oxford University Press: Oxford, UK, 2010. [Google Scholar] [CrossRef]
  10. Erdős, P.; Rényi, A. On Random Graphs I. Publ. Math. 1959, 6, 290–297. [Google Scholar]
  11. Krapivsky, P.L.; Redner, S.; Ben-Naim, E. A Kinetic View of Statisrical Physics; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar] [CrossRef] [Green Version]
  12. Balister, P.; Bollobás, B.; Sarkar, A. Percolation, connectivity, coverage and colouring of random geometric graphs. In Handbook of Large-Scale Random Networks; János Bolyai Mathematics Society Studies; Bollobás, B., Kozma, R., Miklós, D., Eds.; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar] [CrossRef]
  13. Barabási, A.L.; Albert, R. Emergence of scaling in random networks. Science 1999, 286, 509–512. [Google Scholar] [CrossRef] [Green Version]
  14. Watts, D.; Strogatz, S. Collective dynamics of ‘small-world’ networks. Nature 1998, 393, 440–442. [Google Scholar] [CrossRef]
  15. Krioukov, D.; Papadopoulos, F.; Kitsak, M.; Vahdat, A.; Boguna, M. Hyperbolic geometry of complex networks. Phys. Rev. E 2010, 82, 036106. [Google Scholar] [CrossRef] [Green Version]
  16. Boguna, M.; Papadopoulos, F.; Krioukov, D. Sustaining the internet with hyperbolic mapping. Nat. Commun. 2010, 1, 62. [Google Scholar] [CrossRef] [Green Version]
  17. Papadopoulos, F.; Kitsak, M.; Serrano, M.Á.; Boguná, M.; Krioukov, D. Popularity versus similarity in growing networks. Nature 2012, 489, 537–540. [Google Scholar] [CrossRef] [Green Version]
  18. Andrade, J.S., Jr.; Herrmann, H.J.; Andrade, R.F.S.; da Silva, L. Apollonian networks: Simultaneously scale-free, small world, Euclidean, space filling, and with matching graphs. Phys. Rev. Lett. 2005, 94, 018702. [Google Scholar] [CrossRef] [Green Version]
  19. Doye, J.P.K.; Massen, C.P. Self-similar disk packings as model spatial scale-free networks. Phys. Rev. E 2005, 71, 016128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Kasner, E.; Supnick, F. The Apollonian packing of circles. Proc. Natl. Acad. Sci. USA 1943, 29, 378–384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Boyd, D.W. The residual set dimension of the Apollonian packing. Mathematika 1973, 20, 170–174. [Google Scholar] [CrossRef]
  22. Zhang, Z.; Chen, L.; Zhou, S.; Fang, L.; Guan, J.; Zou, T. Analytical solution of average path length for Apollonian networks. Phys. Rev. E 2008, 77, 017102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Bianconi, G.; Ziff, R.M. Topological percolation on hyperbolic simplicial complexes. Phys. Rev. E 2018, 98, 052308. [Google Scholar] [CrossRef] [Green Version]
  24. Andrade, R.F.S.; Hermann, H.J. Magnetic models on Apollonian networks. Phys. Rev. E 2005, 71, 056131. [Google Scholar] [CrossRef] [Green Version]
  25. Lind, P.G.; da Silva, L.R.; Andrade, J.S., Jr.; Hermann, H.J. Spreading gossip in social networks. Phys. Rev. E 2007, 76, 036117. [Google Scholar] [CrossRef] [Green Version]
  26. Zhang, Z.; Rong, L.; Zhou, S. Evolving Apollonian networks with small-world scale-free topologies. Phys. Rev. E 2006, 74, 046105. [Google Scholar] [CrossRef] [Green Version]
  27. Mendes, G.A.; da Silva, L.R.; Hermann, H.J. Traffic gridlock on complex networks. Phys. A 2012, 391, 362. [Google Scholar] [CrossRef] [Green Version]
  28. Huang, Z.-G.; Xu, X.-J.; Wu, Z.-X.; Wang, Y.-H. Walks on Apollonian networks. Eur. Phys. J. B 2006, 51, 549–553. [Google Scholar] [CrossRef] [Green Version]
  29. Zhang, Z.; Guan, J.; Xie, W.; Qi, Y.; Zhou, S. Random walks on the Apollonian network with a single trap. Europhys. Lett. 2009, 86, 10006. [Google Scholar] [CrossRef]
  30. Newberry, M.G.; Savage, V.M. Self-similar processes follow a power law in discrete logarithmic space. Phys. Rev. Lett. 2019, 122, 158303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Zhou, T.; Yan, G.; Zhou, P.-L.; Fu, Z.-Q.; Wang, B.-H. Random apollonian networks. arXiv 2004, arXiv:0409414. [Google Scholar]
  32. Zhou, T.; Yan, G.; Wang, B.-H. Maximal planar networks with large clustering coefficient and power-law degree distribution. Phys. Rev. E 2005, 71, 046141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Kolossvary, I.; Komjathy, J.; Vago, L. Degrees and distances in random and evolving Apollonian networks. Adv. Appl. Prob. 2016, 48, 865–902. [Google Scholar] [CrossRef] [Green Version]
  34. Zhang, Z.; Comellas, F.; Fertin, G.; Rong, L. High-dimensional Apollonian networks. J. Phys. A 2006, 39, 1811. [Google Scholar] [CrossRef] [Green Version]
  35. Zhang, Z.; Rong, L.; Comellas, F. High-dimensional random Apollonian networks. Phys. A 2006, 364, 610–618. [Google Scholar] [CrossRef] [Green Version]
  36. Setna, J.P. Statistical Mechanics: Entropy, Order Parameters, and Complexity; Oxford University Press: Cambridge, UK, 2006. [Google Scholar] [CrossRef]
  37. Oron, G.; Herrmann, H.J. Generalization of space-filling bearings to arbitrary loop size. J. Phys. A 2000, 33, 1417–1434. [Google Scholar] [CrossRef] [Green Version]
  38. Barthelemy, M. Spatial networks. Phys. Rep. 2011, 499, 1–101. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Apollonian network: (a) first (black) and second (black and red) generations of the Apollonian network; (b) first generation of the same Apollonian network represented as a tetrahedron in 3-dimensional space; (c) Apollonian networks of higher generations can be thought of as rooted trees constructed from adjacent tetrahedra. Here, second generation is shown; the shaded face functions as a root of the tree.
Figure 1. Apollonian network: (a) first (black) and second (black and red) generations of the Apollonian network; (b) first generation of the same Apollonian network represented as a tetrahedron in 3-dimensional space; (c) Apollonian networks of higher generations can be thought of as rooted trees constructed from adjacent tetrahedra. Here, second generation is shown; the shaded face functions as a root of the tree.
Physics 03 00063 g001
Figure 2. Constructionof a tetragon-based network: (A) representatives of the tetragon-networks up to the 4th generation; (B) three possible topologically different realizations of the second-generation network; (C) random triangulation of the fourth-generation network.
Figure 2. Constructionof a tetragon-based network: (A) representatives of the tetragon-networks up to the 4th generation; (B) three possible topologically different realizations of the second-generation network; (C) random triangulation of the fourth-generation network.
Physics 03 00063 g002
Figure 3. (a) Degree distributions of the tetragon-based networks of generations 5–14 as indicated. The results shown are obtained after averaging over 2 × 10 5 realizations and logarithmic binning with step 1.1. (b) Same distributions but rescaled as the axes show.
Figure 3. (a) Degree distributions of the tetragon-based networks of generations 5–14 as indicated. The results shown are obtained after averaging over 2 × 10 5 realizations and logarithmic binning with step 1.1. (b) Same distributions but rescaled as the axes show.
Physics 03 00063 g003
Figure 4. Construction of the hexagon-based network (up to 4th generation).
Figure 4. Construction of the hexagon-based network (up to 4th generation).
Physics 03 00063 g004
Figure 5. (a) Degree distributions of the hexagon-based networks of generations 4–12 as indicated. The results shown are obtained after averaging over 2 × 10 5 realizations and logarithmic binning with step 1.1. (b) Same distributions but rescaled as the axes show.
Figure 5. (a) Degree distributions of the hexagon-based networks of generations 4–12 as indicated. The results shown are obtained after averaging over 2 × 10 5 realizations and logarithmic binning with step 1.1. (b) Same distributions but rescaled as the axes show.
Physics 03 00063 g005
Figure 6. (a) Degree distributions of n = 10 th-generation mixed tetragon–hexagon networks with varying p, rainbow changing color from p = 0.9 (red) to p = 0.1 (violet). The results shown are obtained after averaging over 10 5 realizations and logarithmic binning with step 1.05. (b) Same distributions but rescaled: P ( k ) is rescaled by its theoretical behavior k α , with α given by Equation (66), k is rescaled by a factor k 0 ( p , n ) = ( ( 4 p ) / ( 3 p ) ) n , which approximates the growth of the maximal accessible degree with the number of generations n.
Figure 6. (a) Degree distributions of n = 10 th-generation mixed tetragon–hexagon networks with varying p, rainbow changing color from p = 0.9 (red) to p = 0.1 (violet). The results shown are obtained after averaging over 10 5 realizations and logarithmic binning with step 1.05. (b) Same distributions but rescaled: P ( k ) is rescaled by its theoretical behavior k α , with α given by Equation (66), k is rescaled by a factor k 0 ( p , n ) = ( ( 4 p ) / ( 3 p ) ) n , which approximates the growth of the maximal accessible degree with the number of generations n.
Physics 03 00063 g006
Figure 7. (a) Tetragon-based networks are bipartite. (b) An example of a particular realization of a randon tetragon-based graph. (c) An example of a 2nd-generation deterministic pentagon-based graph.
Figure 7. (a) Tetragon-based networks are bipartite. (b) An example of a particular realization of a randon tetragon-based graph. (c) An example of a 2nd-generation deterministic pentagon-based graph.
Physics 03 00063 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tamm, M.V.; Koval, D.G.; Stadnichuk, V.I. Polygon-Based Hierarchical Planar Networks Based on Generalized Apollonian Construction. Physics 2021, 3, 998-1014. https://doi.org/10.3390/physics3040063

AMA Style

Tamm MV, Koval DG, Stadnichuk VI. Polygon-Based Hierarchical Planar Networks Based on Generalized Apollonian Construction. Physics. 2021; 3(4):998-1014. https://doi.org/10.3390/physics3040063

Chicago/Turabian Style

Tamm, Mikhail V., Dmitry G. Koval, and Vladimir I. Stadnichuk. 2021. "Polygon-Based Hierarchical Planar Networks Based on Generalized Apollonian Construction" Physics 3, no. 4: 998-1014. https://doi.org/10.3390/physics3040063

APA Style

Tamm, M. V., Koval, D. G., & Stadnichuk, V. I. (2021). Polygon-Based Hierarchical Planar Networks Based on Generalized Apollonian Construction. Physics, 3(4), 998-1014. https://doi.org/10.3390/physics3040063

Article Metrics

Back to TopTop