
This paper examines if clean energy stocks help investors in managing carbon risk. We use the price of the European Union Allowance (EUA) and European clean energy index (ERIX) for the three phases of the EU-Emission Trading Scheme. Analyzing the time-varying correlation and volatility of EUA stock and ERIX through generalized orthogonal GO-GARCH model, the empirical results reveal relative independence of the European renewable energy market from the carbon market providing diversification benefits and value addition by including carbon assets in clean energy stock portfolio. Furthermore, three portfolios with different weight allocation strategies reveal that the carbon asset provides risk and downside risk benefits when mixed with a clean energy stock portfolio. These results are useful for investors who enter the market for value maximization and the regulators striving to make strategies for managing carbon risk.
Citation: Memoona Kanwal, Hashim Khan. Does carbon asset add value to clean energy market? Evidence from EU[J]. Green Finance, 2021, 3(4): 495-507. doi: 10.3934/GF.2021023
[1] | Ze-Miao Dai, Jia-Bao Liu, Kang Wang . Analyzing the normalized Laplacian spectrum and spanning tree of the cross of the derivative of linear networks. AIMS Mathematics, 2024, 9(6): 14594-14617. doi: 10.3934/math.2024710 |
[2] | Zhi-Yu Shi, Jia-Bao Liu . Topological indices of linear crossed phenylenes with respect to their Laplacian and normalized Laplacian spectrum. AIMS Mathematics, 2024, 9(3): 5431-5450. doi: 10.3934/math.2024262 |
[3] | Yinzhen Mei, Chengxiao Guo . The minimal degree Kirchhoff index of bicyclic graphs. AIMS Mathematics, 2024, 9(7): 19822-19842. doi: 10.3934/math.2024968 |
[4] | Ali Al Khabyah . Mathematical aspects and topological properties of two chemical networks. AIMS Mathematics, 2023, 8(2): 4666-4681. doi: 10.3934/math.2023230 |
[5] | Yuan Shan, Baoqing Liu . Existence and multiplicity of solutions for generalized asymptotically linear Schrödinger-Kirchhoff equations. AIMS Mathematics, 2021, 6(6): 6160-6170. doi: 10.3934/math.2021361 |
[6] | Giovanni G. Soares, Ernesto Estrada . Navigational bottlenecks in nonconservative diffusion dynamics on networks. AIMS Mathematics, 2024, 9(9): 24297-24325. doi: 10.3934/math.20241182 |
[7] | Zenan Du, Lihua You, Hechao Liu, Yufei Huang . The Sombor index and coindex of two-trees. AIMS Mathematics, 2023, 8(8): 18982-18994. doi: 10.3934/math.2023967 |
[8] | Ali N. A. Koam, Ali Ahmad, Yasir Ahmad . Computation of reverse degree-based topological indices of hex-derived networks. AIMS Mathematics, 2021, 6(10): 11330-11345. doi: 10.3934/math.2021658 |
[9] | Jahfar T K, Chithra A V . Central vertex join and central edge join of two graphs. AIMS Mathematics, 2020, 5(6): 7214-7233. doi: 10.3934/math.2020461 |
[10] | Yun-Ho Kim, Hyeon Yeol Na . Multiplicity of solutions to non-local problems of Kirchhoff type involving Hardy potential. AIMS Mathematics, 2023, 8(11): 26896-26921. doi: 10.3934/math.20231377 |
This paper examines if clean energy stocks help investors in managing carbon risk. We use the price of the European Union Allowance (EUA) and European clean energy index (ERIX) for the three phases of the EU-Emission Trading Scheme. Analyzing the time-varying correlation and volatility of EUA stock and ERIX through generalized orthogonal GO-GARCH model, the empirical results reveal relative independence of the European renewable energy market from the carbon market providing diversification benefits and value addition by including carbon assets in clean energy stock portfolio. Furthermore, three portfolios with different weight allocation strategies reveal that the carbon asset provides risk and downside risk benefits when mixed with a clean energy stock portfolio. These results are useful for investors who enter the market for value maximization and the regulators striving to make strategies for managing carbon risk.
Graph categories considered in this study are simple, finite, and linked. Allow the graph G to be made up of VG and edge set EG, i.e., G=(VG,EG). For more graph notations, readers should refer to [1].
If and only if two neighbouring verticies i and j of G, the adjacency matrix A(G)=(aij) is a (0, 1)-matrix, aij=1. The diagonal degree matrix of G is
DG=diag(d1,d2,⋯,dn), |
where di represents the degree of vertex i in G. The difference between the degree matrix DG and adjacency matrix AG of G gives rise to the Laplace matrix, denoted as LG. The normalized Laplacian [2] is defined as
L(G)=I−D(G)12(D(G)−1A(G))D(G)−12=D(G)−12L(G)D(G)−12. |
The (m,n)th-entry of L(G), which is designated
(L(G))mn={1,m=n;−1√dmdn,m≠n,vmisadjacenttovn;0,otherwise. | (1.1) |
The distance between both vi and vj, known as dij=dG(vi,vj), represents the length of the smallest path in question. Wiener and Dobrynin [3,4] introduced the Wiener index for the first time in 1947. In addition, the Wiener index is denoted as
W(G)=∑i<jdij. |
For further information on the Wiener index, please refer to [5,6,7,8,9].
The Gutman index of a simple graph G is introduced [10] and denoted as
Gut(G)=∑i<jdidjdij, |
taking into account the degree di of vertex vi.
The Kirchhoff index[11,12] characterizes graph G by summing the resistance distances between every pair of vertices, similar to the Wiener index, namely
Kf(G)=∑i<jrij. |
The multiplicative degree-Kirchhoff index, initially introduced by Chen and Zhang[13] in 2007. It is an extension of the traditional Kirchhoff index, which is expressed as
Kf∗(G)=∑i<jdidjrij. |
The techniques of the Kirchhoff index and multiplication degree-Kirchhoff index can be found in [14,15,16,17,18]. The multiplication degree-Kirchhoff index has garnered significant attention due to its remarkable contributions in academia and practical applications in computer network science, epidemiology, social economics, and other fields. Further research results on the Kirchhoff index and multiplication degree-Kirchhoff index can be explored through [19,20,21,22,23].
The spanning tree of a graph G, also known as complexity, denoted as τ(G), refers to the number of subgraphs that encompass all vertices in G. This measure serves as a crucial indicator for network stability and plays a significant role in assessing the structural characteristics of graphs. For further insights into related topics such as the count of spanning trees, interested readers are encouraged to consult [24,25,26].
With the rapid advancement of scientific research and the successful application of topology in practical scenarios, topological theory has gained increasing recognition worldwide. The calculation problem concerning the phenylene Wiener index has been effectively resolved by Pavlović and Gutman [27]. Chen and Zhang [28] have developed a precise equation for predicting the Wiener index of random phenylene chains. Additionally, Liu et al. [29] have identified both the degree-Kirchhoff index and the number of spanning trees for Ln dicyclobutadieno derivatives of [n] phenylenes.
Given two automorphic graphs S and K, we define the symbol S⊠K to represent the strong product of these two graphs with V(S)×V(K), which is commonly referred to as the strong product in graph theory literature. Readers can refer to [30] for more comprehensive definitions and concepts. Recently, Pan et al.[25] utilized the resistance distance of a strong prism formed by Pn and Cn to determine the Kirchhoff index. Similarly, Li et al.[31] derived graph invariants and spanning trees from the strong prism of a star Sn. Motivated by [30,31,32,33], we obtain the pentagonal network Rn and its strong product P2n. The pentagonal network consists of numerous adjacent pentagons and quadrilaterals with each quadrangle having a maximum of two non-adjacent pentagons, as shown in Figure 1. The P2n is the strong product of Rn, as depicted in Figure 2. It obviously that
|V(P2n)|=14nand|E(P2n)|=47n−8. |
In this paper, we focus on the strong product of pentagonal networks, specifically examining the graph P2n with n≥1. The subsequent sections are organized as follows: Section 2 provides a comprehensive review of relevant research materials, presenting illustrations, concepts and lemmas. In Section 3, we derive the normalized Laplacian spectrum and present an explicit closed formula for the multiplicative degree-Kirchhoff index. Additionally, we calculate the complexity of P2n. In Section 4, we conclude the paper.
In this section, let Rn represent the penagonal-quadrilateral networks, as illustrated in Figure 1. P2n is composed of Rn and its copy R′n, positioned one in front and one behind, as shown in Figure 2. Moreover,
ΦA(x)=det(xIn−A) |
represents the characteristic polynomial of matrix A.
The fact that
π=(1o,2o,⋯,no)(1,1′)(2,2′)⋯((3n),(3n)′) |
is an automorphism deserves attention. Let
V1={1o,2o,⋯,no,u1,u2,⋯,u3n,v1,⋯,v3n}, |
V2={1′o,2′o,⋯,n′o,u′1,u′2,⋯,u′3n,v′1,⋯,v′3n}, |
|V(P2n)|=14nand|E(P2n)|=47n−8. |
Subsequently, the normalized Laplacian matrix can be represented as a block matrix, that is
L(P2n)=(LV1V1LV1V2LV2V1LV2V2), |
in which
LV1V1=LV2V2,LV1V2=LV2V1. |
Let
W=(1√2I6n1√2I6n1√2I6n−1√2I6n), |
then,
WL(P2n)W′=(LA00LS), |
where
LA=LV1V1+LV1V2andLS=LV1V1−LV1V2. |
Observe that W′ and W are transpose matrices of each other.
The characteristic polynomial of the matrix R, is denoted as
Φ(R):=det(xI−R). |
The decomposition theorem process for P2n is obtained in a similar manner to Pan and Li[34], thus we omit this proof and present it as follows:
Lemma 2.1. [35] Assuming that the determination of LA and LS has been previously described, then
ΦL(Ln)(x)=ΦLA(x)⋅ΦLS(x). |
Lemma 2.2. [13] The graph G is an undirected connected graph with n vertices and m edges. Then
Kf∗(G)=2mn∑k=21λk. |
Lemma 2.3. [2] The number of spanning trees in G, referred to as the graph's complexity, can be considered a fundamental measure in graph theory. Then
τ(G)=12mn∏i=1di⋅n∏j=2λj. |
In this section, we explore the methodology for deriving an explicit analytical expression of the multiplicative Kirchhoff index by traversing the normalized Laplacian matrix. Meanwhile, we determine the computational complexity of P2n. Subsequently, employing the normalized Laplacian, we derive matrices of order 6n as
LV1V1=(1−1√350⋯00−1500⋯00−1√351−17⋯000−170⋯000−171⋯0000−17⋯00⋮⋮⋮⋱⋮⋮⋮⋮⋮⋱⋮⋮000⋯1−1√35000⋯−170000⋯−1√351000⋯0−15−1500⋯001−1√350⋯000−170⋯00−1√351−17⋯0000−17⋯000−171⋯00⋮⋮⋮⋱⋮⋮⋮⋮⋮⋱⋮⋮000⋯−170000⋯1−1√35000⋯0−15000⋯−1√351) |
and
LV1V2=(−15−1√350⋯00−1500⋯00−1√35−17−17⋯000−170⋯000−17−17⋯0000−17⋯00⋮⋮⋮⋱⋮⋮⋮⋮⋮⋱⋮⋮000⋯−17−1√35000⋯−170000⋯−1√35−15000⋯0−15−1500⋯00−15−1√350⋯000−170⋯00−1√35−17−17⋯0000−17⋯000−17−17⋯00⋮⋮⋮⋱⋮⋮⋮⋮⋮⋱⋮⋮000⋯−170000⋯−17−1√35000⋯0−15000⋯−1√35−15). |
Due to
LA=LV1V1(P2n)+LV1V2(P2n) |
and
LS=LV1V1(P2n)−LV1V2(P2n), |
it can be convincingly argued that
LA=2(25−1√350⋯00−1500⋯00−1√3537−17⋯000−170⋯000−1737⋯00000⋯00⋮⋮⋮⋱⋮⋮⋮⋮⋮⋱⋮⋮000⋯37−1√35000⋯−170000⋯−1√3525000⋯0−15−1500⋯0025−1√350⋯000−170⋯00−1√3537−17⋯0000−17⋯000−1737⋯00⋮⋮⋮⋱⋮⋮⋮⋮⋮⋱⋮⋮000⋯−170000⋯37−1√35000⋯0−15000⋯−1√3525) |
and
LS=diag(65,87,87,⋅⋅⋅,87,65,65,87,87,⋅⋅⋅,87,65). |
Utilizing Lemma 2.1, it is revealed that the P2n normalized Laplacian spectrum consists of the eigenvalues from LA and LS. It is established that the LS possesses eigenvalues 65 and 87 with multiplicities of 4 and (6n−4), respectively.
Let
M=(25−1√3500⋯000−1√3537−170⋯0000−1737−17⋯00000−1737⋯000⋮⋮⋮⋮⋱⋮⋮⋮0000⋯37−1700000⋯−1737−1√350000⋯0−1√3525)(3n)×(3n) |
and
N=diag(−15,−17,−17,−17,−17,⋯,−17,−17,−15), |
where the matrices M and N are both of order 3n.
The matrices M and N are combined to form a block matrix, denoted as 12LA, in the following manner:
12LA=(MNNM). |
Suppose that
W=(1√2I3n1√2I3n1√2I3n−1√2I3n) |
is a block matrix. Hence, we can obtain
W(12LA)W′=(M+N00M−N). |
Let J=M+N and K=M−N. Then,
J=(15−1√3500⋯000−1√3527−170⋯0000−1727−17⋯00000−1727⋯000⋮⋮⋮⋮⋱⋮⋮⋮0000⋯27−1700000⋯−1747−1√350000⋯0−1√3535)(3n)×(3n) |
and
K=(35−1√3500⋯000−1√3547−170⋯0000−1747−17⋯00000−1747⋯000⋮⋮⋮⋮⋱⋮⋮⋮0000⋯47−1700000⋯−1747−1√350000⋯0−1√3535)(3n)×(3n), |
in which the diagonal elements are
(35,47,47,47,⋅⋅⋅47,47,35). |
Based on Lemma 2.1, it is evident and demonstrable that the eigenvalues of 12LA are identical to those of J and K. Assume that the eigenvalues of J and K are σi and ςj(i,j=1,2,⋯,3n), respectively, with
σ1≤σ2≤σ3≤⋯≤σ3n,ς1≤ς2≤ς3≤⋯≤ς3n. |
We verify σ1≥0 and ς1≥0. In addition, it is easy to know that the normalized Laplacian spectrum of P2n is {2σ1,2σ2,⋯,2σ3n,2ς1,2ς2,⋯,2ς3n}. Note that
|E(P2n)|=47n−8, |
we can obtain Lemma 3.1 according to Lemma 2.2.
Lemma 3.1. Assume that P2n is the strong product of the pentagonal network. Then,
Kf∗(P2n)=2(47n−8)(2×56+(6n−2)78+123n∑i=21σi+123n∑j=11ςj)=(47n−8)(63n−16+3n∑i=21σi+3n∑j=11ςj). |
Subsequently, we partition the computation of the aforementioned equation into two distinct components and prioritize the initial calculation of ∑3ni=21σi.
Lemma 3.2. Suppose that σi(i=1,2,⋯,3n) is defined as described previously.
3n∑i=21σi=1035n3+142n2+617n2(81n+490). |
Proof. Suppose that
Φ(J)=x3n+a1x3n−1+⋯+a3nx2+a3n+1x=x(x3n−1+a1x3n−2+⋯+a3n−1x+a3n−2). |
Then σ2,σ3,⋯,σ3n fulfil the following equation
x3n−1+a1x3n−2+⋯+a3n−2x+a3n−1=0, |
and we observe that 1σ2,1σ3,⋯,1σ3n are the roots of the following equation
a3n−1x3n−1+a3n−2x3n−2+⋯+a1x+1=0. |
By Vieta's Theorem, one has
3n∑i=21σi=(−1)3n−2a3n−2(−1)3n−1a3n−1. | (3.1) |
For each value of i from 1 to 3n+1, we consider Ji and assign ji as the determinant of Ji. We will derive the formula for ji, which can be utilized to calculate (−1)3n−2a3n−2 and (−1)3n−1a3n−1. Then one has
j1=15,j2=135,j3=1245,j4=11715,j5=112005,j6=184035,j7=1588245,j8=14117715, |
and
{j3i=27j3i−1−149j3i−2,1≤i≤n;j3i+1=27j3i−149j3i−1,0≤i≤n−1;j3i+2=27j3i+1−149j3i,0≤i≤n−1. |
Through a straightforward computation, one can derive the following general formulas:
{j3i=75⋅(1343)i,1≤i≤n;j3i+1=15⋅(1343)i,0≤i≤n−1;j3i+2=135⋅(1343)i,0≤i≤n−1. | (3.2) |
According to Eq (3.1), we divide the numerator and denominator into two facts and reveal them later. For the sake of convenience, we represent the diagonal elements of J as lii in a simplified manner.
Fact 3.3.
(−1)3n−1a3n−1=490+81n25(1343)n. |
Proof. Given that J is a left-right symmetric matrix, the sum of all principal minors can be represented by the number (−1)3n−1a3n−1, where these minors correspond to rows and columns with indices equal to (3n−1).
(−1)3n−1a3n−1=3n∑i=1detLA[i]=3n∑i=1det(Ji−100J3n−i)=3n∑i=1ji−1⋅j3n−i, | (3.3) |
where
J3n−i=(li+1,i+1⋯00⋮⋮⋱⋮0⋯l3n−1,3n−1−1√350⋯−1√35l3n,3n). |
By Eqs (3.2) and (3.3), we have
(−1)3n−1a3n−1=2j3n+1+n∑l=1j3(l−1)+2⋅j3(n−l)+2+n∑l=1j3l⋅j3(n−l)+1+n−1∑l=0j3l+1⋅j3(n−l)=985⋅(1343)n+343n352⋅(1343)n+7n25⋅(1343)n+n⋅(1343)n=490+81n25(1343)n. |
This is the completion of the proof.
Fact 3.4.
(−1)3n−2a3n−2=1035n3+142n2−617n50(1343)n. |
Proof. We note that the sum of all principal minors of order 3n−2 in J can be expressed as (−1)3n−2a3n−2. After that,
(−1)3n−2a3n−2=3n∑1≤i<j|Ji−1000Z000J3n−j|,1≤i<j≤3n−2, |
where
Z=(li+1,i+1⋯0⋮⋱⋮0⋯lj−1,j−1) |
and
J3n−j=(lj+1,j+1⋯00⋮⋱⋮⋮0⋯l3n−1,3n−1−1√350⋯−1√35l3n,3n). |
Note that
(−1)3n−2a3n−2=3n∑1≤i<jdetJi−1⋅detZ⋅detJ3n−j=3n∑1≤i<jdetZ⋅si−1⋅s3n−j. | (3.4) |
According to Eq (3.4), the determinant of Z varies depending on the values of i and j, as well as s and t. Consequently, we can categorize the primary scenarios into six distinct classifications.
Case 1. i=3s, j=3t for 1≤s<t≤n, and
detZ1=|27−1700⋯00−1727−170⋯000−1727−17⋯0000−1727⋯00⋮⋮⋮⋮⋱⋮⋮0000⋯27−170000⋯−1727|(3s−3t−1)=21(s−t)(1343)s−t. |
Case 2. i=3s, j=3t+1 for 1≤s≤t≤n, and
detZ2=|27−1700⋯00−1727−170⋯000−1727−17⋯0000−1727⋯00⋮⋮⋮⋮⋱⋮⋮0000⋯27−170000⋯−1727|(3s−3t)=[3(s−t)+1](1343)s−t, |
or i=3s+2, j=3t for 0≤s<t≤n, and
detZ3=|27−1700⋯00−1727−170⋯000−1727−17⋯0000−1727⋯00⋮⋮⋮⋮⋱⋮⋮0000⋯27−170000⋯−1727|(3s−3t−3)=(3s−3t−2)(1245)s−t−1. |
Case 3. In the same way, i=3s, j=3t+2 for 1≤s≤t≤n, and
detZ4=|27−1700⋯00−1727−170⋯000−1727−17⋯0000−1727⋯00⋮⋮⋮⋮⋱⋮⋮0000⋯27−1√350000⋯−1√3515|(3s−3t+1)=49(3s−3t+2)(1343)s−t+1, |
or i=3s+1, j=3t for 0≤s<t≤n, and
detZ5=|27−1700⋯00−1727−170⋯000−1727−17⋯0000−1727⋯00⋮⋮⋮⋮⋱⋮⋮0000⋯27−170000⋯−1727|(3s−3t−2)=49(3s−3t−1)(1343)s−t−1. |
Case 4. Similarly, i=3s+1, j=3t+1 for 0≤s<t≤n, and
detZ6=|27−1700⋯00−1727−170⋯000−1727−17⋯0000−1727⋯00⋮⋮⋮⋮⋱⋮⋮0000⋯27−170000⋯−1727|(3s−3t−1)=21(s−t)(1343)s−t, |
or i=3s+2, j=3t+2 for 0≤k<l≤n, and
detZ7=|27−1700⋯00−1727−170⋯000−1727−17⋯0000−1727⋯00⋮⋮⋮⋮⋱⋮⋮0000⋯27−1√350000⋯−1√3515|(3s−3t−1)=detZ6. |
Case 5. i=3s+1, j=3t+2 for 0≤s≤t≤n, and
detZ8=|27−1700⋯00−1727−170⋯000−1727−17⋯0000−1727⋯00⋮⋮⋮⋮⋱⋮⋮0000⋯25−1√350000⋯−1√3527|(3s−3t)=(3s−3t+1)(1343)s−t. |
Case 6. i=3s+2, j=3t+1 for 0≤k<l≤n, and
detZ9=|27−1700⋯00−1727−170⋯000−1727−17⋯0000−1727⋯00⋮⋮⋮⋮⋱⋮⋮0000⋯27−170000⋯−1727|(3s−3t−2)=49(3s−3t−1)(1343)s−t. |
Therefore, we can obtain
(−1)3n−2a3n−2=∑1≤i<j≤3ndetZ⋅jk−1⋅j3n−l=℘1+℘2+℘3, | (3.5) |
where
℘1=∑1≤s<t≤ndetZ1⋅j3s−1⋅j3n−3t+∑1≤s≤t≤ndetZ2⋅j3s−1⋅j3n−3t−1+∑1≤s≤t≤n−1detZ4⋅j3s−1⋅j3n−3t+∑1≤s≤ndetJ[3s,3n+2]⋅j3s−1=7n(n2−1)50(1343)n+(n2+2n)(n+1)2450(1343)n+n2(n−1)2450(1343)n+n(3n+1)490(1343)n=345n3+17n2−336n50(1343)n, |
℘2=∑1≤s<t≤ndetZ5⋅j3s⋅j3n−3t+2+∑1≤s<t≤ndetZ6⋅j3s⋅j3n−3t+1+∑1≤s≤t≤n−1detZ8⋅j3s⋅j3n−3t+∑1≤s≤ndetJ[3s+1,3n+2]⋅j3s+∑1≤s≤ndetJ[1,3t]⋅j3n−3t+2+∑1≤t≤ndetS[1,3t+1]⋅s3n−3t+1+∑0≤t≤n−1detJ[1,3t+2]⋅j3n−3t+detJ[1,3n+2]=4950(n3−n2−2n+2)(1343)n−1+49n(n2−1)50(1343)n+49n(n2−1)50(1343)n+7n(3n−1)10(1343)n+7n(3n+1)10(1343)n−1+21n(n+1)10(1343)n+7n(3n−1)10(1343)n+n(3n+1)(1343)n=345n3+73n2−92n50(1343)n |
and
℘3=∑0≤s<t≤ndetZ3⋅j3s+1⋅j3n−3t+2+∑0≤s<t≤n−1detZ7⋅j3s+1⋅j3n−3t+∑0≤s<t≤ndetZ9⋅j3s+1⋅j3n−3t+1+∑0≤s≤n−1detJ[3s+2,3n+2]⋅j3s+1=49n(n2+n−4)50(1343)n−1+n(n2−1)50(1343)n+49n(n2+2n−1)50(1343)n+21n(n+1)10(1343)n−1=345n3+52n2−189n50(1343)n. |
By substituting ℘1, ℘2, and ℘3 into Eq (3.5), the desired outcome can be deduced.
(−1)3n−2a3n−2=℘1+℘2+℘3=1035n3+142n2−617n50(1343)n. |
This is the completion of the proof.
0=ς1<ς2≤ς3≤⋯≤ς3n |
represent the eigenvalues of J. By Facts 3.3 and 3.4, we can further investigate Lemma 3.2. According to Eq (3.1), it is evident that
3n∑i=21σi=(−1)3n−2a3n−2(−1)3n−1a3n−1=1035n3+142n2+617n2(490+81n). |
Considering Lemma 3.1, we will focus on the calculations of 3n∑j=11ςj. Let
δ(n)=10290n(11+2√30)+3600+12430√3060000 |
and
ξ(n)=10290n(11−2√30)+3600−12430√3060000. |
Lemma 3.5. The variable ςj (where j ranges from 1 to 3n+2) is assumed to be defined as previously described. One has
3n∑j=11ςj=(−1)3n−1b3n−1detK, |
where
detK=45+11√30125(11+4√30343)n+45−11√30125(11−4√30343)n |
and
(−1)3n−1b3n−1=δ(n)(11+4√30343)n−ξ(n)(11−4√30343)n. |
Proof. The representation of Φ(K) can be expressed as
y3n+b1y3n−1+⋅⋅⋅+b3n−2y2+b3n−1y=y(y3n+1+b1y3n+⋅⋅⋅+b3n−2y+b3n−1), |
where ς1,ς2,⋯,ς3n represent the roots of the equation.
y3n−1+b1y3n−2+⋯+b3n−2y+b3n−1=0, |
and the equation is determined to possess 1ς1,1ς2,⋯,1ς3n as its solutions
b3n−1y3n−1+b3n−2y3n−2+⋯+b1y+1=0. |
By Vieta's Theorem, one holds that
3n∑j=11ςj=(−1)3n−1b3n−1detK. | (3.6) |
To simplify the analysis, let Rp denote the p-th order principal minors of matrix K and kp represent the determinant of Rp. We will derive an equation for kp that can be utilized to calculate (−1)3n−1b3n−1 and detK for values of p ranging from 1 to 3n−1. Subsequently, we arrive at
k1=15,k2=135,k3=1245,k4=11715,k5=112005,k6=184035,k7=1588245,k8=14117715 |
and
{k3p=27k3p−1−149k3p−2,1≤p≤n;k3p+1=27k3p−149k3p−1,1≤p≤n;k3p+2=27k3p+1−149k3p,0≤p≤n−1. |
After a straightforward computation, the following general formulas can be derived
{k3p=105+14√30150⋅(11+4√30343)p+105−14√30150⋅(11−4√30343)p,1≤p≤n;k3p+1=45+8√30150⋅(11+4√30343)p+45−8√30150⋅(11−4√30343)p,1≤p≤n;k3p+2=11+2√3070⋅(11+4√30343)p+11−2√3070⋅(11−4√30343)p,0≤p≤n−1. |
Fact 3.6.
detK=45+11√30125(11+4√30343)n+45−11√30125(11−4√30343)n. |
Proof. By expanding detK with respect to its last row, we obtain
detK=35k3n+1−135k3n=45+11√30125(11+4√30343)n+45−11√30125(11−4√30343)n. |
The desired outcome has been achieved.
Fact 3.7.
(−1)3n−1b3n−1=δ(n)(11+4√30343)n−ξ(n)(11−4√30343)n, |
where
δ(n)=10290n(11+4√30)+3600+12430√3060000 |
and
ξ(n)=10290n(11−4√30)+3600−12430√3060000. |
Proof. Considering that K has a (3n−1)-row and (3n−1)-column structure, the sum of all principal minors can be expressed as (−1)3n−1b3n−1. Here, lii represents the diagonal entries of K. It is noteworthy that K exhibits bilateral symmetry, which enables us to derive specific information.
(−1)3n−1b3n−1=3n∑i=1detK[i]=3n∑i=1det(Ri−100R3n−i)=3n∑i=1ri−1⋅r3n−i, | (3.7) |
where
R3n−i=(gi+1,i+1⋯00⋮⋮⋱⋮0⋯g3n−1,3n−1−1√350⋯−1√35g3n,3n). |
In line with Eq (3.7), we have
(−1)3n−1b3n−1=2k3n−1+n∑l=1detK[3l]+n−1∑l=0detK[3l+1]+n−1∑l=0detK[3l+2]=2k3n−1+n∑l=1k3(l−1)+2⋅k3(n−l)+2+n−1∑l=0k3l⋅k3(n−l)+1+n−1∑l=0k3l+1⋅k3(n−l). | (3.8) |
Immediately, we can get
n∑l=1k3(l−1)+2⋅k3(n−l)+2=245n20(11+4√30343)n+1−(11−4√30343)n+1+√30600(11+4√30343)n−√30600(11−4√30343)n, | (3.9) |
n∑l=0k3l⋅k3(n−l)+1=2391n300(11+4√30343)n+1−(11−4√30343)n+1+137√3060000(11+4√30343)n−137√3060000(11−4√30343)n, | (3.10) |
n−1∑l=0k3l+1⋅k3(n−l)=2391n300(11+4√30343)n+1−(11−4√30343)n+1+1271√3060000(11+4√30343)n−1271√3060000(11−4√30343)n | (3.11) |
and
2k3n−1=90+16√3075(11+4√30343)n+90−16√3075(11−4√30343)n. | (3.12) |
By incorporating Eqs (3.9)–(3.12) into Eq (3.8), we can achieve Lemma 3.7. Utilizing Eq (3.6), in conjunction with Facts 3.6 and 3.7, Lemma 3.5 can be promptly derived.
By integrating Lemmas 3.1, 3.2 and 3.5, we can readily deduce the Theorems 3.8 and 3.9.
Theorem 3.8. Assume that P2n is the strong product of the pentagonal network. One has
Kf∗(P2n)=1029n3+5736n2+4275n+8503+(57n+30)(−1)3n−1b3n−1detK, |
where
(−1)3n−1b3n−1=δ(n)(11+4√30343)n−ξ(n)(11−4√30343)n,detK=45+11√30125(11+4√30343)n+45−11√30125(11−4√30343)n, |
additionally,
δ(n)=10290n(11+2√30)+3600+12430√3060000 |
and
ξ(n)=10290n(11−2√30)+3600−12430√3060000. |
Theorem 3.9. Let P2n be the strong product of pentagonal network. Then
τ(P2n)=35⋅232n+75((45+11√30)(11+4√30)n+(45−11√30)(11−4√30)n). |
Proof. Based on the proof of Lemma 2.2, it is evident that σ1,σ2,⋯σ2n+1 constitute the roots of the equation
x2n+a1x2n−1+⋯+a2n−1x+a2n=0. |
Accordingly, one has
3n∏i=2σi=(−1)3n−1a3n−1. |
By Fact 3.3, we have
3n∏i=2σi=12+23n25(1343)n. |
By the same method,
3n∏j=1ςj=detK=45+11√30375(11+4√30343)n+45−11√30375((11−4√30)343)n. |
Note that
∏v∈VP2nd(P2n)=54⋅716n−4 |
and
|E(P2n)|=48n−7. |
In conjunction with Lemma 2.3, we get
τ(P2n)=12|E(P2n)|((65)2⋅(87)6n−2⋅3n∏i=22σi⋅3n∏j=12ςj⋅∏v∈VP2nd(P2n))=35⋅232n+75((45+11√30)(11+4√30)n+(45−11√30)(11−4√30)n). |
This is the completion of the proof.
In this study, we have derived explicit expressions for the multiplicative degree-Kirchhoff index and complexity of P2n based on the spectrum of the Laplacian matrix, where P2n=Rn⊠R′n. These two fundamental calculations serve as simple yet reliable graph invariants that effectively capture the stability of diverse networks. Future research should focus on applying our methodology to determine spectra for strong products of automorphic and symmetric networks.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors would like to express their sincere gratitude to the editor and anonymous referees for providing valuable suggestions, which significantly contributed to the substantial improvement of the original manuscript. This work was supported by Anhui Jianzhu University Funding under Grant KJ212005, and by Natural Science Fund of Education Department of Anhui Province under Grant KJ2020A0478.
No potential confilcts of interest were reported by the authors.
[1] | Alexander C (2000) A primer on the orthogonal GARCH model. Manuscript ISMA Centre, University of Reading, UK, 2. |
[2] |
Andersson M, Bolton P, Samama F (2016) Hedging climate risk. Financ Anal J 72: 13-32. doi: 10.2469/faj.v72.n3.4
![]() |
[3] |
Arouri MEH, Lahiani A, Nguyen DK (2015) Cross-market dynamics and optimal portfolio strategies in Latin American equity markets. Eur Bus Rev 27: 161-181. doi: 10.1108/EBR-04-2013-0069
![]() |
[4] | Bai L, Liu Y, Wang Q, et al. (2019) Improving portfolio performance of renewable energy stocks using robust portfolio approach: Evidence from China. Phys A 533: 122059. |
[5] |
Balcılar M, Demirer R, Hammoudeh S, et al. (2016) Risk spillovers across the energy and carbon markets and hedging strategies for carbon risk. Energ Econ 54: 159-172. doi: 10.1016/j.eneco.2015.11.003
![]() |
[6] |
Basher SA, Sadorsky P (2016) Hedging emerging market stock prices with oil, gold, VIX, and bonds: A comparison between DCC, ADCC and GO-GARCH. Energ Econ 54: 235-247. doi: 10.1016/j.eneco.2015.11.022
![]() |
[7] | Benz E, Trück S (2006) CO2 emission allowances trading in Europe—specifying a new class of assets. Prob Perspect Manage 4: 30-40. |
[8] |
Bing X, Bloemhof-Ruwaard J, Chaabane A, et al. (2015) Global reverse supply chain redesign for household plastic waste under the emission trading scheme. J Clean Prod 103: 28-39. doi: 10.1016/j.jclepro.2015.02.019
![]() |
[9] | Bunn D, Fezzi C (2007) Interaction of European carbon trading and energy prices. |
[10] |
Chkili W (2016) Dynamic correlations and hedging effectiveness between gold and stock markets: Evidence for BRICS countries. Res Int Bus Financ 38: 22-34. doi: 10.1016/j.ribaf.2016.03.005
![]() |
[11] | Convery FJ, Redmond L (2007) Market and price developments in the European Union emissions trading scheme. |
[12] | Creti A, Jouvet PA, Mignon V (2012) Carbon price drivers: Phase I versus Phase II equilibrium? Energ Econ 34: 327-334. |
[13] | DeMiguel V, Garlappi L, Uppal R (2009) Optimal versus naive diversification: How inefficient is the 1/N portfolio strategy? Rev Financ Stud 22: 1915-1953. |
[14] | Ding Z (1994) Time series analysis of speculative returns, University of California, San Diego, Department of Economics. |
[15] |
Dutta A (2018) Modeling and forecasting the volatility of carbon emission market: The role of outliers, time-varying jumps and oil price risk. J Clean Prod 172: 2773-2781. doi: 10.1016/j.jclepro.2017.11.135
![]() |
[16] | Dutta A (2019) Impact of carbon emission trading on the European Union biodiesel feedstock market. Biomass Bioenergy 128: 105328. |
[17] |
Dutta A, Bouri E, Noor MH (2018) Return and volatility linkages between CO2 emission and clean energy stock prices. Energy 164: 803-810. doi: 10.1016/j.energy.2018.09.055
![]() |
[18] |
Engle R (2002) Dynamic conditional correlation: A simple class of multivariate generalized autoregressive conditional heteroskedasticity models. J Bus Econ Stat 20: 339-350. doi: 10.1198/073500102288618487
![]() |
[19] | Ghalanos A (2014) The rmgarch models: Background and properties. |
[20] |
Gronwald M, Ketterer J, Trück S (2011) The relationship between carbon, commodity and financial markets: A copula analysis. Econ Record 87: 105-124. doi: 10.1111/j.1475-4932.2011.00748.x
![]() |
[21] |
Jammazi R, Reboredo JC (2016) Dependence and risk management in oil and stock markets. A wavelet-copula analysis. Energy 107: 866-888. doi: 10.1016/j.energy.2016.02.093
![]() |
[22] |
Kara M, Syri S, Lehtilä A, et al. (2008) The impacts of EU CO2 emissions trading on electricity markets and electricity consumers in Finland. Energy Econ 30: 193-211. doi: 10.1016/j.eneco.2006.04.001
![]() |
[23] |
Kazemilari M, Mardani A, Streimikiene D, et al. (2017) An overview of renewable energy companies in stock exchange: Evidence from minimal spanning tree approach. Renewable Energy 102: 107-117. doi: 10.1016/j.renene.2016.10.029
![]() |
[24] |
Keppler JH, Mansanet-Bataller M (2010) Causalities between CO2, electricity, and other energy variables during phase I and phase II of the EU ETS. Energy Policy 38: 3329-3341. doi: 10.1016/j.enpol.2010.02.004
![]() |
[25] |
Koch N, Fuss S, Grosjean G, et al. (2014) Causes of the EU ETS price drop: Recession, CDM, renewable policies or a bit of everything?—New evidence. Energy Policy 73: 676-685. doi: 10.1016/j.enpol.2014.06.024
![]() |
[26] |
Kroner KF, Ng VK (1998) Modeling asymmetric movement of asset prices. Rev Financ Stud 11: 844-871. doi: 10.1093/rfs/11.4.817
![]() |
[27] |
Kumar S, Managi S, Matsuda A (2012) Stock prices of clean energy firms, oil and carbon markets: A vector autoregressive analysis. Energ Econ 34: 215-226. doi: 10.1016/j.eneco.2011.03.002
![]() |
[28] | Lam L, Fung L, Yu IW (2009) Forecasting a large dimensional covariance matrix of a portfolio of different asset classes. Hong Kong Monetary Authority Working Paper, (01). |
[29] |
Lashof DA, Ahuja DR (1990) Relative contributions of greenhouse gas emissions to global warming. Nature 344: 529-531. doi: 10.1038/344529a0
![]() |
[30] |
Leitao J, Ferreira J, Santibanez‐Gonzalez E (2021) Green bonds, sustainable development and environmental policy in the European Union carbon market. Bus Strat Environ 30: 2077-2090. doi: 10.1002/bse.2733
![]() |
[31] |
Lin B, Chen Y (2019) Dynamic linkages and spillover effects between CET market, coal market and stock market of new energy companies: A case of Beijing CET market in China. Energy 172: 1198-1210. doi: 10.1016/j.energy.2019.02.029
![]() |
[32] | Lin B, Jia Z (2020) Is emission trading scheme an opportunity for renewable energy in China? A perspective of ETS revenue redistributions. Appl Energy 263: 114605. |
[33] |
Liu L, Chen C, Zhao Y, et al. (2015) China's carbon-emissions trading: Overview, challenges and future. Renewable Sustainable Energy Rev 49: 254-266. doi: 10.1016/j.rser.2015.04.076
![]() |
[34] |
Luo C, Wu D (2016) Environment and economic risk: An analysis of carbon emission market and portfolio management. Environ Res 149: 297-301. doi: 10.1016/j.envres.2016.02.007
![]() |
[35] |
Marimoutou V, Soury M (2015) Energy markets and CO2 emissions: Analysis by stochastic copula autoregressive model. Energy 88: 417-429. doi: 10.1016/j.energy.2015.05.060
![]() |
[36] |
O'Mahony T (2020) State of the art in carbon taxes: a review of the global conclusions. Green Financ 2: 409-423. doi: 10.3934/GF.2020022
![]() |
[37] |
Papantonis I (2016) Volatility risk premium implications of GARCH option pricing models. Econ Model 58: 104-115. doi: 10.1016/j.econmod.2016.05.004
![]() |
[38] |
Reboredo JC (2013) Modeling EU allowances and oil market interdependence. Implications for portfolio management. Energ Econ 36: 471-480. doi: 10.1016/j.eneco.2012.10.004
![]() |
[39] |
Reboredo JC (2014) Volatility spillovers between the oil market and the European Union carbon emission market. Econ Model 36: 229-234. doi: 10.1016/j.econmod.2013.09.039
![]() |
[40] | Reboredo JC (2015) Is there dependence and systemic risk between oil and renewable energy stock prices? Energ Econ 48: 32-45. |
[41] | Reboredo JC, Ugolini A, Chen Y (2019) Interdependence Between Renewable-Energy and Low-Carbon Stock Prices. Energies 12: 4461. |
[42] |
Sadorsky P (2012) Correlations and volatility spillovers between oil prices and the stock prices of clean energy and technology companies. Energ Econ 34: 248-255. doi: 10.1016/j.eneco.2011.03.006
![]() |
[43] |
Tian Y, Akimov A, Roca E, et al. (2016) Does the carbon market help or hurt the stock price of electricity companies? Further evidence from the European context. J Clean Prod 112: 1619-1626. doi: 10.1016/j.jclepro.2015.07.028
![]() |
[44] |
Uddin N, Holtedahl P (2013) Emission trading schemes-avenues for unified accounting practices. J Clean Prod 52: 46-52. doi: 10.1016/j.jclepro.2013.03.017
![]() |
[45] |
Van der Weide R (2002) GO‐GARCH: a multivariate generalized orthogonal GARCH model. J Appl Econometrics 17: 549-564. doi: 10.1002/jae.688
![]() |
[46] | Wen X, Bouri E, Roubaud D (2017) Can energy commodity futures add to the value of carbon assets? Econ Model 62: 194-206. |
[47] |
Zhang G, Du Z (2017) Co-movements among the stock prices of new energy, high-technology and fossil fuel companies in China. Energy 135: 249-256. doi: 10.1016/j.energy.2017.06.103
![]() |
[48] |
Zhang YJ, Sun YF (2016) The dynamic volatility spillover between European carbon trading market and fossil energy market. J Clean Prod 112: 2654-2663. doi: 10.1016/j.jclepro.2015.09.118
![]() |
[49] | Zhang Y, Liu Z, Yu X (2017) The diversification benefits of including carbon assets in financial portfolios. Sustainability 9: 437. |
1. | Muhammad Shoaib Sardar, Shou-Jun Xu, Resistance Distance and Kirchhoff Index in Windmill Graphs, 2025, 22, 15701794, 159, 10.2174/0115701794299562240606054510 | |
2. | Md. Abdus Sahir, Sk. Md. Abu Nayeem, Characterizing the Degree-Kirchhoff, Gutman, and Schultz Indices in Pentagonal Cylinders and Möbius Chains, 2025, 0278-081X, 10.1007/s00034-025-03130-9 |