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.

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 a molecular/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

Vukičević and Furtula [14] introduced the geometric arithmetic index as follows:

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 , Zagreb indices are defined as follows:

The first and second Zagreb coindices defined in [15] are as follows:

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:

Wang et al. [18], inspired by the ABC index, devised a new modification, which is known as augmented Zagreb index and described as follows:

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 the heat 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 , , , or .

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 . Figure 1 demonstrates the crystal structure of InP [21].

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 three-neck flask attached with Schlenk line and having reflux condenser.

Then this mixture was heated up to for 12 hours. Then 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 and , 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.

Theorem 1. Consider the graph of InP [s, t], which has , and an atom bond connectivity index is corresponding to

Proof. Let denote InP [s, t] crystallographic structure. The following is the atom bond connectivity index result:

Theorem 2. Consider the graph of InP [s, t], which has , and geometric arithmetic index is corresponding to the following:

Proof. Let denote InP [s, t] crystallographic structure. The following is the geometric arithmetic index result:

Theorem 3. Consider the graph of that has , and the first and second Zagreb index are corresponding to

Proof. Let denote InP [s, t] crystallographic structure. The following is the first Zagreb index result:

The following is the second Zagreb index:

The Zagreb coindices for InP [s, t] are computed in Theorem 4.

Theorem 4. The first and second Zagreb coindices for the graph of InP [s, t] with are corresponding to

Proof. Let denote InP [s, t] crystallographic structure. The following is the first Zagreb coindex result:

The second Zagreb coindex is calculated in the following result:

Theorem 5. The forgotten index for the graph of InP [s, t] with is corresponding to

Proof. Let denote InP [s, t] crystallographic structure. The following is the forgotten index result:

Theorem 6. The forgotten index for the graph of InP [s, t] with is corresponding to

Proof. Let denote InP [s, t] crystallographic structure. The following is the augmented index result:

4. Heat of Formation and Entropy of Indium Phosphide

In Section 2, different topological indices , , , , , , , and 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 and number of unit cells are calculated in Table 4. Standard molar enthalpy of indium phosphide is , 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 , 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 =  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.

While in the case of enthalpy

implies

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.

The graphical appearance of formula units and enthalpy/entropy is shown in Figure 3.

The graphical appearance of and indices is shown in Figure 4.

The graphical appearance of , , , and indices is shown in Figure 5.

The graphical appearance of and indices is shown in Figure 6.

4.1. Graphical Models
4.1.1. Heat of Formation versus Indices

A mathematical connection between each index and heat of formation 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 are shown in Figures 714.(i)General model between and where is normalized by mean 279.8 and standard deviation (std) 239.5.Coefficients: .Goodness of fit: .(ii)General model between and where is normalized by mean 256 and std 220.7.Coefficients: .Goodness of fit:.(iii)General model between and where is normalized by mean 1836 and std 1607.Coefficients: .Goodness of fit:.(iv)General model between and where is normalized by mean 3056 and std 2762.Coefficients: .Goodness of fit:.(v)General model between and where is normalized by mean and std .Coefficients: .Goodness of fit:.(vi)General model between and where is normalized by mean and std .Coefficients: .Goodness of fit:.(vii)General model between and where is normalized by mean 6812 and std 6026.Coefficients: .Goodness of fit:.(viii)General model between and where is normalized by mean 3390 and std 3103.Coefficients: .Goodness of fit:.

4.1.2. Entropy versus Indices

The general models between entropy and different indices are given in equations (9)–(16), while their graphical representation for all indices versus are shown in Figures 1522.(i)General model between and where is normalized by mean 279.8 and std 239.5.Coefficients: .Goodness of fit:.(ii)General model between and where is normalized by mean 256 and std 220.7.Coefficients: .Goodness of fit:.(iii)General model between and where is normalized by mean 1836 and std 1607.Coefficients: .Goodness of fit:.(iv)General model between and where is normalized by mean 3056 and std 2762.Coefficients: .Goodness of fit:.(v)General model between and where is normalized by mean and std .Coefficients: .Goodness of fit:.(vi)General model between and where is normalized by mean and std .Coefficients: .Goodness of fit:.(vii)General model between and where is normalized by mean 6812 and std 6026.Coefficients: .Goodness of fit:.(viii)General model between and where is normalized by mean 3390 and std 3103.Coefficients: .Goodness of fit:.

5. Conclusion

In this paper, we have focused to investigate the relationship between underlying graphical properties and the thermodynamical properties of indium phosphide. Figure 2 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].

Data Availability

The data used to support the findings of this study are cited at relevant places within the text as references.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors contributed equally to this study.

Acknowledgments

This research project is supported by the Natural Science Foundation of Anhui Province Higher School (KJ2021A1154).