Next Article in Journal
Reduced Complexity Detection of Narrowband Secondary Synchronization Signal for NB-IoT Communication Systems
Next Article in Special Issue
Extremal Trees with Respect to the Difference between Atom-Bond Connectivity Index and Randić Index
Previous Article in Journal
Subjective Product Evaluation System Based on Kansei Engineering and Analytic Hierarchy Process
Previous Article in Special Issue
Computing Analysis of Connection-Based Indices and Coindices for Product of Molecular Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Properties of General Indices on Random Networks

by
R. Aguilar-Sánchez
1,†,
I. F. Herrera-González
2,†,
J. A. Méndez-Bermúdez
3,4,*,† and
José M. Sigarreta
5,†
1
Facultad de Ciencias Químicas, Benemérita Universidad Autónoma de Puebla, Puebla 72570, Mexico
2
Departamento de Ingeniería, Universidad Popular Autónoma del Estado de Puebla, Puebla 72410, Mexico
3
Departamento de Matemática Aplicada e Estatística, Instituto de Ciências Matemáticas e de Computação, Universidade de São Paulo—Campus de São Carlos, Caixa Postal 668, São Carlos 13560-970, Brazil
4
Instituto de Física, Benemérita Universidad Autónoma de Puebla, Apartado Postal J-48, Puebla 72570, Mexico
5
Facultad de Matemáticas, Universidad Autónoma de Guerrero, Carlos E. Adame No.54 Col. Garita, Acapulco 39650, Mexico
*
Author to whom correspondence should be addressed.
The authors contributed equally to this work.
Symmetry 2020, 12(8), 1341; https://doi.org/10.3390/sym12081341
Submission received: 30 June 2020 / Revised: 31 July 2020 / Accepted: 7 August 2020 / Published: 11 August 2020
(This article belongs to the Special Issue Analytical and Computational Properties of Topological Indices)

Abstract

:
We perform a detailed (computational) scaling study of well-known general indices (the first and second variable Zagreb indices, M 1 α ( G ) and M 2 α ( G ) , and the general sum-connectivity index, χ α ( G ) ) as well as of general versions of indices of interest: the general inverse sum indeg index I S I α ( G ) and the general first geometric-arithmetic index G A α ( G ) (with α R ). We apply these indices on two models of random networks: Erdös–Rényi (ER) random networks G ER ( n ER , p ) and random geometric (RG) graphs G RG ( n RG , r ) . The ER random networks are formed by n ER vertices connected independently with probability p [ 0 , 1 ] ; while the RG graphs consist of n RG vertices uniformly and independently distributed on the unit square, where two vertices are connected by an edge if their Euclidean distance is less or equal than the connection radius r [ 0 , 2 ] . Within a statistical random matrix theory approach, we show that the average values of the indices normalized to the network size scale with the average degree k of the corresponding random network models, where k ER = ( n ER 1 ) p and k RG = ( n RG 1 ) ( π r 2 8 r 3 / 3 + r 4 / 2 ) . That is, X ( G ER ) / n ER X ( G RG ) / n RG if k ER = k RG , with X representing any of the general indices listed above. With this work, we give a step forward in the scaling of topological indices since we have found a scaling law that covers different network models. Moreover, taking into account the symmetries of the topological indices we study here, we propose to establish their statistical analysis as a generic tool for studying average properties of random networks. In addition, we discuss the application of specific topological indices as complexity measures for random networks.

1. Introduction

Harold Wiener can be considered the pioneer in the study of the so-called topological indices. His first investigations in this subject appeared in 1947 [1]. There, he introduced the today known as Wiener index in order tu analyze and correlate physicochemical properties of alkanes. The topological indices associate a numerical value to topological properties of discrete structures, based on its invariants. The initial goal of the study of topological indices was to establish quantitative relationships between different structures and properties of molecules. Winner’s work did not produce an immediate impact. Note that it was about 30 years later, in 1971, that Haruo Hosoya continued the research on topological indices by introducing the Z index [2]; currently known as the Hosoya Z index, which has been successfully applied in establishing quantitative structure–property relationships (QSPRs) and quantitative structure–activity relationships (QSARs).
It is also fair to recognize the relevance of the contributions made by two important groups in the study and development of topological indices. On the one hand, the extinct group of the Ruder Boskovic Institute in Zagreb, where the mathematical and computational properties of the well-known Zagreb indices were first studied. In fact, the first and second Zagreb indices appeared for the first time in 1972, a year later than the Hosoya Z index, in a paper by Gutman and Trinajstić [3,4]. On the other hand, the group leaded by Milan Randić, whose index, the Randić index (introduced in 1975 [5]), is presumably the most studied and applied topological index.
Later on, the concept of variable molecular descriptors was proposed as a way of characterizing heteroatoms in molecules (see [6,7]). The essential idea is that the variables are determined during the regression; this allows making the standard error of the estimate for a particular property (targeted in the study) as small as possible (see, e.g., [8]).
Given a network G = ( V ( G ) , E ( G ) ) , the first and second variable Zagreb indices are defined as [8,9,10]
M 1 α ( G ) = u V ( G ) ( k u ) α
and
M 2 α ( G ) = u v E ( G ) ( k u k v ) α ,
respectively, with α R , k u is the degree of the vertex u, and u v denotes the edge of the network G connecting the vertices u and v. Note that M 1 2 ( G ) is the first Zagreb index M 1 ( G ) and M 1 1 ( G ) is the inverse degree index I D ( G ) . Furthermore, notice that M 2 1 / 2 ( G ) is the usual Randić index R ( G ) and M 2 1 ( G ) is the second Zagreb index M 2 ( G ) . The second variable Zagreb index has found applications in the structure-boiling point modeling of benzenoid hydrocarbons [11].
The general sum-connectivity index was defined in [12] as
χ α ( G ) = u v E ( G ) ( k u + k v ) α .
Some relations of this and other indices are reported in [13]. Note that χ 1 ( G ) is the first Zagreb index M 1 ( G ) , 2 χ 1 ( G ) is the harmonic index H ( G ) and χ 1 / 2 ( G ) is the sum-connectivity index χ ( G ) .
Moreover, here we introduce two new general indices: the general inverse sum indeg index and the general first geometric-arithmetic index, that we define as
I S I α ( G ) = u v E ( G ) k u k v k u + k v α
and
G A α ( G ) = u v E ( G ) k u k v 2 ( k u + k v ) α ,
respectively. Evidently, I S I 1 ( G ) and G A 1 ( G ) are the standard inverse sum indeg index I S I ( G ) and first geometric-arithmetic index G A ( G ) , respectively. We note that a general inverse sum indeg index similar to ours was recently proposed in [14]. Furthermore, for reviews on geometric-arithmetic indices and a new variant of the geometric-arithmetic index (to get a better estimates of the boiling point of alkanes) see Refs. [15,16] and [17], respectively. We stress that the goal of introducing these two general indices is to show that the statistical approach described below is so general that can be applied to any general index, even these new ones with still unknown properties. Thus, we will not focus on their mathematical properties here. Instead we discuss the application of specific topological indices as complexity measures for random networks (see Section 4) and leave the mathematical study and other possible applications for future research.
Here we consider two models of random networks G: Erdös–Rényi (ER) random networks G ER ( n ER , p ) and random geometric (RG) graphs G RG ( n RG , r ) . The ER random networks, introduced by Solomonoff and Rapoport [18] and investigated later in great detail by Erdös and Rényi [19,20,21], are formed by n ER vertices connected independently with probability p [ 0 , 1 ] . While the RG graphs [22,23], often used to study the structure and dynamics of spatially embedded complex systems, consist of n RG vertices uniformly and independently distributed on the unit square, where two vertices are connected by an edge if their Euclidean distance is less or equal than the connection radius r [ 0 , 2 ] . We stress that the statistical study of topological indices we perform here is justified by the random nature of the network models we want to explore. Since a given parameter pair ( n ER , p ) [ ( n RG , r ) ] represents an infinite-size ensemble of ER random networks (RG graphs), the computation of a topological index on a single network is irrelevant. In contrast, the computation of a given topological index on a large ensemble of random networks, all characterized by the same parameter pair ( n ER , p ) [ ( n RG , r ) ], may provide useful average information about the full ensemble. This statistical approach, well known in random matrix theory studies, is not widespread in studies involving topological indices, mainly because topological indices are not commonly applied to random networks; for very recent exceptions see [24,25].
Therefore, the purpose of this work is twofold. First, we propose to establish the statistical analysis of topological indices as a generic tool for studying average properties of random networks, and second, we perform for the first time (to our knowledge), a scaling study of general topological indices on random networks.
Our results are presented separately for ER random networks and RG graphs in Section 2 and Section 3, respectively. Finally, the discussion and conclusions are drawn in Section 4.

2. Erdös–Rényi Random Networks

In what follows, we present the average values of the indices listed in Equations (1)–(5). All averages are computed over ensembles of 10 7 / n ER ER random networks characterized by the parameter pair ( n ER , p ) .

2.1. Computational Properties of General Indices on Erdös–Rényi Random Networks

In Figure 1a, Figure 2a, Figure 3a, Figure 4a and Figure 5a we present, respectively, the average first variable Zagreb index M 1 α ( G ER ) , the average second variable Zagreb index M 2 α ( G ER ) , the average general sum-connectivity index χ α ( G ER ) , the average general inverse sum indeg index I S I α ( G ER ) and the average general first geometric-arithmetic index G A α ( G ER ) as a function of the probability p of ER random networks of size n ER = 1000 . There, we show curves for α [ 2 , 2 ] . As a reference, red curves correspond to α = 0 while blue curves indicate α = 1 .
From these figures we observe that:
(i) 
The curves of M 1 α ( G ER ) , M 2 α ( G ER ) , χ α ( G ER ) and I S I α ( G ER ) show three different behaviors as a function of p depending on the value of α : for α < α 0 , they grow for small p, approach a maximum value and then decrease when p is further increased. For α > α 0 , they are monotonically increasing functions of p. For α = α 0 the curves saturate above a given value of p. See Figure 1a, Figure 2a, Figure 3a and Figure 4a.
(ii) 
α 0 = 0 for M 1 α ( G ER ) , where M 1 α 0 ( G ER ) is the average number of non-isolated vertices V × ( G ER ) ; α 0 = 1 / 2 for M 2 α ( G ER ) , where M 2 α 0 ( G ER ) is the average Randic index R ( G ER ) ; α 0 = 1 for χ α ( G ER ) , where 2 χ 1 ( G ER ) is the average Harmonic index H ( G ER ) ; and α 0 = 1 for I S I α ( G ER ) , so I S I 1 ( G ER ) = V × ( G ER ) .
(iii) 
All curves of G A α ( G ER ) grow linearly with p for all α and n ER ; see the magenta dashed line in Figure 5a, plotted to guide the eye.
(iv) 
For α = 0 , see the red curves in Figure 2a, Figure 3a, Figure 4a and Figure 5a, M 2 0 ( G ER ) , χ 0 ( G ER ) , I S I 0 ( G ER ) and G A 0 ( G ER ) give the average number of edges of the ER random network. That is
M 2 0 ( G ER ) = χ 0 ( G ER ) = I S I 0 ( G ER ) = G A 0 ( G ER ) = n ER ( n ER 1 ) 2 p .
(v) 
When n ER p 1 , we can write k u k v k ER in Equations (1)–(5), with
k ER = ( n ER 1 ) p .
Therefore, for n ER p 1 , the average values of the indices we are computing here are well approximated by:
M 1 α ( G ER ) = n ER [ ( n ER 1 ) p ] α ,
M 2 α ( G ER ) = n ER 2 [ ( n ER 1 ) p ] 1 + 2 α ,
χ α ( G ER ) = n ER 2 1 α [ ( n ER 1 ) p ] 1 + α ,
I S I α ( G ER ) = n ER ( n E R 1 ) p 2 1 + α ,
G A α ( G ER ) = n ER ( n ER 1 ) p 2 1 + 2 α .
In Figure 1b, Figure 2b, Figure 3b, Figure 4b and Figure 5b, we show that Equations (8)–(12) (red-dashed lines) indeed describe well the data (thick black curves) for k ER 10 .

2.2. Scaling Properties of General Indices on Erdös–Rényi Random Networks

Now in Figure 1c, Figure 2c, Figure 3c, Figure 4c and Figure 5c we show X ( G ER ) as a function of the probability p of ER random networks of four different sizes n ER . Here and below, X represents M 1 α , M 2 α , χ α , I S I α and G A α . It is quite clear from these figures that the blocks of curves, characterized by the different network sizes, display similar curves but displaced on both axes. Thus, it is the goal of this section to find the scaling parameters that make the blocks of curves to coincide.
First, we recall that the average degree k ER , see Equation (7), is known to scale both topological and spectral measures applied to ER random networks. In particular, k ER was shown to scale both the normalized Randic index [24] and the normalized Harmonic [25] index on ER random networks.
Second, we observe (not shown here) that the effect of increasing the network size is to displace the blocks of curves X ( G ER ) vs. k ER , characterized by the different network sizes, upwards in the y axis. Moreover, the fact that these blocks of curves, plotted in semi-log scale, are shifted the same amount on the y axis when doubling n ER is a clear signature of scalings of the form X ( G ER ) n ER β k ER . Furthermore, since from Equations (8)–(12) we have that X ( G ER ) n ER k ER , we anticipate that β = 1 for all α .
Therefore, in Figure 1d, Figure 2d, Figure 3d, Figure 4d and Figure 5d we plot X ( G ER ) / n ER as a function of k ER showing that all indices are now properly scaled; i.e., the blocks of curves painted in different colors for different network sizes fall on top of each other. Moreover, following Equations (8)–(12), for k ER 10 we have that
M 1 α ( G ER ) n ER k ER α ,
M 2 α ( G ER ) n ER k ER 1 + 2 α 2 ,
χ α ( G ER ) n ER k ER 1 + α 2 1 α ,
and
I S I α ( G ER ) n ER k ER 1 + α 2 1 + α .
While
G A α ( G ER ) n ER k ER 2 1 + 2 α
is indeed valid for all k ER .

3. Random Geometric Graphs

As in the previous section, here we present the average values of the indices listed in Equations (1)–(5). Again, all averages are computed over ensembles of 10 7 / n RG RG graphs, each ensemble characterized by a fixed parameter pair ( n RG , r ) .

3.1. The Average Degree of Random Geometric Graphs

We want to stress that, in the previous section, the average degree was shown to be the most important quantity in the scaling of the general indices on ER random networks. Thus, since the average degree may also play a prominent role in the scaling of general indices on RG graphs, we focus on it before starting the study of X ( G RG ) . Then, in Figure 6a we present the average index M 1 1 ( G RG ) normalized to the network size (symbols), which in fact corresponds to the average degree k RG :
M 1 1 ( G RG ) n RG = k RG .
Moreover, an expression for k RG was derived in Ref. [26] and reads as
k RG = ( n RG 1 ) f ( r ) ,
where f ( r ) is a quite non-trivial function of r:
f ( r ) = π r 2 8 3 r 3 + 1 2 r 4 .
Indeed, as shown in Figure 6a, Equation (19) (full lines) coincides perfectly with M 1 1 ( G RG ) / n RG (symbols). This fact is of paramount importance to our study because given the similarity of Equation (18) with Equations (13)–(17), we can already interpret k RG as the scaling parameter of the index M 1 1 ( G RG ) . This trivial scaling is shown in Figure 6b where we plot M 1 1 ( G RG ) / n RG versus k RG , as given in Equation (19); there, the curves corresponding to different network sizes collapse on top of the identity function. Moreover, since M 1 1 ( G RG ) is a special case of the average first variable Zagreb index, k RG becomes the perfect candidate to scale the general index M 1 α ( G RG ) and, based on the results of the previous section, also the rest of the general indices we study here.

3.2. Computational Properties of General Indices on Random Geometric Graphs

Now, in Figure 7a, Figure 8a, Figure 9a, Figure 10a and Figure 11a we present, respectively, the average first variable Zagreb index M 1 α ( G RG ) , the average second variable Zagreb index M 2 α ( G RG ) , the average general sum-connectivity index χ α ( G RG ) , the average general inverse sum indeg index I S I α ( G RG ) and the average general first geometric-arithmetic index G A α ( G ER ) as a function of the connection radius r of RG graphs of size n RG = 1000 .
For comparison purposes, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 are equivalent to Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5. In fact, most of the statements (i–v) made in the previous section for ER random networks are also valid for RG graphs. Even though most statements are applicable to RG graphs by just replacing G ER G RG and p f ( r ) , given the fact that this is the first study (to our knowledge) of average topological indices on RG graphs, we want to explicitly write three statements we consider relevant:
(iii’) 
The curves G A α ( G RG ) grow with r, for all α and n RG , as G A α ( G RG ) f ( r ) ; see the magenta dashed line in Figure 11a, plotted to guide the eye.
(iv’) 
For α = 0 , see the red curves in Figure 8a, Figure 9a, Figure 10a and Figure 11a, M 2 0 ( G RG ) , χ 0 ( G RG ) , I S I 0 ( G RG ) and G A 0 ( G RG ) give the average number of edges of the RG graph. That is
M 2 0 ( G RG ) = χ 0 ( G RG ) = I S I 0 ( G RG ) = G A 0 ( G RG ) = n RG ( n RG 1 ) 2 f ( r ) .
(v’) 
When n RG r 1 , we can write k u k v k RG = ( n RG 1 ) f ( r ) in Equations (1)–(5). Therefore, for n RG r 1 , the average values of the indices we are computing here are well approximated by:
M 1 α ( G RG ) = n RG [ ( n RG 1 ) f ( r ) ] α ,
M 2 α ( G RG ) = n RG 2 [ ( n RG 1 ) f ( r ) ] 1 + 2 α ,
χ α ( G RG ) = n RG 2 1 α [ ( n RG 1 ) f ( r ) ] 1 + α ,
I S I α ( G RG ) = n RG ( n R G 1 ) f ( r ) 2 1 + α ,
G A α ( G RG ) = n RG ( n RG 1 ) f ( r ) 2 1 + 2 α .
In Figure 7b, Figure 8b, Figure 9b, Figure 10b and Figure 11b, we show that Equations (21)–(25) (red-dashed lines) indeed describe very well the data (thick black curves) for k RG 10 .

3.3. Scaling Properties of General Indices on Random Geometric Graphs

In Figure 7c, Figure 8c, Figure 9c, Figure 10c and Figure 11c we show X ( G RG ) as a function of the connection radius r of RG graphs of four different sizes n RG . The panorama is equivalent to that observed in the previous section for ER random graphs: the blocks of curves, characterized by the different graph sizes, display similar curves but displaced on both axes.
Finally, in Figure 7d, Figure 8d, Figure 9d, Figure 10d and Figure 11d we plot X ( G RG ) / n RG as a function of k RG showing that all indices are now properly scaled, as anticipated in Section 3.1. Moreover, following Equations (21)–(25), for k RG 10 we have that
M 1 α ( G RG ) n RG k RG α ,
M 2 α ( G RG ) n RG k RG 1 + 2 α 2 ,
χ α ( G RG ) n RG k RG 1 + α 2 1 α ,
and
I S I α ( G RG ) n RG k RG 1 + α 2 1 + α .
While
G A α ( G RG ) n RG k RG 2 1 + 2 α
is indeed valid for all k RG .

4. Discussion and Conclusions

In this paper, within a random matrix theory (RMT) approach, we have performed a thorough statistical study of general indices on two models of random networks. In addition to well-known general indices (the first and second variable Zagreb indices, M 1 α ( G ) and M 2 α ( G ) , and the general sum-connectivity index, χ α ( G ) ) we have introduced general versions of the inverse sum indeg index I S I α ( G ) and the first geometric-arithmetic index G A α ( G ) ; see Equations (3) and (4), respectively. As models of random networks we have used Erdös–Rényi (ER) random networks G ER ( n ER , p ) and random geometric (RG) graphs G RG ( n RG , r ) . It is important to stress that we do not analize the information obtained from the different general indices we use here, instead we focus only on their computational properties. Furthermore, we left for a future investigation the analytical study of the general indices we have introduced in the present work: I S I α ( G ) and G A α ( G ) .
We have found that the average value of the indices, normalized to the network size, scale with the average degree k of the corresponding random network models. This means that k fixes the average value of any index for different combinations of network parameters; i.e., the relevant parameter of both network models is k and not the specific values of the model parameters. This result highlights the relevance of k in random network studies. Moreover, the similarity between Equations (13)–(17) and Equations (26)–(30) allow us to relate the index average values of the two random network models we study here; that is,
X ( G ER ) n ER X ( G RG ) n RG if k ER = k RG .
To verify Equation (31) in Figure 12 we show four normalized indices X ( G ) / n as a function of k for both ER random networks and RG graphs; note that we use networks of different sizes. Indeed, we observe that X ( G ER ) / n ER X ( G RG ) / n RG to a good numerical accuracy.
With Equation (31) we give a step forward in the scaling of topological indices since we have found a scaling law that covers different network models. Clearly, an open question we plan to address in a future investigation is whether Equation (31) could be extended to other random network models.
Finally, we want to recall that in complex systems research there is a continuous search of measures that could serve as complexity indicators. In particular, RMT has provided us with a number of measures able to distinguish between (i) integrable and chaotic (i.e. non-integrable) and (ii) ordered and disordered quantum systems [27,28]. Such measures are computed from the eigenvalues and eigenvectors of quantum Hamiltonian matrices. Examples of eigenvalue-based measures are the distribution of consecutive eigenvalue spacings, the spectrum rigidity and the ratios between consecutive eigenvalue spacings; while the inverse participation ratios and Shannon entropies are popular eigenvector-based measures [27,28]. It is interesting to notice that all these RMT measures have also been successfully applied to study networks and graphs since they can be computed from the eigenvalues and eigenvectors of adjacency matrices; see e.g., [29,30,31] and the references therein. Therefore, these measures are able to distinguish between networks composed by mostly isolated vertices and mostly connected networks. Furthermore, through scaling studies of RMT measures it has been possible to locate the percolation transition point of random graphs models [29,30]. It is worth mentioning that the scaling study of average topological indices performed in this paper has followed a statistical RMT approach; that is, from a detailed numerical study we have been able to identify the average degree as the universal parameter of our random network models: i.e. the parameter that fixes the average values of the topological indices.
Moreover, recently, it has been shown for RG graphs that there is a a huge correlation between the average-scaled Shannon entropy (of the adjacency matrix eigenvectors) and two average-scaled topological indices [32]: the Randić index R ( G ) and the harmonic index H ( G ) . We believe that this is a remarkable result because it validates the use of average topological indices as RMT complexity measures; already suggested in Refs. [24,25] for ER random networks. Now, it is important to stress that not every index could be used as a complexity measure. From our experience, we conclude that good candidates should fulfill a particular requirement: they should be almost constant in the trivial regimes (just as RMT measures are). For example, a useful complexity measure for random networks should be constant in the regime of mostly isolated vertices as well as above the percolation transition. Indeed, this is a property that both R ( G ) and H ( G ) have: R ( G ) H ( G ) 0 for mostly isolated vertices while R ( G ) / n H ( G ) / n 1 / 2 once the network is well above the percolation transition.
Therefore, a straightforward application of our study on general topological indices is the identification of specific indices as complexity measure candidates. Recall that we particularly require, for an average-scaled topological index to work as complexity measure, that X ( G ) / n const . for large enough k . In fact, from Equations (13)–(17) or Equations (26)–(30) we can see that the above condition is fulfilled for M 1 α ( G ) , M 2 α ( G ) , χ α ( G ) and I S I α ( G ) for specific values of α : α = 0 , α = 1 / 2 , α = 1 and α = 1 , respectively. In contrast, note that G A α ( G ) does not saturate for any α , so we discard it as possible complexity measure. Thus, in Figure 13 we plot M 1 0 ( G ) / n , M 2 1 / 2 ( G ) / n , χ 1 ( G ) / n and I S I 1 ( G ) / n as a function of the average degree k for ER random networks and RG graphs. From the behavior of the average-scaled indices reported in Figure 13 we can identify three regimes: (i) a regime of mostly isolated vertices when k < 1 / 10 , where X ( G ) / n 0 , (ii) a regime corresponding to mostly connected networks when k > 10 , where X ( G ) / n const . > 0 , and (iii) a transition regime in the interval 1 / 10 < k < 10 , which is logarithmically symmetric around the percolation transition point k 1 . Accordingly, we propose the use of M 1 0 ( G ) / n , M 2 1 / 2 ( G ) / n , χ 1 ( G ) / n and I S I 1 ( G ) / n as complexity measures for random network models. (We should remark that the behavior reported in Figure 13 for M 2 1 / 2 ( G ) / n and χ 1 ( G ) / n was already expected since they coincide, respectively, with R ( G ) / n and H ( G ) / 2 n ; already studied in Refs. [24,25,32] and used there as measures equivalent to RMT complexity measures.)
We hope that our work may motivate further analytical as well as numerical studies of general topological indices on random networks.

Author Contributions

The authors contributed equally to this work. R.A.-S., I.F.H.-G., J.A.M.-B. and J.M.S. conceived, designed and performed the numerical experiments; R.A.-S., I.F.H.-G., J.A.M.-B. and J.M.S. analyzed the data and wrote the paper. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

J.A.M.-B. acknowledges financial support from FAPESP (Grant No. 2019/06931-2), Brazil, CONACyT (Grant No. 2019-000009-01EXTV-00067) and PRODEP-SEP (Grant No. 511-6/2019.-11821), Mexico. J.M.S. was supported by a grant from Agencia Estatal de Investigación (PID2019-106433GB-I00/AEI/ 10.13039/501100011033), Spain.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Wiener, H. Structural determination of paraffin boiling points. J. Am. Chem. Soc. 1947, 69, 17–20. [Google Scholar] [CrossRef] [PubMed]
  2. Hosoya, H. Topological index. A newly proposed quantity characterizing the topological nature of structural isomers of saturated hydrocarbons. Bull. Chem. Soc. Jpn. 1971, 44, 2331–2339. [Google Scholar] [CrossRef] [Green Version]
  3. Gutman, I.; Trinajstić, N. Graph theory and molecular orbitals. Total ϕ-electron energy of alternant hydrocarbons. Chem. Phys. Lett. 1972, 17, 535–538. [Google Scholar] [CrossRef]
  4. Nikolic, S.; Kovacevic, G.; Milicevic, A.; Trinajstic, N. The Zagreb indices 30 years after. Croat. Chem. Acta 2003, 76, 113–124. [Google Scholar]
  5. Randić, M. On characterization of molecular branching. J. Am. Chem. Soc. 1975, 97, 6609–6615. [Google Scholar] [CrossRef]
  6. Randić, M. Novel graph theoretical approach to heteroatoms in quantitative structure-activity relationships. Chemom. Int. Lab. Syst. 1991, 10, 213–227. [Google Scholar] [CrossRef]
  7. Randić, M. On computation of optimal parameters for multivariate analysis of structure-property relationship. J. Comput. Chem. 1991, 12, 970–980. [Google Scholar] [CrossRef]
  8. Miličević, A.; Nikolić, S. On variable Zagreb indices. Croat. Chem. Acta 2004, 77, 97–101. [Google Scholar]
  9. Li, X.; Zheng, J. A unified approach to the extremal trees for different indices. MATCH Commun. Math. Comput. Chem. 2005, 54, 195–208. [Google Scholar]
  10. Li, X.; Zhao, H. Trees with the first smallest and largest generalized topological indices. MATCH Commun. Math. Comput. Chem. 2004, 50, 57–62. [Google Scholar]
  11. Nikolić, S.; Miličević, A.; Trinajstić, N.; Jurić, A. On use of the variable Zagreb Index in QSPR: Boiling points of Benzenoid hydrocarbons. Molecules 2004, 9, 1208–1221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Zhou, B.; Trinajstić, N. On general sum-connectivity index. J. Math. Chem. 2010, 47, 210–218. [Google Scholar] [CrossRef]
  13. Rodríguez, J.M.; Sigarreta, J.M. New results on the Harmonic index and its generalizations. MATCH Commun. Math. Comput. Chem. 2017, 78, 387–404. [Google Scholar]
  14. Hafeez, S.; Farooq, R. On generalized inverse sum indeg index and energy of graphs. AIMS Math. 2020, 5, 2388–2411. [Google Scholar] [CrossRef]
  15. Das, K.C.; Gutman, I.; Furtula, B. Survey on geometric–arithmetic indices of graphs. MATCH Commun. Math. Comput. Chem. 2011, 65, 595–644. [Google Scholar]
  16. Wilczek, P. New geometric–arithmetic indices. MATCH Commun. Math. Comput. Chem. 2018, 79, 5–54. [Google Scholar]
  17. Aouchiche, M.; Ganesan, V. Adjusting geometric–arithmetic index to estimate boiling point. MATCH Commun. Math. Comput. Chem. 2020, 84, 483–497. [Google Scholar]
  18. Solomonoff, R.; Rapoport, A. Connectivity of random nets. Bull. Math. Biophys. 1951, 13, 107–117. [Google Scholar] [CrossRef]
  19. Erdös, P.; Rényi, A. On random graphs. Publ. Math. 1959, 6, 290–297. [Google Scholar]
  20. Erdös, P.; Rényi, A. On the evolution of random graphs. Inst. Hung. Acad. Sci. 1960, 5, 17–61. [Google Scholar]
  21. Erdös, P.; Rényi, A. On the strength of connectedness of a random graph. Acta Math. Hung. 1961, 12, 261–267. [Google Scholar] [CrossRef]
  22. Dall, J.; Christensen, M. Random geometric graphs. Phys. Rev. E 2002, 66, 016121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Penrose, M. Random Geometric Graphs; Oxford University Press: Oxford, UK, 2003. [Google Scholar]
  24. Martínez-Martínez, C.T.; Méndez-Bermúdez, J.A.; Rodríguez, J.M.; Sigarreta, J.M. Computational and analytical studies of the Randić index in Erdös–Rényi models. Appl. Math. Comput. 2020, 377, 125137. [Google Scholar]
  25. Martínez-Martínez, C.T.; Méndez-Bermúdez, J.A.; Rodríguez, J.M.; Sigarreta, J.M. Computational and analytical studies of the harmonic index in Erdös–Rényi models. MATCH Commun. Math. Comput. Chem. 2021, 85. in press. [Google Scholar]
  26. Estrada, E.; Sheerin, M. Random rectangular graphs. Phys. Rev. E 2015, 91, 042805. [Google Scholar] [CrossRef] [Green Version]
  27. Metha, M.L. Random Matrices; Elsevier: Amsterdam, The Netherlands, 2004. [Google Scholar]
  28. Haake, F. Quantum Signatures of Chaos; Springer: Berlin, Germany, 2010. [Google Scholar]
  29. Mendez-Bermudez, J.A.; Alcazar-Lopez, A.; Martinez-Mendoza, A.J.; Rodrigues, F.A.; Peron, T.K.D.M. Universality in the spectral and eigenfunction properties of random networks. Phys. Rev. E 2015, 91, 032122. [Google Scholar] [CrossRef] [Green Version]
  30. Alonso, L.; Mendez-Bermudez, J.A.; Gonzalez-Melendrez, A.; Moreno, Y. Weighted random-geometric and random-rectangular graphs: Spectral and eigenfunction properties of the adjacency matrix. J. Complex Netw. 2018, 6, 753. [Google Scholar] [CrossRef]
  31. Torres-Vargas, G.; Fossion, R.; Mendez-Bermudez, J.A. Normal mode analysis of spectra of random networks. Physica A 2020, 545, 123298. [Google Scholar] [CrossRef] [Green Version]
  32. Aguilar-Sanchez, R.; Mendez-Bermudez, J.A.; Rodrigues, F.A.; Sigarreta-Almira, J.M. Topological versus spectral properties of random geometric graphs. arXiv 2020, arXiv:2007.02453. [Google Scholar]
Figure 1. (a,b) Average first variable Zagreb index M 1 α ( G ER ) as a function of the probability p of ER random networks of size n ER = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (8). The blue dashed line in (b) marks k ER = 10 . (c) M 1 α ( G ER ) as a function of the probability p of ER random networks of four different sizes n ER . (d) M 1 α ( G ER ) / n ER as a function of the average degree k ER . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Figure 1. (a,b) Average first variable Zagreb index M 1 α ( G ER ) as a function of the probability p of ER random networks of size n ER = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (8). The blue dashed line in (b) marks k ER = 10 . (c) M 1 α ( G ER ) as a function of the probability p of ER random networks of four different sizes n ER . (d) M 1 α ( G ER ) / n ER as a function of the average degree k ER . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Symmetry 12 01341 g001
Figure 2. (a,b) Average second variable Zagreb index M 2 α ( G ER ) as a function of the probability p of ER random networks of size n ER = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (9). The blue dashed line in (b) marks k ER = 10 . (c) M 2 α ( G ER ) as a function of the probability p of ER random networks of four different sizes n ER . (d) M 2 α ( G ER ) / n ER as a function of the average degree k ER . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Figure 2. (a,b) Average second variable Zagreb index M 2 α ( G ER ) as a function of the probability p of ER random networks of size n ER = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (9). The blue dashed line in (b) marks k ER = 10 . (c) M 2 α ( G ER ) as a function of the probability p of ER random networks of four different sizes n ER . (d) M 2 α ( G ER ) / n ER as a function of the average degree k ER . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Symmetry 12 01341 g002
Figure 3. (a,b) Average general sum-connectivity index χ α ( G ER ) as a function of the probability p of ER random networks of size n ER = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (10). The blue dashed line in (b) marks k ER = 10 . (c) χ α ( G ER ) as a function of the probability p of ER random networks of four different sizes n ER . (d) χ α ( G ER ) / n ER as a function of the average degree k ER . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Figure 3. (a,b) Average general sum-connectivity index χ α ( G ER ) as a function of the probability p of ER random networks of size n ER = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (10). The blue dashed line in (b) marks k ER = 10 . (c) χ α ( G ER ) as a function of the probability p of ER random networks of four different sizes n ER . (d) χ α ( G ER ) / n ER as a function of the average degree k ER . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Symmetry 12 01341 g003
Figure 4. (a,b) Average general inverse sum indeg index I S I α ( G ER ) as a function of the probability p of ER random networks of size n ER = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (11). The blue dashed line in (b) marks k ER = 10 . (c) I S I α ( G ER ) as a function of the probability p of ER random networks of four different sizes n ER . (d) I S I α ( G ER ) / n ER as a function of the average degree k ER . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Figure 4. (a,b) Average general inverse sum indeg index I S I α ( G ER ) as a function of the probability p of ER random networks of size n ER = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (11). The blue dashed line in (b) marks k ER = 10 . (c) I S I α ( G ER ) as a function of the probability p of ER random networks of four different sizes n ER . (d) I S I α ( G ER ) / n ER as a function of the average degree k ER . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Symmetry 12 01341 g004
Figure 5. (a,b) Average general first geometric-arithmetic index G A α ( G ER ) as a function of the probability p of ER random networks of size n ER = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.4 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (12). The blue dashed line in (b) marks k ER = 10 . The magenta dashed line in (a) proportional to p is shown to guide the eye. (c) G A α ( G ER ) as a function of the probability p of ER random networks of four different sizes n ER . (d) G A α ( G ER ) / n ER as a function of the average degree k ER . Same curves as in panel (c).
Figure 5. (a,b) Average general first geometric-arithmetic index G A α ( G ER ) as a function of the probability p of ER random networks of size n ER = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.4 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (12). The blue dashed line in (b) marks k ER = 10 . The magenta dashed line in (a) proportional to p is shown to guide the eye. (c) G A α ( G ER ) as a function of the probability p of ER random networks of four different sizes n ER . (d) G A α ( G ER ) / n ER as a function of the average degree k ER . Same curves as in panel (c).
Symmetry 12 01341 g005
Figure 6. (a) Average index M 1 1 ( G RG ) / n RG as a function of the connection radius r of RG graphs of four different sizes n RG (symbols). Full lines correspond to Equation (19). (b) M 1 1 ( G RG ) / n RG as a function of the average degree k RG . Same data as in panel (a).
Figure 6. (a) Average index M 1 1 ( G RG ) / n RG as a function of the connection radius r of RG graphs of four different sizes n RG (symbols). Full lines correspond to Equation (19). (b) M 1 1 ( G RG ) / n RG as a function of the average degree k RG . Same data as in panel (a).
Symmetry 12 01341 g006
Figure 7. (a,b) Average first variable Zagreb index M 1 α ( G RG ) as a function of the connection radius r of RG graphs of size n RG = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (21). The blue dashed line in (b) marks k RG 10 . (c) M 1 α ( G RG ) as a function of the connection radius r of RG graphs of four different sizes n RG . (d) M 1 α ( G RG ) / n RG as a function of the average degree k RG . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Figure 7. (a,b) Average first variable Zagreb index M 1 α ( G RG ) as a function of the connection radius r of RG graphs of size n RG = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (21). The blue dashed line in (b) marks k RG 10 . (c) M 1 α ( G RG ) as a function of the connection radius r of RG graphs of four different sizes n RG . (d) M 1 α ( G RG ) / n RG as a function of the average degree k RG . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Symmetry 12 01341 g007
Figure 8. (a,b) Average second variable Zagreb index M 2 α ( G RG ) as a function of the connection radius r of RG graphs of size n RG = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (22). The blue dashed line in (b) marks k RG 10 . (c) M 2 α ( G RG ) as a function of the connection radius r of RG graphs of four different sizes n RG . (d) M 2 α ( G RG ) / n RG as a function of the average degree k RG . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Figure 8. (a,b) Average second variable Zagreb index M 2 α ( G RG ) as a function of the connection radius r of RG graphs of size n RG = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (22). The blue dashed line in (b) marks k RG 10 . (c) M 2 α ( G RG ) as a function of the connection radius r of RG graphs of four different sizes n RG . (d) M 2 α ( G RG ) / n RG as a function of the average degree k RG . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Symmetry 12 01341 g008
Figure 9. (a,b) Average general sum-connectivity index χ α ( G RG ) as a function of the connection radius r of RG graphs of size n RG = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (23). The blue dashed line in (b) marks k RG 10 . (c) χ α ( G RG ) as a function of the connection radius r of RG graphs of four different sizes n RG . (d) χ α ( G RG ) / n RG as a function of the average degree k RG . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Figure 9. (a,b) Average general sum-connectivity index χ α ( G RG ) as a function of the connection radius r of RG graphs of size n RG = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (23). The blue dashed line in (b) marks k RG 10 . (c) χ α ( G RG ) as a function of the connection radius r of RG graphs of four different sizes n RG . (d) χ α ( G RG ) / n RG as a function of the average degree k RG . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Symmetry 12 01341 g009
Figure 10. (a,b) Average general inverse sum indeg index I S I α ( G RG ) as a function of the connection radius r of RG graphs of size n RG = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (24). The blue dashed line in (b) marks k RG 10 . (c) I S I α ( G RG ) as a function of the connection radius r of RG graphs of four different sizes n RG . (d) I S I α ( G RG ) / n RG as a function of the average degree k RG . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Figure 10. (a,b) Average general inverse sum indeg index I S I α ( G RG ) as a function of the connection radius r of RG graphs of size n RG = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.2 (from bottom to top). The red (blue) curve in (a) corresponds to α = 0 [ α = 1 ]. The red dashed lines in (b) are Equation (24). The blue dashed line in (b) marks k RG 10 . (c) I S I α ( G RG ) as a function of the connection radius r of RG graphs of four different sizes n RG . (d) I S I α ( G RG ) / n RG as a function of the average degree k RG . Same curves as in panel (c). The inset in (d) is the enlargement of the cyan rectangle.
Symmetry 12 01341 g010
Figure 11. (a,b) Average general first geometric-arithmetic index G A α ( G RG ) as a function of the connection radius r of RG graphs of size n RG = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.4 (from bottom to top). The red curve in (a) corresponds to α = 0 . The red dashed lines in (b) are Equation (25). The blue dashed line in (b) marks k RG 10 . The magenta dashed line in (a) proportional to f ( r ) is shown to guide the eye. (c) G A α ( G RG ) as a function of the connection radius r of RG graphs of four different sizes n RG . (d) G A α ( G RG ) / n RG as a function of the average degree k RG . Same curves as in panel (c).
Figure 11. (a,b) Average general first geometric-arithmetic index G A α ( G RG ) as a function of the connection radius r of RG graphs of size n RG = 1000 . Here we show curves for α [ 2 , 2 ] in steps of 0.4 (from bottom to top). The red curve in (a) corresponds to α = 0 . The red dashed lines in (b) are Equation (25). The blue dashed line in (b) marks k RG 10 . The magenta dashed line in (a) proportional to f ( r ) is shown to guide the eye. (c) G A α ( G RG ) as a function of the connection radius r of RG graphs of four different sizes n RG . (d) G A α ( G RG ) / n RG as a function of the average degree k RG . Same curves as in panel (c).
Symmetry 12 01341 g011
Figure 12. (a) M 1 α ( G ) / n , (b) M 2 α ( G ) / n , (c) χ α ( G ) / n and (d) I S I α ( G ) / n as a function of the average degree k for ER random networks (red lines) and RG graphs (black lines). We choose n = 1000 for ER random networks and n = 250 for RG graphs.
Figure 12. (a) M 1 α ( G ) / n , (b) M 2 α ( G ) / n , (c) χ α ( G ) / n and (d) I S I α ( G ) / n as a function of the average degree k for ER random networks (red lines) and RG graphs (black lines). We choose n = 1000 for ER random networks and n = 250 for RG graphs.
Symmetry 12 01341 g012
Figure 13. (a) M 1 0 ( G ) / n , (b) M 2 1 / 2 ( G ) / n , (c) χ 1 ( G ) / n and (d) I S I 1 ( G ) / n as a function of the average degree k for ER random networks (red lines) and RG graphs (black lines). We choose n = 1000 for both ER random networks and RG graphs. Magenta dashed lines indicate k = 1 / 10 and k = 10 .
Figure 13. (a) M 1 0 ( G ) / n , (b) M 2 1 / 2 ( G ) / n , (c) χ 1 ( G ) / n and (d) I S I 1 ( G ) / n as a function of the average degree k for ER random networks (red lines) and RG graphs (black lines). We choose n = 1000 for both ER random networks and RG graphs. Magenta dashed lines indicate k = 1 / 10 and k = 10 .
Symmetry 12 01341 g013

Share and Cite

MDPI and ACS Style

Aguilar-Sánchez, R.; Herrera-González, I.F.; Méndez-Bermúdez, J.A.; Sigarreta, J.M. Computational Properties of General Indices on Random Networks. Symmetry 2020, 12, 1341. https://doi.org/10.3390/sym12081341

AMA Style

Aguilar-Sánchez R, Herrera-González IF, Méndez-Bermúdez JA, Sigarreta JM. Computational Properties of General Indices on Random Networks. Symmetry. 2020; 12(8):1341. https://doi.org/10.3390/sym12081341

Chicago/Turabian Style

Aguilar-Sánchez, R., I. F. Herrera-González, J. A. Méndez-Bermúdez, and José M. Sigarreta. 2020. "Computational Properties of General Indices on Random Networks" Symmetry 12, no. 8: 1341. https://doi.org/10.3390/sym12081341

APA Style

Aguilar-Sánchez, R., Herrera-González, I. F., Méndez-Bermúdez, J. A., & Sigarreta, J. M. (2020). Computational Properties of General Indices on Random Networks. Symmetry, 12(8), 1341. https://doi.org/10.3390/sym12081341

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop