# On Topological Indices for Complex Indium Phosphate Network and Their Applications

Authors: Wang Hui; Lubna Sherin; Sana Javed; Sadia Khalid; Waqar Asghar; Samuel Asefa Fufa
Journal: Complexity
 (2022)

---

## Abstract

A chemical compound in the form of graph terminology is known as a chemical graph. Molecules are usually represented as vertices, while their bonding or interaction is shown by edges in a molecular graph. In this paper, we computed various connectivity indices based on degrees of vertices of a chemical graph of indium phosphide (InP). Afterward, we found the physical measures like entropy and heat of formation of InP. Then, we fitted curves between different indices and the thermodynamical properties, namely, heat of formation and entropy. Curve fitting was done in MATLAB through different methods based on linearity and nonlinearity. Furthermore, we depicted our results numerically and graphically. These numerical systems may give an approach to concentrate on the thermodynamical properties of the compound design of InP at an exceptional level that will help understand the connection between framework measurement and these actions.


---

## Body

## 1. Introduction

An application of graph theory in the field of chemistry is to study the molecular structures of chemical compounds. Graph theory tools are implemented to classify fundamental features entailed in structure-property activity interactions of molecules. The term used to illustrate a molecule/chemical compound in the form of a graph is known as amolecular/chemical graph [1, 2]. Molecules are shown as vertices, while their bondings or interactions are shown by edges in a molecular graph. Mostly, molecule graphs are simple graphs, and the measure of the topological index is invariant under graph isomorphisms. Mostly, the degree or distance measure is used to capture the topology of a graphical structure, so the most common indices are based on either degree or distance between the vertices. Indices comprising degree measurements perform a vigorous part in molecular graph theory. Two isomorphic graphs have the same connectivity index, and the cardinalities of vertex and edge sets of a graph are considered as the topological/connectivity indices as well. A connectivity index explains some helpful details about the structure and analysis of the molecular graph.Many chemical compounds have been analyzed through topological indices in the past few decades. Topological/graphical index is a numeric measure related to chemical compositions asserting the association of chemical structures through numerous physicochemical properties or chemical reactivity [1]. Gao et al. [2] are originators of the field of chemical graph theory. It was acknowledged in 1988 that a lot of researchers had worked in this field and published roughly 500 publications per year. Chemical graph theory, a two-volume comprehensive treatise by Trinajsti, that conveyed the research up to the mid of 1980s, is one of many monographs in the discipline [3].Chemical compounds are represented by graphs in chemical graph theory, and mathematical tools are employed to address chemistry issues [4]. A topological/connectivity descriptor is a numerical measure that describes the topology of a graph [5]. Connectivity descriptors have gained considerable popularity recently, due to their simple nature. Chemical graph theory relies considerably on these graph descriptors. As a consequence, a topological index may be quantitatively characterized a chemical network that is topologically invariant for labelling as well as distinguishing between isomers [6]. In theoretical chemistry and nanotechnology, there are several graphs associated with numerical descriptors that are significant. Degree-based, distance-based, and counting-related graph descriptors are among the most common types [7]. The degree-based graph descriptors have a prominent place among these descriptors and may be used to characterize chemical substances, either organic or inorganic, and forecast their specific physiochemical characteristics [8].In [9], the authors suggested a way to analyze the metal-organic chemical compounds from topological perspectives based on mathematical frameworks. Such frameworks were further utilized in [10] to identify the molecular hydrogen stages within few seconds. In [11], authors studied the metal-organic structure by using different topological descriptors including Zagreb indices and Zagreb coindices. A study on the molecular chain in dialkyltin of complex-2, complex-3, and complex-4 of N-salicylidene-L-valine was conducted in [12], and the authors also analyzed these structures by providing different examples and established a comparison among topological indices.Caporossi et al. [13] established a graph descriptor that is known as atom bond connectivity index and is denoted by(1)ABCG=∑uv∈ΕGςu+ςv−2ςu×ςv.Vukičević and Furtula [14] introduced the geometric arithmetic index as follows:(2)GAG=∑uv∈ΕG2ςu×ςvςu+ςv.The first two indices introduced in the field of chemical graph theory were Zagreb indices that were used to show the connection between chemical structures withπ electrons [3]. For graph G, Zagreb indices are defined as follows:(3)M1G=∑uv∈ΕGςu+ςv,M2G=∑uv∈ΕGςu×ςv.The first and second Zagreb coindices defined in [15] are as follows:(4)M1¯G=∑uv∉ΕGςu+ςv,(5)M2¯G=∑uv∉ΕGςu×ςv.The most seasoned topological index that was introduced by Bollobas and Erdos is the first Zagreb index in light of the degree of vertices of a graph in [16]. Afterward, Furtula and Gutman [17] presented a forgotten connectivity index that was characterized as follows:(6)FG=∑uv∈ΕGςu2+ςv2.Wang et al. [18], inspired by the ABC index, devised a new modification, which is known as augmented Zagreb index and described as follows:(7)AZIG=∑uv∈ΕGςu×ςvςu+ςv−23.To construct a thermodynamical structure, we need to measure some physical quantities, and entropy and heat of formation are two of them.Entropy measure tells us how much heat energy we need to produce more in order to perform some valued work. Since this measure is describing the lack of energy due to which performing valuable work is not possible, it is also termed as a measure of the disorder [19]. The entropy of an isolated system has the highest entropy, according to the second law of thermodynamics. Nonisolated systems can lose entropy if they enhance the entropy of their surroundings by at least the same amount. Because entropy is a state function, every process that moves a system from one state to another, whether reversible or irreversible, will change its entropy.During per unit formation, the heat absorbed or retained is referred to as theheat of formation provided all the elements persist in a normal state. Kilojoule per mole (kJ/mol) is the unit for the measure of the heat of formation. The term enthalpy is also used for the heat of formation. Defining a system in the form of a mathematical framework provides us with an efficient approach to analyze the dynamics of the system. Experimental work is most expensive and very time-consuming so transforming the system into a set of mathematical forms makes this study very coherent. Many software like MATLAB or Python are easily available that provide a very friendly environment to construct mathematical models and study them. As we may fit many mathematical models to the same set of data, it is difficult to decide which one is best suitable for us. There are several statistical tests that might help us decide which mathematical model or framework is the best fit for our data, but we will just consider MSE, RMSE, SSE, or R2.


## 2. Crystal Structure of Indium Phosphide (InP) [s, t]

Indium phosphide, a semiconductor, belongs to group IIIA-VA. Research of Vasil'Ev and Gachon [20] showed that on decomposing InP in bomb calorimeter, two forms of phosphorous (black and red) were found along with indium. Indium phosphide crystallizes in cubic form with lattice parameter value a:5.861. Figure 1 demonstrates the crystal structure of InP [21].Figure 1
Crystal structure of unit cell of indium phosphide: (a) front view and (b) general view.
(a)(b)InP with bandgap 1.35 eV value [22] having low toxic nature, perfect composition, shape, and size as well as solution processability, is being used commercially for many years. Quantum dots (QDs) have got extraordinary importance because of having narrow emission line width, highest photoluminescence stability, reproducibility, and yield [23]. Figure 2 represents the synthesis method of InP QDs by adding indium acetate and palmitic acid in 50ml three-neck flask attached with Schlenk line and having reflux condenser.Figure 2
(a) Synthesis of InP QDs based on the method adopted from [23]. With applications (b) graph representing highest PL intensity [23]. (c) Graph representing highest PL/EL intensity [24]. (d) Graph representing the intensity of fluorescence and absorbance [25].Then this mixture was heated up to120°C for 12 hours. Then N2 gas was filled in the flask, and the temperature was decreased up to 25°C until the white turbid solution appeared. After this, a solution of (TMS)3P and TOP was quickly injected into the flask. Again heating was done up to 305°C for two minutes and then cool down to room temperature. Precipitated QDs were obtained that undergo centrifugation and washing for obtaining pure InP QDs [26]. InP QDs show highly luminescent power with a far-red absorption peak at 631 nm that makes InP QDs attractive for future applications [24]. Indium phosphide nanocrystal is a cadmium-free and low-toxic nanocrystal with extraordinary properties like the highest thermal conductivity, stokes shift emission [25], elongated electro and photoluminescence (PL) generation, near IR emitter, and highest PL quantum yield [27]. InP has got attention because of its high sensitivity for oxidation forming hydroxides and oxides surfaces. Gallium indium phosphide (GaInP) is a semiconductor that has enormous applications in water splitting by using solar energy. But, as its surface is relatively unstable because of corrosion in acidic media, its surface is modified via MoS2 nanoparticles. MoS2 NPs act as a protective layer and enhance catalysis of GaInP because of their anticorrosion properties and are highly recommended as water splitting photocathodes [28]. Indium phosphide is being used in photoelectrochemistry because of having a thin bandgap (1.35 eV) that can absorb solar radiations by decreasing the reflection of photons, resulting in the solar splitting of water. However, use of InP nanowire arrays on Si provides exceptional behavior in this field because of low cost, highest surface area, and catalytic properties [29]. Indium phosphide a promising candidate has proved its applications in optoelectronic as well as in electronic devices. However, nanostructured InP has demonstrated its future applications in the nanophotonic field and optoelectronic devices because of its extraordinary luminescence efficiency harvesting and energy scavenging properties [30]. A scalable artificial synapse based on InP was demonstrated for future applications with neuromorphic processors [31]. These devices when placed to neural spikes reveal the capability of learning without using any external circuit.


## 3. Degree-Based Indices for Indium Phosphide

The number of vertices and edges of InP [s, t] are 10st+3s+3t+2 and 16st, respectively. Furthermore, Table 1 gives details about the edge partition, whereas the comparisons of the indices for InP [s, t] are given in Tables 2 and 3.Table 1
Edge partition of InP [s, t] based on degrees of end vertices.
ςu,ςvFrequency (total edges)Set of edges1,44s+4t−4E12,44st+4s+4tE24,48st−4s−4tE3Table 2
Comparison of indicesABCG and GAG for InP [s, t].
s,tABCGGAG1,116.8414.512,262.4155.763,3133.23120.564,4229.30208.905,5350.62320.796,6497.20456.217,7669.03615.18Table 3
Comparison ofM1G, M2G, M¯1G, and M¯2G indices for InP [s, t].
s,tM1GM2GM¯1GM¯2G1,1921124523542,23805606,4047,4423,38441,32830,54839,7224,41,4842,41693,236127,9145,52,3003,824222,500315,0266,63,2925,552454,052656,3547,74,4607,600831,284121,948Theorem 1.
Consider the graph ofG≅ InP [s, t], which has s,t≥1, and an atom bond connectivity index is corresponding to(8)ABCG=12.626386st+3.843039s+3.843039t−3.464102.Proof. LetG denote InP [s, t] crystallographic structure. The following is the atom bond connectivity index result:(9)ABCG=∑uv∈ΕGςu+ςv−2ςu×ςv=∑uv∈Ε1Gςu+ςv−2ςu×ςv+∑uv∈Ε2Gςu+ςv−2ςu×ςv+∑uv∈Ε3Gςu+ςv−2ςu×ςv=4s+4t−41+4−21×3+4st+4s+4t2+4−22×4+8st−4s−4t4+4−24×4=344s+4t−4+124st+4s+4t+64st−s−t=12.62638st+3.8430939s+3.8433039t−3.464102.Theorem 2.
Consider the graph ofG≅ InP [s, t], which has s,t≥1, and geometric arithmetic index is corresponding to the following:(10)GAG=11.771236st+2.971236s+2.971236t−3.2.Proof. LetG denote InP [s, t] crystallographic structure. The following is the geometric arithmetic index result:(11)GAG=∑uv∈EG2ςu×ςvςu+ςv=∑uv∈E1G2ςu×ςvςu+ςv+∑uv∈E2G2ςu×ςvςu+ςv+∑uv∈E3G2ςu×ςvςu+ςv=4s+4t−421×41+4+4st+4s+4t22×42+4+8st−4s−4t24×44+4=84st+4s+4t3+84s+84t−845+8st−4s−4t=11.77123st+2.971236s+2.971236t−3.2.Theorem 3.
Consider the graph ofG≅InPs,t that has s,t≥1, and the first and second Zagreb index are corresponding to(12)M1G=88st+12s+12t−20,M2G=160st−16s−16t−16.Proof. LetG denote InP [s, t] crystallographic structure. The following is the first Zagreb index result:(13)M1G=∑uv∈EGςu+ςv=∑uv∈E1Gςu+ςv+∑uv∈E2Gςu+ςv+∑uv∈E3Gςu+ςv=4s+4t−41+4+4st+4s+4t2+4+4st+4s+4t4+4=88st+12s+12t−20.The following is the second Zagreb index:(14)M2G=∑uv∈EGςu×ςv=∑uv∈E1Gςu×ςv+∑uv∈E2Gςu×ςv+∑uv∈E3Gςu×ςv=4s+4t−41×4+4st+4s+4t2×4+4st+4s+4t4×4=160st−16s−16t−16.The Zagreb coindices for InP [s, t] are computed in Theorem 4.Theorem 4.
The first and second Zagreb coindices for the graph ofG≅ InP [s, t] with s,t≥1 are corresponding to(15)M1¯G=320s2t2+96s2t+96st2−56st−12s−12t+20,M2¯G=512s2t2−204st+10s+10t+26.Proof. LetG denote InP [s, t] crystallographic structure. The following is the first Zagreb coindex result:(16)M1¯G=∑uv∉EGςu+ςvM1¯G=2EGVG−1−M1G=216st10st+3s+3t+2−1−88st+12s+12t−20=320s2t2+96s2t+96st2−56st−12s−12t+20.The second Zagreb coindex is calculated in the following result:(17)M2¯G=∑uv∉EGςu×ςv=2|EG|2−12M1G−M2G=216st2−1288st+12s+12t−20−160st−16s−16t−16=512s2t2−204st+10s+10t+26.Theorem 5.
The forgotten index for the graph ofG≅ InP [s, t] with s,t≥1 is corresponding to(18)FG=336st+20s+20t−68.Proof. LetG denote InP [s, t] crystallographic structure. The following is the forgotten index result:(19)FG=∑uv∈EGςu2+ςv2=∑uv∈E1Gςu2+ςv2+∑uv∈E2Gςu2+ςv2+∑uv∈E3Gςu2+ςv2=174s+4t−4+204st+4s+4st+328st−4s−4tFG=336st+20s+20t−68.Theorem 6.
The forgotten index for the graph ofG≅ InP [s, t] with s,t≥1 is corresponding to(20)AZIG=496027st−92827s−92827t−25627.Proof. LetG denote InP [s, t] crystallographic structure. The following is the augmented index result:(21)AZIG=∑uv∈EGςu·ςvςu+ςv−23=∑uv∈E1Gςu·ςvςu+ςv−23+∑uv∈E2Gςu·ςvςu+ςu−23+∑uv∈E3Gςu·ςvςu+ςv−23=64274s+4t−4+4st+4s+4t+512278st−4s−4t=496027st−92827s−92827t−25627.


## 4. Heat of Formation and Entropy of Indium Phosphide

In Section2, different topological indices ABCG, GAG, M1G, M2G, M¯1G, M¯2G, FG, and AZIG are calculated for indium phosphide against a different number of unit cells [s, t]. These indices are interrelated with the thermodynamics of InP such as enthalpy/heat of formation and entropy. The values of heat of formation and entropy of InP for 1≤p≤7 and 1≤q≤7 number of unit cells are calculated in Table 4. Standard molar enthalpy of indium phosphide is −69.3kJmol−14, and standard enthalpy for one formula unit was calculated by dividing standard molar enthalpy to Avogadro’s number. Then by multiplying this obtained value with a number of unit cells, enthalpy for that unit cell can be calculated. Standard molar entropy of indium phosphide is 59.75Jk−1mol−1, and standard entropy for one formula unit was calculated by dividing standard molar entropy by Avogadro’s number. Then by multiplying this obtained value with a number of unit cells, entropy for that unit cell can be calculated. The result shows that both enthalpy and entropy increase with an increase in a number of unit cells. From the knowledge of equation S = kB lnW, where W represents microstates (locations/energies of molecules or atoms). As a number of microstates increases, it results in an increase of the entropy of system. See [32] for further explanation.Table 4
Heat of formation and entropy for all formula unites of InP [s, t].
s,tFormula unitesHeat of formationEntropy1,14−4.60443.9732,216−0.188410.158 93,332−0.41430.357 64,464−0.73670.63575,5100−0.01510.993 36,6144−0.01650.014 37,7196−0.02250.019 4While in the case of enthalpy(22)dH=dE+PdV,(23)PV=nRT,implies(24)dH=dE+nRdT.So, according to this equation as a number of moles increases, enthalpy also increases. See [33, 34] for further studies. As chemical graph theory has got a lot of importance in the field of thermochemistry [35], in this study, we demonstrated an explanation of thermodynamic parameters (enthalpy and entropy) of InP using graph theory for easy and better understanding. Tables 2, 3, and 5 are constructed against degree-based indices for InP, and Table 4 is presenting the numerical relation of enthalpy and entropy versus a number of unit cells.Table 5
Comparison ofFG and AZIG indices for InP [s, t].
s,tFGAZIG1,1308105.482,21,356587.853,33,0761,437.624,45,4682,654.815,58,5324,239.406,612, 2686,191.407,716 ,6768,510.81The graphical appearance of formula units and enthalpy/entropy is shown in Figure3.Figure 3
Graph between formula units and enthalpy/entropy.The graphical appearance ofABCG and GAG indices is shown in Figure 4.Figure 4
Comparison ofABCG and GAG indices for InPs,t.The graphical appearance ofM1G, M2G, M1¯G, and M2¯G indices is shown in Figure 5.Figure 5
Comparison ofM1G, M2G, M1¯G, and M2¯G indices for InPs,t.The graphical appearance ofFG and AZIG indices is shown in Figure 6.Figure 6
Comparison ofFG and AZIG indices for InPs,t.


### 4.1. Graphical Models

#### 4.1.1. Heat of Formation versus Indices

A mathematical connection between each index and heat of formationHoF is provided in the following along with the parametric values and the goodness of fit. The general models for all the indices are given in equations (1)–(8), while graphical representations for all the indices versus HoF are shown in Figures 7–14.(i)
General model betweenABC and HoF(25)HoFABC=s1ABC3+t1ABC2+t2ABC+t3,whereABC is normalized by mean 279.8 and standard deviation (std) 239.5.Coefficients:s1=−0.01039−0.01401,−0.006777,t1=1.8371.802,1.871,t2=0.9360.895,0.977,t3=0.13920.1352,0.1431.Goodness of fit:SSE:0.0006455,R2:1,RMSE:0.01467.(ii)
General model betweenGA and HoF(26)HoFGA=s1GA2+s2GA+s3GA3+t1GA2+t2GA+t3,whereGA is normalized by mean 256 and std 220.7.Coefficients:s1=−0.05613−0.2658,0.1535,s2=−0.03786−0.1777,0.102,s3=0.00644−0.07428,0.08716,t1=1.9921.4,2.584,t2=1.1050.4798,1.731,t3=0.13860.09797,0.1793.Goodness of fit:SSE:0.0006455,R2:1,RMSE:0.01467.(iii)
General model betweenM1 and HoF(27)HoFM1=s1M1+s2M13+t1M12+t2M1+t3,whereM1 is normalized by mean 1836 and std 1607.Coefficients:s1=−0.00246−0.03719,0.03227,s2=−0.01174−0.04567,0.0222,t1=1.8391.76,1.917,t2=0.94970.8772,1.022,t3=0.14550.1137,0.1773.Goodness of fit:SSE:0.0006259,R2:1,RMSE:0.01769.(iv)
General model betweenM2 and HoF(28)HoFM2=s1M22+s2M2+s3M23+t1M22+t2M2+t3,whereM2 is normalized by mean 3056 and std 2762.Coefficients:s1=−0.0556−0.2698,0.1586,s2=−0.03933−0.1834,0.1047,s3=0.005922−0.0726,0.08445,t1=2.0051.382,2.627,t2=1.1410.4866,1.795,t3=0.15280.116,0.1897.Goodness of fit:SSE:4.575e−05,R2:1,RMSE:0.006764.(v)
General model betweenM¯1 and HoF(29)HoFM¯1=s1M¯15+t1M¯14+t2M¯13+t3M¯12+t4M¯1+t5,whereM¯1 is normalized by mean 2.341e+05 and std 3.09e+05.Coefficients:s1=−0.0001032−0.0003456,0.0001393,t1=1.1410.1373,2.145,t2=−0.1995−1.352,0.9535,t3=−0.5754−0.7905,−0.3603,t4=−0.152−0.7392,0.4353,t5=0.001936−0.1587,0.1626.Goodness of fit:SSE:0.0005061,R2:1,RMSE:0.0225.(vi)
General model betweenM¯2 and HoF(30)HoFM¯2=s1M¯25+t1M¯24+t2M¯23+t3M¯22+t4M¯2+t5,whereM¯2 is normalized by mean 1.813e+05 and std 2.356e+05.Coefficients:s1=−3.703e−05−0.0003494,0.0002754,t1=1.7061.2,2.212,t2=0.56520.06695,1.064,t3=−0.4214−0.5262,−0.3167,t4=−0.2845−0.4633,−0.1056,t5=−0.04008−0.07154,−0.008617.Goodness of fit:SSE:0.0002722,R2:1,RMSE:0.0165.(vii)
General model betweenF and HoF(31)HoFF=s1F5+t1F4+t2F3+t3F2+t4F+t5,whereF is normalized by mean 6812 and std 6026.Coefficients:s1=0.06101−1.386,1.508,t1=3.727−56.79,64.24,t2=−1.568−8.373,5.238,t3=−9.128−144.2,125.9,t4=−5.501−97.13,86.13,t5=−0.8817−16.64,14.87.Goodness of fit:SSE:0.0001397,R2:1,RMSE:0.01182.(viii)
General model betweenAZI and HoF(32)HoFAZI=s1AZI5+t1AZI4+t2AZI3+t3AZI2+t4AZI+t5,whereAZI is normalized by mean 3390 and std 3103.Coefficients:s1=0.01197−0.4677,0.4917,t1=2.682−27.74,33.11,t2=0.7507−19.85,21.35,t3=−2.333−40.24,35.57,t4=−1.639−32.87,29.59,t5=−0.2712−6.037,5.494.Goodness of fit:SSE:0.0005286,R2:1,RMSE:0.02299.Figure 7
ABC (x-axis) versus HoF (y-axis).Figure 8
GA versus HoF.Figure 9
M1 versus HoF.Figure 10
M2 versus HoF.Figure 11
M¯1 versus HoF.Figure 12
M¯2 versus HoF.Figure 13
F versus HoF.Figure 14
AZI versus HoF.


#### 4.1.2. Entropy versus Indices

The general models between entropyEnt and different indices are given in equations (9)–(16), while their graphical representation for all indices versus Ent are shown in Figures 15–22.(i)
General model betweenABC and Ent(33)EntABC=s1ABC5+t1ABC4+t2ABC3+t3ABC2+t4ABC+t5,whereABC is normalized by mean 279.8 and std 239.5.Coefficients:s1=0.01564−0.05774,0.08902,t1=1.142−4.961,7.245,t2=−0.1919−9.587,9.203,t3=−0.2868−2.726,2.152,t4=0.002129−0.8564,0.8607,t5=0.03417−0.2806,0.3489.Goodness of fit:SSE:0.0003434,R2:1,RMSE:0.01853.(ii)
General model betweenGA and Ent(34)EntGA=s1GA3+t1GA2+t2GA+t3,whereGA is normalized by mean 256 and std 220.7.Coefficients:s1=0.05098−0.01436,0.1163,t1=0.91970.7571,1.082,t2=−0.198−0.3045,−0.09157,t3=0.004903−0.09728,0.1071.Goodness of fit:SSE:0.03282,R2:0.9972,RMSE:0.1046.(iii)
General model betweenM1 and Ent(35)EntM1=s1M15+t1M14+t2M13+t3M12+t4M1+t5,whereM1 is normalized by mean 1836 and std 1607.Coefficients:s1=0.01432−0.04845,0.07708,t1=1.109−4.663,6.881,t2=−0.243−9.205,8.719,t3=−0.311−2.773,2.151,t4=0.007164−0.7696,0.7839,t5=0.03441−0.2578,0.3267.Goodness of fit:SSE:0.0003459,R2:1,RMSE:0.0186.(iv)
General model betweenM2 and Ent(36)EntM2=s1M25+t1M24+t2M23+t3M22+t4M2+t5,whereM2 is normalized by mean 3056 and std 2762.Coefficients:s1=0.01231−0.03602,0.06064,t1=1.061−4.145,6.267,t2=−0.3195−8.509,7.87,t3=−0.3508−2.786,2.084,t4=0.01284−0.6343,0.66,t5=0.03482−0.2217,0.2913.Goodness of fit:SSE:0.0003498,R2:1,RMSE:0.0187.(v)
General model betweenM¯1 and Ent(37)EntM¯1=s1M¯12+s2M¯1+s3M¯13+t1M¯12+t2M¯1+t3,whereM¯1 is normalized by mean 2.341e+05 and std 3.09e+05.Coefficients:s1=−0.08312−0.6134,0.4472,s2=0.01853−0.3932,0.4303,s3=0.06008−0.2452,0.3654,t1=1.019−0.2919,2.329,t2=0.2906−0.4761,1.057,t3=0.06921−0.219,0.3575.Goodness of fit:SSE:0.001714,R2:0.9999,RMSE:0.0414.(vi)
General model betweenM¯2 and Ent(38)EntM¯2=s1M¯25+t1M¯24+t2M¯23+t3M¯22+t4M¯2+t5,whereM¯2 is normalized by mean 1.813e+05 and std 2.356e+05.Coefficients:s1=3.188e−05−0.0002381,0.0003018,t1=1.7031.17,2.236,t2=0.56020.02008,1.1,t3=−0.4242−0.4983,−0.35,t4=−0.285−0.4599,−0.1101,t5=−0.0401−0.07146,−0.008752.Goodness of fit:SSE:0.0002045,R2:1,RMSE:0.0143.(vii)
General model betweenF and Ent(39)EntF=s1F+s2F3+t1F2+t2F+t3,whereF is normalized by mean 6812 and std 6026.Coefficients:s1=0.04603−0.3185,0.4105,s2=0.08438−0.2464,0.4151,t1=0.90480.5207,1.289,t2=−0.1543−0.4216,0.113,t3=0.04557−0.3887,0.4799.Goodness of fit:SSE:0.0295,R2:0.9975,RMSE:0.1215.(viii)
General model betweenAZI and Ent(40)EntAZI=s1AZI+s2AZI4+t1AZI3+t2AZI2+t3AZI+t4,whereAZI is normalized by mean 3390 and std 3103.Coefficients:s1=0.02359−0.2106,0.2578,s2=0.01545−0.2138,0.2447,t1=1.6081.255,1.961,t2=0.48860.3316,0.6456,t3=−0.1166−0.1464,−0.08689,t4=−0.02129−0.2961,0.2536.Goodness of fit:SSE:0.0002636,R2:1,RMSE:0.01624.Figure 15
ABC (x-axis) versus Ent (y-axis).Figure 16
GA versus Ent.Figure 17
M1 versus Ent.Figure 18
M2 versus Ent.Figure 19
M¯1 versus Ent.Figure 20
M¯2 versus Ent.Figure 21
F versus Ent.Figure 22
AZI versus Ent.


## 4.1. Graphical Models

### 4.1.1. Heat of Formation versus Indices

A mathematical connection between each index and heat of formationHoF is provided in the following along with the parametric values and the goodness of fit. The general models for all the indices are given in equations (1)–(8), while graphical representations for all the indices versus HoF are shown in Figures 7–14.(i)
General model betweenABC and HoF(25)HoFABC=s1ABC3+t1ABC2+t2ABC+t3,whereABC is normalized by mean 279.8 and standard deviation (std) 239.5.Coefficients:s1=−0.01039−0.01401,−0.006777,t1=1.8371.802,1.871,t2=0.9360.895,0.977,t3=0.13920.1352,0.1431.Goodness of fit:SSE:0.0006455,R2:1,RMSE:0.01467.(ii)
General model betweenGA and HoF(26)HoFGA=s1GA2+s2GA+s3GA3+t1GA2+t2GA+t3,whereGA is normalized by mean 256 and std 220.7.Coefficients:s1=−0.05613−0.2658,0.1535,s2=−0.03786−0.1777,0.102,s3=0.00644−0.07428,0.08716,t1=1.9921.4,2.584,t2=1.1050.4798,1.731,t3=0.13860.09797,0.1793.Goodness of fit:SSE:0.0006455,R2:1,RMSE:0.01467.(iii)
General model betweenM1 and HoF(27)HoFM1=s1M1+s2M13+t1M12+t2M1+t3,whereM1 is normalized by mean 1836 and std 1607.Coefficients:s1=−0.00246−0.03719,0.03227,s2=−0.01174−0.04567,0.0222,t1=1.8391.76,1.917,t2=0.94970.8772,1.022,t3=0.14550.1137,0.1773.Goodness of fit:SSE:0.0006259,R2:1,RMSE:0.01769.(iv)
General model betweenM2 and HoF(28)HoFM2=s1M22+s2M2+s3M23+t1M22+t2M2+t3,whereM2 is normalized by mean 3056 and std 2762.Coefficients:s1=−0.0556−0.2698,0.1586,s2=−0.03933−0.1834,0.1047,s3=0.005922−0.0726,0.08445,t1=2.0051.382,2.627,t2=1.1410.4866,1.795,t3=0.15280.116,0.1897.Goodness of fit:SSE:4.575e−05,R2:1,RMSE:0.006764.(v)
General model betweenM¯1 and HoF(29)HoFM¯1=s1M¯15+t1M¯14+t2M¯13+t3M¯12+t4M¯1+t5,whereM¯1 is normalized by mean 2.341e+05 and std 3.09e+05.Coefficients:s1=−0.0001032−0.0003456,0.0001393,t1=1.1410.1373,2.145,t2=−0.1995−1.352,0.9535,t3=−0.5754−0.7905,−0.3603,t4=−0.152−0.7392,0.4353,t5=0.001936−0.1587,0.1626.Goodness of fit:SSE:0.0005061,R2:1,RMSE:0.0225.(vi)
General model betweenM¯2 and HoF(30)HoFM¯2=s1M¯25+t1M¯24+t2M¯23+t3M¯22+t4M¯2+t5,whereM¯2 is normalized by mean 1.813e+05 and std 2.356e+05.Coefficients:s1=−3.703e−05−0.0003494,0.0002754,t1=1.7061.2,2.212,t2=0.56520.06695,1.064,t3=−0.4214−0.5262,−0.3167,t4=−0.2845−0.4633,−0.1056,t5=−0.04008−0.07154,−0.008617.Goodness of fit:SSE:0.0002722,R2:1,RMSE:0.0165.(vii)
General model betweenF and HoF(31)HoFF=s1F5+t1F4+t2F3+t3F2+t4F+t5,whereF is normalized by mean 6812 and std 6026.Coefficients:s1=0.06101−1.386,1.508,t1=3.727−56.79,64.24,t2=−1.568−8.373,5.238,t3=−9.128−144.2,125.9,t4=−5.501−97.13,86.13,t5=−0.8817−16.64,14.87.Goodness of fit:SSE:0.0001397,R2:1,RMSE:0.01182.(viii)
General model betweenAZI and HoF(32)HoFAZI=s1AZI5+t1AZI4+t2AZI3+t3AZI2+t4AZI+t5,whereAZI is normalized by mean 3390 and std 3103.Coefficients:s1=0.01197−0.4677,0.4917,t1=2.682−27.74,33.11,t2=0.7507−19.85,21.35,t3=−2.333−40.24,35.57,t4=−1.639−32.87,29.59,t5=−0.2712−6.037,5.494.Goodness of fit:SSE:0.0005286,R2:1,RMSE:0.02299.Figure 7
ABC (x-axis) versus HoF (y-axis).Figure 8
GA versus HoF.Figure 9
M1 versus HoF.Figure 10
M2 versus HoF.Figure 11
M¯1 versus HoF.Figure 12
M¯2 versus HoF.Figure 13
F versus HoF.Figure 14
AZI versus HoF.


### 4.1.2. Entropy versus Indices

The general models between entropyEnt and different indices are given in equations (9)–(16), while their graphical representation for all indices versus Ent are shown in Figures 15–22.(i)
General model betweenABC and Ent(33)EntABC=s1ABC5+t1ABC4+t2ABC3+t3ABC2+t4ABC+t5,whereABC is normalized by mean 279.8 and std 239.5.Coefficients:s1=0.01564−0.05774,0.08902,t1=1.142−4.961,7.245,t2=−0.1919−9.587,9.203,t3=−0.2868−2.726,2.152,t4=0.002129−0.8564,0.8607,t5=0.03417−0.2806,0.3489.Goodness of fit:SSE:0.0003434,R2:1,RMSE:0.01853.(ii)
General model betweenGA and Ent(34)EntGA=s1GA3+t1GA2+t2GA+t3,whereGA is normalized by mean 256 and std 220.7.Coefficients:s1=0.05098−0.01436,0.1163,t1=0.91970.7571,1.082,t2=−0.198−0.3045,−0.09157,t3=0.004903−0.09728,0.1071.Goodness of fit:SSE:0.03282,R2:0.9972,RMSE:0.1046.(iii)
General model betweenM1 and Ent(35)EntM1=s1M15+t1M14+t2M13+t3M12+t4M1+t5,whereM1 is normalized by mean 1836 and std 1607.Coefficients:s1=0.01432−0.04845,0.07708,t1=1.109−4.663,6.881,t2=−0.243−9.205,8.719,t3=−0.311−2.773,2.151,t4=0.007164−0.7696,0.7839,t5=0.03441−0.2578,0.3267.Goodness of fit:SSE:0.0003459,R2:1,RMSE:0.0186.(iv)
General model betweenM2 and Ent(36)EntM2=s1M25+t1M24+t2M23+t3M22+t4M2+t5,whereM2 is normalized by mean 3056 and std 2762.Coefficients:s1=0.01231−0.03602,0.06064,t1=1.061−4.145,6.267,t2=−0.3195−8.509,7.87,t3=−0.3508−2.786,2.084,t4=0.01284−0.6343,0.66,t5=0.03482−0.2217,0.2913.Goodness of fit:SSE:0.0003498,R2:1,RMSE:0.0187.(v)
General model betweenM¯1 and Ent(37)EntM¯1=s1M¯12+s2M¯1+s3M¯13+t1M¯12+t2M¯1+t3,whereM¯1 is normalized by mean 2.341e+05 and std 3.09e+05.Coefficients:s1=−0.08312−0.6134,0.4472,s2=0.01853−0.3932,0.4303,s3=0.06008−0.2452,0.3654,t1=1.019−0.2919,2.329,t2=0.2906−0.4761,1.057,t3=0.06921−0.219,0.3575.Goodness of fit:SSE:0.001714,R2:0.9999,RMSE:0.0414.(vi)
General model betweenM¯2 and Ent(38)EntM¯2=s1M¯25+t1M¯24+t2M¯23+t3M¯22+t4M¯2+t5,whereM¯2 is normalized by mean 1.813e+05 and std 2.356e+05.Coefficients:s1=3.188e−05−0.0002381,0.0003018,t1=1.7031.17,2.236,t2=0.56020.02008,1.1,t3=−0.4242−0.4983,−0.35,t4=−0.285−0.4599,−0.1101,t5=−0.0401−0.07146,−0.008752.Goodness of fit:SSE:0.0002045,R2:1,RMSE:0.0143.(vii)
General model betweenF and Ent(39)EntF=s1F+s2F3+t1F2+t2F+t3,whereF is normalized by mean 6812 and std 6026.Coefficients:s1=0.04603−0.3185,0.4105,s2=0.08438−0.2464,0.4151,t1=0.90480.5207,1.289,t2=−0.1543−0.4216,0.113,t3=0.04557−0.3887,0.4799.Goodness of fit:SSE:0.0295,R2:0.9975,RMSE:0.1215.(viii)
General model betweenAZI and Ent(40)EntAZI=s1AZI+s2AZI4+t1AZI3+t2AZI2+t3AZI+t4,whereAZI is normalized by mean 3390 and std 3103.Coefficients:s1=0.02359−0.2106,0.2578,s2=0.01545−0.2138,0.2447,t1=1.6081.255,1.961,t2=0.48860.3316,0.6456,t3=−0.1166−0.1464,−0.08689,t4=−0.02129−0.2961,0.2536.Goodness of fit:SSE:0.0002636,R2:1,RMSE:0.01624.Figure 15
ABC (x-axis) versus Ent (y-axis).Figure 16
GA versus Ent.Figure 17
M1 versus Ent.Figure 18
M2 versus Ent.Figure 19
M¯1 versus Ent.Figure 20
M¯2 versus Ent.Figure 21
F versus Ent.Figure 22
AZI versus Ent.


## 4.1.1. Heat of Formation versus Indices

A mathematical connection between each index and heat of formationHoF is provided in the following along with the parametric values and the goodness of fit. The general models for all the indices are given in equations (1)–(8), while graphical representations for all the indices versus HoF are shown in Figures 7–14.(i)
General model betweenABC and HoF(25)HoFABC=s1ABC3+t1ABC2+t2ABC+t3,whereABC is normalized by mean 279.8 and standard deviation (std) 239.5.Coefficients:s1=−0.01039−0.01401,−0.006777,t1=1.8371.802,1.871,t2=0.9360.895,0.977,t3=0.13920.1352,0.1431.Goodness of fit:SSE:0.0006455,R2:1,RMSE:0.01467.(ii)
General model betweenGA and HoF(26)HoFGA=s1GA2+s2GA+s3GA3+t1GA2+t2GA+t3,whereGA is normalized by mean 256 and std 220.7.Coefficients:s1=−0.05613−0.2658,0.1535,s2=−0.03786−0.1777,0.102,s3=0.00644−0.07428,0.08716,t1=1.9921.4,2.584,t2=1.1050.4798,1.731,t3=0.13860.09797,0.1793.Goodness of fit:SSE:0.0006455,R2:1,RMSE:0.01467.(iii)
General model betweenM1 and HoF(27)HoFM1=s1M1+s2M13+t1M12+t2M1+t3,whereM1 is normalized by mean 1836 and std 1607.Coefficients:s1=−0.00246−0.03719,0.03227,s2=−0.01174−0.04567,0.0222,t1=1.8391.76,1.917,t2=0.94970.8772,1.022,t3=0.14550.1137,0.1773.Goodness of fit:SSE:0.0006259,R2:1,RMSE:0.01769.(iv)
General model betweenM2 and HoF(28)HoFM2=s1M22+s2M2+s3M23+t1M22+t2M2+t3,whereM2 is normalized by mean 3056 and std 2762.Coefficients:s1=−0.0556−0.2698,0.1586,s2=−0.03933−0.1834,0.1047,s3=0.005922−0.0726,0.08445,t1=2.0051.382,2.627,t2=1.1410.4866,1.795,t3=0.15280.116,0.1897.Goodness of fit:SSE:4.575e−05,R2:1,RMSE:0.006764.(v)
General model betweenM¯1 and HoF(29)HoFM¯1=s1M¯15+t1M¯14+t2M¯13+t3M¯12+t4M¯1+t5,whereM¯1 is normalized by mean 2.341e+05 and std 3.09e+05.Coefficients:s1=−0.0001032−0.0003456,0.0001393,t1=1.1410.1373,2.145,t2=−0.1995−1.352,0.9535,t3=−0.5754−0.7905,−0.3603,t4=−0.152−0.7392,0.4353,t5=0.001936−0.1587,0.1626.Goodness of fit:SSE:0.0005061,R2:1,RMSE:0.0225.(vi)
General model betweenM¯2 and HoF(30)HoFM¯2=s1M¯25+t1M¯24+t2M¯23+t3M¯22+t4M¯2+t5,whereM¯2 is normalized by mean 1.813e+05 and std 2.356e+05.Coefficients:s1=−3.703e−05−0.0003494,0.0002754,t1=1.7061.2,2.212,t2=0.56520.06695,1.064,t3=−0.4214−0.5262,−0.3167,t4=−0.2845−0.4633,−0.1056,t5=−0.04008−0.07154,−0.008617.Goodness of fit:SSE:0.0002722,R2:1,RMSE:0.0165.(vii)
General model betweenF and HoF(31)HoFF=s1F5+t1F4+t2F3+t3F2+t4F+t5,whereF is normalized by mean 6812 and std 6026.Coefficients:s1=0.06101−1.386,1.508,t1=3.727−56.79,64.24,t2=−1.568−8.373,5.238,t3=−9.128−144.2,125.9,t4=−5.501−97.13,86.13,t5=−0.8817−16.64,14.87.Goodness of fit:SSE:0.0001397,R2:1,RMSE:0.01182.(viii)
General model betweenAZI and HoF(32)HoFAZI=s1AZI5+t1AZI4+t2AZI3+t3AZI2+t4AZI+t5,whereAZI is normalized by mean 3390 and std 3103.Coefficients:s1=0.01197−0.4677,0.4917,t1=2.682−27.74,33.11,t2=0.7507−19.85,21.35,t3=−2.333−40.24,35.57,t4=−1.639−32.87,29.59,t5=−0.2712−6.037,5.494.Goodness of fit:SSE:0.0005286,R2:1,RMSE:0.02299.Figure 7
ABC (x-axis) versus HoF (y-axis).Figure 8
GA versus HoF.Figure 9
M1 versus HoF.Figure 10
M2 versus HoF.Figure 11
M¯1 versus HoF.Figure 12
M¯2 versus HoF.Figure 13
F versus HoF.Figure 14
AZI versus HoF.


## 4.1.2. Entropy versus Indices

The general models between entropyEnt and different indices are given in equations (9)–(16), while their graphical representation for all indices versus Ent are shown in Figures 15–22.(i)
General model betweenABC and Ent(33)EntABC=s1ABC5+t1ABC4+t2ABC3+t3ABC2+t4ABC+t5,whereABC is normalized by mean 279.8 and std 239.5.Coefficients:s1=0.01564−0.05774,0.08902,t1=1.142−4.961,7.245,t2=−0.1919−9.587,9.203,t3=−0.2868−2.726,2.152,t4=0.002129−0.8564,0.8607,t5=0.03417−0.2806,0.3489.Goodness of fit:SSE:0.0003434,R2:1,RMSE:0.01853.(ii)
General model betweenGA and Ent(34)EntGA=s1GA3+t1GA2+t2GA+t3,whereGA is normalized by mean 256 and std 220.7.Coefficients:s1=0.05098−0.01436,0.1163,t1=0.91970.7571,1.082,t2=−0.198−0.3045,−0.09157,t3=0.004903−0.09728,0.1071.Goodness of fit:SSE:0.03282,R2:0.9972,RMSE:0.1046.(iii)
General model betweenM1 and Ent(35)EntM1=s1M15+t1M14+t2M13+t3M12+t4M1+t5,whereM1 is normalized by mean 1836 and std 1607.Coefficients:s1=0.01432−0.04845,0.07708,t1=1.109−4.663,6.881,t2=−0.243−9.205,8.719,t3=−0.311−2.773,2.151,t4=0.007164−0.7696,0.7839,t5=0.03441−0.2578,0.3267.Goodness of fit:SSE:0.0003459,R2:1,RMSE:0.0186.(iv)
General model betweenM2 and Ent(36)EntM2=s1M25+t1M24+t2M23+t3M22+t4M2+t5,whereM2 is normalized by mean 3056 and std 2762.Coefficients:s1=0.01231−0.03602,0.06064,t1=1.061−4.145,6.267,t2=−0.3195−8.509,7.87,t3=−0.3508−2.786,2.084,t4=0.01284−0.6343,0.66,t5=0.03482−0.2217,0.2913.Goodness of fit:SSE:0.0003498,R2:1,RMSE:0.0187.(v)
General model betweenM¯1 and Ent(37)EntM¯1=s1M¯12+s2M¯1+s3M¯13+t1M¯12+t2M¯1+t3,whereM¯1 is normalized by mean 2.341e+05 and std 3.09e+05.Coefficients:s1=−0.08312−0.6134,0.4472,s2=0.01853−0.3932,0.4303,s3=0.06008−0.2452,0.3654,t1=1.019−0.2919,2.329,t2=0.2906−0.4761,1.057,t3=0.06921−0.219,0.3575.Goodness of fit:SSE:0.001714,R2:0.9999,RMSE:0.0414.(vi)
General model betweenM¯2 and Ent(38)EntM¯2=s1M¯25+t1M¯24+t2M¯23+t3M¯22+t4M¯2+t5,whereM¯2 is normalized by mean 1.813e+05 and std 2.356e+05.Coefficients:s1=3.188e−05−0.0002381,0.0003018,t1=1.7031.17,2.236,t2=0.56020.02008,1.1,t3=−0.4242−0.4983,−0.35,t4=−0.285−0.4599,−0.1101,t5=−0.0401−0.07146,−0.008752.Goodness of fit:SSE:0.0002045,R2:1,RMSE:0.0143.(vii)
General model betweenF and Ent(39)EntF=s1F+s2F3+t1F2+t2F+t3,whereF is normalized by mean 6812 and std 6026.Coefficients:s1=0.04603−0.3185,0.4105,s2=0.08438−0.2464,0.4151,t1=0.90480.5207,1.289,t2=−0.1543−0.4216,0.113,t3=0.04557−0.3887,0.4799.Goodness of fit:SSE:0.0295,R2:0.9975,RMSE:0.1215.(viii)
General model betweenAZI and Ent(40)EntAZI=s1AZI+s2AZI4+t1AZI3+t2AZI2+t3AZI+t4,whereAZI is normalized by mean 3390 and std 3103.Coefficients:s1=0.02359−0.2106,0.2578,s2=0.01545−0.2138,0.2447,t1=1.6081.255,1.961,t2=0.48860.3316,0.6456,t3=−0.1166−0.1464,−0.08689,t4=−0.02129−0.2961,0.2536.Goodness of fit:SSE:0.0002636,R2:1,RMSE:0.01624.Figure 15
ABC (x-axis) versus Ent (y-axis).Figure 16
GA versus Ent.Figure 17
M1 versus Ent.Figure 18
M2 versus Ent.Figure 19
M¯1 versus Ent.Figure 20
M¯2 versus Ent.Figure 21
F versus Ent.Figure 22
AZI versus Ent.


## 5. Conclusion

In this paper, we have focused to investigate the relationship between underlying graphical properties and the thermodynamical properties of indium phosphide. Figure2 clearly depicts the synthesis of indium phosphide in a systematic way as explained before. At first, the topological degree-based indices were calculated that were lately integrated with thermodynamical properties of indium phosphide. Such a link had been established by fitting the curve between each index and each thermodynamical property. Entropy and heat of formation are two types of thermodynamical properties that were indulged in this study. The rational fitting method was applied using MATLAB software as this method was providing the least mean squared error among all built-in methods. Such a multidisciplinary approach would provide a comprehensive insight into the structural properties of indium phosphide in more detail and depth. [36].



---

*Source: 9707749-2022-01-30.xml*