Symmetry-Adapted Domination Indices: The Enhanced Domination Sigma Index and Their Applications in QSPR Studies of Octane and Its Isomers

1. Introduction

In recent years, mathematical modeling using graphs with parameterized theories or invariants has become increasingly popular across various disciplines in the physical sciences. Disciplines such as computer science, physics, and chemistry have utilized these models to solve complex problems. A critical component of these models is the topological index, which provides valuable information about a graph under graph isomorphism conditions. Topological indices are invariant numbers that convey information about the size, symmetry, branching degree, and cyclicity of a graph. In particular, chemical graph theory has emerged as a leading field that combines graph theory with chemistry to study molecular structures. Among molecular descriptors, topological indices have emerged as the most important ones, providing crucial information about graphs that represent chemical compounds. By utilizing these topological indices, researchers can gain valuable insights into the properties and characteristics of various chemical compounds, allowing them to develop more accurate models and predictions. The integration of graph theory and chemistry has opened up new avenues for research in the physical sciences, and it is an exciting area that promises to yield even more insights and breakthroughs in the years to come. (for more details, [1,2,3,4,5,6]). The study of graph theory and its application to chemistry has become a rapidly growing field of research in recent years. In particular, the use of topological graph indexes has proven to be a valuable tool for understanding the structure-property relationships of chemical compounds. The Chemical Data Bases, which contain over 3000 topological graph indexes, demonstrate the vast amount of information that can be extracted from graphs representing chemical structures. Furthermore, the extensive study of dominating problems within graph theory has led to a wealth of knowledge and research in the field. As evidenced by the 1222 papers listed in the 1998 book, dominating problems have been extensively studied and analyzed, providing a solid foundation for the development of new and innovative approaches to solving problems related to chemical structures. With the continued growth and development of graph theory and its applications in chemistry, we can expect to see even more exciting discoveries and advancements in the future [7,8,9,10].

Highlighted here are some key research studies that demonstrate the significant role of Graph Theory in solving real-world problems. In [11] the authors present a novel approach to solve fractional boundary value problems on the methylpropane graph. The technique involves utilizing fixed point theorems to demonstrate the existence of solutions, as well as studying stability. To demonstrate the efficacy of the method, an example is provided. Turab et al. [12] introduce the isobutane graph and investigates the existence of solutions to fractional boundary value problems using fixed point theory. The study includes two examples to support the findings. While [13], which is Part I of a series of papers, demonstrates that graph theory can be used to provide solutions for problems in distance geometry, potential theory, and theory of metric spaces.

Domination indices are mathematical parameters used in graph theory to study the structural properties of graphs. The concept of domination indices was first introduced in the 1960s, and since then, several types of domination indices have been proposed and investigated [14]. One of the most commonly used domination indices is the total domination number ($\gamma \left(t\right)$), [15] which is defined as the minimum number of vertices in a dominating set of a graph. The total domination number has been extensively studied in the literature, with several properties and bounds known. For example, the total domination number of a tree is at most $\lceil n/2\rceil$, where n is the number of vertices in the tree. Another domination index that has received considerable attention is the connected domination number ($\gamma \left(c\right)$), [16] which is defined as the minimum number of vertices in a dominating set that induces a connected subgraph. The connected domination number has been studied in various contexts, including its relationship with other graph parameters, such as the vertex cover number and the independent domination number. The domination polynomial [17] is another important parameter that has been extensively studied in the literature. The domination polynomial of a graph is a polynomial that encodes the number of dominating sets of each size in the graph. Several properties of the domination polynomial have been established, including its relationship with other graph polynomials, such as the chromatic polynomial and the independence polynomial. Other domination indices that have been investigated in the literature include the independent domination number, the bondage number, the game domination number, and the strong domination number [18]. These parameters have been studied in various contexts, including their relationship with other graph parameters, their computational complexity, and their applicability in real-world problems.

The Sigma index is a topological index introduced by Ivan Gutman in 1978 [19] and [20]. It is defined as the square of the difference between the degrees for all pairs of adjacent vertices in a graph. More specifically, let $\zeta =(V,E)$ be a graph with vertex set V and edge set E, and let $d_u$ be the degree of the vertex u in $\zeta$. Then, the Sigma index is defined as:

$$\sigma \left(\zeta \right)=\sum _{uv\in E\left(\zeta \right)}{(d_u-d_v)}^2$$

(1)

The Sigma index has been extensively studied in the literature due to its applicability in various areas of chemistry, including the prediction of physicochemical properties of molecules [21]. Several properties and bounds of the Sigma index have been established, including its relationship with other topological indices such as the Wiener index [22]. Various modifications and generalizations of the Sigma index have also been proposed, such as the modified Sigma index, which takes into account the number of vertices at a given distance from a central vertex [23]. The degree-based Sigma index has also been introduced, which weights the contributions of each pair of vertices by their respective degrees [24]. Overall, the Sigma index and its variations have been shown to be valuable tools in the study of molecular structure and properties, as well as in the analysis of various types of networks. Detailed discussions of Sigma index applications can be found in [25,26,27,28,29]

Hanan Ahmed introduced the concept of domination topological indices in 2021 [30]. The domination topological index (DTI) is defined as the sum of the distances between each vertex and its nearest dominating vertex in a graph. The domination number of a graph is the minimum number of vertices required to dominate the graph, and it is a special case of the DTI. The DTI has been shown to be a useful tool for predicting various physicochemical properties of organic compounds. Several variations and extensions have been proposed in the literature. For example, Hosamani et al. proposed a modified version of the DTI called the modified domination topological index (MDTI) in their paper [31]. The MDTI is defined as the sum of the distances between each vertex and its nearest dominating vertex, where the dominating set is restricted to a subset of vertices with a fixed size. Another variation of the DTI is the connected domination topological index (CDTI), which was introduced by Merrick et al. [32]. The CDTI is defined as the sum of the distances between each vertex and its nearest dominating vertex in a connected subgraph of the graph. The DTI and its variations have been applied in various fields, including chemistry, biology, and computer science. For example, Yousefi et al. [33] used the DTI to develop quantitative structure-property relationship (QSPR) models for predicting the boiling points of organic compounds. In another study, Yang et al. used the CDTI to predict the toxicity of polycyclic aromatic hydrocarbons (PAHs) [34]. Overall, the DTI and its variations have proven to be useful tools for studying the structural properties of graphs and predicting various physicochemical properties of organic compounds. Extensive explanations regarding the applications of topological domination indices can be found in [35,36,37,38,39,40,41].

The combination of the Gutman topological Sigma index and Hanan et al. domination topological index into a single concept, the domination Sigma index, holds great promise in enhancing our understanding of the structural properties of graphs. This new index combines the mathematical principles of domination topological indices with the concept of Sigma index, creating a distinct index that offers valuable information about the molecular properties. The topological Sigma index provides important information about molecular size, symmetry, branching degree, and cyclicity, while the domination topological index provides insights into dominating sets, subsets of the vertex set such that every vertex outside the set is adjacent to at least one vertex inside the set. The combination of these two indices allows for a more comprehensive analysis of the properties of graphs and has the potential to improve our ability to predict the physicochemical properties of organic compounds. This article presents a novel approach to combining these two indices, providing a framework for further research in the field of mathematical chemistry. Our study presents a novel approach that combines the core principles of the Sigma index with the domination degrees of vertices in a graph to formulate a new composite index. This index is termed as the domination Sigma index and is defined as follows:

$$D\sigma \left(\zeta \right)=\sum _{uv\in E\left(\zeta \right)}{(d_{p_u}-d_{p_v})}^2$$

(2)

In [30] Hanan Ahmed et al. have introduced new degree-based topological indices called domination topological indices, which are based on the domination degree set defined as: For each vertex $v\in V\left(\zeta \right)$, the domination degree of the vertex v is denoted by $d_{p_v}$ and defined as the number of minimal dominating set of $\zeta$ which contains v. The first and second domination Zagreb indices and modified first Zagreb domination indices are defined respectively as:

$$\begin{array}{cc}\hfill D{M}_1\left(\zeta \right)& =\sum _{v\in V\left(\zeta \right)}{d}_{p_v}^2,\hfill \end{array}$$

(3)

$$\begin{array}{cc}\hfill D{M}_2\left(\zeta \right)& =\sum _{uv\in E\left(\zeta \right)}\left[d_{p_u}·d_{p_v}\right],\hfill \end{array}$$

(4)

$$\begin{array}{cc}\hfill D{M}_1^{*}\left(\zeta \right)& =\sum _{uv\in E\left(\zeta \right)}[d_{p_u}+d_{p_v}].\hfill \end{array}$$

(5)

The forgotten domination, hyper domination, and modified forgotten domination indices of graphs are defined respectively as:

$$\begin{array}{cc}\hfill DF\left(\zeta \right)& =\sum _{v\in V\left(\zeta \right)}{d}_{p_v}^3,\hfill \end{array}$$

(6)

$$\begin{array}{cc}\hfill DH\left(\zeta \right)& =\sum _{uv\in E\left(\zeta \right)}{[d_{p_u}+d_{p_v}]}^2,\hfill \end{array}$$

(7)

$$\begin{array}{cc}\hfill D{F}^{*}\left(\zeta \right)& =\sum _{uv\in E\left(\zeta \right)}\left[d_{p_u}^2+d_{p_v}^2\right].\hfill \end{array}$$

(8)

In this study, we aim to investigate the potential of the domination Sigma index and other domination topological indices in predicting the properties of Octanes and its isomers through quantitative structure-property relationship (QSPR) analysis. To achieve this goal, we calculate the domination Sigma index for various families of graphs, including book graphs, compositions of graphs, and special graph classes and find some sharp bounds. By analyzing the domination indices and the new topological index, we hope to gain insights into the properties and structures of these molecules, and apply this knowledge in the design of new chemical compounds with desired properties. However, it is important to acknowledge that the domination Sigma index may not always produce satisfactory results in QSPR analysis. As newly proposed topological indices may not always capture the key features of the chemical structure under investigation, it is not uncommon for a new index to fail to meet expectations. The lack of satisfactory results obtained from the newly introduced topological index may be attributed to its limitations and the nature of the specific property being studied. Nonetheless, it is worth noting that the failure of a new topological index to provide satisfactory results in QSPR analysis does not necessarily imply that the index is of no value. In fact, each new index provides valuable insights into the complex relationship between molecular structure and physical properties. However, the other domination topological indices used in this study have demonstrated good correlation coefficients with the properties of Octanes and its isomers, indicating their effectiveness in predicting such properties. Thus, further research and analysis may be needed to explore the potential of the domination Sigma index in predicting other properties or to modify it to better suit the studied property. Overall, this study contributes to the ongoing process of developing new topological indices and enhancing our understanding of the relationship between molecular structure and physical properties.

2. Preliminaries

In this section, we introduce key concepts, definitions, and assumptions that underpin our research, as well as outline our research questions, methodology, and the structure of our paper.

Let $\zeta$ be connected simple graph with $V\left(\zeta \right)$ a set of vertices and $E\left(\zeta \right)$ a set of edges. A set $D\subseteq V$ is said to bea dominating set of a graph $\zeta$, if for any vertex $v\in V-D$, there is a vertex $u\in D$ such that u and v are adjacent. A dominating set $D=\{v_1,v_2,...,v_r\}$ is minimal if $D-v_i$ is not a dominating set. A dominating set of $\zeta$ of minimum cardinality is said to be a minimum dominating set. Define $▵\left(\zeta \right)=$$min\left|\{D:D\mathrm{is}\mathrm{minimal}\mathrm{dominating}\mathrm{set}\}\right|$, and $\Gamma \left(\zeta \right)=max\left|\left\{D:D\mathrm{is}\mathrm{minimal}\mathrm{dominating}\mathrm{set}\right\}\right|.$ Notation $T_m\left(\zeta \right)$ which indicates the total number of minimal dominating set of $\zeta$ and $\rho$$\left(\zeta \right)={\sum }_{v\in V\left(\zeta \right)}{d}_p\left(v\right).$ A graph is said to be full degree graph if the vertices are all of full degree. The $m-$book graph is defined as the graph Cartesian product $B_m=S_{(m+1)}\times P_2,$ where $S_m$ is a star graph and $P_2$ is the path graph on two vertices.The Windmill graph $W{d}_r^s$ is an undirected graph constructed for $r\ge 2$ and $s\ge 2$ by s copies of the complete graph $K_r$ at a shared universal vertex. Examples of these two graphs are desiplies in Figure 1 and Figure 2. A join of two graphs ${\zeta }_1$ and ${\zeta }_2$ is denoted by ${\zeta }_1+{\zeta }_2$, with disjoint vertex sets $V_1$ and $V_2$ is the graph on the vertex set $V_1\cup V_2$ and the edge set $E_1\cup E_2\cup \{v_1{v}_2:v_1\in V_1,v_2\in V_2\}.$ The corona product of two graphs ${\zeta }_1$ and ${\zeta }_2$ is defined as the graph obtained by taking one copy of ${\zeta }_1$ and $V\left({\zeta }_1\right)$ copies of ${\zeta }_2$ and joining the $i^{-th}$ vertex of ${\zeta }_1$ to every vertex in the $i^{-th}$ copy of ${\zeta }_2$ an example in Figure 3.

Definition 1.

A $k-$domination regular graph is a type of graph in which every vertex has exactly k neighbors that are also adjacent to each other. In other words, every vertex is adjacent to exactly k vertices that are themselves mutually adjacent.

In this work, the concept of $k-$domination regular graph plays a crucial role, and therefore it is important to provide some literature review on this terminology. In the following, we will provide an overview of the relevant literature on $k-$domination regular graphs, highlighting their key properties and applications. By understanding the background and context of this concept, readers will be better equipped to comprehend the significance of our research and its contribution to the field. $K-$domination regular graphs have been studied in several papers in graph theory. For example, in [42], the authors investigate the structure of $k-$domination regular graphs and show that these graphs have many interesting properties. They also provide several examples of $k-$domination regular graphs and use them to study the domination number of these graphs. In [43] the authors studies the $k-$domination number of graphs and provides a characterization of $k-$domination regular graphs. The authors show that a connected graph is $k-$domination regular if and only if it is regular and satisfies a certain condition related to the $k-$domination number. They also investigate the relationship between the $k-$domination number and the total domination number of a graph. In [44] Hansberg et al. investigate some properties of $k-$domination regular graphs and provide some examples of these graphs. They also study the relationship between the $k-$domination number and the independence number of a graph, and show that for certain families of graphs, the $k-$domination number and the independence number are equal. Finally, in [45], the authors provide several methods for constructing $k-$domination regular graphs. They show that $k-$domination regular graphs can be constructed from other $k-$domination regular graphs by several operations, including Cartesian product and composition. They also provide some open problems related to $k-$domination regular graphs, such as finding the smallest k for which a $k-$domination regular graph exists. In summary, $k-$domination regular graphs have been studied extensively in the literature, and they have many interesting properties and applications in graph theory. Further research on these graphs can lead to new insights and solutions in various fields, including computer science, chemistry, and physics.

Our paper aims to address several research questions related to the use of the dominating sigma index in predicting the physicochemical properties of organic compounds. Firstly, we provide an explanation of the dominating sigma index and its relationship with topological indices in mathematical chemistry. Next, we develop linear and non-linear models using the QSPR approach to predict the properties of interest and conduct a comparative analysis to evaluate the effectiveness of the domination sigma index and establish an appropriate domination index that correlates with the physicochemical properties of octane and its isomers. Finally, we explore the relationship between topological indices and molecular properties. The paper is structured as follows: Section 3 presents the main results of the domination sigma index, while Section 4 describes the comparative analysis conducted to evaluate its effectiveness and establish an appropriate domination index. Finally, in Section 5, we provide our conclusions and suggestions for further research into QSAR/QSPR domination indices.

3. Main Results

In this Section, we give results of domination Sigma index for the star, complete bipartite and its complement, book graphs, and for the Windmill graph.

Proposition 1.

• Let $\zeta \simeq S_{r+1}$ the star graph of $r+1$ vertices, then $D\sigma \left(S_{r+1}\right)=0.$
• Let $\zeta \simeq K_r$ the complete graph of r vertices, then $D\sigma \left(K_r\right)=0.$
• Let $\zeta \simeq S_{r+1,s+1}$ the double star graph, then $D\sigma \left(S_{r+1.s+1}\right)=0.$
• Let $\zeta \simeq K_{r,s}$ the complete bipartite graph where $r,s\ge 2,$ then $D\sigma \left(K_{r,s}\right)={\sum }_{uv\in E\left(\zeta \right)}{(s-r)}^2={(s-r)}^{2}rs.$

Proof. 

• Suppose $\zeta \simeq S_{r+1}$ the star graph of r vertices, we have $T_m\left(S_{r+1}\right)=2$ and $d_{p_u}=1$ for all $u\in V\left(S_{r+1}\right).$ Then $D\sigma \left(S_{r+1}\right)=0.$
• Suppose $\zeta \simeq K_r$ the complete graph of r vertices, we have $T_m\left(K_r\right)=r$ and $d_{p_u}=1$ for all $u\in V\left(K_r\right).$ Then $D\sigma \left(K_r\right)=0.$
• Suppose $\zeta \simeq S_{r+1,s+1}$ the double star graph, we have $T_m\left(S_{r+1.s+1}\right)=4$ and $d_{p_u}=2$ for all $u\in V\left(S_{r+1.s+1}\right).$ Then $D\sigma \left(S_{r+1.s+1}\right)=0.$
• Suppose $\zeta \simeq K_{r,s}$ the complete bipartite graph where $r,s\ge 2$, we have $T_m\left(K_{r,s}\right)=rs+2$ and

  $$d_{p_u}=\left\{\begin{array}{c}r+1\mathrm{if}u\in \mathrm{the}\mathrm{set}\mathrm{of}\mathrm{vertices}\mathrm{that}\mathrm{contains}s\mathrm{vertices}\\ s+1\mathrm{if}u\in \mathrm{the}\mathrm{set}\mathrm{of}\mathrm{vertices}\mathrm{that}\mathrm{contains}r\mathrm{vertices}\end{array}\right.$$
  Then

  $$D\sigma \left(K_{r,s}\right)=\sum _{uv\in E\left(\zeta \right)}{(\left(s+1\right)-\left(r+1\right))}^2=\sum _{uv\in E\left(\zeta \right)}{(s-r)}^2={(s-r)}^{2}rs.$$
  This completes the proof.

□

Corollary 1.

Any graph ζ is $k-$domination regular graph if and only if $D\sigma \left(\zeta \right)=0.$

Proof. 

In this case, the sufficient condition is clear. If we consider the necessity condition, $D\sigma \left(\zeta \right)=0,$

$$\begin{array}{cc}\hfill & ⟹\sum _{uv\in E\left(\zeta \right)}{(d_{p_u}-d_{p_v})}^2=0\hfill \\ \hfill & ⟹{(d_{p_u}-d_{p_v})}^2=0\hfill \\ \hfill & ⟹d_{p_u}-d_{p_v}=0\hfill \\ \hfill & ⟹d_{p_u}=d_{p_v}\hfill \end{array}$$

for all $u,v\in V\left(\zeta \right)$, so if $d_{p_u}=d_{p_v}=k$ then $\zeta$ is $k-$domination regular graph. □

Corollary 2.

Let ζ be the complete bipartite graph $K_{r,s}.$ Then $D\sigma \left(K_{r,s}\right)=\sigma \left(K_{r,s}\right).$

Proof. 

We have

$$\begin{array}{cc}\hfill D\sigma \left(K_{r,s}\right)& =\sum _{uv\in E\left(\zeta \right)}{(d_{p_u}-d_{p_v})}^2=\sum _{uv\in E\left(K_{r,s}\right)}{(\left(d_u+1\right)-\left(d_v+1\right))}^2\hfill \\ \hfill & =\sum _{uv\in E\left(K_{r,s}\right)}{(d_u-d_v)}^2=\sigma \left(K_{r,s}\right).\hfill \end{array}$$

□

Proposition 2.

If $\zeta \simeq K_{r,s},$ then $D\sigma \left(\overline{K_{r,s}}\right)=\overline{\sigma \left(K_{r,s}\right)}$

Proof. 

Since $d_{p_{\overline{K_{r,s}}}}=d_{K_{r,s}},$ hence we have

$$\begin{array}{cc}\hfill D\sigma \left(\overline{K_{r,s}}\right)& =\sum _{uv\in E\left(\overline{K_{r,s}}\right)}{(d_{p_u}-d_{p_v})}^2\hfill \\ \hfill & =\sum _{uv\notin E\left(K_{r,s}\right)}{(d_u-d_v)}^2=\overline{\sigma \left(K_{r,s}\right)}.\hfill \end{array}$$

□

Theorem 1.

If $\zeta \simeq W{d}_r^s$ where $s,r\ge 2.$ Then $D\sigma \left(W{d}_r^s\right)=s(r-1){(1-{(r-1)}^{s-1})}^2.$

Proof. 

If $\zeta \simeq W{d}_r^s$ where $s,r\ge 2,$ we have $T_m\left(\zeta \right)={\left(r-1\right)}^s+1$ and

$$d_{p_u}=\left\{\begin{array}{c}1\mathrm{if}u\mathrm{is}\mathrm{the}\mathrm{center}\mathrm{vertex},\\ {\left(r-1\right)}^{s-1}\mathrm{otherwise}.\end{array}\right.$$

Suppose

$$\begin{array}{cc}\hfill E_1& =\left\{e:e\in E\left(W{d}_r^s\right),\mathrm{all}\mathrm{edges}\mathrm{which}\mathrm{are}\mathrm{incident}\mathrm{with}\mathrm{the}\mathrm{center}\mathrm{vertex}\right\}.\hfill \\ \hfill E_2& =\left\{e:e\in E\left(W{d}_r^s\right),\mathrm{all}\mathrm{edges}\mathrm{of}\mathrm{the}\mathrm{complete}\mathrm{graph}\right\}.\hfill \end{array}$$

Then we obtain

$$\begin{array}{cc}\hfill D\sigma \left(W{d}_r^s\right)& =\sum _{uv\in E_1\left(\zeta \right)}{(1-{\left(r-1\right)}^{s-1})}^2+\sum _{uv\in E_2\left(\zeta \right)}{({\left(r-1\right)}^{s-1}-{\left(r-1\right)}^{s-1})}^2\hfill \\ \hfill & ={(1-{\left(r-1\right)}^{s-1})}^2\left|E_1\right|\hfill \\ \hfill & ={(1-{\left(r-1\right)}^{s-1})}^2\left(s\left(r-1\right)\right).\hfill \end{array}$$

Which complete the proof. □

Theorem 2.

If $\zeta \simeq B_r$ the book graph for $r\ge 3.$ Then $D\sigma \left(B_r\right)=r\left(8-2^{r+2}+2^{2r-1}\right).$

Proof. 

If $\zeta \simeq B_r$ the book graph for $r\ge 3,$ we have $T_m\left(B_r\right)=2^r+3$ and

$$d_{p_u}=\left\{\begin{array}{c}3\mathrm{if}u\mathrm{is}\mathrm{the}\mathrm{center}\mathrm{vertex},\\ 2^{r-1}+1\mathrm{otherwise}.\end{array}\right.$$

Suppose $E_1$ denote the set of r edges $u_i{v}_i$ with initial and terminal vertices of the same domination degree $\left(2^{r-1}+1\right).$ Let $E_2$ denote the set containing only one edge $uv$ with initial and terminal vertices of the same domination degree 3. Let $E_3$ denote the set of $2r$ edges of initial vertices of the domination degree 3 and terminal vertices of domination degree $\left(2^{r-1}+1\right).$ Hence

$$\begin{array}{cc}\hfill D\sigma \left(B_r\right)& =\sum _{uv\in E_1\left(\zeta \right)}{(\left(\left(2^{r-1}+1\right)\right)-\left(2^{r-1}+1\right))}^2+\hfill \\ \hfill & \sum _{uv\in E_2\left(\zeta \right)}{(3-3)}^2+\sum _{uv\in E_3\left(\zeta \right)}{(3-\left(2^{r-1}+1\right))}^2\hfill \\ \hfill & =2r{\left(2-2^{r-1}\right)}^2=r\left(8-2^{r+2}+2^{2r-1}\right).\hfill \end{array}$$

Which complete the proof. □

Theorem 3.

If $\zeta \simeq K_{n_1,n_2,...,n_k}$ where $n_1,n_2,...,n_k\ge 2$. Then $D\sigma \left(K_{n_1,n_2,...,n_k}\right)=\sigma \left(K_{n_1,n_2,...,n_k}\right).$

Proof. 

Suppose that $\zeta \simeq K_{n_1,n_2,...,n_k}$ where $n_1,n_2,...,n_k\ge 2$, note that for any vertex $v\in K_{n_1,n_2,...,n_k}$ we have $d_{p_v}=d_v+1$ and $\left|E\left(K_{n_1,n_2,...,n_k}\right)\right|=T_m\left(K_{n_1,n_2,...,n_k}\right)-k,$ where $T_m\left(K_{n_1,n_2,...,n_k}\right)=\sum _{i=2}^k{n}_2{n}_i+...+n_{k-1}{n}_k+k.$ Hence

$$\begin{array}{cc}\hfill D\sigma \left(K_{n_1,n_2,...,n_k}\right)& =\sum _{uv\in E\left(\zeta \right)}{\left((d_u+1)-(d_v+1)\right)}^2\hfill \\ \hfill & =\sum _{uv\in E\left(\zeta \right)}{(d_u-d_v)}^2\hfill \\ \hfill & =\sigma \left(K_{n_1,n_2,...,n_k}\right).\hfill \end{array}$$

□

3.1. Domination Sigma Index of Some Graph Operations

This part, provides the domination Sigma index values for some graph operations such as corona product and join of graphs. In the following Theorem we calculate domination Sigma index for the corona graph of any graph $\zeta$ and the complete graph and its complement.

Theorem 4.

• For any connected graph ζ of $n_1$ vertices and $m_1$ edges. Then $D\sigma \left(\zeta \circ k_{n_2}\right)=0.$
• For any connected graph ζ of $n_1$ vertices and $m_1$ edges. Then $D\sigma \left(\zeta \circ \overline{k_{n_2}}\right)=0.$

Proof. 

(1)

Let $\xi \simeq \zeta \circ k_{n_2},$ we note that there are ${\left(n_2+1\right)}^{n_1}$ minimal dominating sets in $\xi$ and $d_{p_v}={\left(n_2+1\right)}^{n_1-1}.$ Also, there are three types of edges in $\xi .$ All edges of $\zeta ,$ all edges of $k_{n_2},$ and $E_1$ denote the set of all edges that connect vertex from $\zeta$ and a vertex from $k_{n_2}.$ So, we have

$$\begin{array}{cc}\hfill D\sigma \left(\zeta \circ k_{n_2}\right)& =\sum _{uv\in E\left(\zeta \right)}{(d_{p_u}-d_{p_v})}^2+\sum _{uv\in E\left(k_{n_2}\right)}{(d_{p_u}-d_{p_v})}^2+\sum _{uv\in E_1}{(d_{p_u}-d_{p_v})}^2\hfill \\ \hfill & =\sum _{uv\in E\left(\zeta \right)}{({\left(n_2+1\right)}^{n_1-1}-{\left(n_2+1\right)}^{n_1-1})}^2\hfill \\ \hfill & +\sum _{uv\in E\left(k_{n_2}\right)}{({\left(n_2+1\right)}^{n_1-1}-{\left(n_2+1\right)}^{n_1-1})}^2\hfill \\ \hfill & +\sum _{uv\in E_1}{({\left(n_2+1\right)}^{n_1-1}-{\left(n_2+1\right)}^{n_1-1})}^2=0+0+0=0.\hfill \end{array}$$

(2)

Let $\xi \simeq \zeta \circ \overline{k_{n_2}},$ for any vertex $v\in V\left(\xi \right)$ we have $T_m\left(\xi \right)=\sum _{i=1}^{n_1}\left(\begin{array}{c}{n}_1\\ i\end{array}\right)$ and $d_{p_v}=T_m\left(\xi \right)-2^{n_1-1}.$ Hence $\xi$ is $k-$ domination regular graph, where $k=T_m\left(\xi \right)-2^{n_1-1}$ which implies that $D\sigma \left(\zeta \circ \overline{k_{n_2}}\right)=0.$

□

In the following Theorems we calculate domination Sigma index for different cases of join of two graphs. But first we need the following Lemma

Lemma 1.

[35] Let ${\zeta }_1$ and ${\zeta }_2$ be any non complete graphs of $n_1$ and $n_2$ vertices respectively, and there is no vertex in ${\zeta }_1$ or ${\zeta }_2$ of full degree. Then $T_m\left({\zeta }_1+{\zeta }_2\right)=T_m\left({\zeta }_1)+T_m({\zeta }_2\right)+n_1{n}_2$ and

$$d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right)=\left\{\begin{array}{c}{d}_{p_{{\zeta }_1}}\left(v\right)+n_2\mathrm{if}v\in V\left({\zeta }_1\right)\\ d_{p_{{\zeta }_2}}\left(v\right)+n_1\mathrm{if}v\in V\left({\zeta }_2\right).\end{array}\right.$$

Theorem 5.

Let ${\zeta }_1$ and ${\zeta }_2$ be any non complete graphs of $n_1$ and $n_2$ vertices respectively, and there is no vertex in ${\zeta }_1$ or ${\zeta }_2$ of full degree. Then

$$\begin{array}{cc}\hfill D\sigma \left({\zeta }_1+{\zeta }_2\right)& =D\sigma \left({\zeta }_1)+D\sigma ({\zeta }_2\right)+n_2\sum _{u\in V\left({\zeta }_1\right)}{(d_{p_{{\zeta }_1}}\left(u\right)+n_2)}^2\hfill \\ \hfill & +n_1\sum _{v\in V\left({\zeta }_2\right)}{(d_{p_{{\zeta }_2}}\left(v\right)+n_1)}^2-2\rho \left({\zeta }_1\right)\rho \left({\zeta }_2\right)\hfill \\ \hfill & -2n_1{n}_2\rho \left({\zeta }_1\right)-2n_1{n}_2\rho \left({\zeta }_2\right)-2n_1^2{n}_2^2.\hfill \end{array}$$

Proof. 

Let ${\zeta }_1$ and ${\zeta }_2$ be any non complete graphs of $n_1$ and $n_2$ vertices respectively, and there is no vertex in ${\zeta }_1$ or ${\zeta }_2$ of full degree. Then

$$\begin{array}{cc}\hfill D\sigma \left({\zeta }_1+{\zeta }_2\right)& =\stackrel{\left(1\right)}{\overbrace{{\sum }_{uv\in E\left({\zeta }_1\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}+\stackrel{\left(2\right)}{\overbrace{{\sum }_{uv\in E\left({\zeta }_2\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}\hfill \\ \hfill & +\stackrel{\left(3\right)}{\overbrace{{\sum }_{\begin{array}{c}u\in V\left({\zeta }_1\right)\\ v\in V\left({\zeta }_2\right)\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}\hfill \end{array}$$

For (1):

$$\begin{array}{cc}\hfill \sum _{uv\in E\left({\zeta }_1\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{uv\in E\left({\zeta }_1\right)}(d_{p_{{\zeta }_1}}\left(u\right)+n_2)-(d_{p_{{\zeta }_1}}\left(v\right)+n_2){)}^2\hfill \\ \hfill & =\sum _{uv\in E\left({\zeta }_1\right)}{(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_1}}\left(v\right))}^2=D\sigma \left({\zeta }_1\right)\hfill \end{array}$$

For (2):

$$\begin{array}{cc}\hfill \sum _{uv\in E\left({\zeta }_2\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{uv\in E\left({\zeta }_2\right)}(d_{p_{{\zeta }_2}}\left(u\right)+n_1)-(d_{p_{{\zeta }_2}}\left(v\right)+n_1){)}^2\hfill \\ \hfill & =\sum _{uv\in E\left({\zeta }_2\right)}{(d_{p_{{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right))}^2=D\sigma \left({\zeta }_2\right)\hfill \end{array}$$

For (3):

$$\begin{array}{cc}\hfill & \sum _{{}_{\begin{array}{c}u\in V\left({\zeta }_1\right)\\ v\in V\left({\zeta }_2\right)\end{array}}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2\hfill \\ \hfill & ={\left[\left(d_{p_{{\zeta }_1}}\left(u_1\right)+n_2\right)-\left(d_{p_{{\zeta }_2}}\left(v_1\right)+n_1\right)\right]}^2+\hfill \\ \hfill & {\left[\left(d_{p_{{\zeta }_1}}\left(u_1\right)+n_2\right)-\left(d_{p_{{\zeta }_2}}\left(v_2\right)+n_1\right)\right]}^2+...+\hfill \\ \hfill & {\left[\left(d_{p_{{\zeta }_1}}\left(u_1\right)+n_2\right)-\left(d_{p_{{\zeta }_2}}\left(v_{n_2}\right)+n_1\right)\right]}^2+\hfill \\ \hfill & {\left[\left(d_{p_{{\zeta }_1}}\left(u_2\right)+n_2\right)-\left(d_{p_{{\zeta }_2}}\left(v_1\right)+n_1\right)\right]}^2+\hfill \\ \hfill & {\left[\left(d_{p_{{\zeta }_1}}\left(u_2\right)+n_2\right)-\left(d_{p_{{\zeta }_2}}\left(v_2\right)+n_1\right)\right]}^2+...+\hfill \\ \hfill & {\left[\left(d_{p_{{\zeta }_1}}\left(u_2\right)+n_2\right)-\left(d_{p_{{\zeta }_2}}\left(v_{n_2}\right)+n_1\right)\right]}^2+...+\hfill \\ \hfill & {\left[\left(d_{p_{{\zeta }_1}}\left(u_{n_1}\right)+n_2\right)-\left(d_{p_{{\zeta }_2}}\left(v_1\right)+n_1\right)\right]}^2+\hfill \\ \hfill & {\left[\left(d_{p_{{\zeta }_1}}\left(u_{n_1}\right)+n_2\right)-\left(d_{p_{{\zeta }_2}}\left(v_2\right)+n_1\right)\right]}^2+...+\hfill \\ \hfill & {\left[\left(d_{p_{{\zeta }_1}}\left(u_{n_1}\right)+n_2\right)-\left(d_{p_{{\zeta }_2}}\left(v_{n_2}\right)+n_1\right)\right]}^2\hfill \\ \hfill & =n_2{\left(d_{p_{{\zeta }_1}}\left(u_1\right)+n_2\right)}^2+\sum _{v\in V\left({\zeta }_2\right)}{\left(d_{p_{{\zeta }_2}}\left(v\right)+n_1\right)}^2-\hfill \\ \hfill & 2d_{p_{{\zeta }_1}}\left(u_1\right)\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)-2n_1{n}_2{d}_{p_{{\zeta }_1}}\left(u_1\right)-\hfill \\ \hfill & 2n_2\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)-2n_1{n}_2^2+n_2{\left(d_{p_{{\zeta }_1}}\left(u_2\right)+n_2\right)}^2+\hfill \\ \hfill & \sum _{v\in V\left({\zeta }_2\right)}{\left(d_{p_{{\zeta }_2}}\left(v\right)+n_1\right)}^2-\hfill \\ \hfill & 2d_{p_{{\zeta }_1}}\left(u_2\right)\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)-2n_1{n}_2{d}_{p_{{\zeta }_1}}\left(u_2\right)-\hfill \\ \hfill & 2n_1{n}_2{d}_{p_{{\zeta }_1}}\left(u_2\right)-2n_2\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)-2n_1{n}_2^2+...+\hfill \\ \hfill & n_2{\left(d_{p_{{\zeta }_1}}\left(u_{n_1}\right)+n_2\right)}^2+\sum _{v\in V\left({\zeta }_2\right)}{\left(d_{p_{{\zeta }_2}}\left(v\right)+n_1\right)}^2\hfill \\ \hfill & -2d_{p_{{\zeta }_1}}\left(u_{n_1}\right)\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)-2n_1{n}_2{d}_{p_{\zeta }}\left(u_{n_1}\right)-\hfill \\ \hfill & 2n_2\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)-2n_1{n}_2^2\hfill \\ \hfill & =n_2\sum _{u\in V\left({\zeta }_1\right)}{\left(d_{p_{{\zeta }_1}}\left(u\right)+n_2\right)}^2+n_1\sum _{v\in V\left({\zeta }_2\right)}{\left(d_{p_{{\zeta }_2}}\left(v\right)+n_1\right)}^2\hfill \\ \hfill & -2\sum _{u\in V\left({\zeta }_1\right)}{d}_{p_{{\zeta }_1}}\left(u\right)\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)-\hfill \\ \hfill & 2n_1{n}_2\sum _{u\in V\left({\zeta }_1\right)}{d}_{p_{{\zeta }_1}}\left(u\right)-2n_1{n}_2\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)-2n_1^2{n}_2^2\hfill \\ \hfill & =n_2\sum _{u\in V\left({\zeta }_1\right)}{\left(d_{p_{{\zeta }_1}}\left(u\right)+n_2\right)}^2+n_1\sum _{v\in V\left({\zeta }_2\right)}{\left(d_{p_{{\zeta }_2}}\left(v\right)+n_1\right)}^2-2\rho \left({\zeta }_1\right)\rho \left({\zeta }_2\right)-2n_1{n}_2\rho \left({\zeta }_1\right)-\hfill \\ \hfill & 2n_1{n}_2\rho \left({\zeta }_2\right)-2n_1^2{n}_2^2.\hfill \end{array}$$

Now adding (1), (2) and (3) we conclude our result. □

Using the fact [35] that for any graph $\zeta ,$ we have $\rho$$\left(\zeta \right)\le T_m\left({\zeta }_1\right)\Gamma \left(\zeta \right).$ We conclude the following

Corollary 3.

Let ${\zeta }_1$ and ${\zeta }_2$ be any non complete graphs of $n_1$ and $n_2$ vertices respectively, and there is no vertex in ${\zeta }_1$ or ${\zeta }_2$ of full degree. Then

$$\begin{array}{cc}\hfill D\sigma \left({\zeta }_1+{\zeta }_2\right)& \le D\sigma \left({\zeta }_1)+D\sigma ({\zeta }_2\right)+n_2\sum _{u\in V\left({\zeta }_1\right)}{(d_{p_{{\zeta }_1}}\left(u\right)-n_2)}^2+n_1\sum _{v\in V\left({\zeta }_2\right)}{(d_{p_{{\zeta }_2}}\left(v\right)-n_1)}^2-\hfill \\ \hfill & 2T_m\left({\zeta }_1\right)\Gamma \left({\zeta }_1\right)T_m\left({\zeta }_2\right)\Gamma \left({\zeta }_2\right)-2n_1{n}_2{T}_m\left({\zeta }_1\right)\Gamma \left({\zeta }_1\right)-2n_1{n}_2{T}_m\left({\zeta }_2\right)\Gamma \left({\zeta }_2\right)-2n_1^2{n}_2^2.\hfill \end{array}$$

Proposition 3.

Let ${\zeta }_1$ and ${\zeta }_2$ are both complete graphs of $n_1$ and $n_2$ vertices respectively. Then $D\sigma \left({\zeta }_1+{\zeta }_2\right)=0.$

Proof. 

Let ${\zeta }_1$ and ${\zeta }_2$ are both complete graphs of $n_1$ and $n_2$ vertices respectively, then $T_m\left({\zeta }_1+{\zeta }_2\right)=n_1+n_2$ and $d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right)=1$ for all $v\in V\left({\zeta }_1+{\zeta }_2\right)$ which implies that ${\zeta }_1+{\zeta }_2$ is $1-$domination regular graph. Hence $D\sigma \left({\zeta }_1+{\zeta }_2\right)=0.$ □

Theorem 6.

If ${\zeta }_1$ is a complete graph and ${\zeta }_2$ is not a complete graphs of $n_1$ and $n_2$ vertices respectively. Then

$$D\sigma \left({\zeta }_1+{\zeta }_2\right)=D\sigma \left({\zeta }_2\right)+n_1{n}_2+n_{1}D{M}_1\left({\zeta }_2\right)-2n_1\rho \left({\zeta }_2\right).$$

Proof. 

Let ${\zeta }_1$ be a complete graph and ${\zeta }_2$ is not a complete graphs of $n_1$ and $n_2$ vertices respectively, then $T_m\left({\zeta }_1+{\zeta }_2\right)=n_1+T_m\left({\zeta }_2\right)$ and

$$d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right)=\left\{\begin{array}{c}1\mathrm{if}v\in V\left({\zeta }_1\right)\\ d_{p_{{\zeta }_2}}\left(v\right)\mathrm{if}v\in V\left({\zeta }_2\right).\end{array}\right.$$

Hence,

$$\begin{array}{cc}\hfill D\sigma \left({\zeta }_1+{\zeta }_2\right)& =\stackrel{\left(1\right)}{\overbrace{{\sum }_{uv\in E\left({\zeta }_1\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}+\hfill \\ \hfill & \stackrel{\left(2\right)}{\overbrace{{\sum }_{uv\in E\left({\zeta }_2\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}+\hfill \\ \hfill & \stackrel{\left(3\right)}{\overbrace{{\sum }_{\begin{array}{c}u\in V\left({\zeta }_1\right)\\ v\in V\left({\zeta }_2\right)\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}\hfill \end{array}$$

For (1):

$$\sum _{uv\in E\left({\zeta }_1\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2=\sum _{uv\in E\left({\zeta }_1\right)}{(1-1)}^2=0.$$

For (2):

$$\sum _{uv\in E\left({\zeta }_2\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2=\sum _{uv\in E\left({\zeta }_2\right)}{(d_{p_{{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right))}^2=D\sigma \left({\zeta }_2\right).$$

For (3):

$$\begin{array}{cc}\hfill \sum _{\begin{array}{c}u\in V\left({\zeta }_1\right)\\ v\in V\left({\zeta }_2\right)\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& ={\left(1-d_{p_{{\zeta }_2}}\left(v_1\right)\right)}^2+{\left(1-d_{p_{{\zeta }_2}}\left(v_2\right)\right)}^2+...+{\left(1-d_{p_{{\zeta }_2}}\left(v_{n_2}\right)\right)}^2+\hfill \\ \hfill & {\left(1-d_{p_{{\zeta }_2}}\left(v_1\right)\right)}^2+{\left(1-d_{p_{{\zeta }_2}}\left(v_2\right)\right)}^2\hfill \\ \hfill & +...+{\left(1-d_{p_{{\zeta }_2}}\left(v_{n_2}\right)\right)}^2+...+\hfill \\ \hfill & {\left(1-d_{p_{{\zeta }_2}}\left(v_1\right)\right)}^2+{\left(1-d_{p_{{\zeta }_2}}\left(v_2\right)\right)}^2\hfill \\ \hfill & +...+{\left(1-d_{p_{{\zeta }_2}}\left(v_{n_2}\right)\right)}^2\hfill \\ \hfill & =(1+d_{p_{{\zeta }_2}}^2\left(v_1\right)-2d_{p_{{\zeta }_2}}\left(v_1\right))+\hfill \\ \hfill & (1+d_{p_{{\zeta }_2}}^2\left(v_2\right)-2d_{p_{{\zeta }_2}}\left(v_2\right))+...+\hfill \\ \hfill & (1+d_{p_{{\zeta }_2}}^2\left(v_{n_2}\right)-2d_{p_{{\zeta }_2}}\left(v_{n_2}\right))+\hfill \\ \hfill & (1+d_{p_{{\zeta }_2}}^2\left(v_1\right)-2d_{p_{{\zeta }_2}}\left(v_1\right))+\hfill \\ \hfill & (1+d_{p_{{\zeta }_2}}^2\left(v_2\right)-2d_{p_{{\zeta }_2}}\left(v_2\right))+...+\hfill \\ \hfill & (1+d_{p_{{\zeta }_2}}^2\left(v_{n_2}\right)-2d_{p_{{\zeta }_2}}\left(v_{n_2}\right))+...+\hfill \\ \hfill & (1+d_{p_{{\zeta }_2}}^2\left(v_1\right)-2d_{p_{{\zeta }_2}}\left(v_1\right))+\hfill \\ \hfill & (1+d_{p_{{\zeta }_2}}^2\left(v_2\right)-2d_{p_{{\zeta }_2}}\left(v_2\right))+...+\hfill \\ \hfill & (1+d_{p_{{\zeta }_2}}^2\left(v_{n_2}\right)-2d_{p_{{\zeta }_2}}\left(v_{n_2}\right))\hfill \\ \hfill & =\left(n_2+\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}^2\left(v\right)-2\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)\right)+\hfill \\ \hfill & \left(n_2+\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}^2\left(v\right)-2\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)\right)+...+\hfill \\ \hfill & \left(n_2+\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}^2\left(v\right)-2\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)\right)\hfill \\ \hfill & =n_1{n}_2+n_1\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}^2\left(v\right)-2n_1\sum _{v\in V\left({\zeta }_2\right)}{d}_{p_{{\zeta }_2}}\left(v\right)\hfill \\ \hfill & =n_1{n}_2+n_{1}DM\left({\zeta }_2\right)-2n_1\rho \left({\zeta }_2\right).\hfill \end{array}$$

Adding $\left(1\right),\left(2\right)$, and $\left(3\right)$ we complete the proof. □

Theorem 7.

If ${\zeta }_1$ is not a complete graph and ${\zeta }_2$ is a complete graph of $n_1$ and $n_2$ vertices respectively. Then

$$D\sigma \left({\zeta }_1+{\zeta }_2\right)=D\sigma \left({\zeta }_1\right)+n_1{n}_2+n_{2}D{M}_1\left({\zeta }_1\right)-2n_1\rho \left({\zeta }_1\right).$$

Proof. 

Let ${\zeta }_1$ be not complete graph and ${\zeta }_2$ is a complete graphs of $n_1$ and $n_2$ vertices respectively, then $T_m\left({\zeta }_1+{\zeta }_2\right)=n_2+T_m\left({\zeta }_1\right)$ and

$$d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right)=\left\{\begin{array}{c}1\mathrm{if}v\in V\left({\zeta }_2\right)\\ d_{p_{{\zeta }_1}}\left(v\right)\mathrm{if}v\in V\left({\zeta }_1\right).\end{array}\right.$$

Hence,

$$\begin{array}{cc}\hfill D\sigma \left({\zeta }_1+{\zeta }_2\right)& =\stackrel{\left(1\right)}{\overbrace{{\sum }_{uv\in E\left({\zeta }_1\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}+\stackrel{\left(2\right)}{\overbrace{{\sum }_{uv\in E\left({\zeta }_2\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}+\hfill \\ \hfill & \stackrel{\left(3\right)}{\overbrace{{\sum }_{\begin{array}{c}u\in V\left({\zeta }_1\right)\\ v\in V\left({\zeta }_2\right)\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}\hfill \end{array}$$

For (1):

$$\sum _{uv\in E\left({\zeta }_1\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2=\sum _{uv\in E\left({\zeta }_2\right)}{(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_1}}\left(v\right))}^2=D\sigma \left({\zeta }_1\right).$$

For (2):

$$\sum _{uv\in E\left({\zeta }_2\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2=\sum _{uv\in E\left({\zeta }_2\right)}{(1-1)}^2=0.$$

For (3):

$$\begin{array}{cc}\hfill \sum _{\begin{array}{c}u\in V\left({\zeta }_1\right)\\ v\in V\left({\zeta }_2\right)\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& ={\left(d_{p_{{\zeta }_1}}\left(u_1\right)-1\right)}^2+{\left(d_{p_{{\zeta }_1}}\left(u_1\right)-1\right)}^2+...+{\left(d_{p_{{\zeta }_1}}\left(v_1\right)-1\right)}^2+\hfill \\ \hfill & {\left(d_{p_{{\zeta }_1}}\left(u_2\right)-1\right)}^2+{\left(d_{p_{{\zeta }_1}}\left(u_2\right)-1\right)}^2+...\hfill \\ \hfill & +{\left(d_{p_{{\zeta }_1}}\left(v_{n_2}\right)-1\right)}^2+...+\hfill \\ \hfill & {\left(d_{p_{{\zeta }_1}}\left(u_{n_1}\right)-1\right)}^2+{\left(d_{p_{{\zeta }_1}}\left(u_{n_1}\right)-1\right)}^2\hfill \\ \hfill & +...+{\left(d_{p_{{\zeta }_1}}\left(u_{n_1}\right)-1\right)}^2\hfill \\ \hfill & =n_2{\left(d_{p_{{\zeta }_1}}\left(u_1\right)-1\right)}^2+n_2{\left(d_{p_{{\zeta }_1}}\left(u_2\right)-1\right)}^2\hfill \\ \hfill & +...+n_2{\left(d_{p_{{\zeta }_1}}\left(u_{n_1}\right)-1\right)}^2\hfill \\ \hfill & =n_2\sum _{u\in V\left({\zeta }_1\right)}{\left(d_{p_{{\zeta }_1}}\left(u\right)-1\right)}^2\hfill \\ \hfill & =n_2\sum _{u\in V\left({\zeta }_1\right)}\left(d_{p_{{\zeta }_1}}^2\left(u\right)-2d_{p_{{\zeta }_1}}\left(u\right)+1\right)\hfill \\ \hfill & =n_2\sum _{u\in V\left({\zeta }_1\right)}{d}_{p_{{\zeta }_1}}^2\left(u\right)-2n_2\sum _{u\in V\left({\zeta }_1\right)}{d}_{p_{{\zeta }_1}}\left(u\right)+\sum _{u\in V\left({\zeta }_1\right)}1\hfill \\ \hfill & =n_2{n}_1+n_{2}DM\left({\zeta }_1\right)-2n_1\rho \left({\zeta }_1\right).\hfill \end{array}$$

Adding $\left(1\right),\left(2\right)$, and $\left(3\right)$ we complete the proof. □

Lemma 2.

If ${\zeta }_1$ and ${\zeta }_2$ are any non complete graphs of $\left(n_1,m_1\right)$ and $\left(n_2,m_2\right)$ vertices and edges respectively, such that $A\subseteq {\zeta }_1$ contains all vertices of full degrees and $B\subseteq {\zeta }_2$ contains all vertices of full degrees. Then $T_m\left({\zeta }_1+{\zeta }_2\right)=T_m\left({\zeta }_1\right)+T_m\left({\zeta }_2\right)+\left(n_1-A\right)\left(n_2-B\right)$ and

$$d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right)=\left\{\begin{array}{c}{d}_{p_{{\zeta }_1}}\left(v\right)+\left(n_2-\left|B\right|\right)\mathrm{if}v\in V\left({\zeta }_1\right)\mathrm{and}v\notin A,\\ d_{p_{{\zeta }_2}}\left(v\right)\left(n_1-\left|A\right|\right)\mathrm{if}v\in V\left({\zeta }_2\right)\mathrm{and}v\notin B,\\ d_{p_{{\zeta }_1}}\left(v\right)\mathrm{if}v\in V\left({\zeta }_1\right)\mathrm{and}v\in A,\\ d_{p_{{\zeta }_2}}\left(v\right)\mathrm{if}v\in V\left({\zeta }_2\right)\mathrm{and}v\in B.\end{array}\right.$$

Theorem 8.

If ${\zeta }_1$ and ${\zeta }_2$ are any non complete graphs of $\left(n_1,m_1\right)$ and $\left(n_2,m_2\right)$ vertices and edges respectively, such that $A\subseteq V\left({\zeta }_1\right)$ contains all vertices of full degrees and $B\subseteq V\left({\zeta }_2\right)$ contains all vertices of full degrees. Then

$$\begin{array}{ccc}\hfill D\sigma \left({\zeta }_1+{\zeta }_2\right)& =& D\sigma \left({\zeta }_1\right)+D\sigma \left({\zeta }_2\right)+{\left(n_2-\left|B\right|\right)}^2\left|E_5\right|+{\left(n_1-\left|A\right|\right)}^2\left|E_6\right|-\hfill \\ & & 2\left[\begin{array}{c}\left(n_2-\left|B\right|\right)\sum _{\begin{array}{c}uv\in E_5\\ u\notin A\\ v\in A\end{array}}\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_1}}\left(v\right)\right)+\\ \left(n_1-\left|A\right|\right)\sum _{\begin{array}{c}uv\in E_6\\ u\notin B\\ v\in B\end{array}}\left(d_{p_{{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)\end{array}\right]+\hfill \\ & & \sum _{\begin{array}{c}uv\in E_8\\ u\in V\left({\zeta }_1\right),u\in A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}{\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)}^2+{\left(n_1-\left|A\right|\right)}^2\left|E_8\right|-\hfill \\ & & 2\left(n_1-\left|A\right|\right)\sum _{\begin{array}{c}uv\in E_8\\ u\in V\left({\zeta }_1\right),u\in A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)+\hfill \\ & & \sum _{\begin{array}{c}uv\in E_9\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\in B\end{array}}{\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)}^2+{\left(n_2-\left|B\right|\right)}^2\left|E_9\right|+\hfill \\ & & 2\left(n_2-\left|B\right|\right)\sum _{\begin{array}{c}uv\in E_9\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\in B\end{array}}\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)+\hfill \\ & & \sum _{\begin{array}{c}uv\in E_{10}\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}{\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)}^2+{\left(n_2-n_1+\left|A\right|-\left|B\right|\right)}^2\left|E_{10}\right|+\hfill \\ & & 2\left(n_2-n_1+\left|A\right|-\left|B\right|\right)\sum _{\begin{array}{c}uv\in E_{10}\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right).\hfill \end{array}$$

Proof. 

Suppose ${\zeta }_1$ and ${\zeta }_2$ are any non complete graphs of $\left(n_1,m_1\right)$ and $\left(n_2,m_2\right)$ vertices and edges respectively, such that $A\subseteq V\left({\zeta }_1\right)$ contains all vertices of full degrees and $B\subseteq V\left({\zeta }_2\right)$ contains all vertices of full degrees. Then $T_m\left({\zeta }_1+{\zeta }_2\right)=T_m\left({\zeta }_1\right)+T_m\left({\zeta }_2\right)+\left(n_1-A\right)\left(n_2-B\right)$ and

$$d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right)=\left\{\begin{array}{c}{d}_{p_{{\zeta }_1}}\left(v\right)+\left(n_2-\left|B\right|\right)\mathrm{if}v\in V\left({\zeta }_1\right)\mathrm{and}v\notin A,\\ d_{p_{{\zeta }_1}}\left(v\right)\mathrm{if}v\in V\left({\zeta }_1\right)\mathrm{and}v\in A,\\ d_{p_{{\zeta }_2}}\left(v\right)+\left(n_1-\left|A\right|\right)\mathrm{if}v\in V\left({\zeta }_2\right)\mathrm{and}v\notin B,\\ d_{p_{{\zeta }_2}}\left(v\right)\mathrm{if}v\in V\left({\zeta }_2\right)\mathrm{and}v\in B.\end{array}\right.$$

We divide the edge set of $\left({\zeta }_1+{\zeta }_2\right)$ according to the degree domination of the vertices as follows:

$$\begin{array}{cc}\hfill E_1& =\{e=uv:u,v\in A\},E_2=\{e=uv:u,v\in B\},\hfill \\ \hfill E_3& =\{e=uv:u,v\notin A\},E_4=\{e=uv:u,v\notin B\},\hfill \\ \hfill E_5& =\{e=uv:u\in A,v\in \left({\zeta }_1-A\right)\},E_6=\{e=uv:u\in B,v\in \left({\zeta }_2-B\right)\},\hfill \\ \hfill E_7& =\{e=uv:u\in A,v\in B\},E_8=\{e=uv:u\in A,v\in \left({\zeta }_2-B\right)\},\hfill \\ \hfill E_9& =\{e=uv:u\in \left({\zeta }_1-A\right),v\in B\},E_{10}=\{e=uv:u\in \left({\zeta }_1-A\right),v\in \left({\zeta }_2-B\right)\}.\hfill \end{array}$$

Hence,

$$\begin{array}{cc}\hfill D\sigma \left({\zeta }_1+{\zeta }_2\right)& =\stackrel{\left(1\right)}{\overbrace{{\sum }_{uv\in E\left({\zeta }_1\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}+\hfill \\ \hfill & \stackrel{\left(2\right)}{\overbrace{{\sum }_{uv\in E\left({\zeta }_2\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}+\hfill \\ \hfill & \stackrel{\left(3\right)}{\overbrace{{\sum }_{\begin{array}{c}u\in V\left({\zeta }_1\right)\\ v\in V\left({\zeta }_2\right)\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}\hfill \\ \hfill & =\sum _{uv\in E_1}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2+...+\sum _{uv\in E_{10}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2.\hfill \end{array}$$

Now we calculate the summation for each term as follows:

$$\begin{array}{cc}\hfill \sum _{uv\in E_1}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{uv\in E_1}{\left(1-1\right)}^2=0\hfill \\ \hfill \sum _{uv\in E_2}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{uv\in E_2}{\left(1-1\right)}^2=0\hfill \\ \hfill \sum _{uv\in E_3}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{uv\in E_3}{\left[\left(d_{p_{{\zeta }_1}}\left(u\right)+\left(n_2-\left|B\right|\right)\right)-\left(d_{p_{{\zeta }_1}}\left(v\right)+\left(n_2-\left|B\right|\right)\right)\right]}^2\hfill \\ \hfill & =\sum _{uv\in E_3}{(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_1}}\left(v\right))}^2\hfill \\ \hfill \sum _{uv\in E_4}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{uv\in E_4}{\left[\left(d_{p_{{\zeta }_2}}\left(u\right)+\left(n_1-\left|A\right|\right)\right)-\left(d_{p_{{\zeta }_2}}\left(v\right)+\left(n_1-\left|A\right|\right)\right)\right]}^2\hfill \\ \hfill & =\sum _{uv\in E_4}{(d_{p_{{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right))}^2\hfill \\ \hfill \sum _{\begin{array}{c}uv\in E_5\\ u\notin A\\ v\in A\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{\begin{array}{c}uv\in E_5\\ u\notin A\\ v\in A\end{array}}{\left[\left(d_{p_{{\zeta }_1}}\left(u\right)+\left(n_2-\left|B\right|\right)\right)-d_{p_{{\zeta }_1}}\left(v\right)\right]}^2\hfill \\ \hfill & =\sum _{\begin{array}{c}uv\in E_5\\ u\notin A\\ v\in A\end{array}}\left[\begin{array}{c}{d}_{p_{{\zeta }_1}}^2\left(u\right)+\left(n_2-\left|B\right|\right)d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_1}}\left(u\right)d_{p_{{\zeta }_1}}\left(v\right)\\ +\left(n_2-\left|B\right|\right)d_{p_{{\zeta }_1}}\left(u\right)+{\left(n_2-\left|B\right|\right)}^2-(n_2-\left|B\right|)d_{p_{{\zeta }_1}}\left(v\right)\\ -d_{p_{{\zeta }_1}}\left(u\right)d_{p_{{\zeta }_1}}\left(v\right)-(n_2-\left|B\right|)d_{p_{{\zeta }_1}}\left(v\right)+d_{p_{{\zeta }_1}}^2\left(v\right)\end{array}\right]\hfill \\ \hfill & =\sum _{\begin{array}{c}uv\in E_5\\ u\notin A\\ v\in A\end{array}}{(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_1}}\left(v\right))}^2+{(n_2-\left|B\right|)}^2\left|E_5\right|-\hfill \\ \hfill & 2(n_2-\left|B\right|)\sum _{\begin{array}{c}uv\in E_5\\ u\notin A\\ v\in A\end{array}}\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_1}}\left(v\right)\right).\hfill \end{array}$$

We notice that the edges of ${\zeta }_1$ are $E_1,E_3,$ and $E_5$. Hence

$$\begin{array}{cc}\hfill \stackrel{\left(1\right)}{\overbrace{{\sum }_{uv\in E\left({\zeta }_1\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}& =D\sigma \left({\zeta }_1\right)+{\left(n_2-\left|B\right|\right)}^2\left|E_5\right|-\hfill \\ \hfill & 2\left(n_2-\left|B\right|\right)\sum _{\begin{array}{c}uv\in E_5\\ u\notin A\\ v\in A\end{array}}\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_1}}\left(v\right)\right).\hfill \end{array}$$

Now

$$\begin{array}{cc}\hfill \sum _{\begin{array}{c}uv\in E_6\\ u\notin B\\ v\in B\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{\begin{array}{c}uv\in E_6\\ u\notin B\\ v\in B\end{array}}{\left[\left(d_{p_{{\zeta }_2}}\left(u\right)+\left(n_1-\left|A\right|\right)\right)-d_{p_{{\zeta }_2}}\left(v\right)\right]}^2\hfill \\ \hfill & =\sum _{\begin{array}{c}uv\in E_6\\ u\notin B\\ v\in B\end{array}}[{\left(d_{p_{{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)}^2+{\left(n_1-\left|A\right|\right)}^2-\hfill \\ \hfill & 2\left(n_1-\left|A\right|\right)\sum _{\begin{array}{c}uv\in E_6\\ u\notin B\\ v\in B\end{array}}\left(d_{p_{{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)]\hfill \\ \hfill & =\sum _{\begin{array}{c}uv\in E_6\\ u\notin B\\ v\in B\end{array}}[{\left(d_{p_{{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)}^2+{\left(n_1-\left|A\right|\right)}^2\left|E_6\right|-\hfill \\ \hfill & 2\left(n_1-\left|A\right|\right)\sum _{\begin{array}{c}uv\in E_6\\ u\notin B\\ v\in B\end{array}}\left(d_{p_{{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)].\hfill \end{array}$$

We notice that the edges of the graph ${\zeta }_2$ are the edges of $E_2,E_4,$ and $E_6.$ Hence

$$\begin{array}{cc}\hfill \stackrel{\left(2\right)}{\overbrace{{\sum }_{uv\in E\left({\zeta }_2\right)}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}& =D\sigma \left({\zeta }_2\right)+{\left(n_1-\left|A\right|\right)}^2\left|E_6\right|-\hfill \\ \hfill & 2\left(n_1-\left|A\right|\right)\sum _{\begin{array}{c}uv\in E_6\\ u\notin B\\ v\in B\end{array}}\left(d_{p_{{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right).\hfill \end{array}$$

Now we move to calculate the results on edges between vertices of ${\zeta }_1$ and vertices of ${\zeta }_2$

$$\begin{array}{cc}\hfill \sum _{\begin{array}{c}uv\in E_7\\ u\in V\left({\zeta }_1\right),u\in A\\ v\in V\left({\zeta }_2\right),v\in B\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{\begin{array}{c}uv\in E_7\\ u\in V\left({\zeta }_1\right),u\in A\\ v\in V\left({\zeta }_2\right),v\in B\end{array}}{\left(1-1\right)}^2=0.\hfill \\ \hfill \sum _{\begin{array}{c}uv\in E_8\\ u\in V\left({\zeta }_1\right),u\in A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{\begin{array}{c}uv\in E_8\\ u\in V\left({\zeta }_1\right),u\in A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}{(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right))}^2+{\left(n_1-\left|A\right|\right)}^2\left|E_8\right|-\hfill \\ \hfill & 2\left(n_1-\left|A\right|\right)\sum _{\begin{array}{c}uv\in E_8\\ u\in V\left({\zeta }_1\right),u\in A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)).\hfill \\ \hfill \sum _{\begin{array}{c}uv\in E_9\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\in B\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{\begin{array}{c}uv\in E_9\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\in B\end{array}}{(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right))}^2+{\left(n_2-\left|B\right|\right)}^2\left|E_9\right|+\hfill \\ \hfill & 2\left(n_2-\left|B\right|\right)\sum _{\begin{array}{c}uv\in E_9\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\in B\end{array}}(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)).\hfill \\ \hfill \sum _{\begin{array}{c}uv\in E_{10}\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2& =\sum _{\begin{array}{c}uv\in E_{10}\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}{[(d_{p_{{\zeta }_1}}\left(u\right)+\left(n_2-\left|B\right|\right))-(d_{p_{{\zeta }_2}}\left(v\right)+\left(n_1-\left|A\right|\right))]}^2\hfill \\ \hfill & =\sum _{\begin{array}{c}uv\in E_{10}\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}\left[\begin{array}{c}{\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)}^2+\\ {\left(n_2-n_1+\left|A\right|-\left|B\right|\right)}^2+\\ 2\left(n_2-n_1+\left|A\right|-\left|B\right|\right)\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)\end{array}\right]\hfill \\ \hfill & =\sum _{\begin{array}{c}uv\in E_{10}\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}{\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)}^2+{\left(n_2-n_1+\left|A\right|-\left|B\right|\right)}^2\left|E_{10}\right|+\hfill \\ \hfill & 2\left(n_2-n_1+\left|A\right|-\left|B\right|\right)\sum _{\begin{array}{c}uv\in E_{10}\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right).\hfill \end{array}$$

Hence we have

$$\begin{array}{cc}\hfill \stackrel{\left(3\right)}{\overbrace{{\sum }_{\begin{array}{c}u\in V\left({\zeta }_1\right)\\ v\in V\left({\zeta }_2\right)\end{array}}{(d_{p_{{\zeta }_1+{\zeta }_2}}\left(u\right)-d_{p_{{\zeta }_1+{\zeta }_2}}\left(v\right))}^2}}& =\sum _{\begin{array}{c}uv\in E_8\\ u\in V\left({\zeta }_1\right),u\in A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}{(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right))}^2+{\left(n_1-\left|A\right|\right)}^2\left|E_8\right|-\hfill \\ \hfill & 2\left(n_1-\left|A\right|\right)\sum _{\begin{array}{c}uv\in E_8\\ u\in V\left({\zeta }_1\right),u\in A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right))+\hfill \\ \hfill & \sum _{\begin{array}{c}uv\in E_9\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\in B\end{array}}{(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right))}^2+{\left(n_2-\left|B\right|\right)}^2\left|E_9\right|+\hfill \\ \hfill & 2\left(n_2-\left|B\right|\right)\sum _{\begin{array}{c}uv\in E_9\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\in B\end{array}}(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right))+\hfill \\ \hfill & \sum _{\begin{array}{c}uv\in E_{10}\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}{\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right)}^2+{\left(n_2-n_1+\left|A\right|-\left|B\right|\right)}^2\left|E_{10}\right|+\hfill \\ \hfill & 2\left(n_2-n_1+\left|A\right|-\left|B\right|\right)\sum _{\begin{array}{c}uv\in E_{10}\\ u\in V\left({\zeta }_1\right),u\notin A\\ v\in V\left({\zeta }_2\right),v\notin B\end{array}}\left(d_{p_{{\zeta }_1}}\left(u\right)-d_{p_{{\zeta }_2}}\left(v\right)\right).\hfill \end{array}$$

Now adding (1),(2), and (3) we obtain our result. □

3.2. Bounds for Domination Sigma Index

To make accurate predictions using topological indices, it is crucial to comprehend the potential range of values that these indices can assume. Determining the bounds of topological indices is a significant task when it comes to predicting the properties of chemical compounds. The bounds refer to the minimum and maximum values that an index can attain. Being aware of the bounds of a topological index can assist researchers in identifying the range of values that are physically relevant and in gaining an understanding of the connection between the topological characteristics of a molecular graph and the properties of a chemical compound. In this section, we will present some lower and upper bounds for the Domination Sigma index.

Theorem 9.

If $\zeta =\left(V\left(\zeta \right),E\left(\zeta \right)\right)$ with order n and size m is not a $K-$domination regular graph then

$$n-2\sum _{u\in V\left(\zeta \right)}{d}_{p\left(u\right)}+D{M}_1\left(\zeta \right)\le D\sigma \left(\zeta \right).$$

Proof. 

We have $1\le d_{p\left(v\right)},$ then

$$\begin{array}{cc}\hfill {\left(1-d_{p\left(u_1\right)}\right)}^2& \le {\left(d_{p\left(v_1\right)}-d_{p\left(u_1\right)}\right)}^2\hfill \\ \hfill {\left(1-d_{p\left(u_2\right)}\right)}^2& \le {\left(d_{p\left(v_2\right)}-d_{p\left(u_2\right)}\right)}^2\hfill \\ \hfill & ⋮\hfill \\ \hfill {\left(1-d_{p\left(u_n\right)}\right)}^2& \le {\left(d_{p\left(v_n\right)}-d_{p\left(u_n\right)}\right)}^2\hfill \end{array}$$

In the left side we take the summation over all vertices of $\zeta ,$ and in the right side we take the summation over all edges of $\zeta .$ Hence

$$\begin{array}{cc}\hfill 1-2d_{p\left(u_n\right)}+d_{p\left(u_n\right)}^2& \le {\left(d_{p\left(v_n\right)}-d_{p\left(u_n\right)}\right)}^2\hfill \\ \hfill n-2\sum _{u\in V\left(\zeta \right)}{d}_{p\left(u\right)}+D{M}_1\left(\zeta \right)& \le D\sigma \left(\zeta \right).\hfill \end{array}$$

Which complete the proof. □

Theorem 10.

If $\zeta =\left(V\left(\zeta \right),E\left(\zeta \right)\right)$ with order n and size m is not a $K-$domination regular graph then

$$D\sigma \left(\zeta \right)<n{T}_m^2\left(\zeta \right)-2T_m\left(\zeta \right)\sum _{u\in V\left(\zeta \right)}{d}_{p\left(u\right)}+D{M}_1\left(\zeta \right).$$

Proof. 

We have $d_{p\left(u\right)}<T_m\left(\zeta \right),$ hence

$$\begin{array}{cc}\hfill {\left(d_{p\left(v_1\right)}-d_{p\left(u_1\right)}\right)}^2& <{\left(T_m\left(\zeta \right)-d_{p\left(u_1\right)}\right)}^2\hfill \\ \hfill {\left(d_{p\left(v_2\right)}-d_{p\left(u_2\right)}\right)}^2& <{\left(T_m\left(\zeta \right)-d_{p\left(u_2\right)}\right)}^2\hfill \\ \hfill {\left(d_{p\left(v_3\right)}-d_{p\left(u_3\right)}\right)}^2& <{\left(T_m\left(\zeta \right)-d_{p\left(u_3\right)}\right)}^2\hfill \\ \hfill & ⋮\hfill \\ \hfill {\left(d_{p\left(v_n\right)}-d_{p\left(u_n\right)}\right)}^2& <{\left(T_m\left(\zeta \right)-d_{p\left(u_n\right)}\right)}^2\hfill \end{array}$$

Similarly of Theorem 9 we get

$$\sum _{uv\in E\left(\zeta \right)}{\left(d_{p\left(v\right)}-d_{p\left(u\right)}\right)}^2<T_m^2\left(\zeta \right)-T_m\left(\zeta \right)d_{p\left(u\right)}+d_{p\left(u\right)}^2$$

Hence $D\sigma \left(\zeta \right)\le n{T}_m^2\left(\zeta \right)-2T_m\left(\zeta \right){\sum }_{u\in V\left(\zeta \right)}{d}_{p\left(u\right)}+D{M}_1\left(\zeta \right).$ □

Proposition 4.

If $\zeta =\left(V\left(\zeta \right),E\left(\zeta \right)\right)$ with order n and size m is not a $K-$domination regular graph then

• $n-2T_m\left(\zeta \right)\gamma \left(\zeta \right)+D{M}_1\left(\zeta \right)\le D\sigma \left(\zeta \right),$
• $D\sigma \left(\zeta \right)<n{T}_m^2\left(\zeta \right)-T_m^2\left(\zeta \right)\Gamma \left(\zeta \right)+D{M}_1\left(\zeta \right).$

Proof. 

We have $T_m\left(\zeta \right)\gamma \left(\zeta \right)\le {\sum }_{u\in V\left(\zeta \right)}{d}_{p\left(u\right)}\le T_m\left(\zeta \right)\Gamma \left(\zeta \right).$ Hence by Theorems 9 and 10 we obtain our results. □

4. Statistical Validity of Domination Indices

Physicochemical and biological properties of molecules are predicted and modeled through QSPR analysis. To extract maximum information from a data set, chemometrics is a powerful tool that incorporates statistical and mathematical methods. Chemical descriptors of a molecule’s chemical structure are used in QSPR to describe how physicochemical properties vary as a consequence of chemometric methods. Hence, calculated descriptors can replace expensive biological tests or experiments concerning a particular physicochemical property. In turn, these descriptors can be used to predict the properties of interest for upcocompounds. QSPR works by finding a quantitative relationship that can be utilized for the prediction of compounds’ properties, even those that can’t be measured. As a matter of fact, QSPR models are primarily affected by the molecular structure parameters employed. Alternative molecular descriptors have been developed that are derived solely from chemical structure information. Researchers have focused much attention on connecting and composing molecules in order to determine "topological indices" that can be used in QSPR analyses. In addition to having the advantage of simplicity, the topological index is also fast to calculate, so scientists are interested in using it. A compound’s physicochemical and biological properties are influenced by its molecular structure, according to several chemical and medical experiments. Quantitative structure property/activity relationship $(QSAR/QSPR)$ models are generated by employing mathematical/statistical tools to determine this dependence. Regression models are used for relating physicochemical/biological properties to molecular descriptors. The graph-theoretic topological indices generate $QSAR/QSPR$ models by converting compounds into chemical graphs. To be accepted by Milan Randic [46], topological indices must meet certain criteria, the most significant of which is to be positively correlated with at least one physicochemical property. The purpose of this section is to investigate the significance of these newly developed domination topological indices. Octane and some of its isomers are described in Table 1 according to their experimental data [47], and also https://pubchem.ncbi.nlm.nih.gov (accessed on 26 March 2022). Computed domination indices values are shown in Table 2. Our analysis has shown that these indices play role in evaluating Entropy (E), Acentric factor (AF), Enthalpy of vaporization (HVAP), and Standard enthalpy of vaporization (DHVAP). A correlation coefficient $\left(R\right)$ between these indices and some physicochemical properties can be seen in Table 3.

4.1. Regression Model

QSPR analysis of domination topological indices will be discussed in this subsection. Furthermore, we demonstrate a positive correlation between the characteristics and the physicochemical characteristics of octane and its isomers. Here we discuss how topological indices can be used to predict physicochemical properties. We calculated six domination topological indices and one physicochemical property using R-software. Based on the below nonlinear regression model, we can derive different nonlinear models for the topological indices of domination:

$$ln\left(CP\right)=A+Bln\left(DI\right)$$

(9)

For the domination Sigma index we use the linear regression model because we have some values of zeros

$$\left(CP\right)=A+B\left(DI\right)$$

(10)

where $\left(CP\right)$ is the physical and chemical properties of octanes and its isomers and $\left(DI\right)$ represents the domination topological indices. By using Equations 9 and 10, we can obtain different linear and non-linear models for the domination topological indices as follows:

• The domination first Zagreb index $D{M}_1\left(\zeta \right):$

  $$\begin{array}{cc}\hfill ln\left(AF\right)& =-1.6+0.12ln\left(D{M}_1\left(\zeta \right)\right)\hfill \end{array}$$

  (11)

  $$\begin{array}{cc}\hfill ln\left(S\right)& =4.4+0.04ln\left(D{M}_1\left(\zeta \right)\right)\hfill \end{array}$$

  (12)

  $$\begin{array}{cc}\hfill ln\left(HVAP\right)& =4.05+0.03ln\left(D{M}_1\left(\zeta \right)\right)\hfill \end{array}$$

  (13)

  $$\begin{array}{cc}\hfill ln\left(DHVAP\right)& =1.9+0.05ln\left(D{M}_1\left(\zeta \right)\right)\hfill \end{array}$$

  (14)
• The domination second Zagreb index $D{M}_2\left(\zeta \right):$

  $$\begin{array}{cc}\hfill ln\left(AF\right)& =-1.5+0.1ln\left(D{M}_2\left(\zeta \right)\right)\hfill \end{array}$$

  (15)

  $$\begin{array}{cc}\hfill ln\left(S\right)& =4.5+0.04ln\left(D{M}_2\left(\zeta \right)\right)\hfill \end{array}$$

  (16)

  $$\begin{array}{cc}\hfill ln\left(HVAP\right)& =4.07+0.03ln\left(D{M}_2\left(\zeta \right)\right)\hfill \end{array}$$

  (17)

  $$\begin{array}{cc}\hfill ln\left(DHVAP\right)& =1.9+0.05ln\left(D{M}_2\left(\zeta \right)\right)\hfill \end{array}$$

  (18)
• The domination modified first Zagreb index $D{M}_1^{*}\left(\zeta \right):$

  $$\begin{array}{cc}\hfill ln\left(AF\right)& =-1.9+0.2ln\left(D{M}_1^{*}\left(\zeta \right)\right)\hfill \end{array}$$

  (19)

  $$\begin{array}{cc}\hfill ln\left(S\right)& =4.3+0.08ln\left(D{M}_1^{*}\left(\zeta \right)\right)\hfill \end{array}$$

  (20)

  $$\begin{array}{cc}\hfill ln\left(HVAP\right)& =3.9+0.07ln\left(D{M}_1^{*}\left(\zeta \right)\right)\hfill \end{array}$$

  (21)

  $$\begin{array}{cc}\hfill ln\left(DHVAP\right)& =1.7+0.1ln\left(D{M}_1^{*}\left(\zeta \right)\right)\hfill \end{array}$$

  (22)
• The domination forgotten index $DF\left(\zeta \right):$

  $$\begin{array}{cc}\hfill ln\left(AF\right)& =-1.5+0.08ln\left(DF\left(\zeta \right)\right)\hfill \end{array}$$

  (23)

  $$\begin{array}{cc}\hfill ln\left(S\right)& =4.43+0.4ln\left(DF\left(\zeta \right)\right)\hfill \end{array}$$

  (24)

  $$\begin{array}{cc}\hfill ln\left(HVAP\right)& =4.08+0.02ln\left(DF\left(\zeta \right)\right)\hfill \end{array}$$

  (25)

  $$\begin{array}{cc}\hfill ln\left(DHVAP\right)& =1.9+0.03ln\left(DF\left(\zeta \right)\right)\hfill \end{array}$$

  (26)
• The domination hyper index $DH\left(\zeta \right):$

  $$\begin{array}{cc}\hfill ln\left(AF\right)& =-1.7+0.11ln\left(DH\left(\zeta \right)\right)\hfill \end{array}$$

  (27)

  $$\begin{array}{cc}\hfill ln\left(S\right)& =4.39+0.04ln\left(DH\left(\zeta \right)\right)\hfill \end{array}$$

  (28)

  $$\begin{array}{cc}\hfill ln\left(HVAP\right)& =4.02+0.03ln\left(DH\left(\zeta \right)\right)\hfill \end{array}$$

  (29)

  $$\begin{array}{cc}\hfill ln\left(DHVAP\right)& =1.8+0.05ln\left(DH\left(\zeta \right)\right)\hfill \end{array}$$

  (30)
• The domination modified forgotten index $D{F}^{*}\left(\zeta \right):$

  $$\begin{array}{cc}\hfill ln\left(AF\right)& =-1.6+0.1ln\left(D{F}^{*}\left(\zeta \right)\right)\hfill \end{array}$$

  (31)

  $$\begin{array}{cc}\hfill ln\left(S\right)& =4.439+0.04ln\left(D{F}^{*}\left(\zeta \right)\right)\hfill \end{array}$$

  (32)

  $$\begin{array}{cc}\hfill ln\left(HVAP\right)& =4.05+0.03ln\left(D{F}^{*}\left(\zeta \right)\right)\hfill \end{array}$$

  (33)

  $$\begin{array}{cc}\hfill ln\left(DHVAP\right)& =1.9+0.05ln\left(D{F}^{*}\left(\zeta \right)\right)\hfill \end{array}$$

  (34)
• The domination Sigma index $D\sigma \left(\zeta \right):$

  $$\begin{array}{cc}\hfill \left(AF\right)& =0.3+0.003\left(D\sigma \left(\zeta \right)\right)\hfill \end{array}$$

  (35)

  $$\begin{array}{cc}\hfill \left(S\right)& =102.75+0.4\left(D\sigma \left(\zeta \right)\right)\hfill \end{array}$$

  (36)

  $$\begin{array}{cc}\hfill \left(HVAP\right)& =68.03+0.16\left(D\sigma \left(\zeta \right)\right)\hfill \end{array}$$

  (37)

  $$\begin{array}{cc}\hfill \left(DHVAP\right)& =8.9+0.03\left(D\sigma \left(\zeta \right)\right)\hfill \end{array}$$

  (38)

Now, the predicted values of physicochemical properties are given in Table 4, Table 5, Table 6 and Table 7. The experimental values vs the predicted values of all properties are shown in Figure 4, Figure 5, Figure 6 and Figure 7.

Table 4. The (AF) values predicted by domination topological indices

octane and its isomers	$\mathit{D}{\mathit{M}}_1$	$\mathit{D}{\mathit{M}}_2$	$\mathit{D}{\mathit{M}}_1^{*}$	$\mathit{D}\mathit{F}$	$\mathit{D}\mathit{H}$	$\mathit{D}{\mathit{F}}^{*}$
n-octane	$0.3990$	$0.3862$	$0.3611$	$0.4064$	$0.38918$	$0.37479$
2-Methyl-heptane	$0.3545$	$0.3491$	$0.3270$	$0.3617$	$0.34927$	$0.35391$
3-Methyl-heptane	$0.3956$	$0.3809$	$0.3566$	$0.4036$	$0.38389$	$0.37066$
4-Methyl-heptane	$0.3525$	$0.3458$	$0.3244$	$0.3603$	$0.34642$	$0.33859$
3-Ethyl-hexane	$0.3478$	$0.3454$	$0.3230$	$0.3546$	$0.34444$	$0.33559$
2,2-Dimethyl-hexane	$0.3487$	$0.3408$	$0.3217$	$0.3561$	$0.34262$	$0.33186$
2,3-Dimethyl-hexane	$0.3572$	$0.3522$	$0.3296$	$0.3641$	$0.35233$	$0.34295$
2,4-Dimethyl-hexane	$0.3704$	$0.3636$	$0.3403$	$0.3777$	$0.36475$	$0.35378$
2,5-Dimethyl-hexane	$0.3509$	$0.3436$	$0.3217$	$0.3576$	$0.34262$	$0.33410$
3,3-Dimethyl-hexane	$0.3572$	$0.3522$	$0.3296$	$0.3641$	$0.35233$	$0.34295$
3-Ethyl-2-methyl-pentane	$0.3450$	$0.3286$	$0.3795$	$0.3588$	$0.32913$	$0.32449$
2,2,3-Trimethyl-pentane	$0.3614$	$0.3577$	$0.3346$	$0.3675$	$0.35776$	$0.35152$
2,2,4-Trimethyl-pentane	$0.3024$	$0.3066$	$0.2869$	$0.3083$	$0.30251$	$0.29856$
2,3,3-Trimethyl-pentane	$0.3614$	$0.35767$	$0.3346$	$0.3675$	$0.35755$	$0.35152$
2,3,4-Trimethyl-pentane	$0.3614$	$0.3577$	$0.3346$	$0.3675$	$0.35755$	$0.35152$
2,2,3,3-Tetramethylbutane	$0.3060$	$0.3114$	$0.2913$	$0.3112$	$0.30698$	$0.30196$

Table 4. The (AF) values predicted by domination topological indices

octane and its isomers	$\mathit{D}{\mathit{M}}_1$	$\mathit{D}{\mathit{M}}_2$	$\mathit{D}{\mathit{M}}_1^{*}$	$\mathit{D}\mathit{F}$	$\mathit{D}\mathit{H}$	$\mathit{D}{\mathit{F}}^{*}$
n-octane	$0.3990$	$0.3862$	$0.3611$	$0.4064$	$0.38918$	$0.37479$
2-Methyl-heptane	$0.3545$	$0.3491$	$0.3270$	$0.3617$	$0.34927$	$0.35391$
3-Methyl-heptane	$0.3956$	$0.3809$	$0.3566$	$0.4036$	$0.38389$	$0.37066$
4-Methyl-heptane	$0.3525$	$0.3458$	$0.3244$	$0.3603$	$0.34642$	$0.33859$
3-Ethyl-hexane	$0.3478$	$0.3454$	$0.3230$	$0.3546$	$0.34444$	$0.33559$
2,2-Dimethyl-hexane	$0.3487$	$0.3408$	$0.3217$	$0.3561$	$0.34262$	$0.33186$
2,3-Dimethyl-hexane	$0.3572$	$0.3522$	$0.3296$	$0.3641$	$0.35233$	$0.34295$
2,4-Dimethyl-hexane	$0.3704$	$0.3636$	$0.3403$	$0.3777$	$0.36475$	$0.35378$
2,5-Dimethyl-hexane	$0.3509$	$0.3436$	$0.3217$	$0.3576$	$0.34262$	$0.33410$
3,3-Dimethyl-hexane	$0.3572$	$0.3522$	$0.3296$	$0.3641$	$0.35233$	$0.34295$
3-Ethyl-2-methyl-pentane	$0.3450$	$0.3286$	$0.3795$	$0.3588$	$0.32913$	$0.32449$
2,2,3-Trimethyl-pentane	$0.3614$	$0.3577$	$0.3346$	$0.3675$	$0.35776$	$0.35152$
2,2,4-Trimethyl-pentane	$0.3024$	$0.3066$	$0.2869$	$0.3083$	$0.30251$	$0.29856$
2,3,3-Trimethyl-pentane	$0.3614$	$0.35767$	$0.3346$	$0.3675$	$0.35755$	$0.35152$
2,3,4-Trimethyl-pentane	$0.3614$	$0.3577$	$0.3346$	$0.3675$	$0.35755$	$0.35152$
2,2,3,3-Tetramethylbutane	$0.3060$	$0.3114$	$0.2913$	$0.3112$	$0.30698$	$0.30196$

Table 5. The (E) values predicted by domination topological indices

octane and its isomers	$\mathit{D}{\mathit{M}}_1$	$\mathit{D}{\mathit{M}}_2$	$\mathit{D}{\mathit{M}}_1^{*}$	$\mathit{D}\mathit{F}$	$\mathit{D}\mathit{H}$	$\mathit{D}{\mathit{F}}^{*}$
n-octane	$102.214$	$112.1$	$104.85$	$108.29$	$106.16$	$996.73$
2-Methyl-heptane	$98.26$	$107.67$	$100.78$	$103.66$	$102.07$	$792.54$
3-Methyl-heptane	$101.92$	$111.48$	$104.32$	$108.009$	$105.64$	$953.49$
4-Methyl-heptane	$98.079$	$107.26$	$100.45$	$103.51$	$101.76$	$663.94$
3-Ethyl-hexane	$97.64$	$107.20$	$100.28$	$102.89$	$101.55$	$640.7$
2,2-Dimethyl-hexane	$97.72$	$106.62$	$100.11$	$103.05$	$101.36$	$612.73$
2,3-Dimethyl-hexane	$98.508$	$108.04$	$101.09$	$103.92$	$102.39$	$698.77$
2,4-Dimethyl-hexane	$99.708$	$109.43$	$102.39$	$105.35$	$103.69$	$791.38$
2,5-Dimethyl-hexane	$97.925$	$106.98$	$100.11$	$103.22$	$101.36$	$629.40$
3,3-Dimethyl-hexane	$98.508$	$108.04$	$101.09$	$103.92$	$102.39$	$698.77$
3-Ethyl-2-methyl-pentane	$97.381$	$105.09$	$98.38$	$102.81$	$99.89$	$560.02$
2,2,3-Trimethyl-pentane	$98.89$	$108.71$	$101.7$	$104.28$	$102.94$	$771.2$
2,2,4-Trimethyl-pentane	$93.19$	$102.21$	$95.64$	$97.64$	$96.87$	$401.34$
2,3,3-Trimethyl-pentane	$98.89$	$108.71$	$101.7$	$104.28$	$102.94$	$771.29$
2,3,4-Trimethyl-pentane	$98.89$	$108.71$	$101.7$	$104.28$	$102.94$	$771.29$
2,2,3,3-Tetramethylbutane	$93.56$	$102.85$	$96.2$	$97.98$	$97.39$	$419.95$

Table 5. The (E) values predicted by domination topological indices

octane and its isomers	$\mathit{D}{\mathit{M}}_1$	$\mathit{D}{\mathit{M}}_2$	$\mathit{D}{\mathit{M}}_1^{*}$	$\mathit{D}\mathit{F}$	$\mathit{D}\mathit{H}$	$\mathit{D}{\mathit{F}}^{*}$
n-octane	$102.214$	$112.1$	$104.85$	$108.29$	$106.16$	$996.73$
2-Methyl-heptane	$98.26$	$107.67$	$100.78$	$103.66$	$102.07$	$792.54$
3-Methyl-heptane	$101.92$	$111.48$	$104.32$	$108.009$	$105.64$	$953.49$
4-Methyl-heptane	$98.079$	$107.26$	$100.45$	$103.51$	$101.76$	$663.94$
3-Ethyl-hexane	$97.64$	$107.20$	$100.28$	$102.89$	$101.55$	$640.7$
2,2-Dimethyl-hexane	$97.72$	$106.62$	$100.11$	$103.05$	$101.36$	$612.73$
2,3-Dimethyl-hexane	$98.508$	$108.04$	$101.09$	$103.92$	$102.39$	$698.77$
2,4-Dimethyl-hexane	$99.708$	$109.43$	$102.39$	$105.35$	$103.69$	$791.38$
2,5-Dimethyl-hexane	$97.925$	$106.98$	$100.11$	$103.22$	$101.36$	$629.40$
3,3-Dimethyl-hexane	$98.508$	$108.04$	$101.09$	$103.92$	$102.39$	$698.77$
3-Ethyl-2-methyl-pentane	$97.381$	$105.09$	$98.38$	$102.81$	$99.89$	$560.02$
2,2,3-Trimethyl-pentane	$98.89$	$108.71$	$101.7$	$104.28$	$102.94$	$771.2$
2,2,4-Trimethyl-pentane	$93.19$	$102.21$	$95.64$	$97.64$	$96.87$	$401.34$
2,3,3-Trimethyl-pentane	$98.89$	$108.71$	$101.7$	$104.28$	$102.94$	$771.29$
2,3,4-Trimethyl-pentane	$98.89$	$108.71$	$101.7$	$104.28$	$102.94$	$771.29$
2,2,3,3-Tetramethylbutane	$93.56$	$102.85$	$96.2$	$97.98$	$97.39$	$419.95$

Table 6. The (HVAP) values predicted by domination topological indices

octane and its isomers	$\mathit{D}{\mathit{M}}_1$	$\mathit{D}{\mathit{M}}_2$	$\mathit{D}{\mathit{M}}_1^{*}$	$\mathit{D}\mathit{F}$	$\mathit{D}\mathit{H}$	$\mathit{D}{\mathit{F}}^{*}$
n-octane	$68.05$	$69.03$	$67.25$	$68.7$	$68.46$	$69.1$
2-Methyl-heptane	$66.07$	$66.97$	$64.96$	$66.73$	$66.46$	$67.92$
3-Methyl-heptane	$67.90$	$68.74$	$66.95$	$68.58$	$68.20$	$68.87$
4-Methyl-heptane	$65.97$	$66.78$	$64.77$	$66.67$	$66.32$	$67.02$
3-Ethyl-hexane	$65.75$	$66.75$	$64.68$	$66.40$	$66.21$	$66.84$
2,2-Dimethyl-hexane	$65.79$	$66.48$	$64.58$	$66.47$	$66.12$	$66.6$
2,3-Dimethyl-hexane	$66.19$	$67.15$	$65.14$	$66.84$	$66.62$	$67.28$
2,4-Dimethyl-hexane	$66.79$	$67.79$	$65.87$	$67.46$	$67.2$	$67.9$
2,5-Dimethyl-hexane	$65.9$	$66.65$	$64.58$	$66.54$	$66.12$	$66.76$
3,3-Dimethyl-hexane	$66.19$	$67.15$	$65.1$	$66.84$	$66.62$	$67.28$
3-Ethyl-2-methyl-pentane	$65.62$	$65.76$	$63.6$	$66.37$	$65.4$	$66.17$
2,2,3-Trimethyl-pentane	$66.39$	$67.46$	$65.48$	$67.004$	$66.89$	$67.78$
2,2,4-Trimethyl-pentane	$63.49$	$64.41$	$62.05$	$64.12$	$63.9$	$64.54$
2,3,3-Trimethyl-pentane	$66.39$	$67.46$	$65.48$	$67.004$	$66.89$	$67.78$
2,3,4-Trimethyl-pentane	$66.39$	$67.46$	$65.48$	$67.004$	$66.89$	$67.78$
2,2,3,3-Tetramethylbutane	$63.68$	$64.71$	$62.38$	$64.27$	$64.17$	$64.76$

Table 6. The (HVAP) values predicted by domination topological indices

octane and its isomers	$\mathit{D}{\mathit{M}}_1$	$\mathit{D}{\mathit{M}}_2$	$\mathit{D}{\mathit{M}}_1^{*}$	$\mathit{D}\mathit{F}$	$\mathit{D}\mathit{H}$	$\mathit{D}{\mathit{F}}^{*}$
n-octane	$68.05$	$69.03$	$67.25$	$68.7$	$68.46$	$69.1$
2-Methyl-heptane	$66.07$	$66.97$	$64.96$	$66.73$	$66.46$	$67.92$
3-Methyl-heptane	$67.90$	$68.74$	$66.95$	$68.58$	$68.20$	$68.87$
4-Methyl-heptane	$65.97$	$66.78$	$64.77$	$66.67$	$66.32$	$67.02$
3-Ethyl-hexane	$65.75$	$66.75$	$64.68$	$66.40$	$66.21$	$66.84$
2,2-Dimethyl-hexane	$65.79$	$66.48$	$64.58$	$66.47$	$66.12$	$66.6$
2,3-Dimethyl-hexane	$66.19$	$67.15$	$65.14$	$66.84$	$66.62$	$67.28$
2,4-Dimethyl-hexane	$66.79$	$67.79$	$65.87$	$67.46$	$67.2$	$67.9$
2,5-Dimethyl-hexane	$65.9$	$66.65$	$64.58$	$66.54$	$66.12$	$66.76$
3,3-Dimethyl-hexane	$66.19$	$67.15$	$65.1$	$66.84$	$66.62$	$67.28$
3-Ethyl-2-methyl-pentane	$65.62$	$65.76$	$63.6$	$66.37$	$65.4$	$66.17$
2,2,3-Trimethyl-pentane	$66.39$	$67.46$	$65.48$	$67.004$	$66.89$	$67.78$
2,2,4-Trimethyl-pentane	$63.49$	$64.41$	$62.05$	$64.12$	$63.9$	$64.54$
2,3,3-Trimethyl-pentane	$66.39$	$67.46$	$65.48$	$67.004$	$66.89$	$67.78$
2,3,4-Trimethyl-pentane	$66.39$	$67.46$	$65.48$	$67.004$	$66.89$	$67.78$
2,2,3,3-Tetramethylbutane	$63.68$	$64.71$	$62.38$	$64.27$	$64.17$	$64.76$

Table 7. The (DHVAP) values predicted by domination topological indices

octane and its isomers	$\mathit{D}{\mathit{M}}_1$	$\mathit{D}{\mathit{M}}_2$	$\mathit{D}{\mathit{M}}_1^{*}$	$\mathit{D}\mathit{F}$	$\mathit{D}\mathit{H}$	$\mathit{D}{\mathit{F}}^{*}$
n-octane	$8.88$	$8.79$	$8.5$	$8.37$	$8.53$	$9.10$
2-Methyl-heptane	$8.45$	$8.36$	$8.09$	$8.01$	$8.12$	$8.85$
3-Methyl-heptane	$8.84$	$8.73$	$8.45$	$8.34$	$8.47$	$9.05$
4-Methyl-heptane	$8.43$	$8.32$	$8.06$	$8.002$	$8.09$	$8.65$
3-Ethyl-hexane	$8.38$	$8.31$	$8.04$	$7.95$	$8.07$	$8.61$
2,2-Dimethyl-hexane	$8.39$	$8.26$	$8.02$	$7.96$	$8.05$	$8.57$
2,3-Dimethyl-hexane	$8.47$	$8.39$	$8.12$	$8.03$	$8.15$	$8.71$
2,4-Dimethyl-hexane	$8.6$	$8.53$	$8.25$	$8.14$	$8.28$	$8.85$
2,5-Dimethyl-hexane	$8.41$	$8.29$	$8.02$	$7.97$	$8.05$	$8.60$
3,3-Dimethyl-hexane	$8.47$	$8.39$	$8.12$	$8.03$	$8.15$	$8.71$
3-Ethyl-2-methyl-pentane	$8.35$	$8.11$	$7.85$	$7.94$	$7.01$	$8.47$
2,2,3-Trimethyl-pentane	$8.52$	$8.46$	$8.18$	$8.06$	$8.23$	$8.82$
2,2,4-Trimethyl-pentane	$7.91$	$7.83$	$7.58$	$7.54$	$7.61$	$8.13$
2,3,3-Trimethyl-pentane	$8.5$	$8.46$	$8.18$	$8.06$	$8.23$	$8.82$
2,3,4-Trimethyl-pentane	$8.5$	$8.46$	$8.18$	$8.06$	$8.23$	$8.82$
2,2,3,3-Tetramethylbutane	$7.95$	$7.89$	$7.63$	$7.54$	$7.65$	$8.17$

Table 7. The (DHVAP) values predicted by domination topological indices

octane and its isomers	$\mathit{D}{\mathit{M}}_1$	$\mathit{D}{\mathit{M}}_2$	$\mathit{D}{\mathit{M}}_1^{*}$	$\mathit{D}\mathit{F}$	$\mathit{D}\mathit{H}$	$\mathit{D}{\mathit{F}}^{*}$
n-octane	$8.88$	$8.79$	$8.5$	$8.37$	$8.53$	$9.10$
2-Methyl-heptane	$8.45$	$8.36$	$8.09$	$8.01$	$8.12$	$8.85$
3-Methyl-heptane	$8.84$	$8.73$	$8.45$	$8.34$	$8.47$	$9.05$
4-Methyl-heptane	$8.43$	$8.32$	$8.06$	$8.002$	$8.09$	$8.65$
3-Ethyl-hexane	$8.38$	$8.31$	$8.04$	$7.95$	$8.07$	$8.61$
2,2-Dimethyl-hexane	$8.39$	$8.26$	$8.02$	$7.96$	$8.05$	$8.57$
2,3-Dimethyl-hexane	$8.47$	$8.39$	$8.12$	$8.03$	$8.15$	$8.71$
2,4-Dimethyl-hexane	$8.6$	$8.53$	$8.25$	$8.14$	$8.28$	$8.85$
2,5-Dimethyl-hexane	$8.41$	$8.29$	$8.02$	$7.97$	$8.05$	$8.60$
3,3-Dimethyl-hexane	$8.47$	$8.39$	$8.12$	$8.03$	$8.15$	$8.71$
3-Ethyl-2-methyl-pentane	$8.35$	$8.11$	$7.85$	$7.94$	$7.01$	$8.47$
2,2,3-Trimethyl-pentane	$8.52$	$8.46$	$8.18$	$8.06$	$8.23$	$8.82$
2,2,4-Trimethyl-pentane	$7.91$	$7.83$	$7.58$	$7.54$	$7.61$	$8.13$
2,3,3-Trimethyl-pentane	$8.5$	$8.46$	$8.18$	$8.06$	$8.23$	$8.82$
2,3,4-Trimethyl-pentane	$8.5$	$8.46$	$8.18$	$8.06$	$8.23$	$8.82$
2,2,3,3-Tetramethylbutane	$7.95$	$7.89$	$7.63$	$7.54$	$7.65$	$8.17$

Table 8. Statical parameters for the non-linear QSPR model for first Zagreb domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.549$	$0.3040$	$9.75$	$0.029$	$0.00750$
Entropy $\left(E\right)$	$0.541$	$0.4141$	$9.51$	$3.745$	$0.0081$
Enthalpy of vaporization $\left(HVAP\right)$	$0.690$	$0.3236$	$16.01$	$1.512$	$0.0013$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.757$	$0.2774$	$20.31$	$0.271$	$0.0005$

Table 8. Statical parameters for the non-linear QSPR model for first Zagreb domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.549$	$0.3040$	$9.75$	$0.029$	$0.00750$
Entropy $\left(E\right)$	$0.541$	$0.4141$	$9.51$	$3.745$	$0.0081$
Enthalpy of vaporization $\left(HVAP\right)$	$0.690$	$0.3236$	$16.01$	$1.512$	$0.0013$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.757$	$0.2774$	$20.31$	$0.271$	$0.0005$

Table 9. Statical parameters for the non-linear QSPR model for the second domination Zagreb index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.465$	$0.3827$	$7.17$	$0.031$	$0.0180$
Entropy $\left(E\right)$	$0.461$	$0.3911$	$7.07$	$3.955$	$0.0187$
Enthalpy of vaporization $\left(HVAP\right)$	$0.614$	$0.3418$	$12.28$	$1.616$	$0.0035$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.659$	$0.2994$	$14.46$	$0.294$	$0.0019$

Table 9. Statical parameters for the non-linear QSPR model for the second domination Zagreb index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.465$	$0.3827$	$7.17$	$0.031$	$0.0180$
Entropy $\left(E\right)$	$0.461$	$0.3911$	$7.07$	$3.955$	$0.0187$
Enthalpy of vaporization $\left(HVAP\right)$	$0.614$	$0.3418$	$12.28$	$1.616$	$0.0035$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.659$	$0.2994$	$14.46$	$0.294$	$0.0019$

Table 10. Statical parameters for the linear QSPR model for modified zagreb domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.476$	$0.3794$	$4.10$	$0.034$	$0.0423$
Entropy $\left(E\right)$	$0.465$	$0.4188$	$7.19$	$3.944$	$0.0179$
Enthalpy of vaporization $\left(HVAP\right)$	$0.615$	$0.3041$	$12.31$	$1.615$	$0.0035$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.659$	$0.3094$	$14.85$	$0.292$	$0.0018$

Table 10. Statical parameters for the linear QSPR model for modified zagreb domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.476$	$0.3794$	$4.10$	$0.034$	$0.0423$
Entropy $\left(E\right)$	$0.465$	$0.4188$	$7.19$	$3.944$	$0.0179$
Enthalpy of vaporization $\left(HVAP\right)$	$0.615$	$0.3041$	$12.31$	$1.615$	$0.0035$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.659$	$0.3094$	$14.85$	$0.292$	$0.0018$

Table 11. Statical parameters for the non-linear QSPR model for forgotten domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.569$	$0.2949$	$10.45$	$0.029$	$0.0060$
Entropy $\left(E\right)$	$0.561$	$0.3445$	$10.20$	$3.690$	$0.0065$
Enthalpy of vaporization $\left(HVAP\right)$	$0.704$	$0.3236$	$16.82$	$1.492$	$0.0011$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.694$	$0.3127$	$21.08$	$0.265$	$0.0004$

Table 11. Statical parameters for the non-linear QSPR model for forgotten domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.569$	$0.2949$	$10.45$	$0.029$	$0.0060$
Entropy $\left(E\right)$	$0.561$	$0.3445$	$10.20$	$3.690$	$0.0065$
Enthalpy of vaporization $\left(HVAP\right)$	$0.704$	$0.3236$	$16.82$	$1.492$	$0.0011$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.694$	$0.3127$	$21.08$	$0.265$	$0.0004$

Table 12. Statical parameters for the linear QSPR model for hyper domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.489$	$0.3480$	$7.80$	$0.031$	$0.0144$
Entropy $\left(E\right)$	$0.482$	$0.4127$	$7.63$	$3.903$	$0.0153$
Enthalpy of vaporization $\left(HVAP\right)$	$0.635$	$0.3434$	$13.24$	$1.587$	$0.0027$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.691$	$0.2654$	$16.06$	$0.377$	$0.0012$

Table 12. Statical parameters for the linear QSPR model for hyper domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.489$	$0.3480$	$7.80$	$0.031$	$0.0144$
Entropy $\left(E\right)$	$0.482$	$0.4127$	$7.63$	$3.903$	$0.0153$
Enthalpy of vaporization $\left(HVAP\right)$	$0.635$	$0.3434$	$13.24$	$1.587$	$0.0027$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.691$	$0.2654$	$16.06$	$0.377$	$0.0012$

Table 13. Statical parameters for the non-linear QSPR model for modified forgotten domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.486$	$0.3819$	$7.74$	$0.031$	$0.0147$
Entropy $\left(E\right)$	$0.476$	$0.0063$	$7.49$	$3.917$	$0.0161$
Enthalpy of vaporization $\left(HVAP\right)$	$0.618$	$0.3343$	$12.47$	$1.610$	$0.0033$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.681$	$0.2799$	$15.54$	$0.288$	$0.0015$

Table 13. Statical parameters for the non-linear QSPR model for modified forgotten domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.486$	$0.3819$	$7.74$	$0.031$	$0.0147$
Entropy $\left(E\right)$	$0.476$	$0.0063$	$7.49$	$3.917$	$0.0161$
Enthalpy of vaporization $\left(HVAP\right)$	$0.618$	$0.3343$	$12.47$	$1.610$	$0.0033$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.681$	$0.2799$	$15.54$	$0.288$	$0.0015$

Table 14. Statical parameters for the linear QSPR model for sigma domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.320$	$0.033$	$5.111$	$0.0014$	$0.0402$
Entropy $\left(E\right)$	$0.267$	$4.002$	$0.1697$	$0.0012$	$0.0224$
Enthalpy of vaporization $\left(HVAP\right)$	$0.213$	$1.963$	$3.81$	$0.0833$	$0.0714$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.236$	$0.366$	$4.324$	$0.0155$	$0.0564$

Table 14. Statical parameters for the linear QSPR model for sigma domination index.

Physical Properties	${\mathit{R}}^2$	$\mathit{S}\mathit{e}$	F	$\mathit{S}\mathit{F}$	$\mathit{p}-\mathit{v}\mathit{a}\mathit{l}\mathit{u}\mathit{e}$
Acentric factor $\left(AF\right)$	$0.320$	$0.033$	$5.111$	$0.0014$	$0.0402$
Entropy $\left(E\right)$	$0.267$	$4.002$	$0.1697$	$0.0012$	$0.0224$
Enthalpy of vaporization $\left(HVAP\right)$	$0.213$	$1.963$	$3.81$	$0.0833$	$0.0714$
Standard enthalpy of vaporization $\left(DHVAP\right)$	$0.236$	$0.366$	$4.324$	$0.0155$	$0.0564$

4.2. Results and Discussion

To begin with, it was discovered see Figure 8, that any structure-property relationship could be achieved using the domination forgotten index $DF\left(\zeta \right)$. Based on Table 3 we can determine that the domination forgotten index $DF\left(\zeta \right)$ is the most appropriate index to model Standard enthalpy of vaporization $\left(DHVAP\right)$, Enthalpy of vaporization $\left(HVAP\right)$, Entropy $\left(E\right)$, and Acentric factor $\left(AF\right)$ with $R=0.879,0.839,0.749$, and $0.754$, respectively. We have found that our approach using the forgotten domination index $DF\left(\zeta \right)$ has provided a significant improvement in predicting the physicochemical properties of octane and its isomers compared to some of the recent studies. For example, in a study by Xie et al. [48], they used machine learning techniques to predict the physicochemical properties of organic compounds, including octane and its isomers, and reported correlation coefficients between $0.648$ and $0.873$ for various properties. In contrast, our approach using the $DF\left(\zeta \right)$ index achieved higher correlation coefficients of $0.749$ to $0.879$ for the same properties. It has been demonstrated that the Standard Enthalpy of Vaporization $\left(DHVAP\right)$ with correlation coefficient range between $0.812-0.870$ (except the correlation coefficient for the dominion Sigma index $D\sigma \left(\zeta \right)$ which was $0.486$) is best physicochemical property predicted by the new domination indices. We can see from Table 3 that for $\left(E\right),\left(AF\right),\left(HVAP\right),$ and $\left(DHVAP\right)$ the domination first Zagreb index $D{M}_1\left(\zeta \right)$ gives the second highest correlation coefficients $R=0.736,0.741,0.830,$ and $0.870$ respectively. Furthermore, in a recent study by Zhang et al. [49], they used graph theory-based topological indices to predict the boiling points of organic compounds, including octane and its isomers, and reported correlation coefficients of $0.827$ to $0.873$. In our study, we have also used graph theory-based topological indices, and our approach using the domination first Zagreb index $D{M}_1\left(\zeta \right)$ achieved a correlation coefficient of $0.870$ for predicting the standard enthalpy of vaporization $\left(DHVAP\right)$ of octane and its isomers. Therefore, we believe that our work provides a significant contribution to the field of predicting the physicochemical properties of organic compounds and offers a promising alternative to the existing approaches. It is worth noting that our finding that the domination first Zagreb index $D{M}_1\left(\zeta \right)$ provides a high degree of correlation for predicting various physicochemical properties of octane and its isomers is consistent with several recent studies in the field. For example, a study by Ghorbani et al. [50] also found that the $M_1\left(\zeta \right)$ index was a reliable predictor of various physicochemical properties of hydrocarbons, including octane. Similarly, a study by Moosavi et al. [51] also reported that the $M_1\left(\zeta \right)$ index was a useful predictor of the thermodynamic properties of various hydrocarbons, including isomers of octane. Overall, our findings, along with those of other studies, highlight the potential usefulness of the $M_1\left(\zeta \right)$ index in predicting the physicochemical properties of octane and its isomers. However, it is important to note that the effectiveness of this index, like any other index or method, may depend on various factors such as the size and diversity of the dataset, the modeling techniques used, and the specific properties being predicted. Therefore, further investigation is needed to fully understand the potential and limitations of the $D{M}_1\left(\zeta \right)$ index and other related indices in the context of predicting physicochemical properties of octane and its isomers.

It is important to note that the other domination topological indices used in this study have demonstrated good effectiveness in predicting the properties of Octanes and its isomers. However, it is also essential to acknowledge that the domination Sigma index, while providing valuable insights into the complex relationship between molecular structure and physical properties, showed poor correlation in this particular study, since it does not allow modifying the values of it using the non-linear model. The limitations of newly proposed topological indices and the specific property being studied can affect its effectiveness. As such, this study highlights the need for further research to explore modifications to the domination Sigma index for better suitability or to investigate its potential in predicting other properties. Additionally, investigating the effectiveness of the newly proposed index for graphs that are not $K-$domination regular graph can be an exciting avenue for future research.

5. Conclusion

The concept of domination in graph theory is a fundamental problem that involves finding the smallest set of vertices in a graph that collectively cover all the other vertices. The domination Sigma index is a measure that quantifies this problem by assigning a numerical value to each graph. In this paper, the authors focus on computing the domination Sigma index for different types of graphs, such as book graphs and compositions of graphs. One of the key contributions of this paper is the introduction of the k-domination regular graph. This type of graph is particularly interesting because it provides a useful tool for studying domination theory. A $k-$domination regular graph is a graph where every vertex has the same number k of neighbors in each dominating set of the graph. The authors use this concept to compute the domination Sigma index for some special graph classes and also identify areas where further investigation is needed. The authors also apply their findings to the study of octane and its isomers using QSPR analysis. This analysis assesses the significance of domination indices and a newly invented topological index. However, the results for the domination Sigma index of octane and its isomers show a poor correlation coefficient for most molecular structures and properties under study. This finding highlights the need for further investigation into the newly invented domination index, which utilizes graphs that are not k-domination regular graphs. In conclusion, this paper makes a valuable contribution to the field of graph theory by exploring the domination Sigma index and $k-$domination regular graphs. The application of QSPR analysis to assess the significance of these indices also highlights the interdisciplinary nature of graph theory and its applications. Further investigation into the newly invented domination index and its application to real-world problems could potentially lead to new insights and solutions in various fields. As future work, researchers could investigate the domination Sigma index for other types of graphs and explore their properties such as Mycielski [52,53], Thorn [54], Sierpiński [55,56], Helm [57] graphs. This would allow for a deeper understanding of the index’s capabilities and limitations in various contexts. Additionally, exploring the properties of non $k-$domination regular graphs could also yield valuable insights into the index’s behavior. For example, investigating the relationship between the domination index and other graph properties, such as connectivity or edge density, could help determine the usefulness of the index in different graph-theoretic applications [58,59,60,61,62]. The interdisciplinary nature of graph theory and its applications presents numerous opportunities for future research. Finally, this paper lays a foundation for future researchers to explore the domination Sigma index and its applications in greater depth.