Next Article in Journal
On the Extremal Weighted Mostar Index of Bicyclic Graphs
Previous Article in Journal
Transfer Learning for Logistic Regression with Differential Privacy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Limiting Behaviors of the Gutman and Schultz Indices in Random 2k-Sided Chains

School of Mathematics and Big Data, Anhui University of Science & Technology, Huainan 232001, China
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(8), 518; https://doi.org/10.3390/axioms13080518
Submission received: 26 May 2024 / Revised: 5 July 2024 / Accepted: 17 July 2024 / Published: 30 July 2024

Abstract

:
The study of complex networks with topological indices has flourished in recent years. The aim of this paper is to study the limiting behaviors of Gutman and Schultz indices in random polygonal chains, whose graph-theoretic mathematical properties and their future applications have attracted the interest of scientists. By applying the concepts of symmetry and asymptotics as well as the knowledge of probability theory, we obtain explicit analytic expressions for the Gutman and Schultz indices of n random 2k-vertex chains and prove that they converge to a normal distribution, which contributes to a deeper understanding of the structural features of random polygonal chains and plays a crucial role in the study of the limiting behavior of topological indices and their applications.

1. Introduction

In this paper, we only discuss finite, simple, and symmetrical connected graphs, which are mentioned in the work of Bondy and Murty [1] and the references they cite. Graph theory is closely intertwined with chemistry, and its applications within the field are extensive. The description of chemical compounds is based on chemical graph theory, with vertices representing atoms and edges representing the covalent bonds between atoms, whose fundamental technique involves creating a model of a compound’s molecular structure [2,3].
The molecular formula of a compound can give an indication of different molecular structures and properties, but theoretical chemists are concerned with the physical and chemical properties of compounds and their relationship to their molecular formula [4,5].
We will use the following basic notations, as defined in [1]. Consider a graph G represented by the tuple G = ( V G , E G ) , where V G is the set of vertices and E G is the set of edges. Then, let d G ( u , v ) (simply d ( u , v ) ) denote the distance between two vertices u and v in G, and the distance is defined as the length of the shortest path between u and v in G. The Wiener index (also known as the Wiener transmission) of a graph G with edge set E G is defined as the sum of distances between all possible pairs of vertices in G. The index was introduced by H. Wiener in 1947 [6] and can be expressed mathematically as follows:
W ( G ) = 1 2 { u , v } V G d G ( u , v ) .
The Wiener index is a well-researched, comprehensively understood, and extensively employed molecular shape descriptor rooted in graph theory [7,8]. Furthermore, its usage and study have been progressively growing, as evidenced by [9,10].
A weighted graph ( G , ω ) is composed of a graph G = ( V G , E G ) and a weight function ω : V G N + . Suppose ⊕ represents one of the four arithmetic operations +, −, ×, and ÷. Consequently, the weighted Wiener index W ( G , ω ) can be defined as follows:
W ( G , ω ) = 1 2 u V G v V G ( ω ( u ) ω ( v ) ) d G ( u , v ) .
It is evident that when ω 1 and ⊕ represent the operation ×, W ( G , ω ) = W ( G ) .
Suppose ⊕ denotes the operation × and ω ( · ) d G ( · ) . In such a scenario, Equation (2) is identical to
G u t ( G ) = 1 2 u V G v V G ( d G ( u ) d G ( v ) ) d G ( u , v ) = { u . v } V G ( d G ( u ) d G ( v ) ) d G ( u , v ) ,
which is the Gutman index [11]. For more information on the prospective chemical uses of the Gutman index, alongside related metrics, as well as their theoretical analyses, wherein polycyclic molecules present more complex cases, please refer to [12].
Suppose ⊕ denotes the operation + and ω ( · ) d G ( · ) . In such a scenario, Equation (2) is identical to
S ( G ) = 1 2 u V G v V G ( d G ( u ) + d G ( v ) ) d G ( u , v ) = { u . v } V G ( d G ( u ) + d G ( v ) ) d G ( u , v ) ,
which is the Schultz index. For further developments in topology indexes for [13,14,15,16], encompassing mathematical properties, discrimination, and applications, please refer to [17] for additional articles.
In this paper, we mainly talk about the polygonal chain G n with n regular 2k-gon graphs, which are constructed in the following way. Firstly, G 1 is a 2k-sided polygon, and G 2 is the graph with two 2k-sided polygons connected by an edge, as shown in Figure 1. Secondly, to construct the random polygon chain G n + 1 , a new terminal polygon with 2 k sides, denoted as H n + 1 , is attached to the existing polygon chain G n which contains n 2k-gon polygons. This process is illustrated in Figure 2. As demonstrated in Figure 3, we can attach the terminal polygon with 2 k sides H n + 1 to G n with k methods and denote the resulting figures with G n + 1 1 , G n + 1 2 , G n + 1 3 , ⋯, G n + 1 k , respectively. At each step, one of the following possible constructions is randomly selected:
  • G n G n + 1 1 with probability p 1 ;
  • G n G n + 1 2 with probability p 2 ;
  • G n G n + 1 3 with probability p 3 ;
    ⋮   ⋮        ⋮
  • G n G n + 1 k 1 with probability p k 1 ;
  • G n G n + 1 k with probability p k = 1 p 1 p 2 p 3 p k 1 .
We define k random variables Z n 1 , Z n 2 , Z n 3 , ⋯, Z n k to represent our options. If we choose G n + 1 i , we set Z n i = 1 ; otherwise, Z n i = 0 for i = 1 , 2 , 3 , , k . It is easy to see that
P ( Z n i = 1 ) = p i , P ( Z n i = 0 ) = 1 p i , i = 1 , 2 , 3 , , k
and Z n 1 + Z n 2 + Z n 3 + + Z n k = 1 .
Through the aforementioned procedures, we generate a random polygonal chain G n ( p 1 , p 2 , p 3 , , p k ) whose total number of sides is 2 k . That is to say, each cycle has 2 k sides and it is called a chain of 2k cycles. We consistently abbreviate G n ( p 1 , p 2 , p 3 , , p k ) as G n .
Noting that G n is a random graph, both G u t ( G n ) and S ( G n ) are random variables in probability. In the view of probability, one problem appears naturally: when n is big enough, the distribution of G u t ( G n ) and S ( G n ) will either look or not look like Gaussian distribution.
In this paper, we quantify the distributions of values of G u t ( G n ) and S ( G n ) .
Polyphenylenes, cyclooctatetraenes, and their derivatives have been the subject of considerable attention among chemists for many years, closely associated with certain classes of graphs. These graphs exhibit intriguing structures and properties, both of which have corresponding chains of chemical compounds, as detailed in [18,19,20]. We speculate whether the 2k cycle chains exhibit properties similar to those of polyphenylenes, cyclooctatetraenes, and their derivatives. The research methodologies in chemical graph theory assist in analyzing and describing the topological features, connectivity, and possible reaction pathways of these chains. By considering various indices and parameters, we can gain better insights into harnessing the role of 2k-gon chains in chemistry and materials science, paving the way for future innovations and advancements [21,22].
Researchers have investigated the mathematical properties of random polyphenylene chains and random cyclooctatetraene chains [23,24,25,26]. Zhu [27] conducted research on establishing explicit analytical expressions for the simple formulas of the expected values of Gutman and Schultz indices in a random polygon. By applying the concepts of symmetry and probability theory, we aim to derive the exact analytical expressions, as well as their asymptotic behaviors, for the variances of the Schultz index and Gutman index of a random chain composed of n 2k-gon polygons. This research plays a crucial role in the study and application of topological indices [28].
In this paper, we make the following hypothesis.
Hypothesis 1.
Our choice of attaching the new terminal 2k-gon polygon H n + 1 to G n , n = 2 , 3 , is random and independent. In more accurate words, the series of random variables Z n 1 , Z n 2 , Z n 3 , , Z n k n = 2 is independent and has the same law (1.5). For i 1 , 2 , 3 , , k , we have 0 < p i < 1 .
(a)
Under Hypothesis 1, we give analytical expressions for the variances of G u t ( G n ) and S ( G n ) .
(b)
We prove that the random variables G u t ( G n ) and S ( G n ) are asymptotic to normal distributions when n . That is
lim n sup x R P ( ( G u t ( G n ) E ( G u t ( G n ) ) V a r ( G u t ( G n ) ) ) x ) x 1 2 π e t 2 2 d t = 0
where E ( X ) and V a r ( X ) denote the expected value and variance of the random variable X, respectively.

2. The Limiting Behavior of the Gutman Index of a Random Polygon Chain

For a random chain G n , which is composed of n 2k-gon polygons, both the Gutman index and the Schultz index are random variables. In this section, we will explore and compute the variances of G u t ( G n ) and S ( G n ) . As depicted in Figure 1, G n + 1 is obtained by connecting G n to a new terminal hexagon H n + 1 by the edge u n x 1 , with the vertices of H n + 1 labeled x 1 , x 2 , ⋯, x 2 k in clockwise orientation, where n 2 . On the one hand, for all v V G n V G n + 1 , x i H n + 1 , it holds that
d ( x 1 , v ) = d ( u n , v ) + 1 , d ( x 2 , v ) = d ( u n , v ) + 2 , , d ( x k , v ) = d ( u n , v ) + k ,
d ( x k + 1 , v ) = d ( u n , v ) + k + 1 , d ( x k + 2 , v ) = d ( u n , v ) + k , , d ( x 2 k , v ) = d ( u n , v ) + 2 ,
v V G n d G n + 1 ( v ) = ( 4 k + 2 ) n 1 .
On the other hand,
i = 1 2 k d ( x i ) d ( x 1 , x i ) = 2 k 2 , i = 1 2 k d ( x i ) d ( x 2 , x i ) = 2 k 2 + 1 , , i = 1 2 k d ( x i ) d ( x k + 1 , x i ) = 2 k 2 + k ,
i = 1 2 k d ( x i ) d ( x k + 2 , x i ) = 2 k 2 + k + 1 , i = 1 2 k d ( x i ) d ( x k + 3 , x i ) = 2 k 2 + k , , i = 1 2 k d ( x i ) d ( x 2 k , x i ) = 2 k 2 + 1 .
In [27] Theorem 1, the author demonstrates that
E ( G u t ( G n ) ) = { ( 18 k 3 + 16 k 2 + 10 k + 2 ) ( 2 k + 1 ) i = 1 k 1 [ 4 k 2 ( 4 i 2 ) k 2 i ] p i } n 3 3 { ( 4 k 2 + 6 k + 2 ) ( 2 k + 1 ) i = 1 k 1 [ 4 k 2 ( 4 i 2 ) k 2 i ] p i } n 2 + { ( 4 k 3 4 k 2 + 8 k + 7 ) 2 ( 2 k + 1 ) i = 1 k 1 [ 4 k 2 ( 4 i 2 ) k 2 i ] p i } n 3 1 .
This section begins with the introduction of the first main result.
Theorem 1.
The following results obtained are dependent on whether the following Hypothesis 1 is true.
  • (i) The Gutman index of G n , which is a random 2k-sided chain, is G u t ( G n ) , and the variance in it is determined by
V a r ( G u t ( G n ) ) = 1 30 ( σ 2 n 5 5 r n 4 + 10 σ ˜ 2 n 3 + ( 65 r 30 σ 2 45 σ ˜ 2 ) n 2 + ( 120 r + 59 σ 2 + 65 σ ˜ 2 ) n + ( 60 r 30 σ 2 30 σ ˜ 2 ) ) .
where
σ 2 = i = 1 k ( a i 2 p i ) ( i = 1 k ( a i p i ) ) 2 σ ˜ 2 = i = 1 k ( b i 2 p i ) ( i = 1 k ( b i p i ) ) 2 r = i = 1 k ( a i b i p i ) i = 1 k ( a i p i ) i = 1 k ( b i p i )
with
a i = ( i + 1 ) ( 4 k + 2 ) 2 b i = ( 4 k + 2 ) [ ( i + 1 ) ( 4 k + 3 ) ( 2 k 2 + i ) ] i = 1 , 2 , 3 , , k
(ii) As n , the limiting behavior of G u t ( G n ) converges to a normal distribution. In other words,
lim n sup x R P ( ( G u t ( G n ) E ( G u t ( G n ) ) V a r ( G u t ( G n ) ) ) x ) x 1 2 π e t 2 2 d u = 0
Proof. 
(i) As concluded in [27], we can derive that
G u t ( G n + 1 ) = G u t ( G n ) + ( 4 k + 2 ) v V G n d ( v ) d ( u n , v ) + ( 8 k 3 + 20 k 2 + 12 k + 2 ) n + ( 4 k 3 4 k 1 ) .
Let
A n : = ( 4 k + 2 ) v V G n d ( v ) d ( u n , v ) .
and we can then obtain
G u t ( G n + 1 ) = G u t ( G n ) + A n + ( 8 k 3 + 20 k 2 + 12 k + 2 ) n + ( 4 k 3 4 k 1 ) .
 □
With regard to the random variables Z n 1 , Z n 2 , Z n 3 , ⋯, Z n k mentioned in the method of constructing G n + 1 from G n in Section 1, we can establish the following four equalities.
  • Equality 1.
A n Z n 1 = { A n 1 + 2 ( 4 k + 2 ) 2 n ( 4 k + 2 ) [ 2 ( 4 k + 3 ) ( 2 k 2 + 1 ) ] } Z n 1
The above equality is trivial if Z n 1 = 0 , so our focus lies on the case in which Z n 1 = 1 , indicating the transition G n G n + 1 1 . In this scenario, vertex u n in G n coincides with either vertex x 2 or x 2 k in H n . That is to say, u n = x 1 would be replaced by u n = x 2 or u n = x 2 k , as shown in Figure 4. As a result, A n can be expressed as
( 4 k + 2 ) v V G n d ( v ) d ( x 2 , v ) = ( 4 k + 2 ) v V G n 1 d ( v ) d ( x 2 , v ) + ( 4 k + 2 ) v V H n d ( v ) d ( x 2 , v ) = ( 4 k + 2 ) v V G n 1 d ( v ) ( d ( v , u n 1 ) + d ( x 2 , u n 1 ) ) + ( 4 k + 2 ) v V H n d ( v ) d ( x 2 , v ) = ( 4 k + 2 ) v V G n 1 d ( v ) ( d ( v , u n 1 ) + 2 ) + ( 4 k + 2 ) ( 2 k 2 + 1 ) = A n 1 + 2 ( 4 k + 2 ) v V G n 1 d ( v ) + ( 4 k + 2 ) ( 2 k 2 + 1 ) = A n 1 + 2 ( 4 k + 2 ) ( ( 4 k + 2 ) n ( 4 k + 3 ) ) + ( 4 k + 2 ) ( 2 k 2 + 1 ) , = A n 1 + 2 ( 4 k + 2 ) 2 n ( 4 k + 2 ) [ 2 ( 4 k + 3 ) ( 2 k 2 + 1 ) ]
By employing Equations (6)–(8) mentioned above, we are able to establish the desired equality.
  • Equality 2.
A n Z n 2 = { A n 1 + 3 ( 4 k + 2 ) 2 n ( 4 k + 2 ) [ 3 ( 4 k + 3 ) ( 2 k 2 + 2 ) ] } Z n 2
Similar to the proof of Equality 1, we solely focus on the scenario where Z n 2 = 1 , resulting in the transformation G n G n + 1 2 . The proof follows a similar approach, and we have chosen to omit the specific details.
Similarly, when reaching the vertex labeled as x k + 1 in H n , we just only consider the case Z n k = 1 , which corresponds to G n G n + 1 k . By doing so, we obtain the following result:
  • Equality 3.
A n Z n k = { A n 1 + ( k + 1 ) ( 4 k + 2 ) 2 n ( 4 k + 2 ) [ ( k + 1 ) ( 4 k + 3 ) ( 2 k 2 + k ) ] } Z n k
The proof follows the same steps as the previous equalities, and we will omit the details. Therefore, we can derive the general form of the analytic formula for A n Z n i as described below:
  • Equality 4.
A n Z n i = { A n 1 + ( i + 1 ) ( 4 k + 2 ) 2 n ( 4 k + 2 ) [ ( i + 1 ) ( 4 k + 3 ) ( 2 k 2 + i ) ] } Z n i
where i 1 , 2 , , k . We will restrict our analysis to the case where Z n i = 1 , corresponding to G n G n + 1 i . The proof follows a similar approach as for the previously mentioned equalities and we will skip the details.
Taking into account the observation that Z n 1 + Z n 2 + + Z n i + + Z n k = 1 , as mentioned earlier, we can conclude that
A n = A n ( Z n 1 + Z n 2 + + Z n i + + Z n k ) = { A n 1 + 2 ( 4 k + 2 ) 2 n ( 4 k + 2 ) [ 2 ( 4 k + 3 ) ( 2 k 2 + 1 ) ] } Z n 1 + { A n 1 + 3 ( 4 k + 2 ) 2 n ( 4 k + 2 ) [ 3 ( 4 k + 3 ) ( 2 k 2 + 2 ) ] } Z n 2 + + { A n 1 + ( i + 1 ) ( 4 k + 2 ) 2 n ( 4 k + 2 ) [ ( i + 1 ) ( 4 k + 3 ) ( 2 k 2 + i ) ] } Z n i + + { A n 1 + ( k + 1 ) ( 4 k + 2 ) 2 n ( 4 k + 2 ) [ ( k + 1 ) ( 4 k + 3 ) ( 2 k 2 + k ) ] } Z n k = A n 1 ( Z n 1 + Z n 2 + + Z n i + + Z n k ) + [ 2 ( 4 k + 2 ) 2 Z n 1 + 3 ( 4 k + 2 ) 2 Z n 2 + + ( i + 1 ) ( 4 k + 2 ) 2 Z n i + + ( k + 1 ) ( 4 k + 2 ) 2 Z n k ] n { ( 4 k + 2 ) [ 2 ( 4 k + 3 ) ( 2 k 2 + 1 ) ] Z n 1 + ( 4 k + 2 ) [ 3 ( 4 k + 3 ) ( 2 k 2 + 2 ) ] Z n 2 + + ( 4 k + 2 ) [ ( i + 1 ) ( 4 k + 3 ) ( 2 k 2 + i ) ] Z n i + + ( 4 k + 2 ) [ ( k + 1 ) ( 4 k + 3 ) ( 2 k 2 + k ) ] Z n k } = A n 1 + i = 1 k ( i + 1 ) ( 4 k + 2 ) 2 Z n i i = 1 k ( 4 k + 2 ) [ ( i + 1 ) ( 4 k + 3 ) ( 2 k 2 + i ) ] } Z n i = A n 1 + i = 1 k a i Z n i i = 1 k b i Z n i = A n 1 + n U n V n
where for each n, i { 1 , 2 , 3 , , k } ,
U n = i = 1 k a i Z n i , V n = i = 1 k b i Z n i .
From Equation (12), we can deduce
G u t ( G n ) = G u t ( G 1 ) + i = 1 n 1 A i + i = 1 n 1 [ ( 8 k 3 + 20 k 2 + 12 k + 2 ) n + ( 4 k 3 4 k 1 ) ] = G u t ( G 1 ) + i = 1 n 1 ( j = 1 i 1 ( A j + 1 A j ) + A 1 ) + i = 1 n 1 [ ( 8 k 3 + 20 k 2 + 12 k + 2 ) n + ( 4 k 3 4 k 1 ) ] = G u t ( G 1 ) + i = 1 n 1 j = 1 i 1 ( A j + 1 A j ) + ( n 1 ) A 1 + i = 1 n 1 [ ( 8 k 3 + 20 k 2 + 12 k + 2 ) n + ( 4 k 3 4 k 1 ) ] = G u t ( G 1 ) + i = 1 n 1 j = 1 i 1 ( ( j + 1 ) U j + 1 V j + 1 ) + ( n 1 ) A 1 + i = 1 n 1 [ ( 8 k 3 + 20 k 2 + 12 k + 2 ) n + ( 4 k 3 4 k 1 ) ]
Since both A 1 and i = 1 n 1 [ ( 8 k 3 + 20 k 2 + 12 k + 2 ) n + ( 4 k 3 4 k 1 ) ] are constant terms and both are independent of i, G u t ( G n ) can be written as
G u t ( G n ) = G u t ( G 1 ) + i = 1 n 1 j = 1 i 1 ( ( j + 1 ) U j + 1 V j + 1 ) + O ( n 2 )
Directly calculating shows that V a r ( U j ) = σ 2 , V a r ( V j ) = σ ˜ 2 , and C o v ( U j , V j ) = r , and C o v ( X , Y ) : = E ( X Y ) E ( X ) E ( Y ) denotes the covariance of any two random variables X and Y. Applying the properties of variance and the rule of interchanging the order of sums yields the following
V a r ( G u t ( G n ) ) = V a r ( i = 1 n 1 j = 1 i 1 ( ( j + 1 ) U j + 1 V j + 1 ) ) = V a r ( j = 1 n 2 i = j + 1 n 1 ( ( j + 1 ) U j + 1 V j + 1 ) ) = V a r ( j = 1 n 2 ( ( j + 1 ) U j + 1 V j + 1 ) ( n j 1 ) ) = j = 1 n 2 ( n j 1 ) 2 V a r ( ( j + 1 ) U j + 1 V j + 1 ) = j = 1 n 2 ( n j 1 ) 2 C o v ( ( j + 1 ) U j + 1 V j + 1 , ( j + 1 ) U j + 1 V j + 1 ) = j = 1 n 2 ( n j 1 ) 2 ( ( j + 1 ) 2 C o v ( U j + 1 , U j + 1 ) 2 ( j + 1 ) C o v ( U j + 1 , V j + 1 ) + C o v ( V j + 1 , V j + 1 ) ) = j = 1 n 2 ( n j 1 ) 2 ( ( j + 1 ) 2 σ 2 2 ( j + 1 ) r + σ ˜ 2 )
By utilizing computational, ad hoc tools (see Appendix A for details), the aforementioned equality suggests the intended outcome of V a r ( G u t ( G n ) ) .
(ii) We are now ready to demonstrate the second part of this theorem. First, for any n N , we set
U n = i = 1 n 1 j = 1 i 1 ( j + 1 ) U j + 1 , V n = i = 1 n 1 j = 1 i 1 V j + 1
and μ = E ( U j ) , ϕ ( t ) = E exp t ( U j μ ) . With these notations, we can explicitly obtain that
E exp t ( U n E ( U n ) ) = E exp t ( i = 1 n 1 j = 1 i 1 ( j + 1 ) ( U j + 1 μ ) ) = E exp t ( j = 1 n 2 i = j + 1 n 1 ( j + 1 ) ( U j + 1 μ ) ) = E exp t ( j = 1 n 2 ( n j 1 ) ( j + 1 ) ( U j + 1 μ ) ) = Π j = 1 n 2 E exp t ( n j 1 ) ( j + 1 ) ( U j + 1 μ ) = Π j = 1 n 2 ϕ ( t ( n j 1 ) ( j + 1 ) )
There exists C > 0 such that
V n = i = 1 n 1 j = 1 i 1 V j + 1 C n 2 .
If without confusion, we consider f ( n ) , g ( n ) two functions of n, we can define f ( n ) g ( n ) if lim n f ( n ) g ( n ) = 1 , and for some C > 0 , we denote f ( n ) = O ( g ( n ) ) if lim n ¯ f ( n ) g ( n ) C . Considering the relationship V a r ( G u t ( G n ) ) 1 30 σ 2 n 5 and the above equations, by using Taylor’s formula, one can deduce that
lim n E exp t G u t ( G n ) E ( G u t ( G n ) ) V a r ( G u t ( G n ) ) = lim n E exp t ( G u t ( G 1 ) + U n V n + O ( n 2 ) ) E ( G u t ( G 1 ) + U n V n + O ( n 2 ) ) σ n 5 2 30 = lim n E exp t 30 ( U n E ( U n ) ) σ n 5 2 = lim n Π j = 1 n ϕ 30 t ( n j 1 ) ( j + 1 ) σ n 5 2 = lim n E exp j = 1 n ln ϕ 30 t ( n j 1 ) ( j + 1 ) σ n 5 2
where we employ ϕ ( t ) = 1 + σ 2 2 t 2 + o ( t 2 ) such that
lim n E exp t G u t ( G n ) E ( G u t ( G n ) ) V a r ( G u t ( G n ) ) = lim n E exp j = 1 n ln 1 + σ 2 2 30 t 2 ( n j 1 ) 2 ( j + 1 ) 2 σ 2 n 5 + o ( 1 n ) = lim n E exp j = 1 n 2 σ 2 2 30 t 2 ( n j 1 ) 2 ( j + 1 ) 2 σ 2 n 5 + o ( 1 n ) = e t 2 2
where we applied j = 1 n 2 ( n j 1 ) 2 ( j + 1 ) 2 n 5 30 n 2 + 59 n 30 30 1 30 n 5 , which can be confirmed by computational methods. When i is in the complex field and satisfies the condition i 2 = 1 , by substituting t with i t in the above equation, we obtain
lim n E exp t i G u t ( G n ) E ( G u t ( G n ) ) V a r ( G u t ( G n ) ) = e t 2 2
Based on the aforementioned equality in Chapter 1 of [29] and the continuity theory in Chapter 15 of [30] for characteristic functions in probability, we have successfully provided proof of Theorem 1 (ii).

3. The Limiting Behavior of the Schultz Index of a Random Polygon Chain

We will now delve into the variance of the Schultz index for a random polygonal chain. In [27], Theorem 2.3 states that:
E ( S ( G n ) ) = { ( 8 k 3 + 12 k 2 + 4 k ) 2 i = 1 k 1 [ 4 k 3 ( 4 i 2 ) k 2 2 i k ] p i } n 3 3 { ( 2 k 2 + 2 k ) 2 i = 1 k 1 [ 4 k 3 ( 4 i 2 ) k 2 2 i k ] p i } n 2 + { ( 4 k 3 6 k 2 + 2 k ) 2 · 2 i = 1 k 1 [ 4 k 3 ( 4 i 2 ) k 2 2 i k ] p i } n 3 .
Next, we present our findings.
Theorem 2.
Given Hypothesis 1, the following results are valid:
  • (i) The Schultz index of G n , which is a random 2 k –sided chain, is S ( G n ) , and the variance in it is determined by
V a r ( S ( G n ) ) = 1 30 ( σ 2 n 5 5 r n 4 + 10 σ ˜ 2 n 3 + ( 65 r 30 σ 2 45 σ ˜ 2 ) n 2 + ( 120 r + 59 σ 2 + 65 σ ˜ 2 ) n + ( 60 r 30 σ 2 30 σ ˜ 2 ) ) .
where
σ 2 = i = 1 k ( c i 2 p i ) ( i = 1 k ( c i p i ) ) 2 σ ˜ 2 = i = 1 k ( d i 2 p i ) ( i = 1 k ( d i p i ) ) 2 r = i = 1 k ( c i d i p i ) i = 1 k ( c i p i ) i = 1 k ( d i p i )
with
c i = ( i + 1 ) ( 16 k 2 + 8 k ) d i = i ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k i = 1 , 2 , 3 , , k
(ii) As n , the limiting behavior of S ( G n ) converges to a normal distribution. In other words,
lim n sup x R P ( ( S ( G n ) E ( S ( G n ) ) V a r ( S ( G n ) ) ) x ) x 1 2 π e t 2 2 d u = 0
Proof. 
(i) As concluded in [27], we can derive that
S ( G n + 1 ) = S ( G n ) + ( 4 k + 2 ) v V G n d ( u n , v ) + 2 k v V G n d ( v ) d ( u n , v ) + ( 8 k 3 + 18 k 2 + 6 k ) n + ( 4 k 3 2 k ) .
Let
B n : = v V G n [ ( 4 k + 2 ) + 2 k d ( v ) ] d ( u n , v ) .
and we can then obtain
S ( G n + 1 ) = S ( G n ) + B n + ( 8 k 3 + 18 k 2 + 6 k ) n + ( 4 k 3 2 k ) .
 □
With regard to the random variables Z n 1 , Z n 2 , Z n 3 , ⋯, Z n k mentioned in Section 1, this signifies the approach we have taken in constructing G n + 1 from G n . Consequently, we can establish the following four equalities:
  • Equality 1.
B n Z n 1 = { B n 1 + 2 ( 16 k 2 + 8 k ) n [ 2 ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] } Z n 1
The above equality is trivial if Z n 1 = 0 , so our focus lies on the case where Z n 1 = 1 , indicating the transition G n G n + 1 1 . In this scenario, vertex u n in G n coincides with either vertex x 2 or x 2 k in H n , as shown in Figure 4. As a result, B n can be expressed as
B n = v V G n [ ( 4 k + 2 ) + 2 k d ( v ) ] d ( x 2 , v ) = v V G n 1 [ ( 4 k + 2 ) + 2 k d ( v ) ] d ( x 2 , v ) + v V H n [ ( 4 k + 2 ) + 2 k d ( v ) ] d ( x 2 , v ) = v V G n 1 [ ( 4 k + 2 ) + 2 k d ( v ) ] ( d ( u n 1 , v ) + 2 ) + ( 4 k + 2 ) k 2 + 2 k ( 2 k 2 + 1 ) = v V G n 1 [ ( 4 k + 2 ) + 2 k d ( v ) ] d ( u n 1 , v ) + 2 v V G n 1 [ ( 4 k + 2 ) + 2 k d ( v ) ] + ( 4 k + 2 ) k 2 + 2 k ( 2 k 2 + 1 ) = v V G n 1 [ ( 4 k + 2 ) + 2 k d ( v ) ] d ( u n 1 , v ) + 2 { ( 4 k + 2 ) 2 k ( n 1 ) + 2 k [ ( 4 k + 2 ) ( n 1 ) 1 ] } + ( 4 k + 2 ) k 2 + 2 k ( 2 k 2 + 1 ) = B n 1 + 2 ( 16 k 2 + 8 k ) n 2 ( 16 k 2 + 8 k ) + 8 k 3 + 2 k 2 + 2 k = B n 1 + 2 ( 16 k 2 + 8 k ) n [ 2 ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ]
By employing Equations (6)–(8) mentioned above, we are able to establish the desired equality.
  • Equality 2.
B n Z n 2 = { B n 1 + 3 ( 16 k 2 + 8 k ) n [ 3 ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] } Z n 2
Similar to the proof of Equality 1, we solely focus on the scenario where Z n 2 = 1 , resulting in the transformation G n G n + 1 2 . The proof follows a similar approach, and we have chosen to omit the specific details.
Similarly, when reaching the vertex labeled as x k + 1 in H n , we just only consider the case Z n k = 1 , which corresponds to G n G n + 1 k . In doing so, we obtain the following result:
  • Equality 3.
B n Z n k = { B n 1 + ( k + 1 ) ( 16 k 2 + 8 k ) n [ k ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] } Z n k
The proof follows the same steps as the previous equalities, and we will omit the details. Therefore, we can derive the general form of the analytic formula for A n Z n i as described below:
  • Equality 4.
B n Z n i = { B n 1 + ( i + 1 ) ( 16 k 2 + 8 k ) n [ i ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] } Z n i
where i 1 , 2 , , k . We will restrict our analysis to the case where Z n i = 1 , corresponding to G n G n + 1 i . The proof follows a similar approach as for the previously mentioned equalities, and we will skip the details.
Taking into account the observation that Z n 1 + Z n 2 + + Z n i + + Z n k = 1 , as mentioned earlier, we can conclude that
B n = B n ( Z n 1 + Z n 2 + + Z n i + + Z n k ) = { B n 1 + 2 ( 16 k 2 + 8 k ) n [ 1 ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] } Z n 1 + { B n 1 + 3 ( 16 k 2 + 8 k ) n [ 2 ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] } Z n 2 + + { B n 1 + ( i + 1 ) ( 16 k 2 + 8 k ) n [ i ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] } Z n i + + { B n 1 + ( k + 1 ) ( 16 k 2 + 8 k ) n [ k ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] } Z n k = B n 1 ( Z n 1 + Z n 2 + + Z n i + + Z n k ) + [ 2 ( 16 k 2 + 8 k ) Z n 1 + 3 ( 16 k 2 + 8 k ) Z n 2 + + ( i + 1 ) ( 16 k 2 + 8 k ) Z n i + + ( k + 1 ) ( 16 k 2 + 8 k ) Z n k ] n { [ 1 ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] Z n 1 + [ 2 ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] Z n 2 + + [ i ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] Z n i + + [ k ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] Z n k } = B n 1 + i = 1 k ( i + 1 ) ( 16 k 2 + 8 k ) Z n i i = 1 k [ i ( 16 k 2 + 8 k ) 8 k 3 + 14 k 2 + 10 k ] Z n i = B n 1 + i = 1 k c i Z n i i = 1 k d i Z n i = B n 1 + n U n V n
where for each n, i { 1 , 2 , 3 , , k } ,
U n = i = 1 k c i Z n i , V n = i = 1 k d i Z n i .
From Equation (12), we can deduce
S ( G n ) = S ( G 1 ) + i = 1 n 1 B i + i = 1 n 1 [ ( 8 k 3 + 18 k 2 + 6 k ) n + ( 4 k 3 2 k ) ] = S ( G 1 ) + i = 1 n 1 ( j = 1 i 1 ( B j + 1 B j ) + B 1 ) + i = 1 n 1 [ ( 8 k 3 + 18 k 2 + 6 k ) n + ( 4 k 3 2 k ) ] = S ( G 1 ) + i = 1 n 1 j = 1 i 1 ( B j + 1 B j ) + ( n 1 ) B 1 + i = 1 n 1 [ ( 8 k 3 + 18 k 2 + 6 k ) n + ( 4 k 3 2 k ) ] = S ( G 1 ) + i = 1 n 1 j = 1 i 1 ( ( j + 1 ) U j + 1 V j + 1 ) + ( n 1 ) B 1 + i = 1 n 1 [ ( 8 k 3 + 18 k 2 + 6 k ) n + ( 4 k 3 2 k ) ] = S ( G 1 ) + i = 1 n 1 j = 1 i 1 ( ( j + 1 ) U j + 1 V j + 1 ) + O ( n 2 )
Directly calculating shows that V a r ( U j ) = σ 2 , V a r ( V j ) = σ ˜ 2 , and C o v ( U j , V j ) = r , and C o v ( X , Y ) : = E ( X Y ) E ( X ) E ( Y ) denote the covariance of any two random variables X and Y. Applying the properties of variance and the rule of interchanging the order of sums yields the result that
V a r ( S ( G n ) ) = V a r ( i = 1 n 1 j = 1 i 1 ( ( j + 1 ) U j + 1 V j + 1 ) ) = V a r ( j = 1 n 2 i = j + 1 n 1 ( ( j + 1 ) U j + 1 V j + 1 ) ) = V a r ( j = 1 n 2 ( ( j + 1 ) U j + 1 V j + 1 ) ( n j 1 ) ) = j = 1 n 2 ( n j 1 ) 2 V a r ( ( j + 1 ) U j + 1 V j + 1 ) = j = 1 n 2 ( n j 1 ) 2 C o v ( ( j + 1 ) U j + 1 V j + 1 , ( j + 1 ) U j + 1 V j + 1 ) = j = 1 n 2 ( n j 1 ) 2 ( ( j + 1 ) 2 C o v ( U j + 1 , U j + 1 ) 2 ( j + 1 ) C o v ( U j + 1 , V j + 1 ) + C o v ( V j + 1 , V j + 1 ) ) = j = 1 n 2 ( n j 1 ) 2 ( ( j + 1 ) 2 σ 2 2 ( j + 1 ) r + σ ˜ 2 )
Using the same procedure as used in Gutman index, we obtain V a r ( S ( G n ) ) .
(ii) Similarly to the aforementioned derivation process of the Gutman index, we can obtain that
lim n E exp t i S ( G n ) E ( S ( G n ) ) V a r ( S ( G n ) ) = e t 2 2

4. Concluding Remarks

In this paper, we obtain explicit analytical expressions for the variances of the Gutman index and the Schultz index of a random polygonal chain composed of n 2k-sided polygons. Furthermore, we proved that the random variables G u t ( G n ) and ( S ( G n ) approach normal distributions asymptotically as n . These findings will significantly contribute to our comprehension of the inherent patterns and trends exhibited by the Gutman and Schultz indices in random-sided graphs. By further elucidating the limiting behaviors, we can gain deeper insights into the structural characteristics and properties of these graphs. Furthermore, the continuous advancement and progress of science are driving forces behind the continuous discovery and creation of new molecules, leading to advancements in various fields and applications. Researchers are extensively studying polygonal chains in chemical graph theory. Establishing exact formulas for the variances of certain indices in a random polygonal chain composed of n regular polygons is within the realm of possibility. Furthermore, we can anticipate obtaining the additional behaviors of the random-sided chain graph and its physicochemical properties in the near future through these studies. We are able not only to obtain precise expressions for the variance of random-sided graph chains but also to study their asymptotic behavior. As a result of these discoveries, future research in this field will undergo remarkable growth and expansion.

Author Contributions

Investigation, S.T.; Writing—original draft, C.T.; Supervision, X.G. All authors have read and agreed to the published version of the manuscript.

Funding

The research is partially supported by National Science Foundation of China (Grant No. 11671164), the Natural Science Foundation of Anhui Province (Grant No. 2008085MA01), and the Anhui University of Science and Technology Graduate Student Innovation Fund Program (Grant No. 2023CX2147).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

matlab codes; syms n k d1 r d2; f = ( n k 1 ) 2 ( ( k + 1 ) 2 d 1 2 2 ( k + 1 ) r + d 2 2 ) ;   f 2 = s i m p l i f y ( s y m s u m ( f , k , 1 , n 2 ) ) ; c o l l e c t ( f 2 , n ) .

References

  1. Bondy, J.A.; Murty, U.S.R. Graph Theory; Springer: New York, NY, USA, 2008. [Google Scholar]
  2. Flower, D.R. On the Properties of Bit String-Based Measures of Chemical Similarity. J. Chem. Inf. Comput. Sci. 1998, 38, C379–C386. [Google Scholar] [CrossRef]
  3. Pavlyuchko, A.I.; Vasiliev, E.V.; Gribov, L.A. Quantum chemical estimation of the overtone contribution to the IR spectra of hydrocarbon halogen derivatives. J. Struct. Chem. 2010, 51, C1045–C1051. [Google Scholar] [CrossRef]
  4. Hosoya, H. A Newly Proposed Quantity Characterizing the Topological Nature of Structural Isomers of Saturated Hydrocarbons. J. Bull. Chem. Soc. Jpn. 2006, 44, 2332–2339. [Google Scholar] [CrossRef]
  5. Estrada, E. Edge Adjacency Relationships and a Novel Topological Index Related to Molecular Volume. J. Chem. Inf. Comput. Sci. 1995, 35, 31–33. [Google Scholar] [CrossRef]
  6. Wiener, H. Structrual determination of paraffin boiling points. J. Am. Chem. Soc. 1947, 69, 17–20. [Google Scholar] [CrossRef] [PubMed]
  7. Buckley, F.; Harary, F. Distance in Graphs; Addison-Wesley: Reading, PA, USA, 1989. [Google Scholar]
  8. Entringer, R.C.; Jackson, D.E.; Snyder, D.A. Distance in graphs. Czechoslov. Math. J. 1976, 26, 283–296. [Google Scholar] [CrossRef]
  9. Chen, A.L.; Zhang, F.J. Wiener index and perfect matchings in random phenylene chains. MATCH Commun. Math. Comput. Chem. 2009, 61, 623–630. [Google Scholar]
  10. Zhou, Q.N.; Wang, L.G.; Lu, Y. Wiener index and Harary index on Hamiltonconnected graphs with large minimum degree. J. Discret. Appl. Math. 2018, 247, 180–185. [Google Scholar] [CrossRef]
  11. Ayache, A.; Alameri, A. Topological indices of the mk-graph. J. Assoc. Arab. Univ. Basic Appl. Sci. 2017, 24, 283–291. [Google Scholar]
  12. Gutman, I. Selected properties of the Schultz molecular topological index. J. Chem. Inf. Comput. Sci. 1994, 34, 1087–1089. [Google Scholar] [CrossRef]
  13. Chen, S. Modified Schultz index of zig-zag polyhex nanotubes. J. Comput. Theor. Nanosci. 2009, 6, 1499–1503. [Google Scholar] [CrossRef]
  14. Farahani, M.R. Hosoya, Schultz, modified Schultz polynomials and their topological indices of benzene molecules: First members of polycyclic aromatic hydrocarbons (PAHs). Int. J. Theor. Chem. 2013, 1, 9–16. [Google Scholar]
  15. Heydari, A. On the modified Schultz index of C4C8(S) nanotubes and nanotorus. Digest. J. Nanomater. Biostruct. 2010, 5, 51–56. [Google Scholar]
  16. Xiao, Z.; Chen, S. The modified Schultz index of armchair polyhex nanotubes. J. Comput. Theor. Nanosci. 2009, 6, 1109–1114. [Google Scholar] [CrossRef]
  17. Mukwembi, S.; Munyira, S. MunyiraDegree distance and minimum degree. Bull. Aust. Math. Soc. 2013, 87, 255–271. [Google Scholar] [CrossRef]
  18. Person, W.B.; Pimentel, G.C.; Pitzer, K.S. The Structure of Cyclooctatetraene. J. Am. Chem. Soc. 1952, 74, 3437–3438. [Google Scholar] [CrossRef]
  19. Cope, A.C.; Pike, R.M.; Polyolefins, C. Cyclooctatetraene Derivatives from Copolymerization and Side Chain Modification. J. Am. Chem. Soc. 1953, 75, 3220–3223. [Google Scholar] [CrossRef]
  20. Donald, H.L.; Whitehead, M.A. Molecular geometry and bond energy. III. Cyclooctatetraene and related compounds. J. Am. Chem. Soc. 1969, 91, 238–242. [Google Scholar]
  21. Serbezeanu, D.; Popa, A.M.; Sava, I.; Carja, I.D.; Amberg, M.; Rossi, R.M.; Fortunato, G. Design and synthesis of polyimide—Gold nanofibers with tunable optical properties. J. Eur. Polym. J. 2015, 64, 10–20. [Google Scholar] [CrossRef]
  22. Bonchev, D.; Rouvray, D.H. Chemical Graph Theory: Introduction and Fundamentals. J. Math. Chem. 1991, 7, 5971–5987. [Google Scholar]
  23. Zhang, L.L.; Li, Q.S.; Li, S.C.; Zhang, M.J. The expected values for the Schultz index, Gutman index, multiplicative degree-Kirchhoff index and additive degree-Kirchhoff index of a random polyphenylene chain. J. Discret. Applited Math. 2020, 282, 243–256. [Google Scholar] [CrossRef]
  24. Qi, J.; Ni, J.; Geng, X. The expected values for the Kirchhoff indices in the random cyclooctatetraene and spiro chains. J. Discret. Applited Math. 2022, 321, 240–249. [Google Scholar] [CrossRef]
  25. Zhang, J.L.; Peng, X.H.; Chen, H.L. The limiting behaviours for the Gutman index, Schultz index, multiplicative degree-Kirchhoff index and additive degree-Kirchhoff index of a random polyphenylene chain. Discrete Appl. Math. 2021, 299, C62–C73. [Google Scholar] [CrossRef]
  26. Tao, C.; Tang, S.J.; Geng, X.Y. Statistical Analyses of a Class of Random Cyclooctatetraene Chain Networks with Respect to Several Topological Properties. Symmetry 2023, 15, 1971. [Google Scholar] [CrossRef]
  27. Zhu, W.L.; Fang, M.L.; Geng, X.Y. Enumeration of the Gutman and Schultz indices in the random polygonal chains. Math. Biosci. Eng. 2022, 19, 10826–10845. [Google Scholar] [CrossRef] [PubMed]
  28. Gutman, I. Degree-Based Topological Indices. Croat. Chem. Acta 2013, 86, 351–361. [Google Scholar] [CrossRef]
  29. Petrov, V.V. Limit Theorems of Probability Theory Sequences of Independent Random Variables; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
  30. Feller, W. An Introduction to Probability Theory and Its Applications; Wiley Series in Probability and Statistics; John Wiley & Sons: Hoboken, NJ, USA, 1971. [Google Scholar]
Figure 1. Graph G 2 .
Figure 1. Graph G 2 .
Axioms 13 00518 g001
Figure 2. The formation of G n + 1 by combining G n and H n + 1 .
Figure 2. The formation of G n + 1 by combining G n and H n + 1 .
Axioms 13 00518 g002
Figure 3. k methods for attaching the new terminal 2k-gon polygon H n + 1 to G n .
Figure 3. k methods for attaching the new terminal 2k-gon polygon H n + 1 to G n .
Axioms 13 00518 g003
Figure 4. G n G n + 1 1 .
Figure 4. G n G n + 1 1 .
Axioms 13 00518 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tao, C.; Tang, S.; Geng, X. The Limiting Behaviors of the Gutman and Schultz Indices in Random 2k-Sided Chains. Axioms 2024, 13, 518. https://doi.org/10.3390/axioms13080518

AMA Style

Tao C, Tang S, Geng X. The Limiting Behaviors of the Gutman and Schultz Indices in Random 2k-Sided Chains. Axioms. 2024; 13(8):518. https://doi.org/10.3390/axioms13080518

Chicago/Turabian Style

Tao, Chen, Shengjun Tang, and Xianya Geng. 2024. "The Limiting Behaviors of the Gutman and Schultz Indices in Random 2k-Sided Chains" Axioms 13, no. 8: 518. https://doi.org/10.3390/axioms13080518

APA Style

Tao, C., Tang, S., & Geng, X. (2024). The Limiting Behaviors of the Gutman and Schultz Indices in Random 2k-Sided Chains. Axioms, 13(8), 518. https://doi.org/10.3390/axioms13080518

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