1.
Introduction
Every graph taken into consideration here is a simple graph Γ=(V,E), where |V|=ν and |E|=m. The count of edges connected to a vertex v∈V, denoted by dΓ(v) is the degree of that vertex. The degree of an edge e=uv is defined by dΓ(e)=dΓ(u)+dΓ(v)−2. The parameters Δ (resp. δ) and Δ′(δ′) represent the maximum (resp. minimum) vertex and edge degree of a graph, respectively. If ΔΓ=δΓ=ν−1 then the graph Γ is said to be a complete graph denoted by Kν. The graph Γ is called bipartite if V(Γ) is partitioned into two sets say, M and N (partite sets), such that every edge in Γ has one endpoint in M and the other in N. If every vertex of M is adjacent to every vertex of N with |M|=r and |N|=s, then the graph is called the complete bipartite graph, denoted as Kr,s. The graph K1,ν−1 denoted by Sν is a star graph, and the graph Kr,r is called the equi-bipartite graph. The graph ¯Γ represents the complement of a graph Γ, which is defined on the same vertex set as Γ such that if two vertices are adjacent in ¯Γ, then they are not adjacent in Γ. For more graph-theoretic terminologies, we refer to the book by Harary [19]. For any real number x, the floor function is the greatest integer less than or equal to x, denoted as ⌊x⌋.
The first degree-based molecular descriptor, the Zagreb index, was developed by Gutman and Trinajstić [9]. It first emerged in the topological formula for conjugated molecules regarding their total π-electron energy. The first Zagreb index is defined as:
The Zagreb indices were reformulated in 2004 by Miliˊceviˊc et al. [28] in terms of edge degree, where the edge degree is given by dΓ(e)=dΓ(u)+dΓ(v)−2 for e=uv∈E(Γ). Thus, the first reformulated Zagreb index is given by
The reader is referred to [9] regarding applications of Zagreb indices.
The sum of the absolute values of eigenvalues of the adjacency matrix of Γ gives us the energy E(Γ) of a graph Γ. This quantity is introduced in [10]. Suppose λ1≥λ2≥…≥λν are the eigenvalues of the adjacency matrix A(Γ), then the energy of the graph Γ is given by
Other graph energies were also introduced and studied. Indulal et al. [20] presented some results on the distance energy of graphs. Gutman and Zhou [8] investigated the Laplcaian energy of graphs and derived some extremal results from it.
The extended adjacency matrix of graph Γ was proposed by Yang et al. [36] in 1994, denoted by Aex(Γ). Its (i,j)-entry is defined to be equal to 12(djdi+didj) if vi∼vj and 0 otherwise. Since Aex is a real symmetric matrix of order ν, all its eigenvalues are real, which are denoted as η1≥η2≥…≥ην. Yang et al. [36] also investigated the extended graph energy by summing the absolute values of the eigenvalues of the Aex-matrix, defined as
For recent studies on graph energy, we refer to [5,6,23,24,33,35].
1.1. Edge-weighted graph energy
For a given graph Γ, we define the following terminologies:
Definition 1.1. Let V={v1,v2,…,vν} be the vertex set and {w1,w2,…,wν} be the vertex-weights, then the vertex weight for vi∈V is defined as w(vi)=dΓ(vi). The range for a vertex-weight of vi∈V is 0≤w(vi)≤ΔΓ (here w(vi)=0 if Γ is disconnected).
Definition 1.2. Let EΓ={e1,e2,…,em} be the edge set, then the edge-weight of ei=uivi is defined by w(ei)=dΓ(ui)+dΓ(vi)−2. The range for an edge-weight of e∈E is 0≤w(e)≤M1(Γ)2−m (where w(e)=0 if and only if Γ=K2.
Definition 1.3. The weighted degree of a vertex v∈V(Γ) is defined as
It is very clear that w(e)=dΓ(e) for every e∈E.
Observe that
Motivated by the extended adjacency matrix of graph Γ, we introduce a new edge-weighted adjacency matrix of graph Γ denoted by Aw(Γ). It is defined in such a way that, for any vi that is adjacent to vj, its (i,j)-entry equals dΓ(vi)+dΓ(vj)−2; otherwise, it equals to 0. In fact, Aw(Γ) is a real symmetric matrix of order ν. Hence, all its eigenvalues are real and can be arranged as λw1≥λw2≥…≥λwν, where the largest eigenvalue λw1 is called the spectral radius of Aw(Γ). The edge-weighted graph energy of Γ is given by
Following the types of adjacencies in [34], the type of adjacency employed by the edge-weighted matrix is the vertex-based adjacency.
1.2. Some useful identities
This subsection presents some basic properties of edge-weighted eigenvalues of graphs.
(1) ∑νi=1λwi=0.
(2) ∑νi=1(λwi)2=2∑mi=1[w(ei)]2=2∑mi=1dΓ(ei)2=2EM1(Γ).
(3) ∑0≤i≤jλwiλwj=−∑mi=1[w(ei)]2=−∑mi=1dΓ(ei)2=−EM1(Γ).
Moreover, observe that,
(1) ∑νi=1.(λwi)2=2EM1(Γ).
(2) ∑0≤i≤jλwiλwj=−EM1(Γ).
The following example delivers the edge-weighted energy of standard graph families.
Example 1.4. (1) The edge-weighted energy Ew(Kν) of Kν is 4(ν−1)(ν−2).
(2) The edge-weighted energy Ew(Kr,s) of Kr,s is 2(r+s−2)√rs.
(3) The edge-weighted energy Ew(Sν) of Sν is 2(ν−2)√(ν−1).
(4) The edge-weighted energy Ew(Kr,r) of Kr,r is 4r(r−1).
Proof. The edge-weighted energy Ew(Kr,s) of Kr,s is 2(r+s−2)√rs. Replacing s by r in Ew(Kr,s) of Kr,s, we get Ew(Kr,r) of Kr,r as 2(r+r−2)√r2=4r(r−1). □
2.
Auxiliary results
In subsequent sections, we need the following already established results:
Lemma 2.1. [35] Let C=(cij) and D=(dij) be real symmetric non-negative matrices of order ν. If C≥D, i.e., cij≥dij for all i,j, then λ1(C)≥λ1(D), whereas λ1 is the largest eigenvalue.
Lemma 2.2. [35] Let Γ be a connected graph of order ν with m edges. Then
with equality if and only if Γ≅K1,ν−1 or Γ≅Kν.
Here we deliver the well-known Cauchy–Schwartz inequality.
Lemma 2.3. (Cauchy–Schwartz inequality)[3] Let ri and si, 1≤i≤ν be any real numbers, then
The Ozeki inequality is frequently used in the spectral analysis of graphs.
Lemma 2.4. (Ozeki inequality)[30] If ri and si (1≤i≤ν) are non-negative real numbers, then
where P1=max1≤i≤ν{ri}, P2=max1≤i≤ν{si}, p1=min1≤i≤ν{ri}, p2=min1≤i≤ν{si}.
The following inequality has been retrieved from Dragomir [7].
Lemma 2.5. [7] Let pi, qi, ri and si be sequences of real numbers, and mi, νi are non-negative for i=1,2,…,ν. Then the following inequality is valid:
Jog and Gurjar [23] used the following inequality while studying bounds on the distance energy of graphs:
Lemma 2.6. Let ri, 1≤i≤ν be any real numbers, then
The following inequality was shown by Pólya and Szegó in their book [31].
Lemma 2.7. [31] Suppose ri and si, 1≤i≤ν are positive real numbers, then
where P1=max1≤i≤ν{ri}, P2=max1≤i≤ν{si}, p1=min1≤i≤ν{ri}, p2=min1≤i≤ν{si}.
The following classical inequality was proven by Biernacki et al. [2].
Lemma 2.8. [2] Suppose ri and si, 1≤i≤ν are positive real numbers, then
where r, s, R, and S are real constants such that for each i, 1≤i≤ν, r≤ri≤R, and s≤si≤S. Further, α(ν)=ν⌊ν2⌋(1−1ν⌊ν2⌋).
Diaz and Metcalf [4] delivered a proof of the following inequality:
Lemma 2.9. [4] Let ri and si, 1≤i≤ν are nonnegative real numbers, then
where p and P are real constants, so that for each i, 1≤i≤ν, holds, pri≤si≤Pri.
Next, we prove the following inequality on the edge-connected eigenvalues:
Lemma 2.10. Let Γ be a connected graph of order ν≥2. Then λw1>λw2.
Proof. Let us assume, for the sake of contradiction, that λw1=λw2. Since Γ is connected, Aw(Γ) is an irreducible non-negative ν×ν matrix. By the Perron–Frobenius theorem, the eigenvector x corresponding to λw1 has all components positive. Let y be an eigenvector corresponding to λw2. Since λw1=λw2, any linear combination of x and y would be an eigenvector corresponding to λw1. This implies that it would be possible to construct an eigenvector with some zero components, which contradicts the fact that all components of x are positive. Hence, we must have λw1>λw2. □
Next, we deliver a characterization for an ν-vertex satisfying |λw1|=|λw2|=…=|λwν|.
Proposition 2.11. Let Γ be a graph of order ν. Then |λw1|=|λw2|=…=|λwν| if and only if Γ≅¯Kν or Γ≅ν2K2.
Proof. First, assume that |λw1|=|λw2|=…=|λwν|. Let S be the number of isolated vertices in Γ. If S≥1, then λw1=λw2=…=λwν=0, hence Γ≅¯Kν or Γ≅ν2K2. Otherwise, if the maximum degree Δ≥2, then Γ contains a connected component H with at least 3 vertices. If H=Kν, ν≥3, then by Lemma 2.10, |λw1|=2(ν−1)(ν−2) and |λw2|=2(ν−2), clearly |λw1|>|λw2|, a contradiction. Otherwise, if H is not a complete graph, then by Lemma 2.10, λw1>λw2, a contradiction.
Conversely, one can easily check that |λw1|=|λw2|=…=|λwν| holds for ¯Kν and ν2K2. □
3.
Main results
This section delivers various upper/lower extremal values for Ew(Γ) and λw1 for an (m,ν)-graph with ν-vertices and m-edges.
A sharp upper bound on λw1 (i.e., edge-weighted spectral radius) is being proven in the following result.
Theorem 3.1. Let Γ be an (m,ν)-graph possessing the maximum degree Δ. Then
with λw1=2(Δ−1)√1−ν+2m⟺ Γ≅Kν, ν≥2.
Proof. Since Aw(Γ)≤2(Δ−1)A(Γ) and λw1 be its spectral radius. Employing Lemma 2.1, one has
By Lemma 2.2, one obtains
Also, λw1=2(Δ−1)√1−ν+2m in (3.1) ⟺Γ≅Kν, ν≥2. □
The next theorem delivers an upper extremal value considering the first–formulated Zagreb EM1 index of Γ.
Theorem 3.2. If Γ is an ν-vertex graph having λw1 as its spectral radius (the largest eigenvalue), then
Proof. Since ∑νk=1λwk=0 it can be rewritten as ∑νk=2λwk=−λw1. Further, (∑νk=1(λwk)2)=2EM1(Γ), (∑νk=2(λwk)2)=(∑νk=1(λwk)2−(λw1)2)=(2EM1(Γ)−(λw1)2) and (∑νk=21)=(ν−1).
Put rk=1 and sk=λwk in Lemma 2.3 and we obtain
Hence, we have furnished the proof. □
Some lower and upper extremal values on Ew i.e., edge-weighted energy of Γ.
Theorem 3.3. For a connected ν-vertex graph Γ, Ew satisfies
Proof. For the upper bound, consider Lemma 2.3, i.e.,
put rk=1 and sk=|λwk|2 in Lemma 2.3, we obtain
since (∑νk=1|λwk|)=Ew(Γ), (∑νk=112)=1 and (∑νk=1(|λwk|)2)=2EM1(Γ), we have
Similarly, for a lower bound, consider Lemma 2.6, i.e.,
put rk=|λwk| in Lemma 2.6, we obtan
Hence, by combining the upper bound and the lower bound, we obtain the required result. □
In terms of the EM1 index, our next result delivers another lower extremal value on Ew.
Theorem 3.4. Any (m,ν)-graph Γ satisfies
where |λwν| (resp. |λw1|) are minimum (resp. maximum) values of |λwk|.
Proof. Let |λw1|≥|λw2|≥…≥|λwν| be the eigenvalues of Aw(Γ). By putting rk=1, sk=|λwk| P1=1, P2=|λw1|, p1=1 and p2=|λwν| in Lemma 2.4, one gets
since (∑νk=1(λwk)2)=2EM1(Γ) we have
Hence, the proof has been furnished. □
The following theorem further refines the bound in Theorem 3.4.
Theorem 3.5. For an (ν,m)-graph Γ, assume λw1≥λw2≥…≥λwν are eigenvalues of Aw(Γ). This implies that,
where α(ν)=ν⌊ν2⌋(1−1ν⌊ν2⌋).
Proof. Let |λw1|≥|λw2|≥…≥|λwν| be the eigenvalues of Aw(Γ). By putting rk=|λwk|=sk, R=|λw1|=S, and r=|λwν|=s in Lemma 2.8, one obtains
Hence the proof has been furnished. □
For a non-zero eigenvalue of a graph, the following theorem delivers yet another lower bound on the edge-weighted energy Ew of graphs.
Theorem 3.6. If the eigenvalues of Aw(Γ) are non-zero, then
where |λw1| and |λwν| are the maximum and minimum of |λwk|.
Proof. Let |λw1|≥|λw2|≥…≥|λwν| be the eigenvalues of Aw(Γ). By putting rk=|λk| and sk=1 in Lemma 2.7, we obtain
Hence the proof is completed. □
The following two results deliver lower bounds on Ew in terms of EM1, the smallest and largest eigenvalues of graphs.
Theorem 3.7. Assuming λw1≥λw2≥…≥λwν to be the eigenvalues of Aw(Γ), where Γ is an (ν,m)-graph. Then
Proof. Let λw1≥λw2≥…≥λwν be the eigenvalues of Aw(Γ). By putting rk=|λwk|, sk=1, p=λwν and P=λw1, in Lemma 2.9, we obtain
Hence, the proof has been completed. □
Theorem 3.8. Assuming λw1≥λw2≥…≥λwν to be the eigenvalues of Aw(Γ), where Γ is an (ν,m)-graph, then
Proof. Let λw1≥λw2≥…≥λwν be the eigenvalues of Aw(Γ). By putting rk=|λwk|, sk=1, p=λwν and P=λw1, in Lemma 2.9, we obtain
Hence the proof has been furnished. □
Next, a sharp upper extremal value on Ew is proven.
Theorem 3.9. If Γ is a non-empty graph of order ν. Then
Proof. Let λw1≥λw2≥…≥λwν be the eigenvalues of Aw(Γ). Substituting pk=|λwk|=qk and rk=sk=mk=νk=1 in Lemma 2.5, we obtain
□
4.
Applications on the edge-weighted energy of graphs
Until now, many researchers have studied the predictive potential of molecular descriptors (mainly degree-, distance-, and eigenvalue-based) for estimating the π-electron energy (Eπ) of benzenoid hydrocarbons (BHs) and also for estimating the enthalpy of formation ΔHof and boiling point Bp of BHs. For instance, Hayat et al. [16] (resp. Hayat and coauthors [14,27]) investigated the efficiency of degree-dependent (resp. eigenvalues-dependent) graphical descriptors for estimating Eπ of BHs. Similar studies were conducted by Hayat et al. [15] (Hayat and Liu [12]) for distance-related and temperature-based graphical descriptors. For predicting physicochemical properties such as Bp and ΔHof of BHs, comparative studies for distance-dependent, temperature-related, degree-related, and eigenvalue-related descriptors were conducted in [13,17,18,26], respectively. For more studies on QSPR, we refer to [1,11,29,32].
In this section, we calculate energy and edge-weighted energy for molecular graphs of 22 benzenoid hydrocarbons which are listed in Table 1. The adjacency matrix is defined, and eigenvalues, energy, and edge-weighted energy are calculated using the Python programming language. The correlation and regression models are obtained for the physicochemical properties of 22 BHs, for which the data is taken from [25] (refer to Table 2), and the predictive ability is tested using energy E(Γ) and edge-weighted energy Ew(Γ). In Table 3, the values of 11 molecular descriptors for 22 BHs are listed, which are from [25]. Note that we consider carbon–carbon structures as chemical graphs. One can consider other construction, of chemical graph based on molecular alignment [22].
This subsection records all the data sets that we employ for our structure-property models.
5.
Results and discussions
The intercorrelation between the physicochemical properties of polycyclic aromatic hydrocarbons, such as Kovats retention index (RI), acentric factor (ω), octanol-water partition coefficient (logP), boiling point (BP), enthalpy of formation (ΔHf) and entropy (S), with graph energy E(Γ) and edge-weighted energy Ew(Γ), is analysed in Table 4. Also, the intercorrelation between 11 molecular descriptors such as atom bond connectivity index (ABC), forgotten index (F), 1st and 2nd Zagreb invariants (M1 and M2), sum division degree index (SDD), reciprocal Randiˊc index (RR), classical Randić index (R), redefined Zagreb index (ReZM) with graph energy E(Γ) and edge-weighted energy Ew(Γ) is listed in Table 5. We observe that the 11 molecular descriptors are highly intercorrelated with edge-weighted energy Ew(Γ) with r>0.97, which is highlighted in Table 5.
The value of r for Ew(Γ) ranges from 0.972 to 0.980.
5.1. Regression models
The quadratic regression models for physico-chemical properties (PPs) such as Kovats retention index (RI), acentric factor (ω), octanol–water partition coefficient (logP), boiling point (BP), enthalpy of formation (ΔHf), and entropy (S) are derived with respect to graph energy E(Γ) and edge-weighted energy Ew(Γ). The symbols ν, r, F, and SE are used to represent population, correlation coefficient, F-values, and the standard error of the estimate, respectively. Note that, in general, quadratic models have very bad predictive power, even if they have good estimating power. The reader is referred to [21] for diversity in detailed regression analysis.
The quadratic regression model is defined as
The quadratic regression of PP with E(Γ) is as follows:
The quadratic regression of PP with Ew(Γ) is as follows:
5.2. Analysis
The following analysis can be made from the quadratic regression models:
● The correlation coefficient r for quadratic regression models gives high predictability for physicochemical properties with respect to graph energy E(Γ) and edge-weighted energy Ew(Γ).
● The quadratic regression models for E(Γ) give high intercorrelation with a correlation value of r=0.9823 for the acentric factor.
● The degree-2 polynomial regression model for E(Γ) gives appreciable intercorrelation with a correlation value of r=0.9230 for the boiling point, r=0.9241 for the log P, r=0.9428 for the retention index.
● The quadratic regression model for E(Γ) is weakly intercorrelation with correlation value r=0.8491 for the entropy, and r=0.8831 for the enthalpy.
● The quadratic regression models for Ew(Γ) give high intercorrelation with a correlation value of r=0.9705 for the boiling point and r=0.9741 for retention index.
● The quadratic regression model for Ew(Γ) gives appreciable intercorrelation with a correlation value of r=0.9208 for entropy, r=0.9544 for acentfac, r=0.9695 for log P, and r=0.9386 for the enthalpy.
● From all the 12 quadratic regression models, it has been observed that the significance F is 0.000.
The scattered curve diagram of E(Γ) with physiochemical properties are depicted in Figures 1–6.
The scattered curve diagram of Ew(Γ) with physicochemical characteristics are depicted in Figures 7–12.
6.
Conclusions
This paper puts forward the edge-weighted adjacency matrix Aw(Γ) of a graphical structure Γ. The energy Ew(Γ) as well as the spectral radius λw1 of the Aw(Γ) have been studied, and lower and upper extremes are derived for λw1 and Ew(Γ) in terms of other graphical parameters. Further, we calculated the graph energy and edge-weighted energy of 22 BHs by drawing their molecular graphs to check the predictive potential of the physicochemical characteristics of BHs. Polynomials of degree-2 regression models were generated for Kovats retention index (RI), acentric factor (ω), octanol-water partition coefficient (logP), boiling point (BP), enthalpy of formation (ΔHf), and entropy (S) using theses two graph energies. We also found correlation coefficients of the physicochemical properties and molecular descriptors of BHa corresponding to the two graph energies E(Γ) and Ew(Γ).
Author contributions
All authors contributed equally to this paper. All authors have read and approved the final version of the manuscript for publication.
Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
Acknowledgments
S. Hayat is supported by UBD Faculty Research Grants under Grant Number UBD/RSCH/1.4/FICBF(b)/2022/053 and the National Natural Science Foundation of China (Grant No. 622260-101). A. Khan was sponsored by the National Natural Science Foundation of China grant Nos. 622260-101 and 12250410247, and also by the Ministry of Science and Technology of China, grant No. WGXZ2023054L. Mohammed J. F. Alenazi extends his appreciation to Researcher Supporting Project number (RSPD2024R582), King Saud University, Riyadh, Saudi Arabia.
Conflict of interest
The authors declare that they have no known competing financial interests.