Next Article in Journal
Complete Quantum Information in the DNA Genetic Code
Previous Article in Journal
Statistical Validation Verifies That Enantiomorphic States of Chiral Cells Are Determinant Dictating the Left- or Right-Handed Direction of the Hindgut Rotation in Drosophila
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Topological and Thermodynamic Entropy Measures for COVID-19 Pandemic through Graph Theory

by
G. Kirithiga Nandini
1,
R. Sundara Rajan
2,
A. Arul Shantrinal
2,
T. M. Rajalaxmi
3,
Indra Rajasingh
4 and
Krishnan Balasubramanian
5,*
1
Department of Computer Science and Engineering, Hindustan Institute of Technology and Science, Chennai 603 103, India
2
Department of Mathematics, Hindustan Institute of Technology and Science, Chennai 603 103, India
3
Department of Mathematics, Sri Sivasubramaniya Nadar College of Engineering, Chennai 603 110, India
4
School of Advanced Sciences, Vellore Institute of Technology, Chennai 600 127, India
5
School of Molecular Sciences, Arizona State University, Tempe, AZ 85287-1604, USA
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(12), 1992; https://doi.org/10.3390/sym12121992
Submission received: 14 October 2020 / Revised: 27 November 2020 / Accepted: 30 November 2020 / Published: 2 December 2020
(This article belongs to the Section Life Sciences)

Abstract

:
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the global pandemic, coronavirus disease-2019 (COVID-19) which has resulted in 60.4 million infections and 1.42 million deaths worldwide. Mathematical models as an integral part of artificial intelligence are designed for contact tracing, genetic network analysis for uncovering the biological evolution of the virus, understanding the underlying mechanisms of the observed disease dynamics, evaluating mitigation strategies, and predicting the COVID-19 pandemic dynamics. This paper describes mathematical techniques to exploit and understand the progression of the pandemic through a topological characterization of underlying graphs. We have obtained several topological indices for various graphs of biological interest such as pandemic trees, Cayley trees, Christmas trees, and the corona product of Christmas trees and paths. We have also obtained an analytical expression for the thermodynamic entropies of pandemic trees as a function of R0, the reproduction number, and the level of spread, using the nested wreath product groups. Our plots of entropy and logarithms of topological indices of pandemic trees accentuate the underlying severity of COVID-19 over the 1918 Spanish flu pandemic.

1. Introduction

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), commonly called the novel coronavirus (COVID-19), has resulted in 1.42 million deaths worldwide up to now [1,2]. The first case was detected in the city of Wuhan, China on 31 December 2019, and it was reported by Zhou et al. [3] as a pneumonia outbreak associated with a new coronavirus of bat origin. It is presently understood that the disease spreads by human-to-human or aerosol transmission via either droplet in the air or through direct or indirect contacts. An underestimated mean incubation period of 2.24–3.58 days [4] has now been revised to a median R0 value close to 5.7 in some regions of a susceptible population. This dramatically contrasts with the previous Spanish flu virus which has a R0 value close to 1.8. There are seven types of coronaviruses that infect humans [5], including the newly discovered SARS-CoV-2 virus. The level and duration of infectious virus replication are important factors in the risk assessments [6]. A study on interspecies transmission of the virus and its genetic diversity can provide mitigation strategies against the infection [7].
Biological models can simulate the dynamics of the spread and genetic evolution of a disease as to how cells interact in a single patient and how genetic mutations spread across people living in different geographical locations [2,8]. Forster et al. [9] recently considered phylogenetic networks as powerful tools for understanding the genetic evolution of the coronavirus 2019. Graph models also enhance our understanding of the effect of introducing an affected individual in a pool of unaffected individuals. Recursive trees constructed from phylogenetics have provided graphical representation of microbiomes and have been proven to be powerful in the perturbations induced to genomes and proteomes by environment or toxins [10,11,12,13], and hence, graph techniques are especially suited to understand the dynamics of the COVID-19 epidemic. Furthermore, it has been shown that tree pruning methods and other efficient topological techniques can be applied to characterize various biological trees such as phylogenetic trees, epidemic trees, Cayley trees, and the lattices through graph theoretical entities [14,15,16,17]. The combinatorial theory of genetic mutations and the applications of graph theory to the computation of evolutionary genomics have been pioneered by Sellers and others [12,13,18,19,20]. Yun et al. [21] demonstrated the graph mining approaches in identifying the patterns in complex and large data sets such as the chemical and genome data sets.
Network models help to capture the transmission, forecast the future of the pandemic, and estimate the sequel of ongoing interventions [22]. A pandemic tree provides a pictorial representation of the epidemic dynamics where each node is connected to k other nodes, where k is an integer rounded from the R0 value, an epidemiological measure of the degree of infection, suggesting that an infected individual in turn infects R0 others in a susceptible population pool. Consequently, quantitative measures of a pandemic network can provide significant new insights into the epidemic dynamics. Topological indices have been developed that aid in the prediction of chemical, physical, pharmaceutical, and biological properties [23,24,25,26,27,28,29]. The eccentric connectivity index is one such topological index that is currently being used for modeling biological activities of chemical compounds. In anticonvulsant, anti-inflammatory, and diuretic activities, this index exhibits high degree of predictability [30]. Ghorbani et al. [31] obtained various eccentricity–entropy-based topological measures of fullerenes.
We consider several trees of biological and phylogenetic interests and corona products of graphs in order to obtain the topological measures of the associated networks. In particular, we consider several classes of recursive networks of biological importance, namely pandemic trees, Cayley trees, Christmas trees, and the corona products of Christmas trees and paths. In Section 2, basic definitions are stated, while Section 3 computes the eccentricity-based topological indices. Section 4 highlights various applications of topological indices to the COVID-19 pandemic, while in Section 5, we study the thermodynamic entropy of pandemic trees. Section 6 discusses stochasticity in pandemic tree generation. The conclusion is given in Section 7.

2. Basic Concepts

Let G = (V(G), E(G)) be a graph with V(G) as the set of vertices and E(G) as the set of edges. A tree is a connected, acyclic graph. In a binary tree, each node can have at most two descendants. It is said to be a complete binary tree if each internal node has exactly two descendants.
The length of the shortest path between two vertices x and y in a graph G is the topological distance between x and y and it is denoted by d(x,y). The eccentricity value of a vertex u V ( G ) denoted by ec(u) is the longest distance between the vertex u and any other vertex v of G. The neighborhood of a vertex u, denoted by n(u), is the set of all vertices which are adjacent to u. Let d(u) denote the degree of a vertex u. For an edge u v E ( G ) ,   let nu(e) be the number of vertices lying closer to the vertex u than the vertex v and let nv(e) denote the number of vertices lying closer to the vertex v than the vertex u.
The fourth geometric-arithmetic eccentricity index was introduced by Ghorbani and Khaki [32]: as
G A 4 ( G ) = u v E ( G ) 2 e c ( u ) e c ( v ) e c ( u ) + e c ( v ) .  
Furthermore, the fourth Zagreb index, the fourth multiplicative Zagreb index, the sixth Zagreb index, and the sixth multiplicative Zagreb index [33,34] are defined as follows:
Z g 4 ( G ) = u v E ( G ) ( e c ( u ) + e c ( v ) ) ; Π 4 * ( G ) = u v E ( G ) ( e c ( u ) + e c ( v ) ) ; Z g 6 ( G ) = u v E ( G ) e c ( u )   e c ( v ) ; Π 6 * ( G ) = u v E ( G ) e c ( u )   e c ( v ) .
The fourth Zagreb polynomial Z g 4 ( G ,   x ) and the sixth Zagreb polynomial Z g 6 ( G ,   x ) [33,34] are defined by:
Z g 4 ( G ,   x ) = u v E ( G ) x e c ( u ) +   e c ( v ) ; Z g 6 ( G ,   x ) = u v E ( G ) x e c ( u )   e c ( v ) .
The fifth multiplicative atom bond connectivity index [33] is defined as:
A B C 5 Π ( G ) = u v E ( G ) e c ( u ) + e c ( v ) 2 e c ( u )   e c ( v ) .
The first Zagreb index M1(G) [35] is defined as:
M 1 ( G ) = v i V ( G ) d i 2 .
The second Zagreb index M2(G) [35] of graph G is defined as:
M 2 ( G ) = v i v j E ( G ) d i d j .
The Atom-bond Connectivity index (ABC) [36] is defined as:
A B C ( G ) = v i v j E ( G )   d i + d j 2 d i d j .
The Padmakar–Ivan (PI) index [37] is defined as:
P I ( G ) = e = u v E ( G )   [ n u ( e ) + n v ( e ) ] .
The Szeged Index [38] is defined as:
S z ( G ) = e = u v E ( G )   [ n u ( e ) × n v ( e ) ] .

3. Main Results for the Topological Indices

In epidemiology, R0, the reproduction number measures the severity of the pandemic, and it is defined as the number of individuals an infected person can in turn infect in a group of susceptible population set. For example, Figure 1 shows a pandemic tree for an epidemic with an R0 value of 4. In this section, we have obtained the eccentricity-related indices for the pandemic tree network (also called a complete k-ary tree network), Cayley tree network, Christmas trees, and the corona product of Christmas tree and a path, all of which are defined in the respective sections.

3.1. Topological Indices of Pandemic Trees

A pandemic tree is a complete R0-ary tree, where R0 is the reproduction number of a pandemic rounded to the nearest integer.
Definition 1 ([39]).
A k -ary tree is a rooted tree in which each node has no more than k children. All children of one node can be described as descendants of this node. In a k -ary tree, the maximum distance l from leaves to the root node is called the height of the k -ary tree. The root node is said to be at level 0. Inductively, the descendants of nodes at level i are at level i + 1 , 0 i l 1 . A k -ary tree is said to be a complete k -ary tree if each internal node has exactly k descendants. Thus, a pandemic tree is a complete k -ary tree, with the rounded epidemiologicalR0 value of k . We denote this tree by T l k , where l is the height of the tree, k , l ≥ 2.
The number of vertices in level i ,   0 i l of   T l k is k i . Hence, the total number of vertices and edges in T l k are k l + 1 1 k 1 and k l + 1 1 k 1 1 , respectively. The degree d ( v ) of a vertex v in T l k   is   1 , if v is a leaf; k , if v is the root and k + 1 , if v is an internal vertex. For illustration, the pandemic tree T 5 3 of level 5, and the pandemic tree T 6 6 of level 6 are shown in Figure 2a,b, respectively.
Theorem 1.
Let T l k denote a pandemic tree on l levels with a reproduction number k.
Then:
(i) 
G A 4 ( T l k ) = i = 1 l k i × 2 ( l + i 1 ) ( l + i ) 2 ( l + i 1 ) + 1 ;
(ii) 
Z g 4 ( T l k ) = k ( k + 2 l 2 l k + 1 ) ( k 1 ) 2 k l + 1 ( k + 4 l 4 l k + 1 ) ( k 1 ) 2 ;
(iii) 
Π 4 * ( T l k ) = i = 1 l ( 2 ( l + i ) 1 ) k i ;
(iv) 
Z g 6 ( T l k ) = 2 k l + 1 ( 2 k 2 l 2 l k 2 4 k l 2 + k + 2 l 2 + l ) ( k 1 ) 3 k ( k 2 l 2 l k 2 2 k l 2 + 2 k + l 2 + l ) ( k 1 ) 3 ;
(v) 
Π 6 * ( T l k ) = i = 1 l ( ( l + i 1 ) ( l + i ) ) k i ;
(vi) 
Z g 4 ( T l k ,   x ) = i = 1 l k i × x 2 ( l + i 1 ) + 1 ;
(vii) 
Z g 6 ( T l k ,   x ) = i = 1 l k i × x ( l + i 1 ) ( l + i ) ;
(viii) 
A B C 5 ( T l k ) = i = 1 l ( 2 ( l + i 1 ) 1 ( l + i 1 ) ( l + i ) ) k i ;
(ix) 
P I ( T l k ) = k ( k l 1 ) ( k l + 1 1 ) ( k 1 ) 2 ;
(x) 
S z ( T l k ) = k l + 1 ( 2 k l 2 k × k l + k l k l × k l + k 2 + l × l ) ( k 1 ) 3 .
Proof. 
We prove the results based on the edge partition and structure analysis of T l k :
  • E l + i 1 ,   l + i = { e = u v E ( T l k ) | e c ( u ) = l + i 1   and   e c ( v ) = l + i } and n l + i 1 , l + i = | E l + i 1 , l + i | = k i , where i { 1 , , l } :
    and
  • E k + 1 , k = { e = u v E ( T l k ) |   d ( u ) = k + 1   and   d ( v ) = k ,   n ( u ) = ( k l 1 k 1 )   and   n ( v ) = ( k l + 1 1 k 1 ) ( k l 1 k 1 ) }   and   n k + 1 , k = | E k + 1 , k | = k ;
  • E 1 , k + 1 = { e = u v E ( T l k ) |   d ( u ) = 1   and   d ( v ) = k + 1 ,   n ( u ) = 1   and   n ( v ) = ( k l + 1 1 k 1 ) 1 }   and   | E 1 , k | = k l ;
  • E k + 1 , k + 1 = { e = u v E ( T l k ) |   d ( u ) = d ( v ) = k + 1 , n ( u ) = i = 1 l 2 ( k i + 1 1 k 1 ) and   n ( v ) = ( k l + 1 1 k 1 ) i = 1 l 2 ( k i + 1 1 k 1 ) }   and   | E k , k | = i = 1 l 2 ( k l i ) .
From the definitions we have:
G A 4 ( T l k ) = i = 1 l k i × 2 ( l + i 1 ) ( l + i ) 2 ( l + i 1 ) + 1 ; Z g 4 ( T l k ) = i = 1 l k i × 2 ( l + i 1 ) + 1 = k ( k + 2 l 2 l k + 1 ) ( k 1 ) 2 k l + 1 ( k + 4 l 4 l k + 1 ) ( k 1 ) 2 ; Π 4 * ( T l k ) = i = 1 l ( 2 ( l + i ) 1 ) k i ; Z g 6 ( T l k ) = i = 1 l k i × ( l + i 1 ) ( l + i ) = 2 k l + 1 ( 2 k 2 l 2 l k 2 4 k l 2 + k + 2 l 2 + l ) ( k 1 ) 3 k ( k 2 l 2 l k 2 2 k l 2 + 2 k + l 2 + l ) ( k 1 ) 3 ; Π 6 * ( T l k ) = i = 1 l ( ( l + i 1 ) ( l + i ) ) k i ; Z g 4 ( T l k ,   x ) = i = 1 l k i × x 2 ( l + i 1 ) + 1 ; Z g 6 ( T l k ,   x ) = i = 1 l k i × x ( l + i 1 ) ( l + i ) ; A B C 5 ( T l k ) = i = 1 l ( 2 ( l + i 1 ) 1 ( l + i 1 ) ( l + i ) ) k i .
Further:
P I ( T l k )   = ( ( k l 1 k 1 ) + ( ( k l + 1 1 k 1 ) ( k l 1 k 1 ) ) ) × k + ( 1 + ( ( k l + 1 1 k 1 ) 1 ) ) × k l + i = 1 l 2 ( k l i ) { ( k i + 1 1 k 1 ) + ( ( k l + 1 1 k 1 ) ( k i + 1 1 k 1 ) ) } = k ( k l 1 ) ( k l + 1 1 ) ( k 1 ) 2   ; S z ( T l k )   = ( ( k l 1 k 1 ) × ( ( k l + 1 1 k 1 ) ( k l 1 k 1 ) ) ) × k + ( 1 × ( ( k l + 1 1 k 1 ) 1 ) ) × k l + i = 1 l 2 ( k l i ) { ( k i + 1 1 k 1 ) × ( ( k l + 1 1 k 1 ) ( k i + 1 1 k 1 ) ) } = k l + 1 ( 2 k l 2 k × k l + k l k l × k l + k 2 + l × l ) ( k 1 ) 3 .
 □
Theorem 2.
The various degree based topological indices of T l k are given by:
(i) 
M 1 ( T l k ) = k 2 + ( k + 1 ) 2 × k l k k 1 + k l ;
(ii) 
M 2 ( T l k ) = k l ( k + 1 ) + k l k 2 k 1 ( k + 1 ) 2 + k 2 ( k + 1 ) ;
(iii) 
A B C ( T l k ) =   k 2 k 1 k ( k + 1 ) + k l k k + 1   + k l k 2 k 1 × 2 k k + 1 .
Proof. 
Let the vertex and edge sets of T l k be partitioned as follows:
  • P 1 = { v |   d ( v ) = k }   and   | P 1 | = 1 ;
  • P 2 = { v |   d ( v ) = k + 1 }   and   | P 2 | = k l k k 1 ;
  • P 3 = { v |   d ( v ) = 1 }   and   | P 3 | = k l ;
and
  • E k , k + 1 = { e = u v E ( T l k ) |   d ( u ) = k   and   d ( v ) = k + 1 }   and   | E k , k + 1 | = k ;
  • E 1 , k + 1 = { e = u v E ( T l k ) |   d ( u ) = 1   and   d ( v ) = k + 1 }   and   | E 1 , k + 1 | = k l ;
  • E k + 1 , k + 1 = { e = u v E ( T l k ) |   d ( u ) = k + 1   and   d ( v ) = k + 1 }   and   | E k + 1 , k + 1 | = i = 2 l 1 k i .
By the definition of first Zagreb index, we have:
M 1 ( T l k )   = v ϵ P 1 ( d ( v ) ) 2 + v ϵ P 2 ( d ( v ) ) 2 + v ϵ P 3 ( d ( v ) ) 2 = k 2 + ( k + 1 ) 2 × k l k k 1 + k l .
From the definitions of M 2 and A B C we obtain:
M 2 ( T l k )   = e ϵ E k , k + 1 d ( u ) d ( v ) + e ϵ E 1 , k + 1 d ( u ) d ( v ) + e ϵ E k + 1 , k + 1 d ( u ) d ( v ) = ( ( k ) ( k + 1 ) ) k + ( ( 1 ) ( k + 1 ) )   k l + ( k + 1 ) ( k + 1 )   k l k 2 2 = k l ( k + 1 ) +   k l k 2 k 1 ( k + 1 ) 2 + k 2 ( k + 1 ) ; A B C ( T l k )   = k ( k ) + ( k + 1 ) 2 ( k ) ( k + 1 ) + k l ( 1 ) + ( k + 1 ) 2 ( 1 ) ( k + 1 ) + i = 2 l 1 k i ( k + 1 ) + ( k + 1 ) 2 ( k + 1 ) ( k + 1 ) = k 2 k 1 k ( k + 1 ) + k l k k + 1   + k l k 2 k 1 × 2 k k + 1 .
 □
Remark 1.
The results obtained from Topo-Chemie-2020 [40] are shown in Table 1 for a comparison with the results obtained from Theorems 1 and 2.

3.2. Topological Indices of Cayley Trees

Definition 2.
A tree is said to be a k-Cayley tree with l levels, denoted by C(k, l), k ,   l 3 if all vertices have the same degree k except the leaves [41]. The number of vertices and edges in C(k, l) are k ( k 1 ) l 2 k 2 and k ( k 1 ) l 2 k 2 1 , respectively. For illustration, the Cayley trees C(3,6) and C(6,6) are shown in Figure 3a,b, respectively.
Theorem 3.
For a Cayley tree C ( k , l ) ,   w i t h   k ,   l 3 ,   w e   o b t a i :
(i) 
G A 4 ( C ( k , l ) ) = i = 1 l k ( k 1 ) i 1 × 2 ( l + i 1 ) ( l + i ) 2 ( l + i 1 ) + 1 ;
(ii) 
Z g 4 ( C ( k , l ) ) = ( k ( 4 l + k 2 l k 8 l ( k 1 ) l k ( k 1 ) l + 4 l k ( k 1 ) l ) ) ( k 2 ) 2 ;
(iii) 
Π 4 * ( C ( k , l ) ) = i = 1 l ( 2 ( l + i ) 1 ) k ( k 1 ) i 1 ;
(iv) 
  Z g 6 ( C ( k , l ) ) = k ( ( 4 l 2 k l 2 k 2 4 l 2 + l k 2 2 l k 2 k + 2 ( k 2 ) 3 ) + ( 2 ( k 1 ) l ( 2 l 2 k 2 8 l 2 k + 8 l 2 l k 2 + 2 l k + k 1 ) ( k 2 ) 3 ) ) ;
(v) 
  Π 6 * ( C ( k , l ) ) = i = 1 l ( ( l + i 1 ) ( l + i ) ) k ( k 1 ) i 1 ;
(vi) 
Z g 4 ( C ( k , l ) ,   x ) = i = 1 l k ( k 1 ) i 1 × x 2 ( l + i 1 ) + 1 ;
(vii) 
Z g 6 ( C ( k , l ) ,   x ) = i = 1 l k ( k 1 ) i 1 × x ( l + i 1 ) ( l + i ) ;
(viii) 
A B C 5 ( C ( k , l ) ) = i = 1 l ( 2 ( l + i 1 ) 1 ( l + i 1 ) ( l + i ) ) k ( k 1 ) i 1 ;
(ix) 
P I ( C ( k , l ) ) = ( k ( k ( k 1 ) l 2 ) ( ( k 1 ) l 1 ) ) ( k 2 ) 2 ;
(x) 
S z ( C ( k , l ) ) = ( k ( 2 k ( k 1 ) l 2 k ( k 1 ) 2 l + ( k 1 ) 2 l 2 k l ( k 1 ) 2 l + l k 2 ( k 1 ) 2 l 1 ) ) ( k 2 ) 3 .
Proof. 
We compute the edge partition with respect to eccentricities and the edge partitions with respect to the degrees of the associated vertices from the edge set of C ( k ,   l ) as follows:
  • E l + i 1 ,   l + i = { e = u v E ( C ( k , l ) ) | e c ( u ) = l + i 1   and   e c ( v ) = l + i } and n l + i 1 ,   l + i = | E l + i 1 ,   l + i | = k ( k 1 ) i 1 , where i { 1 , , l } :
and
  • E k , k = { e = u v E ( C ( k , l ) ) |   d ( u ) = d ( v ) = k ,   n ( u ) = ( ( k 1 ) l 1 k 2 )   and   n ( v ) = ( k ( k 1 ) l ( k 1 ) l 1 k 2 )   }   and | E k , k | = k ;
  • E 1 , k = { e = u v E ( C ( k , l ) ) |   d ( u ) = 1   and   d ( v ) = k ,   n ( u ) = 1   and   n ( v ) = ( k ( k 1 ) l 2 k 2 1 ) }   and   | E 1 , k | = k ( k 1 ) l 1 ;
  • E k , k = { e = u v E ( C ( k , l ) ) |   d ( u ) = d ( v ) = k , n ( u ) = i = 1 l 2 ( ( k 1 ) i + 1 1 k 2 ) and   n ( v ) = i = 1 l 2 ( k ( k 1 ) d 2 k 2 ( k 1 ) i + 1 1 k 2 ) }   and   | E k , k | = i = 1 l 2 k ( k 1 ) l i 1 .
From the definitions of the eccentric version indices, we obtain:
G A 4 ( C ( k , l ) ) = i = 1 l k ( k 1 ) i 1 × 2 ( l + i 1 ) ( l + i ) 2 ( l + i 1 ) + 1 ; Z g 4 ( C ( k , l ) ) = ( k ( 4 l + k 2 l k 8 l ( k 1 ) l k ( k 1 ) l + 4 l k ( k 1 ) l ) ) ( k 2 ) 2 ; Π 4 * ( C ( k , l ) ) = i = 1 l ( 2 ( l + i ) 1 ) k ( k 1 ) i 1 ; Z g 6 ( C ( k , l ) ) = k ( ( 4 l 2 k l 2 k 2 4 l 2 + l k 2 2 l k 2 k + 2 ( k 2 ) 3 ) + ( 2 ( k 1 ) l ( 2 l 2 k 2 8 l 2 k + 8 l 2 l k 2 + 2 l k + k 1 ( k 2 ) 3 ) ) ; Π 6 * ( C ( k , l ) ) = i = 1 l ( ( l + i 1 ) ( l + i ) ) k ( k 1 ) i 1 ; Z g 4 ( C ( k , l ) ,   x ) = i = 1 l k ( k 1 ) i 1 × x 2 ( l + i 1 ) + 1 ; Z g 6 ( C ( k , l ) ,   x ) = i = 1 l k ( k 1 ) i 1 × x ( l + i 1 ) ( l + i ) ; A B C 5 ( C ( k , l ) ) = i = 1 l ( 2 ( l + i 1 ) 1 ( l + i 1 ) ( l + i ) ) k ( k 1 ) i 1 .
By the definitions of the PI and Sz indices we have:
P I ( C ( k , l ) )   = ( ( ( k 1 ) l 1 k 2 ) + ( k ( k 1 ) l ( k 1 ) l 1 k 2 ) ) × k + ( 1 + ( k ( k 1 ) l 2 k 2 1 ) ) × k ( k 1 ) l 1 + i = 1 l 2 k ( k 1 ) l i 1 { ( ( k 1 ) i + 1 1 k 2 ) + ( ( k ( k 1 ) l 2 k 2 ) ( ( k 1 ) i + 1 1 k 2 ) ) } =   ( k ( k ( k 1 ) l 2 ) ( ( k 1 ) l 1 ) ) ( k 2 ) 2 ; S z ( C ( k , l ) )   = ( ( ( k 1 ) l 1 k 2 ) × ( k ( k 1 ) l ( k 1 ) l 1 k 2 ) ) × k + ( 1 × ( k ( k 1 ) l 2 k 2 1 ) ) × k ( k 1 ) l 1 + i = 1 l 2 k ( k 1 ) l i 1 { ( ( k 1 ) i + 1 1 k 2 ) × ( ( k ( k 1 ) l 2 k 2 ) ( ( k 1 ) i + 1 1 k 2 ) ) } = ( k ( 2 k ( k 1 ) l 2 k ( k 1 ) 2 l + ( k 1 ) 2 l 2 k l ( k 1 ) 2 l + l k 2 ( k 1 ) 2 l 1 ) ) ( k 2 ) 3 .
 □
Theorem 4.
The various degree based topological indices of C ( k , l ) are given by:
(i) 
M 1 ( C ( k , l ) ) = k ( k ( k 1 ) l 2 k + 2 ( k 1 ) l k 2 ;
(ii) 
  M 2 ( C ( k , l ) ) = k 2 ( 2 ( k 1 ) l k ) k 2 ;
(iii) 
A B C ( C ( k , l ) ) = ( k ( k 1 ) l k + 1 k 2 3 k + 2 ) 2 ( k 1 ) k + k ( k 1 ) l 1 k 1 k .
Proof. 
Partition the vertex and edge sets of C ( k , l ) as follows:
  • P 1 = { v |   d ( v ) = k }   and   | P 1 | = ( k ( k 1 ) l 1 2 k 2 ) ;
  • P 2 = { v |   d ( v ) = 1 }   and   | P 2 | = k ( k 1 ) l 1 :
and
  • E k , k = { e = u v E ( C ( k , l ) ) |   d ( u ) = k   and   d ( v ) = k }   and   | E k , k | = ( k ( k 1 ) l k + 1 k 2 3 k + 2 ) ;
  • E 1 , k = { e = u v E ( C ( k , l ) ) |   d ( u ) = 1   and   d ( v ) = k }   and   | E 1 , k | = k ( k 1 ) l 1 .
By the definition of M 1 :
M 1 ( C ( k , l ) )   = v ϵ P 1 ( d ( v ) ) 2 + v ϵ P 2 ( d ( v ) ) 2 = k 2 ( k ( k 1 ) l 1 2 k 2 ) + ( k ( k 1 ) l 1 ) × 1 = k ( k ( k 1 ) l 2 k + 2 ( k 1 ) l k 2 .
Again, by the definitions of   M 2 and A B C we obtain:
M 2 ( C ( k , l ) )   = e ϵ E k , k d ( u ) d ( v ) + e ϵ E 1 , k d ( u ) d ( v ) = ( k × k ) ( k ( k 1 ) l k + 1 k 2 3 k + 2 ) + ( k × 1 ) (   k ( k 1 ) l 1 ) = k 2 ( 2 ( k 1 ) l k ) k 2 ; A B C ( C ( k , l ) )   = e ϵ E k , k d ( u ) + d ( v ) 2 d ( u ) d ( v ) + e ϵ E 1 , k d ( u ) + d ( v ) 2 d ( u ) d ( v ) = ( k ( k 1 ) l k + 1 k 2 3 k + 2 ) ( k ) + ( k ) 2 ( k ) ( k ) + k ( k 1 ) l 1 ( 1 ) + ( k ) 2 ( 1 ) ( k ) = ( k ( k 1 ) l k + 1 k 2 3 k + 2 ) 2 ( k 1 ) k + k ( k 1 ) l 1 k 1 k .
 □
Remark 2.
The results obtained from Topo-Chemie-2020 [40] shown in Table 2 corroborate with the results obtained from Theorems 3 and 4.
Remark 3.
Gutman et al. [42] proved that the Szeged index and the Wiener index are equal for trees. Hence, in particular, W and Sz indices are equal for the pandemic and Cayley trees.

3.3. Christmas Tree Network

Definition 3 ([43]).
For s ≥ 2, a Christmas tree CT(s) is composed of an sth slim tree ST(s) = ( V 1 ,   E 1 ,   u 1 ,   l 1 ,   r 1 ) and an (s + 1)th slim tree ST(s + 1) = ( V 2 ,   E 2 ,   u 2 ,   l 2 ,   r 2 ) together with the edges ( u 1 ,   u 2 ), ( l 1 ,   r 2 ), and ( l 2 ,   r 1 ), where ST(s) = (V, E, u, l, r), with V as the node set, E as the edge set, u ∈ V as the root node, l ∈ V as the left node, and r ∈ V as the right node defined below:
1 
ST(2) is the complete graph K3 with its nodes labelled with u, l and r.
2 
The sth slim tree ST(s), with s ≥ 3, is composed of a root node u and two disjoint copies of (s − 1)th slim trees as the left subtree and right subtree, denoted by S T l ( s 1 ) = ( V 1 ,   E 1 ,   u 1 ,   l 1 ,   r 1 ) and S T r ( s 1 ) = ( V 2 ,   E 2 ,   u 2 ,   l 2 ,   r 2 ) , respectively and ST(s) = (V, E, u, l, r) is given by V = V 1 V 2 { u } ,   E = E 1 E 2 { ( u ,   u 1 ) , ( u ,   u 2 ) , ( r 1 ,   l 2 ) } ,   l = l 1 ,   r = r 2 .
For illustration, the Christmas tree CT(3) is shown in Figure 4. The number of nodes and edges of CT(s) are ( 3 × 2 s ) 2 and 9 × 2 s 6 2 , respectively.
Theorem 5.
We have the following results for CT(s):
(i) 
G A 4 ( C T ( s ) ) = i = 1 s 3 ( 2 i 1 ) 2 ( s + i 1 ) ( s + i ) 2 ( s + i 1 ) + 1 + 3 ( 2 s 1 ) ;
(ii) 
Z g 4 ( C T ( s ) ) = 12 s ( 2 s ) 6 s 9 ( 2 s ) + 9 + 12 s ( 2 s 1 ) ;
(iii) 
  Π 4 * ( C T ( s ) ) = ( 4 s ) 3 ( 2 s 1 ) i = 1 s ( 2 ( s + i ) 1 ) 3 ( 2 i 1 ) ;
(iv) 
Z g 6 ( C T ( s ) ) = 9 s + s 2 ( 12 × 2 s 1 3 ) + 3 × 2 s + 1 ( 2 s 2 3 s + 2 ) 12 ;
(v) 
Π 6 * ( C T ( s ) ) = ( 4 s 2 ) 3 ( 2 s 1 ) i = 1 s ( ( s + i 1 ) ( s + i ) ) 3 ( 2 i 1 ) ;
(vi) 
Z g 4 ( C T ( s ) ,   x ) = i = 1 s 3 ( 2 i 1 ) x 2 ( s + i ) 1 + 3 ( 2 s 1 ) x 4 s ;
(vii) 
Z g 6 ( C T ( s ) ,   x ) = i = 1 s 3 ( 2 i 1 ) x ( s + i 1 ) ( s + i ) + 3 ( 2 s 1 ) x 4 s 2 ;
(viii) 
A B C 5 Π ( C T ( s ) ) = i = 1 s ( 2 ( s + i ) 3 ( s + i 1 ) ( s + i ) ) 3 ( 2 i 2 ) × ( 4 s 2 4 s 2 ) 3 ( 2 s 2 ) .
Proof. 
The edge partitions of C T ( s ) are defined below:
  •   E s + i 1 ,   s + i = { e = u v E ( C T ( s ) ) | e c ( u ) = s + i 1   and   e c ( v ) = s + i } and n s + i 1 , s + i = | E s + i 1 , s + i | = 3 ( 2 i 1 )   , where i { 1 , , s } ;
  • E 2 s ,   2 s = { e = u v E ( C T ( s ) ) | e c ( u ) = e c ( v ) = 2 s } and n 2 s ,   2 s = | E 2 s ,   2 s | = 3 ( 2 s 1 ) .
Computations similar to those of Theorems 1 and 3 yield the following results. □
Theorem 6.
The various degree based topological indices of CT(s) are given by:
(i) 
M 1   ( C T ( s ) ) = 9 ( 3 × 2 s 2 ) ;
(ii) 
M 2   ( C T ( s ) ) = 9 ( 9 × 2 s 6 2 ) ;
(iii) 
A B C ( C T ( s ) ) = ( 3 × 2 s 2 ) .
Proof. 
As CT(s) is 3-regular, it is enough to consider the following partitions:
  • P 1 = { v |   d ( v ) = 3 }   and   | P 1 | = ( 3.2 s 2 ) :
and
  • E 3 , 3 = { e = u v E ( C T ( s ) )   | d ( u ) =   d ( v ) = 3 } and | E 3 , 3 | = ( 9 × 2 s 6 ) 2 .
A direct application of the definitions yields the required results. □
Remark 4.
The results obtained from Topo-Chemie-2020 [40] are shown in Table 3 for a comparison with the results obtained from Theorems 5 and 6.

3.4. Corona Product of Christmas Tree and a Path

Graph operations facilitate decomposition of a graph G into two or more isomorphic subgraphs. The corona product G1ʘG2 of two graphs G1 (with n1 vertices and m1 edges) and G2 (with n2 vertices and m2 edges) is defined as the graph obtained by taking a copy of G1 and n1 copies of G2, and then joining the ith vertex of G1 with edges to every vertex in the ith copy of G2. It follows from the definition of the corona product that G1ʘG2 has n1 + n1n2 vertices and m1 + n1m2 + n1n2 edges. It is easy to see that G1ʘG2 is not in general isomorphic to G2ʘG1 [44].
The corona product CT(3)ʘP3 is shown in Figure 5. The number of nodes and edges of CT(sPn are 3 × 2 s 2 + ( 3 × 2 s 2 ) n and 9 × 2 s 6 2 + ( ( 3 × 2 s ) 2 ) ( 2 n 1 ) , respectively. The eccentricity-based topological indices for CT(s) ʘ Pn, s, n ≥ 2 are given in Theorem 7.
Theorem 7.
For s , n 2 , we have the following for C T ( s ) ʘ P n :
(i) 
G A 4 ( C T ( s ) ʘ P n ) = ( 6 + 2 n ) ( s + 1 ) ( s + 2 ) 2 s + 3 + i = 2 s 3 ( 2 i 2 ) ( 2 + n ) 2 ( s + i ) ( s + i + 1 ) 2 ( s + i ) + 1
+ 6 n ( 2 s 1 ) ( 2 s + 1 ) ( 2 s + 2 ) 4 s + 3 + i = 1 s 3 ( 2 i 1 ) ( n 1 ) 2 ( s + i + 2 ) 2 ( s + i + 1 ) + 2
+ 3 ( 2 s 1 ) + n 1 ;
(ii) 
  Z g 4 ( C T ( s ) ʘ P n   ) = n ( 9 × 2 s 2 + 24 s × 2 s ) + s ( 6 × 2 s 8 n 2 ) 6 × 2 s + 5 ;
(iii) 
Π 4 * ( C T ( s ) ʘ P n   ) = ( 2 s + 3 ) 3 + n × i = 2 s ( 2 ( s + i ) + 1 ) ( 6 + 3 n ) 2 i 2 × ( 4 s + 3 ) 3 n ( 2 s 1 )
× i = 1 s ( 2 ( s + i + 2 ) ) ( 3 n 3 ) 2 i 1 × ( 4 s + 2 ) 3 ( 2 s 1 ) × ( 2 s + 4 ) n 1 ;
(iv) 
Z g 6 ( C T ( s ) ʘ P n ) = s 2 ( 6 × 2 s 4 n + 24 n × 2 s 1 ) + s ( 5 2 n 12 × 2 s + 18 n × 2 s ) + 2 s ( 15 n 3 2 ) 9 n 1 ;
(v) 
Π 6 * ( C T ( s ) ʘ P n ) = ( s 2 + 3 s + 2 ) 3 + n × i = 2 s ( ( s + i ) ( s + i + 1 ) ) ( 6 + 3 n ) 2 i 2
× ( 2 ( 2 s 2 + 3 s + 1 ) ) 3 n ( 2 s 1 ) × i = 1 s ( s + i + 2 ) ( 3 n 3 ) 2 i × ( 2 s + 1 ) 6 ( 2 s 1 )
× ( s + 2 ) 2 ( n 1 ) ;
(vi) 
Z g 4 ( C T ( s ) ʘ P n   , x ) = ( 3 + n ) x 2 s + 3 + ( 6 + 3 n ) i = 2 s 2 i 2 x 2 ( s + i ) + 1 + 3 n ( 2 s 1 ) x 4 s + 3
+ 3 ( n 1 ) i = 1 s 2 i 1 x 2 ( s + i + 2 ) + 3 ( 2 s 1 ) x 4 s + 2 + ( n 1 ) x 2 ( s + 2 ) ;
(vii) 
Z g 6 ( C T ( s ) ʘ P n   , x ) = ( 3 + n ) x s 2 + 3 s + 2 + ( 6 + 3 n ) i = 2 s 2 i 2 x ( s + i ) ( s + i + 1 )
+ 3 n ( 2 s 1 ) x ( 2 s + 1 ) ( 2 s + 2 ) + 3 ( n 1 ) i = 1 s 2 i 1 x ( s + i + 2 ) 2
+ 3 ( 2 s 1 ) x ( 2 s + 1 ) 2 + ( n 1 ) x ( s + 2 ) 2 ;
(viii) 
A B C 5 Π ( C T ( s ) ʘ P n   ) = ( 2 s + 1 ( s + 1 ) ( s + 2 ) ) 3 + n × i = 2 s ( 2 ( s + i ) 1 ( s + i + 1 ) ( s + i ) ) ( 6 + 3 n ) ( 2 i 3 )
× ( 4 s + 1 ( 2 s + 1 ) ( 2 s + 2 ) ) 3 n ( 2 s 2 ) × i = 1 s ( 2 ( s + i + 1 ) ( s + i + 2 ) 2 ) ( 3 n 3 ) ( 2 i 2 )
× ( 4 s ( 2 s + 1 ) ) 3 ( 2 s 1 ) × ( 2 ( s + 1 ) s + 2 ) n 1 .
Proof. 
In C T ( s ) ʘ P n   , we have the following partitions:
  • E s + 1 ,   s + 2 = { e = u v E ( C T ( s ) ʘ P n   ) | e c ( u ) = s + 1   and   e c ( v ) = s + 2 } and n s + 1 ,   s + 2 = | E s + 1 ,   s + 2 | = 3 + n ;
  • E s + i ,   s + i + 1 = { e = u v E ( C T ( s ) ʘ P n   ) | e c ( u ) = s + i   and   e c ( v ) = s + i + 1 } and n s + i ,   s + i + 1 = | E s + i ,   s + i + 1 | = 3 ( 2 i 2 ) ( 2 + n ) , where i { 2 , , s } ;
  • E 2 s + 1 ,   2 s + 2 = { e = u v E ( C T ( s ) ʘ P n   ) | e c ( u ) = 2 s + 1   and   e c ( v ) = 2 s + 2 } and n 2 s + 1 ,   2 s + 2 = | E 2 s + 1 ,   2 s + 2 | = 3 n ( 2 s 1 ) ;
  • E s + i + 2 ,   s + i + 2 = { e = u v E ( C T ( s ) ʘ P n   ) | e c ( u ) = e c ( v ) = s + i + 2   } and n s + i + 2 ,   s + i + 2 = | E s + i + 2 ,   s + i + 2 | = 3 ( 2 i 1 ) ( n 1 ) , where i { 1 , , s } ;
  • E 2 s + 1 ,   2 s + 1 = { e = u v E ( C T ( s ) ʘ P n   ) | e c ( u ) = e c ( v ) = 2 s + 1   } and n 2 s + 1 ,   2 s + 1 = | E 2 s + 1 ,   2 s + 1 | = 3 ( 2 s 1 ) ;
  • E s + 2 ,   s + 2 = { e = u v E ( C T ( s ) ʘ P n   ) | e c ( u ) = e c ( v ) = s + 2   } and n s + 2 ,   s + 2 = | E s + 2 ,   s + 2 | = n 1 .
Let the edge colors red, blue, yellow, black, green, and sky blue in Figure 5 represent the edge partitions E 3 + 1 ,   3 + 2 , E s + i ,   s + i + 1 , E 6 + 1 ,   6 + 2 , E s + i + 2 ,   s + i + 2 , E 6 + 1 ,   6 + 1 , and E 3 + 2 ,   3 + 2 , respectively. Hence, the edge partitions of CT(3) ʘ P3 are given as follows:
  • E 3 + 1 ,   3 + 2 = { e = u v E ( C T ( 3 ) ʘ P 3   ) | e c ( u ) = 3 + 1   a n d   e c ( v ) = 3 + 2 } and n 3 + 1 ,   3 + 2 = n 4 , 5 = | E 4 ,   5 | = 3 + 3 = 6 ;
  • E s + i ,   s + i + 1 = { e = u v E ( C T ( 3 ) ʘ P 3   ) | e c ( u ) = s + i   and   e c ( v ) = s + i + 1 } and n s + i ,   s + i + 1 = | E s + i ,   s + i + 1 | = 3 ( 2 i 2 ) ( 2 + 3 ) = 15 ( 2 i 2 )   , where i { 2 , 3 } ;
  • E 6 + 1 ,   6 + 2 = { e = u v E ( C T ( 3 ) ʘ P 3   ) | e c ( u ) = 6 + 1   and   e c ( v ) = 6 + 2 } and n 6 + 1 ,   6 + 2 = n 7 , 8 = | E 7 ,   8 | = 3 ( 3 ( 2 3 1 ) ) = 36 ;
  • E s + i + 2 ,   s + i + 2 = { e = u v E ( C T ( 3 ) ʘ P 3   ) | e c ( u ) = e c ( v ) = s + i + 2   } and n s + i + 2 ,   s + i + 2 = | E s + i + 2 ,   s + i + 2 | = 3 ( 2 i 1 ) ( n 1 ) , where i { 1 , 2 , 3 } ;
  • E 6 + 1 ,   6 + 1 = { e = u v E ( C T ( 3 ) ʘ P 3   ) | e c ( u ) = e c ( v ) = 6 + 1   } and n 7 , 7 = | E 7 , 7 | = 3 ( 2 3 1 ) = 12 ;
  • E 3 + 2 ,   3 + 2 = { e = u v E ( C T ( 3 ) ʘ P 3   ) | e c ( u ) = e c ( v ) = 3 + 2   } and n 3 + 2 , 3 + 2 = n 5 , 5 = | E 5 , 5 | = 3 1 = 2 .
Now,
Z g 4 ( C T ( s ) ʘ P n   ) = ( 3 + n ) ( 2 s + 3 ) + ( 6 + 3 n ) i = 2 s 2 i 2 ( 2 ( s + i ) + 1 ) + n ( 12 s + 9 ) ( 2 s 1 ) + 3 ( n 1 ) i = 1 s 2 i 1 × 2 ( s + i + 2 ) + 6 ( 2 s 1 ) ( 2 s + 1 ) + 2 ( n 1 ) ( s + 2 ) = n ( 9 × 2 s 2 + 24 s × 2 s ) + s ( 6 × 2 s 8 n 2 ) 6 × 2 s + 5 ;
and
Z g 6 ( C T ( s ) ʘ P n   ) = ( s 2 + 3 s ) ( 3 + n ) + 2 ( n + 3 ) + ( 6 + 3 n ) i = 2 s 2 i 2 ( s + i ) ( s + i + 1 ) + 6 n ( 2 s 1 ) ( 2 s 2 + 3 s + 1 ) + 3 ( n 1 ) i = 1 s 2 i 1 ( s + i + 2 ) 2 + 3 ( 2 s 1 ) ( 2 s + 1 ) 2 + ( n 1 ) ( s + 2 ) 2 = s 2 ( 6 × 2 s 4 n + 24 n × 2 s 1 ) + s ( 5 2 n 12 × 2 s + 18 n × 2 s ) + 2 s ( 15 n 3 2 ) 9 n 1 .
 □
All the results follow from the edge partitions defined above and the definitions of the respective indices.
Theorem 8.
The first and second Zagreb, and the atom bond connectivity indices of C T ( s ) ʘ P n are given by:
(i) 
M 1   ( C T ( s ) ʘ P n ) = ( 3 × 2 s 2 ) ( n 2 + 15 n 1 ) ;
(ii) 
M 2   ( C T ( s ) ʘ P n   ) = ( ( 3 × 2 s 2 ) ( 9 n 2 + 50 n 15 ) ) / 2 ;
(iii) 
A B C ( C T ( s ) ʘ P n ) = ( 9 × 2 s 6 ) 2 × 2 ( n + 2 ) ( n + 3 ) + 2 ( 3 × 2 s 2 ) 1 2
+ ( n 2 ) ( 3 × 2 s 2 ) ( n + 4 ) ( 3 ) ( n + 3 ) + ( 3 × 2 s 2 ) 2
+ ( n 3 ) ( 3 × 2 s 2 ) ( 2 / 3 ) .
Proof. 
Let the vertex and edge partitions of C T ( s ) ʘ P n be given by:
  • P 1 = { v |   d ( v ) = 2 }   and   | P 1 | = 2 ( 3 × 2 s 2 ) ;
  • P 2 = { v |   d ( v ) = 3 }   and   | P 2 | = ( n 2 ) ( 3 × 2 s 2 ) ;
  • P 3 = { v |   d ( v ) = n + 3 }   and   | P 3 | = ( 3 × 2 s 2 ) :
and
  • E n + 3 , n + 3 = { e = u v E ( C T ( s ) ʘ P n )   | d ( u ) = n + 3   and   d ( v ) = n + 3 } and | E n + 3 , n + 3 | = ( 9 × 2 s 6 ) / 2 ;
  • E 2 , n + 3 = { e = u v E ( C T ( s ) ʘ P n )   | d ( u ) = 2   and   d ( v ) = n + 3 } and | E 2 , n + 3 | = 2 ( 3 × 2 s 2 ) ;
  • E 3 , n + 3 = { e = u v E ( C T ( s ) ʘ P n )   | d ( u ) = 3   and   d ( v ) = n + 3 } and | E 3 , n + 3 | = ( n 2 ) ( 3 × 2 s 2 ) ;
  • E 2 , 3 = { e = u v E ( C T ( s ) ʘ P n )   | d ( u ) = 2   and   d ( v ) = 3   } and | E 2 , 3 | = 2 ( 3 × 2 s 2 ) ;
  • E 3 , 3 = { e = u v E ( C T ( s ) ʘ P n )   | d ( u ) = d ( v ) = 3   } and | E 3 , 3 | = ( n 3 ) ( 3 × 2 s 2 ) .
Now,
M 1 ( C T ( s ) ʘ P n   )   = v ϵ P 1 ( d ( v ) ) 2 + v ϵ P 2 ( d ( v ) ) 2 + v ϵ P 3 ( d ( v ) ) 2 = 2 ( 3 × 2 s 2 ) ( 2 2 ) + ( n 2 ) ( 3 × 2 s 2 ) ( 3 2 ) +   ( 3 × 2 s 2 ) ( n + 3 ) 2 = ( 3 × 2 s 2 ) ( n 2 + 15 n 1 ) ; M 2   ( C T ( s ) ʘ P n   ) = ( 9 × 2 s 6 ) 2 ( n + 3 ) 2 + 2 ( 3 × 2 s 2 ) + ( ( 2 ) ( n + 3 ) ) + ( n 2 ) ( 3 × 2 s 2 )   ( ( 3 ) ( n + 3 ) ) + 2 ( 3 × 2 s 2 ) ( ( 2 ) ( 3 ) ) + ( n 3 ) ( 3 × 2 s 2 ) ( ( 3 ) ( 3 ) ) = ( ( 3 × 2 s 2 ) ( 9 n 2 + 50 n 15 ) ) / 2 ;
Likewise, we can obtain A B C ( C T ( s ) ʘ P n ) . □
Remark 5.
The results obtained from Topo-Chemie-2020 [40] are shown in Table 4 for a comparison with the results obtained from Theorems 7 and 8.

4. Various Applications of Topological Indices for the COVID-19 Pandemic

There are two types of emerging applications of graph theory and topological indices. One of the more well studied applications of topological indices is the prediction of physicochemical properties and biological activities of chemicals and drugs through an approach called the quantitative structure activity relation (QSAR). This approach relies on the fact that the property of a molecule bears a mathematical relation to the structure or the associated connectivity of the molecule. In this context, several topological indices have been shown to statistically correlate with a number of observable properties such as dermal penetration, protonation constants, chromatographic retention indices, lipophilicity, octanol partition coefficients, and so on [11,45,46]. When topological indices are combined with other quantum chemically derived electronic parameters [46,47,48,49] such as the Highest Occupied Molecular Orbital (HOMO)- Lowest Unoccupied Molecular Orbital (LUMO) energy gaps, hardness, softness, polarizability, molecular electrostatic potentials, Natural Bond Orbital (NBO) analysis, etc., one could obtain quantitative measures of the relative stabilities, reactivities and binding potentials of a drug as a viral inhibitor which can then assist in finding potential cures for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). One could invoke statistical techniques such as the principal component analysis by including several of the indices developed in this study to discover among various topological indices which ones would offer the best correlation and predictive power compared to the observed properties [46,47,48,49].
For example, it has been shown earlier [45] that the logP values of polyacenes have a very good correlation with the Szeged index that we computed in the present study using Theorem 1, as per the following statistical correlation:
logP(Lh) = 1.0875 × 10−4Sz + 9.210,
with an r value of 0.9258. We anticipate a similar correlation for the potential drugs for the coronavirus-2, although the actual additive and multiplier constants would depend on available experimental data for these compounds.
In a recent work, Mondal et al. [47] studied degree-based and neighborhood degree sum-based topological indices in order to provide structure-activity correlations for a set of investigated antiviral drugs for COVID-19 such as remdesivir (GS-5734), chloroquine, hydroxychloroquine, theaflavin, etc., by using graph polynomial approaches. Mondal et al. [47] noted that the results obtained from such topological studies can aid in the development of potentially new drugs for the treatment of COVID-19, although it is skeptical that topological indices alone can result in the discovery of new drugs. It can, however, be stated that the developed indices have the efficacy to predict many properties and activities such as boiling point, entropy, enthalpy, acentric factor, vapor pressure etc., as these properties do not depend on guest–host interactions and exhibit direct structural dependency. However, the discovery of potentially new drugs needs more sophisticated approaches such as quantum chemical docking, molecular dynamics etc., [48,49,50] that focus on the viral protease–drug inhibitor interactions in order to enhance the efficacy of a drug as an inhibitor of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) protease. Consequently, in the area of drug discovery, topological indices play a limited role, that is, only as starting points for molecular similarity and QSAR analysis, and in the selection of various ligands from a library. Moreover, topological indices cannot serve as the end point in silico drug discovery for which sophisticated tools based on quantum mechanics are required to analyze and understand the guest–host interactions involved in potential drug discovery [48,49,50]. Figure 6 and Figure 7 show the 3D plots of log of the Szeged and fourth geometric-arithmetic indices respectively for pandemic trees. As can be seen from Figure 6, the magnitude of log of Szeged index for k = 6 at a given level is several orders larger than that of k = 2 at the same level. This clearly demonstrates the severity of the ongoing COVID-19 which has a R0 value close to 6 in a susceptible population set compared to k = 2 for the 1918 Spanish flu. As seen from Figure 7, we obtain a similar trend for the log of the fourth geometric-arithmetic indices for pandemic trees.
Plots such as the ones obtained in these figures from the computed topological indices provide significant new insights into the severity of the ongoing COVID-19 pandemic, underscoring the importance of intervention steps not only through quarantining but also through several other mitigation measures such as therapeutic interventions, vaccines and improved public awareness of hygiene, social distancing, facial shields, masks, and other measures to prevent the spread.
The second area where our developed topological indices may help us in understanding of the evolutionary dynamics of the COVID-19 pandemic in terms of genomics is in the identification of various mutations of the genes that occur during the epidemic propagation, and in the evolutionary relationship of the mutations as a function of the propagation to the presumed genome of the Wuhan bat origin. Secondly, the graph theoretical network aids in enhancing our understanding of the dynamics of the spread of the pandemic as a function of regions, sectors, states or provinces, and countries.
In relationship to the epidemiological dynamics, as several of the topological indices that we computed here are eccentricity-based indices, they measure the extent of the spread and the magnitude of the effect of the pandemic (see Figure 6, Figure 7 and Figure 8) on a given population set. Another important feature of the topological indices considered here is that the indices can provide insight into the effect of the introduction of an infected individual in a pool of uninfected individuals through vertex coloring to identify the infected person and recomputing the indices with the resulting weighted graphs. As the topological indices computed here vary as a function of level, one obtains quantitative measures in terms of these indices of the epidemic progression at different levels. Graph theoretical tools developed here can facilitate juxtaposing different scenarios by superposing one network over the other, for example, under the conditions that a given individual infects 2 others compared to an individual infecting 4 other individuals and so on. Consequently, the network analysis presented here can aid in the development of artificial intelligence tools in a computer-assisted study of the epidemiological dynamics of the pandemic, and thus in the development of mitigation strategies. There are several such applications of these indices which can become the topic of future studies.
It should be emphasized that all of the topological indices in Section 4 and entropies obtained in Section 5 and comparisons that are made, are based on the assumed R0, which is a result of the illness itself and various interventions. It is a macroscopic rather than a microscopic measure. Thus, the comparison made here is by no means indicative of the total number of infections or deaths, as they depend on multiple factors such as the discovery of vaccines, adherence to social distancing, more advanced technological tools available today for diagnosis, contract tracing, treatments, mitigation measures, and prognosis compared to 1918. Hence, there are several substantial differences between the two pandemics, namely COVID 19 and the Spanish Flu. First, they are both caused by different types of viruses with different modes of infection. We have denser populations and a far more connected world with air travel which can contribute to the spread of COVID-19 globally. Specifically, the comparisons of Covid-19 and Spanish flu provided in this section do not take directly into consideration a number of factors enumerated below which could alter the assumed R0 values dynamically and hence the computed indices:
  • The difference in the awareness among all the countries and the general population on the severity of the disease and the necessary preventive actions that are needed as a result of internet and other forms of communications.
  • The difference in the resources available for healthcare.
  • The difference in number of qualified and trained virologists, doctors, nurses, government officials, government/private institutions, and other important frontline workers.
  • The difference in information content available through research, innovations, and the data from the past pandemic, all of which can contribute to improved predictions. Researchers worldwide are working diligently to find a vaccine against the virus causing the COVID-19 pandemic. The WHO is working in collaboration with scientists, business, and global health organizations to accelerate the vaccine effectiveness and discovery [51].
  • Other social, economic, and medical changes that occurred in the last 100 years.
Therefore, comparisons of entropies computed in Section 5 and other topological indices in this section were made with presently available R0 for COVID-19, which itself could dynamically change with discoveries of vaccine, technological advances to effectively control the epidemic and the role of super spreaders. Furthermore, entropy measures from Section 5 suggest that the ongoing pandemic left to itself without such mitigating measures, poses an enormous and unique challenge globally, and the ongoing battle is nowhere near the end. Hence, significant interventions and countermeasures must be taken such as therapeutic interventions including vaccines, social distancing, vigorous quarantine measures, facial masks, improved oral and hand hygiene measures. The availability of preventive and therapeutic tools against viral infections and their related complications is an important factor for the disaster risk assessment [52].

5. Thermodynamic Entropy of Pandemic Trees

As demonstrated earlier, a complete k-ary tree T l k is the pandemic tree with an R0 value of k. The R0 epidemiological measure is exponential similar to other natural disasters such as an earthquake measured in a logarithmic Richter scale. We introduce here yet another measure of the chaos created by the COVID-19 pandemic which is based on the thermodynamic Boltzmann’s definition of entropy, a measure of the degree of disorder. At the very outset, it is pointed out that the thermodynamic entropy is different from the graph theoretical entropy discussed by Dehmer and co-workers [53,54,55] although in principle there should be a correspondence between statistical thermodynamic entropy and information theoretic entropy. The graph theoretical entropy concept invoked by these authors is based on Shannon’s information theoretic formulation as applied to graph theoretical invariants. On the other hand, the thermodynamic entropy is simply based on the celebrated Boltzmann’s definition and formula for the entropy. Consequently, in the present study, the term entropy refers to the thermodynamic entropy which is simply based on the number of distinct configurations or the number of inequivalent ways to label the vertices of a graph under the action of the automorphism group of the graph. The set RD of functions from the set D of vertices to the set R of colors would form equivalence classes under the action of the automorphism group. Among these classes, the number of classes where each vertex is colored with a unique color corresponds to the number of unique ways to label the graph, and hence, this number measures the number of distinct configurations a graph can take under the group action.
The entropy of a pandemic tree is measured directly by the number of unique ways of labeling the pandemic tree. In general, the number of inequivalent ways to label any graph with n vertices is given by:
N L a b = n ! | A u t ( G ) | ,
where Aut(G) is the automorphism group of the graph G under consideration, and |Aut(G)| is the number of elements in the automorphism group. The automorphism group of a graph with n vertices is defined as edge-preserving permutations of the vertices or it is comprised of those permutations of the vertices with permutation matrices P satisfying the condition:
A = PTAP,
where A is the adjacency matrix of the graph. The automorphism groups of several types of polyhedral graphs were studied recently by Ghorbani et al. [56].
Balasubramanian [15] previously studied the automorphisms of Cayley trees, complete k-ary trees, etc., and showed that the automorphism group of a T l k tree with l levels is given by:
A u t ( T l k ) = S k [ S k [ S k ] ] ,
where the group is a nested wreath product of the symmetric group S k recurring l times in the above product. The wreath product [15] of two groups G and H denoted by G[H] is a group defined by (H1 x H2 x…Hn)^ G’, where x is a direct product of two groups, ^ is a semi-direct product and G′ is a group isomorphic with G acting on the whole set of vertices. The group is comprised of |H|n|G| operations. The order of each Sk group composing the wreath product of a pandemic tree is k!. The order of a nested wreath product group was obtained by Balasubramanian [15] as:
| { [ S k ] } l |   =   ( k ! ) a l ,     a l =   ( k l 1 ) ( k 1 ) .
As it was shown in a previous section, the number of vertices in T l k = a l + 1 = k l + 1 1 k 1 , we arrive at the result for the number of unique ways to label a pandemic tree T l k as:
N L a b   =   ( a l + 1 ) ! ( k ! ) a l ,
where NLab is also referred to as the number of different configurations denoted by in statistical mechanics. The entropy of any system, denoted by S, is defined by the Boltzmann equation and it measures the degree of disorder. It is given by:
S = kB ln(Ω),
where kB is the Boltzmann constant.
We can obtain the entropy S(l, k) of a pandemic tree using the above equation with defined as the number of inequivalent ways of labeling the pandemic tree that we have obtained above. This combined with the Stirling approximation for the factorial of any large number, we arrive at the following expression for S(l, k) of a pandemic tree T l k   as :
S ( l , k ) = k B { a l + 1 ln ( a l + 1 ) a l + 1 + 1 2 ln ( 2 π a l + 1 ) a l ln ( k ! ) } , a l + 1 = k l + 1 1 k 1 .
In Figure 8 we have plotted S(l, k)/kB as function of both k and l, and we can see that the order of S increases as a function of k and l. This shows the severity of the ongoing pandemic which has a k value of 6 and as l goes to 6 nearly the entire susceptible population is infected. Moreover, the entropy measure of the ongoing COVID-19 pandemic is dramatically larger than that of the 1918 Spanish flu virus, underscoring the severity on the novel coronavirus 2. Furthermore, entropy is measured in natural logarithmic scale analogous to the Richter scale of an earthquake and thus the ratios of the entropies of COVID-19 to Spanish flu of 1104.6 at level 6 would suggest about 103 greater severity as measured in a natural log scale. Thus, entropy provides yet another quantitative measure for the epidemiologists to measure the degree of disorder and hence the severity of the ongoing pandemic.

6. Stochasticity in Pandemic Tree Generation

Our study so far has captured the worst case, when an infected individual spreads the infection to R0 other individuals, by studying a complete k-ary tree with k = R0 as a particular case. This may not be the case always except in worst affected locations. In such cases, the pandemic trees will just be R0-ary trees. Stochasticity in the generation of such trees leads to considering the adjacency matrices of the trees for implementing distance-based graph algorithms. Degree and eccentricity-based topological indices can be easily computed from these matrices.
As an illustration, a random tree T of level l = 6 with epicenter of degree R0 = 6 at level 0 and with every vertex of level i having R0-i children, 0 i 5 , is considered in Figure 9.
A few topological indices of this tree T of level l = 6 are listed below:
  • M1(T) = v i V ( T ) d i 2 = 1957l2 – 15,660l + 33,270;
  • M 2 (T) = v i v j E ( T ) d i d j = 1956l2 – 13,710l+25,440;
  • ABC(T) = v i v j E ( T ) d i + d j 2 d i d j
    = ( 10 6 ) l + 30 2 l 3 l 2 l + 120 2 l 5 l 2 3 l + 2 + 360 2 l 7 l 2 5 l + 6 + 120 2 l 9 l 2 7 l + 12 + 720 2 l 11 l 2 9 l + 20 ;
  • PI(T) = e = u v E ( T ) [ n u ( e ) + n v ( e ) ] = 1957l + 3,816,150;
  • Sz(T) = e = u v E ( T ) [ n u ( e ) × n v ( e ) ] = 531,706l + 15,153,306.
Evidently such generation of random k-ary trees needs to be considered through computational techniques based on algorithms such as the Monte Carlo algorithm and subsequent iterative computation of the various topological indices numerically. Such a technique would also encompass the possibility of including a super spreader in the random pandemic tree.

7. Conclusions

In the current study, we have shown that the developed topological indices play an important role in quantifying the network information contained in pandemic trees. Furthermore, the graph theoretical techniques presented here can aid in a number of predictions concerning the dynamics of the ongoing pandemic, for example, how the introduction of an infected person in a pool of uninfected group would alter the spread dynamics. This is tantamount to introducing a new vertex in the tree and associated edges. Then, one can recompute the various topological indices to see how the pandemic evolves as a function of perturbations and quarantines. Such detailed investigations specific to the pandemic can be the topic of future studies. Finally, we computed the entropies of pandemic trees and the related Cayley trees, which provided significant new insights into the severity of the ongoing COVID-19 pandemic as compared to the 1918 Spanish flu which pales into insignificance in terms of entropy and topological measures compared to the severity of COVID-19.

Author Contributions

Conceptualization, K.B., R.S.R. and G.K.N.; methodology, K.B., R.S.R., G.K.N. and A.A.S.; software, K.B.; validation, K.B., T.M.R., G.K.N. and A.A.S.; formal analysis, K.B., I.R.; investigation, R.S.R. and G.K.N.; resources, K.B.; data curation, K.B.; writing—original draft preparation, R.S.R., G.K.N., A.A.S., T.M.R. and I.R.; writing—review and editing, K.B. and R.S.R.; visualization, K.B. and R.S.R.; supervision, R.S.R., K.B. and I.R.; project administration, R.S.R.; funding acquisition, R.S.R. All authors have read and agreed to the published version of the manuscript.

Funding

The work of R. Sundara Rajan is supported by Project No. ECR/2016/1993, Science and Engineering Research Board (SERB), Department of Science and Technology (DST), Government of India.

Acknowledgments

The authors sincerely thank the reviewers for their comments and suggestions which significantly improved the quality of this paper. Further, we thank A. John Sujith, Alagappa University, Karaikudi, and K. Jagadeesh Kumar, Department of Mathematics, Hindustan Institute of Technology and Science, Chennai, India, for their fruitful suggestions.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

References

  1. Whitelaw, S.; A Mamas, M.; Topol, E.; Van Spall, H.G. Applications of digital technology in COVID-19 pandemic planning and response. Lancet Digit. Health 2020, 2, e435–e440. [Google Scholar] [CrossRef]
  2. Blackwood, J.C.; Childs, L.M. An introduction to compartmental modeling for the budding infectious disease modeler. Lett. Biomath. 2018, 5, 195–221. [Google Scholar] [CrossRef]
  3. Zhou, P.; Yang, X.-L.; Wang, X.-G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.-R.; Zhu, Y.; Li, B.; Huang, C.-L.; et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [Green Version]
  4. Lai, C.-C.; Shih, T.-P.; Ko, W.-C.; Tang, H.-J.; Hsueh, P.-R. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): The epidemic and the challenges. Int. J. Antimicrob. Agents 2020, 55, 105924. [Google Scholar] [CrossRef]
  5. Andersen, K.G.; Rambaut, A.; Lipkin, W.I.; Holmes, E.C.; Garry, R.F. The proximal origin of SARS-CoV-2. Nat. Med. 2020, 26, 450–452. [Google Scholar] [CrossRef] [Green Version]
  6. Zhou, F.; Yu, T.; Du, R.; Fan, G.; Liu, Y.; Liu, Z.; Xiang, J.; Wang, Y.; Song, B.; Gu, X.; et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: A retrospective cohort study. Lancet 2020, 395, 1054–1062. [Google Scholar] [CrossRef]
  7. Cui, J.; Li, F.; Shi, Z.-L. Origin and evolution of pathogenic coronaviruses. Nat. Rev. Genet. 2019, 17, 181–192. [Google Scholar] [CrossRef] [Green Version]
  8. Wu, J.T.; Leung, K.; Leung, G.M. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: A modelling study. Lancet 2020, 395, 689–697. [Google Scholar] [CrossRef] [Green Version]
  9. Forster, P.; Forster, L.; Renfrew, C.; Forster, M. Phylogenetic network analysis of SARS-CoV-2 genomes. Proc. Natl. Acad. Sci. USA 2020, 117, 9241–9243. [Google Scholar] [CrossRef] [Green Version]
  10. Balasubramanian, K.; Khokhani, K.; Basak, S.C. Complex Graph Matrix Representations and Characterizations of Proteomic Maps and Chemically Induced Changes to Proteomes. J. Proteome Res. 2006, 5, 1133–1142. [Google Scholar] [CrossRef]
  11. Basak, S.C.; Grunwald, G.D.; Gute, B.D.; Balasubramanian, K.; Opitz, D. Use of statistical and neural net approaches in predicting toxicity of chemicals. J. Chem. Inf. Comput. Sci. 2000, 40, 885–890. [Google Scholar] [CrossRef]
  12. Matsen, F.A. Phylogenetics and the Human Microbiome. Syst. Biol. 2015, 64, e26–e41. [Google Scholar] [CrossRef]
  13. Delucchi, E. Nested set complexes of Dowling lattices and complexes of Dowling trees. J. Algebraic Comb. 2007, 26, 477–494. [Google Scholar] [CrossRef] [Green Version]
  14. Balasubramanian, K. Tree pruning and lattice statistics on Bethe lattices. J. Math. Chem. 1988, 2, 69–82. [Google Scholar] [CrossRef]
  15. Balasubramanian, K. Nested wreath groups and their applications to phylogeny in biology and Cayley trees in chemistry and physics. J. Math. Chem. 2017, 55, 195–222. [Google Scholar] [CrossRef]
  16. Balasubramanian, K. Spectra of chemical trees. Int. J. Quantum Chem. 1982, 21, 581–590. [Google Scholar] [CrossRef]
  17. Balasubramanian, K. Symmetry groups of chemical graphs. Int. J. Quantum Chem. 1982, 21, 411–418. [Google Scholar] [CrossRef]
  18. Sellers, P.H. On the Theory and Computation of Evolutionary Distances. SIAM J. Appl. Math. 1974, 26, 787–793. [Google Scholar] [CrossRef]
  19. Sellers, P.H. An algorithm for the distance between two finite sequences. J. Comb. Theory Ser. A 1974, 16, 253–258. [Google Scholar] [CrossRef] [Green Version]
  20. Fischer, M.; Klaere, S.; Nguyen, M.A.T.; Von Haeseler, A. On the group theoretical background of assigning stepwise mutations onto phylogenies. Algorithms Mol. Biol. 2012, 7, 36. [Google Scholar] [CrossRef] [Green Version]
  21. Yun, U.; Lee, G.; Kim, C.-H. The Smallest Valid Extension-Based Efficient, Rare Graph Pattern Mining, Considering Length-Decreasing Support Constraints and Symmetry Characteristics of Graphs. Symmetry 2016, 8, 32. [Google Scholar] [CrossRef] [Green Version]
  22. Xue, L.; Jing, S.; Miller, J.C.; Sun, W.; Li, H.; Estrada-Franco, J.G.; Hyman, J.M.; Zhu, H. A data-driven network model for the emerging COVID-19 epidemics in Wuhan, Toronto and Italy. Math. Biosci. 2020, 326, 108391. [Google Scholar] [CrossRef]
  23. Basavanagoud, B.; Desai, V.R.; Patil, S. (β,α)−Connectivity Index of Graphs. Appl. Math. Nonlinear Sci. 2017, 2, 21–30. [Google Scholar] [CrossRef] [Green Version]
  24. Gao, W.; Wang, W. The Vertex Version of Weighted Wiener Number for Bicyclic Molecular Structures. Comput. Math. Methods Med. 2015, 2015, 1–10. [Google Scholar] [CrossRef] [Green Version]
  25. Gao, W.; Farahani, M.R.; Shi, L. The forgotten topological index of some drug structures. Acta. Medica. Mediterr. 2016, 32, 579–585. [Google Scholar]
  26. Gao, W.; Wang, W.; Farahani, M.R. Topological Indices Study of Molecular Structure in Anticancer Drugs. J. Chem. 2016, 2016, 1–8. [Google Scholar] [CrossRef] [Green Version]
  27. Martínez-Pérez, Á.; Rodriguez, J.M. New Bounds for Topological Indices on Trees through Generalized Methods. Symmetry 2020, 12, 1097. [Google Scholar] [CrossRef]
  28. Atanasov, R.; Furtula, B.; Škrekovski, R. Trees with Minimum Weighted Szeged Index Are of a Large Diameter. Symmetry 2020, 12, 793. [Google Scholar] [CrossRef]
  29. Knor, M.; Imran, M.; Jamil, M.K.; Škrekovski, R. Remarks on Distance Based Topological Indices for ℓ-Apex Trees. Symmetry 2020, 12, 802. [Google Scholar] [CrossRef]
  30. Liu, J.-B.; Shaker, H.; Nadeem, I.; Farahani, M.R. Eccentric Connectivity Index of t-Polyacenic Nanotubes. Adv. Mater. Sci. Eng. 2019, 2019, 1–9. [Google Scholar] [CrossRef] [Green Version]
  31. Ghorbani, M.; Dehmer, M.; Emmert-Streib, F. Properties of Entropy-Based Topological Measures of Fullerenes. Mathematics 2020, 8, 740. [Google Scholar] [CrossRef]
  32. Ghorbani, M.; Khaki, A. A note on the fourth version of geometric-arithmetic index. Optoelectron. Adv. Mat. 2010, 4, 2212–2215. [Google Scholar]
  33. Gao, W.; Wu, H.; Siddiqui, M.K.; Baig, A.Q. Study of biological networks using graph theory. Saudi J. Biol. Sci. 2018, 25, 1212–1219. [Google Scholar] [CrossRef]
  34. Farahani, M.R.; Kanna, R.M.R. Fourth zagreb index of circumcoronene series of benzenoid. Leonardo Electron. J. Pract. Technol. 2015, 27, 155–161. [Google Scholar]
  35. 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]
  36. Das, K.C.; Gutman, I.; Furtula, B. On atom-bond connectivity index. Chem. Phys. Lett. 2011, 511, 452–454. [Google Scholar] [CrossRef] [Green Version]
  37. Khadikar, P.V.; Karmarkar, S.; Agrawal, V.K. A Novel PI Index and Its Applications to QSPR/QSAR Studies. J. Chem. Inf. Comput. Sci. 2001, 41, 934–949. [Google Scholar] [CrossRef]
  38. Gutman, I. A formula for the Wiener number of trees and its extension to graphs containing cycles. Graph Theory Notes NY 1994, 27, 9–15. [Google Scholar]
  39. Song, C.; Hao, R.-X. Antimagic orientations for the complete k-ary trees. J. Comb. Optim. 2019, 38, 1077–1085. [Google Scholar] [CrossRef]
  40. Balasubramanian, K. Topo-Chemie-2020 is a package of codes that computes numerous degree-based, distance-based, eccentricity-based, neighbour-based topological indices, characteristic polynomials, matching polynomials, distance polynomials, distance degree vector sequences, walks and self-returning walks, and automorphisms of graphs.
  41. Weisstein, E.W. Cayley Tree, MathWorld-A Wolfram Web Resource. Available online: https://mathworld.wolfram.com/CayleyTree.html (accessed on 17 September 2020).
  42. Gutman, I.; Dobrynin, A.A. The Szeged index—A success story. Graph Theory Notes NY 1998, 34, 37–44. [Google Scholar]
  43. Hung, C.-N.; Hsu, L.-H.; Sung, T.-Y. Christmas tree: A versatile 1-fault-tolerant design for token rings. Inf. Process. Lett. 1999, 72, 55–63. [Google Scholar] [CrossRef]
  44. Nada, S.; Elrokh, A.; Elsakhawi, E.; Sabra, D. The corona between cycles and paths. J. Egypt. Math. Soc. 2017, 25, 111–118. [Google Scholar] [CrossRef]
  45. Khadikar, P.; Karmarkar, S.; Agrawal, V.; Singh, J.; Khadikar, P.V.; Lukovits, I.; Diudea, M.V. Szeged Index—Applications for Drug Modeling. Lett. Drug Des. Discov. 2005, 2, 606–624. [Google Scholar] [CrossRef] [Green Version]
  46. Basak, S.C.; Mills, D.; Mumtaz, M.M.; Balasubramanian, K. Use of topological indices in predicting aryl hydrocarbon receptor binding potency of dibenzofurans: A hierarchical QSAR approach. Indian J. Chem. 2003, 42, 1385–1391. [Google Scholar]
  47. Mondal, S.; De, N.; Pal, A. Topological Indices of Some Chemical Structures Applied for the Treatment of COVID-19 Patients. Polycycl. Aromat. Compd. 2020, 1–15. [Google Scholar] [CrossRef]
  48. Balasubramanian, K.; Gupta, S.P. Quantum Molecular Dynamics, Topological, Group Theoretical and Graph Theoretical Studies of Protein-Protein Interactions. Curr. Top. Med. Chem. 2019, 19, 426–443. [Google Scholar] [CrossRef]
  49. Balasubramanian, K. Mathematical and Computational Techniques for Drug Discovery: Promises and Developments. Curr. Top. Med. Chem. 2019, 18, 2774–2799. [Google Scholar] [CrossRef]
  50. Patil, V.M.; Narkhede, R.R.; Masand, N.; Rameshwar, S.; Cheke, R.S.; Balasubramanian, K. Molecular insights into Resveratrol and its analogs as SARS-CoV-2 (COVID-19) protease inhibitors. Coronaviruses 2020, in press. [Google Scholar]
  51. Available online: https://www.who.int/emergencies/diseases/novel-coronavirus-2019/covid-19-vaccines (accessed on 25 November 2020).
  52. Javelle, E.; Raoult, D. COVID-19 pandemic more than a century after the Spanish flu. Lancet Infect. Dis. 2020. [Google Scholar] [CrossRef]
  53. Dehmer, M.; Mowshowitz, A. A history of graph entropy measures. Inf. Sci. 2011, 181, 57–78. [Google Scholar] [CrossRef]
  54. Mowshowitz, A.; Dehmer, M. Entropy and the Complexity of Graphs Revisited. Entropy 2012, 14, 559–570. [Google Scholar] [CrossRef]
  55. Mowshowitz, A. Entropy and the complexity of graphs: I. An index of the relative complexity of a graph. Bull. Math. Biol. 1968, 30, 175–204. [Google Scholar] [CrossRef]
  56. Ghorbani, M.; Dehmer, M.; Rahmani, S.; Rajabi-Parsa, M. A Survey on Symmetry Group of Polyhedral Graphs. Symmetry 2020, 12, 370. [Google Scholar] [CrossRef] [Green Version]
Figure 1. A pandemic tree with an epidemiological R0 value of 4.
Figure 1. A pandemic tree with an epidemiological R0 value of 4.
Symmetry 12 01992 g001
Figure 2. (a): A 5-level pandemic tree T 5 3 , with reproduction number 3. (b): A 6-level pandemic tree T 6 6 , with a reproduction number 6.
Figure 2. (a): A 5-level pandemic tree T 5 3 , with reproduction number 3. (b): A 6-level pandemic tree T 6 6 , with a reproduction number 6.
Symmetry 12 01992 g002
Figure 3. (a): A 3-Cayley tree, C(3,6), of degree 3 with 6 levels. (b): A 6-Cayley tree, C(6,6), of degree 6 with 6 levels.
Figure 3. (a): A 3-Cayley tree, C(3,6), of degree 3 with 6 levels. (b): A 6-Cayley tree, C(6,6), of degree 6 with 6 levels.
Symmetry 12 01992 g003
Figure 4. A Christmas tree CT(3) with 3 levels.
Figure 4. A Christmas tree CT(3) with 3 levels.
Symmetry 12 01992 g004
Figure 5. The corona product CT(3)ʘP3 of a Christmas tree CT(3) and a path P3.
Figure 5. The corona product CT(3)ʘP3 of a Christmas tree CT(3) and a path P3.
Symmetry 12 01992 g005
Figure 6. log(Sz) values obtained for the Szeged index of pandemic trees T l k demonstrate the underlying severity of COVID-19 as compared to the 1918 Spanish flu.
Figure 6. log(Sz) values obtained for the Szeged index of pandemic trees T l k demonstrate the underlying severity of COVID-19 as compared to the 1918 Spanish flu.
Symmetry 12 01992 g006
Figure 7. log(GA4) values obtained for the fourth geometric-arithmetic index of pandemic trees T l k demonstrate the underlying severity of COVID-19 as compared to the 1918 Spanish flu.
Figure 7. log(GA4) values obtained for the fourth geometric-arithmetic index of pandemic trees T l k demonstrate the underlying severity of COVID-19 as compared to the 1918 Spanish flu.
Symmetry 12 01992 g007
Figure 8. Entropy values measured in the units of kB or S(l, k) obtained for pandemic tree T l k accentuates the underlying severity of COVID-19 compared to the 1918 Spanish flu.
Figure 8. Entropy values measured in the units of kB or S(l, k) obtained for pandemic tree T l k accentuates the underlying severity of COVID-19 compared to the 1918 Spanish flu.
Symmetry 12 01992 g008
Figure 9. A random tree T with epicenter of degree 6.
Figure 9. A random tree T with epicenter of degree 6.
Symmetry 12 01992 g009
Table 1. Results obtained for the pandemic tree T l k with the computer code, compared with the results from the expressions in Theorems 1 and 2.
Table 1. Results obtained for the pandemic tree T l k with the computer code, compared with the results from the expressions in Theorems 1 and 2.
IndexDimension l, kFrom Expressions 1 and 2Topo-Chemie-2020 [40]
GA4 ( T l k ) l = 3, k = 338.801638.801700143357856
l = 3, k = 483.594783.59488962413361
l = 4, k = 3119.6838119.6839561076481
l = 4, k = 4339.1496339.1498266447425
Zg4 ( T l k ) l = 3, k = 3399399
l = 3, k = 4876876
l = 4, k = 316921692
l = 4, k = 448844884
Π 4 * ( T l k ) l = 3, k = 31.74212456 × 10391.742124563637115 × 1039
l = 3, k = 41.98337192 × 10851.9833719240465464 × 1085
l = 4, k = 33.75944893 × 101373.759448938495563 × 10137
l = 4, k = 4Infinity+Inf
Zg6 ( T l k ) l = 3, k = 310261026
l = 3, k = 422882288
l = 4, k = 360006000
l = 4, k = 417,58417,584
Π 6 * ( T l k ) l = 3, k = 36.74664061 × 10546.746640616477462 × 1054
l = 3, k = 44.66622007 × 101194.666220065428966 × 10119
l = 4, k = 34.24758724 × 102024.247587242244764 × 10202
l = 4, k = 4Infinity+Inf
ABC5Π ( T l k ) l = 3, k = 32.0850970590 × 10−102.0850970590395853 × 10−10
l = 3, k = 47.2443933360 × 10−227.2443933361227945 × 10−22
l = 4, k = 32.6373891592 × 10−372.6373891591075294 × 10−37
l = 4, k = 41.26583295 × 10−1051.2658329528339406 × 10−105
P I ( T l k ) l = 3, k = 315601560
l = 3, k = 471407140
l = 4, k = 314,52014,520
l = 4, k = 4115,940115,940
S z ( T l k ) l = 3, k = 334023402
l = 3, k = 417,15217,152
l = 4, k = 344,71244,712
l = 4, k = 4389,120389,120
M 1 ( T l k ) l = 3, k = 3228228
l = 3, k = 4580580
l = 4, k = 3714714
l = 4, k = 423722372
M 2 ( T l k ) l = 3, k = 3288288
l = 3, k = 4800800
l = 4, k = 3936936
l = 4, k = 433603360
A B C ( T l k ) l = 3, k = 330.830830.83052949654569
l = 3, k = 468.660768.66073893642222
l = 4, k = 394.1394.12995706469174
l = 4, k = 4276.5956276.59462680515827
Table 2. Results obtained for the Cayley Trees C ( k , l ) with the computer code, compared with the results from the expressions in Theorems 3 and 4.
Table 2. Results obtained for the Cayley Trees C ( k , l ) with the computer code, compared with the results from the expressions in Theorems 3 and 4.
IndexDimension k, lExpressions 3 and 4Topo-Chemie-2020 [40]
GA4 ( C ( k , l ) ) k = 3, l = 320.882388241420.88238824139906
k = 3, l = 444.867630784444.86763078442802
k = 4, l = 351.735600191151.73560019114378
k = 4, l = 4159.578608144159.57860814353083
Zg4 ( C ( k , l ) ) k = 3, l = 3207207
k = 3, l = 4609609
k = 4, l = 3532532
k = 4, l = 422562256
Π 4 * ( C ( k , l ) ) k = 3, l = 35.72086104 × 10205.7208610362887386 × 1020
k = 3, l = 45.06517068 × 10505.065170680143067 × 1050
k = 4, l = 32.09623003 × 10522.0962300336315637 × 1052
k = 4, l = 42.71330692 × 101832.7133069091840815 × 10183
Zg6 ( C ( k , l ) ) k = 3, l = 3516516
k = 3, l = 420882088
k = 4, l = 313681368
k = 4, l = 480008000
Π 6 * ( C ( k , l ) ) k = 3, l = 35.87731231 × 10285.877312307199999 × 1028
k = 3, l = 41.5897237 × 10741.589723697730939 × 1074
k = 4, l = 31.27482362 × 10731.2748236216396078 × 1073
k = 4, l = 41.482028 × 102701.4820280048671849 × 10270
ABC5Π ( C ( k , l ) ) k = 3, l = 30.000008406490.000008406491899037542
k = 3, l = 44.4189082 × 10−144.419672639996076 × 10−14
k = 4, l = 31.23642668 × 10−131.2364266794260217 × 10−13
k = 4, l = 41.66132298 × 10−491.6913453312569212 × 10−49
P I ( C ( k , l ) ) k = 3, l = 3462462
k = 3, l = 420702070
k = 4, l = 327562756
k = 4, l = 425,76025,760
S z ( C ( k , l ) ) k = 3, l = 3909909
k = 3, l = 456615661
k = 4, l = 363046304
k = 4, l = 482,33682,336
M 1 ( C ( k , l ) ) k = 3, l = 3102102
k = 3, l = 4222222
k = 4, l = 3308308
k = 4, l = 4956956
M 2 ( C ( k , l ) ) k = 3, l = 3117117
k = 3, l = 4261261
k = 4, l = 3400400
k = 4, l = 412641264
A B C ( C ( k , l ) ) k = 3, l = 315.797958971415.797958971132715
k = 3, l = 433.595917942233.59591794226542
k = 4, l = 340.974873507440.974873507372514
k = 4, l = 4125.374110265125.37411026490041
Table 3. Results obtained for the Christmas trees, CT(s) with the computer code, compared with the results from the expressions in Theorems 5 and 6.
Table 3. Results obtained for the Christmas trees, CT(s) with the computer code, compared with the results from the expressions in Theorems 5 and 6.
IndexDimensionsFrom Expressions 5 and 6 Topo-Chemie-2020 [40]
GA4 ( C T ( s ) ) s = 332.882332.882388241399056
s = 468.867168.86763078442803
s = 5140.8350140.83501683597572
Zg4 ( C T ( s ) ) s = 3351351
s = 4993993
s = 525712571
  Π 4 * ( C T ( s ) ) s = 35.10073079 × 10335.10077716500643 × 1033
s = 44.01304165 × 10794.0130416580886176 × 1079
s = 51.91126114 × 101771.9112611351919165 × 10177
Zg6 ( C T ( s ) ) s = 3948948
s = 436243624
s = 511,86211,862
  Π 6 * ( C T ( s ) ) s = 32.7848947 × 10472.7848946955924445 × 1047
s = 43.54520226 × 101173.5452023119163584 × 10117
s = 53.37859549 × 102693.3785954904846475 × 10269
ABC5Π ( C T ( s ) ) s = 33.8618890813 × 10−93.861889064427591 × 10−9
s = 45.3059950832 × 10−225.305995082939164 × 10−22
s = 51.7589082206 × 10−501.7589082206095107 × 10−50
M 1 ( C T ( s ) ) s = 3198198
s = 4414414
s = 5846846
M 2 ( C T ( s ) ) s = 3297297
s = 4621621
s = 512691269
A B C ( C T ( s ) ) s = 32222.000000000000004
s = 44645.99999999999997
s = 59494.00000000000011
Table 4. Results obtained for the corona product CT(sPn of a Christmas tree CT(s) and a path Pn with the computer code, compared with the results from the expressions in Theorems 7 and 8.
Table 4. Results obtained for the corona product CT(sPn of a Christmas tree CT(s) and a path Pn with the computer code, compared with the results from the expressions in Theorems 7 and 8.
IndexDimension s, nFrom Expressions 7 and 8 Topo-Chemie-2020 [40]
GA4(CT(sPn)s = 3, n = 3142.7317142.73175770071703
s = 4, n = 4390.5798390.5798602154981
Zg4(CT(sPn)s = 3, n = 319611961
s = 4, n = 468696869
Π 4 * (CT(sPn)s = 3, n = 39.65986697 × 101619.659866968127555 × 10161
s = 4, n = 4Infinity+Inf
Zg6(CT(sPn)s = 3, n = 368246824
s = 4, n = 430,56730,567
Π 6 * (CT(sPn)s = 3, n = 34.3849717 × 102374.384971723725398 × 10237
s = 4, n = 4Infinity+Inf
ABC5Π(CT(sPn)s = 3,n = 31.4118921672 × 10−431.411892167070263 × 10−43
s = 4, n = 46.08702627 × 10−1366.087026262790365 × 10−136
M 1 ( C T ( s ) ʘ P n ) s = 3, n = 311661166
s = 4, n = 434503450
M 2 ( C T ( s ) ʘ P n ) s = 3, n = 323762376
s = 4, n = 475677567
A B C ( C T ( s ) ʘ P n ) s = 3, n = 393.3373342935293.33733429351358
s = 4, n = 4251.70409168311251.70409168311375
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nandini, G.K.; Rajan, R.S.; Shantrinal, A.A.; Rajalaxmi, T.M.; Rajasingh, I.; Balasubramanian, K. Topological and Thermodynamic Entropy Measures for COVID-19 Pandemic through Graph Theory. Symmetry 2020, 12, 1992. https://doi.org/10.3390/sym12121992

AMA Style

Nandini GK, Rajan RS, Shantrinal AA, Rajalaxmi TM, Rajasingh I, Balasubramanian K. Topological and Thermodynamic Entropy Measures for COVID-19 Pandemic through Graph Theory. Symmetry. 2020; 12(12):1992. https://doi.org/10.3390/sym12121992

Chicago/Turabian Style

Nandini, G. Kirithiga, R. Sundara Rajan, A. Arul Shantrinal, T. M. Rajalaxmi, Indra Rajasingh, and Krishnan Balasubramanian. 2020. "Topological and Thermodynamic Entropy Measures for COVID-19 Pandemic through Graph Theory" Symmetry 12, no. 12: 1992. https://doi.org/10.3390/sym12121992

APA Style

Nandini, G. K., Rajan, R. S., Shantrinal, A. A., Rajalaxmi, T. M., Rajasingh, I., & Balasubramanian, K. (2020). Topological and Thermodynamic Entropy Measures for COVID-19 Pandemic through Graph Theory. Symmetry, 12(12), 1992. https://doi.org/10.3390/sym12121992

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