Estrada Index and Laplacian Estrada Index of Random Interdependent Graphs

Abstract

Let G be a simple graph of order n. The Estrada index and Laplacian Estrada index of G are defined by $EE\left(G\right)={\sum }_{i=1}^n{e}^{{\lambda }_i\left(A\left(G\right)\right)}$ and $LEE\left(G\right)={\sum }_{i=1}^n{e}^{{\lambda }_i\left(L\left(G\right)\right)}$, where ${\left\{{\lambda }_i\left(A\left(G\right)\right)\right\}}_{i=1}^n$ and ${\left\{{\lambda }_i\left(L\left(G\right)\right)\right\}}_{i=1}^n$ are the eigenvalues of its adjacency and Laplacian matrices, respectively. In this paper, we establish almost sure upper bounds and lower bounds for random interdependent graph model, which is fairly general encompassing Erdös-Rényi random graph, random multipartite graph, and even stochastic block model. Our results unravel the non-triviality of interdependent edges between different constituting subgraphs in spectral property of interdependent graphs.

1. Introduction

We consider a simple graph $G=(V,E)$ on the vertex set $V=\{1,2,\cdots ,n\}$ with $\left|V\right|=n$ and the edge set E consisting of unordered pairs of vertices. The adjacency matrix of G is a $(0,1)$-matrix denoted by $A\left(G\right)=\left(a_{ij}\right)\in {\mathbb{R}}^{n\times n}$, where $a_{ij}=a_{ji}=1$ when i and j are adjacent, and $a_{ij}=a_{ji}=0$ otherwise. The degree of vertex i is $d_i={\sum }_{j=1}^n{a}_{ij}$, i.e., the number of its incident edges. The Laplacian matrix of G is defined to be $L\left(G\right)=D\left(G\right)-A\left(G\right)\in {\mathbb{R}}^{n\times n}$, where the diagonal matrix $D\left(G\right)=\mathrm{diag}(d_1,d_2,\cdots ,d_n)$ is called the degree matrix of G. Since G is undirected, both $A\left(G\right)$ and $L\left(G\right)$ are symmetric. Moreover, it follows from algebraic graph theory (e.g., [1]) that $A\left(G\right)$ has n real eigenvalues arranged in the non-increasing order ${\lambda }_1\left(A\left(G\right)\right)\ge {\lambda }_2\left(A\left(G\right)\right)\ge \cdots \ge {\lambda }_n\left(A\left(G\right)\right)$, and $L\left(G\right)$ has n real and nonnegative eigenvalues ordered non-increasingly as ${\lambda }_1\left(L\left(G\right)\right)\ge {\lambda }_2\left(L\left(G\right)\right)\ge \cdots \ge {\lambda }_n\left(L\left(G\right)\right)=0$.

The Estrada index of G, given by

$$EE\left(G\right)=\sum _{i=1}^n{e}^{{\lambda }_i\left(A\left(G\right)\right)},$$

is a graph spectral invariant introduced by Estrada [2] in the year 2000. It has notable applications in biochemistry and complex networks including quantifying the degree of folding of long-chain molecules [3,4,5] and network resilience [6,7]. An important variant of Estrada index is the Laplacian Estrada index [8], which is defined by invoking Laplacian eigenvalues as

$$LEE\left(G\right)=\sum _{i=1}^n{e}^{{\lambda }_i\left(L\left(G\right)\right)}.$$

A variety of mathematical properties, including upper and lower bounds, of these metric have been investigated; see e.g., [9,10,11,12,13] and references therein.

In addition to fixed graphs, Estrada and Laplacian Estrada indices have been recently investigated for classical Erdös-Rényi random graph model as well as random multipartite graphs [14,15,16,17]. These results are noteworthy in the sense that they not only contribute to the understanding of spectral theory of random networks but also presenting estimates to $EE$ and $LEE$ for almost all graphs (as the number of vertices goes to infinity), which are typically much sharper than previous bounds for fixed graphs. Important ramifications for distance Estrada index [18] and Gaussian Estrada index [19,20] have also been explored lately.

In this paper, we study the Estrada index and Laplacian Estrada index for the class of random interdependent graphs, which consist of m subgraphs with edges between different subgraphs appearing independently with probability p, where $p\in (0,1)$ is a constant. Formally, an interdependent graph $G_{m,n}=(G_1,G_2,\cdots ,G_m,G)$ with $m=m\left(n\right)\ge 2$ is defined by a family of subgraphs $G_l=(V_l,E_l)$ for $l=1,2,\cdots ,m$ with $|V_l|=n_l$, $n_1+n_2+\cdots +n_m=n$, and an m-partite graph $G=(V_1\cup V_2\cup \cdots \cup V_m,E)$ with E containing edges of unordered pairs between $V_s$ and $V_t$ with $s\ne t$. Here, G and $G_l$$(l=1,\cdots ,m)$ are deterministic. A random interdependent graph, denoted by $G_{m,n}\left(p\right)=(G_1,G_2,\cdots ,G_m,G\left(p\right))$, is a graph from the probability space $({\mathcal{S}}_n,{\mathcal{A}}_n,{\mathbb{P}}_n)$ with ${\mathcal{S}}_n$ containing all possible interdependent graphs $G_{m,n}$, $\sigma$-algebra ${\mathcal{A}}_n$ being the power set of ${\mathcal{S}}_n$, and probability measure ${\mathbb{P}}_n$ defining the probability of each graph $G_{m,n}$ in ${\mathcal{S}}_n$ by assigning probability p for each possible edge of G independently.

Note that the random interdependent graph $G_{m,n}\left(p\right)$ is very general in that the random m-partite graph $G\left(p\right)$ is independent of $G_l$ ($l=1,2,\cdots ,m$) and no assumption is made on the topologies inside these fixed subgraphs ${\left\{G_l\right\}}_{l=1}^m$. It is easy to see that Erdös-Rényi random graph $G_n\left(p\right)$ is a special case with each subgraph $G_l$ being a single vertex, and random multipartite graph is also a special case of $G_{m,n}\left(p\right)$ with each $G_l$ being an empty graph. Moreover, the stochastic block model extensively utilized in statistics and machine learning literature (see e.g., [21,22]) can also be viewed as a special case by setting different edge probabilities in each subgraph.

The rest of the paper is organized as follows. In Section 2, we examine Estrada index and derive some upper and lower bounds for random interdependent graphs. In Section 3, we investigate the upper and lower bounds of Laplacian Estrada index for random interdependent graphs. Under certain conditions, the exact estimates for both $EE\left(G_{m,n}\left(p\right)\right)$ and $LEE\left(G_{m,n}\left(p\right)\right)$ are obtained. We conclude the paper in Section 4.

2. Estrada Index of Random Interdependent Graphs

In this section, we estimate the Estrada index $EE\left(G_{m,n}\left(p\right)\right)$ of random interdependent graphs. Recall that $n_1,n_2,\cdots ,n_m$ are the orders of subgraphs $G_1,G_2,\cdots ,G_m$. We re-label these subgraphs by ${\{n_{\left(l\right)}\}}_{l=1}^m$ so that $n_{\left(1\right)}\ge n_{\left(2\right)}\ge \cdots \ge n_{\left(m\right)}$. Standard Landau asymptotic notations will be used here. For example, for two functions $f\left(n\right)$ and $g\left(n\right)$, $f\left(n\right)=o\left(g\right(n\left)\right)$ implies that ${lim}_{n\to \infty }f\left(n\right)/g\left(n\right)=0$; $f\left(n\right)=O\left(g\right(n\left)\right)$ means that $\left|f\right(n)/g(n\left)\right|\le C$ for some constant C for sufficiently large n. A property $\mathcal{B}$ for a random graph model holds asymptotically almost surely (a.a.s.) if its probability tends to 1, namely, ${\mathbb{P}}_n\left(\mathcal{B}\right|{\mathcal{A}}_n)\to 1$ as $n\to \infty$, where ${\mathcal{A}}_n$ is a flow of algebras, and ${\mathcal{A}}_{k_1}\subseteq {\mathcal{A}}_{k_2}$ for each $k_1<k_2$.

Lemma 1

(Weyl’s inequality [23]). Consider symmetric matrices $X\in {\mathbb{R}}^{n\times n},Y\in {\mathbb{R}}^{n\times n}$ and $Z\in {\mathbb{R}}^{n\times n}$ satisfying $X=Y+Z$. Assume that ${\lambda }_1\left(X\right)\ge {\lambda }_2\left(X\right)\ge \cdots \ge {\lambda }_n\left(X\right)$, ${\lambda }_1\left(Y\right)\ge {\lambda }_2\left(Y\right)\ge \cdots \ge {\lambda }_n\left(Y\right)$, and ${\lambda }_1\left(Z\right)\ge {\lambda }_2\left(Z\right)\ge \cdots \ge {\lambda }_n\left(Z\right)$ are their eigenvalues, respectively. We have

$$\underset{j+k=i+n}{max}\{{\lambda }_j\left(Y\right)+{\lambda }_k\left(Z\right)\}\le {\lambda }_i\left(X\right)\le \underset{j+k=i+1}{min}\{{\lambda }_j\left(Y\right)+{\lambda }_k\left(Z\right)\},$$

for all $i=1,2,\cdots ,n$.

In many applications of Weyl’s inequality, it suffices to use a special result ${\lambda }_i\left(Y\right)+{\lambda }_n\left(Z\right)\le {\lambda }_i\left(X\right)\le {\lambda }_i\left(Y\right)+{\lambda }_1\left(Z\right)$. We will resort to the full power of this lemma in the following sections.

Theorem 1.

Let $G_{m,n}\left(p\right)$ be a random interdependent graph in $({\mathcal{S}}_n,{\mathcal{A}}_n,{\mathbb{P}}_n)$. We have

$$\begin{array}{c}\hfill e^{np-n_{\left(1\right)}p}(e^{O\left(\sqrt{n}\right)}+o\left(1\right))\le EE\left(G_{m,n}\left(p\right)\right)\le e^{np+n_{\left(1\right)}}(e^{O\left(\sqrt{n}\right)}+o\left(1\right))\hspace{1em}a.a.s.\end{array}$$

Proof. 

By the construction of random interdependent graph, the adjacency matrix $A\left(G_{m,n}\left(p\right)\right)$ satisfies the following relations:

$$\begin{array}{c}\hfill A\left(G_{m,n}\left(p\right)\right)=\overline{A}+A\left(G\left(p\right)\right),\end{array}$$

(1)

and

$$\begin{array}{c}\hfill A\left(G_n\left(p\right)\right)=A\left(G\left(p\right)\right)+\tilde{A},\end{array}$$

(2)

where $\overline{A}=\mathrm{diag}(A\left(G_1\right),A\left(G_2\right),\cdots ,A\left(G_m\right))$ and $\tilde{A}=\mathrm{diag}(A\left(G_{n_1}\left(p\right)\right),A\left(G_{n_2}\left(p\right)\right),$$\cdots ,A\left(G_{n_m}\left(p\right)\right))$ are two n-dimensional block diagonal matrices.

We will first estimate $EE\left(G\right(p\left)\right)$ and the eigenvalues of $A\left(G\right(p\left)\right)$. By using (2) and Lemma 1, we have

$${\lambda }_1\left(A\left(G_n\left(p\right)\right)\right)-{\lambda }_1\left(\tilde{A}\right)\le {\lambda }_1\left(A\left(G\left(p\right)\right)\right)\le {\lambda }_1\left(A\left(G_n\left(p\right)\right)\right)-{\lambda }_n\left(\tilde{A}\right).$$

For Erdös-Rényi random graph $G_n\left(p\right)$, it is known that [14] ${\lambda }_1\left(A\left(G_n\left(p\right)\right)\right)=np+O\left(\sqrt{n}\right)$ and ${\lambda }_i\left(A\left(G_n\left(p\right)\right)\right)=O\left(\sqrt{n}\right)$ for $i=2,3,\cdots ,n$ a.a.s. Therefore, we derive

$$np-n_{\left(1\right)}p+O\left(\sqrt{n}\right)\le {\lambda }_1\left(A\left(G\left(p\right)\right)\right)\le np+O\left(\sqrt{n}\right),\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.$$

For ${\lambda }_i\left(A\left(G\left(p\right)\right)\right)$, $i=2,3,\cdots ,n-m+1$, it follows from Lemma 1 that

$$\begin{array}{cc}\hfill {\lambda }_{i+m-1}\left(A\left(G_n\left(p\right)\right)\right)-{\lambda }_m\left(\tilde{A}\right)& ={\lambda }_{i+m-1}\left(A\left(G_n\left(p\right)\right)\right)+{\lambda }_{n-m+1}(-\tilde{A})\hfill \\ \hfill & \le {\lambda }_i\left(A\left(G\left(p\right)\right)\right)\hfill \\ \hfill & \le {\lambda }_i\left(A\left(G_n\left(p\right)\right)\right)+{\lambda }_1(-\tilde{A})={\lambda }_i\left(A\left(G_n\left(p\right)\right)\right)-{\lambda }_n\left(\tilde{A}\right).\hfill \end{array}$$

Similarly, we have ${\lambda }_i\left(A\left(G\left(p\right)\right)\right)=O\left(\sqrt{n}\right)$ a.a.s. for $i=2,3,\cdots ,n-m+1$. A further application of Lemma 1 and decomposition (2) yields

$${\lambda }_i\left(A\left(G_n\left(p\right)\right)\right)-{\lambda }_1\left(\tilde{A}\right)\le {\lambda }_i\left(A\left(G\left(p\right)\right)\right)\le {\lambda }_i\left(A\left(G_n\left(p\right)\right)\right)-{\lambda }_n\left(\tilde{A}\right)$$

for $i=n-m+2,n-m+3,\cdots ,n$. Hence, $-n_{\left(1\right)}p+O\left(\sqrt{n}\right)\le {\lambda }_i\left(A\left(G\left(p\right)\right)\right)\le O\left(\sqrt{n}\right)$ a.a.s. for $i=n-m+2,n-m+3,\cdots ,n$. Combining the above discussion, we arrive at

$$\begin{array}{c}\hfill e^{np-n_{\left(1\right)}p}(e^{O\left(\sqrt{n}\right)}+o\left(1\right))\le EE\left(G\left(p\right)\right)\le e^{np}(e^{O\left(\sqrt{n}\right)}+o\left(1\right))\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.\end{array}$$

(3)

It is well known that Estrada index $EE\left(G\right)$ can be interpreted as weighted sum of closed walks of all lengths in G, and hence it changes increasingly with respect to edge addition. More results on the perturbation of Estrada index have been reported in [24]. In light of (1) and (3), we obtain

$$\begin{array}{c}\hfill e^{np-n_{\left(1\right)}p}(e^{O\left(\sqrt{n}\right)}+o\left(1\right))\le EE\left(G\left(p\right)\right)\le EE\left(G_{m,n}\left(p\right)\right).\end{array}$$

(4)

On the other hand, if $G_l$ ($l=1,2,\cdots ,m$) are complete graphs, the Estrada index $EE\left(G_{m,n}\left(p\right)\right)$ will attain its maximum. In this case, the eigenvalues of $\overline{A}$ satisfy ${\lambda }_1,{\lambda }_2,\cdots ,{\lambda }_m\in \{n_1-1,n_2-1,\cdots ,n_m-1\}$ and ${\lambda }_{m+1}={\lambda }_{m+2}=\cdots ={\lambda }_n=-1$. Therefore, using (1) and Lemma 1, we have

$${\lambda }_1\left(A\left(G_{m,n}\left(p\right)\right)\right)\le {\lambda }_1\left(\overline{A}\right)+{\lambda }_1\left(A\left(G\left(p\right)\right)\right)\le n_{\left(1\right)}-1+np+O\left(\sqrt{n}\right)\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.$$

For ${\lambda }_i\left(A\left(G_{m,n}\left(p\right)\right)\right)$, $i=m+2,m+3,\cdots ,n$, similarly,

$${\lambda }_i\left(A\left(G_{m,n}\left(p\right)\right)\right)\le {\lambda }_{m+1}\left(\overline{A}\right)+{\lambda }_{i-m}\left(A\left(G\left(p\right)\right)\right)\le -1+O\left(\sqrt{n}\right)\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.,$$

and for ${\lambda }_i\left(A\left(G_{m,n}\left(p\right)\right)\right)$, $i=2,3,\cdots ,m+1$, we have

$${\lambda }_i\left(A\left(G_{m,n}\left(p\right)\right)\right)\le {\lambda }_1\left(\overline{A}\right)+{\lambda }_i\left(A\left(G\left(p\right)\right)\right)\le n_{\left(1\right)}-1+O\left(\sqrt{n}\right)\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.$$

It follows from the above comments and the definition of Estrada index that

$$\begin{array}{cc}\hfill EE\left(G_{m,n}\left(p\right)\right)& \le e^{n_{\left(1\right)}+np+O\left(\sqrt{n}\right)}+(n-m-1)e^{O\left(\sqrt{n}\right)}+m{e}^{n_{\left(1\right)}+O\left(\sqrt{n}\right)}\hfill \\ \hfill & =e^{np+n_{\left(1\right)}}(e^{O\left(\sqrt{n}\right)}+o\left(1\right))\hspace{1em}a.a.s.,\hfill \end{array}$$

(5)

where the last equality holds since

$$\begin{array}{c}\hfill \frac{e^{n_{\left(1\right)}+np+O\left(\sqrt{n}\right)}}{e^{np+n_{\left(1\right)}}}+\frac{(n-m-1)e^{O\left(\sqrt{n}\right)}}{e^{np+n_{\left(1\right)}}}+\frac{m{e}^{n_{\left(1\right)}+O\left(\sqrt{n}\right)}}{e^{np+n_{\left(1\right)}}}=e^{O\left(\sqrt{n}\right)}+o\left(1\right)\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.\end{array}$$

The theorem thus follows from (4) and (5). ☐

Corollary 1.

Let $G_{m,n}\left(p\right)$ be a random interdependent graph in $({\mathcal{S}}_n,{\mathcal{A}}_n,{\mathbb{P}}_n)$. If $n_{\left(1\right)}=O\left(\sqrt{n}\right)$, we have

$$\begin{array}{c}\hfill EE\left(G_{m,n}\left(p\right)\right)=e^{np}(e^{O\left(\sqrt{n}\right)}+o\left(1\right))\hspace{1em}a.a.s.\end{array}$$

Proof. 

Recall that p is a constant. If $n_{\left(1\right)}=O\left(\sqrt{n}\right)$, then $EE\left(G_{m,n}\left(p\right)\right)=e^{np}(e^{O\left(\sqrt{n}\right)}+o\left(1\right))$ a.a.s. by using (4) and (5). ☐

Corollary 1 gives an exact estimate for Estrada index of random interdependent graphs. This result reveals that $EE\left(G_{m,n}\left(p\right)\right)$ only depends on the inter-subgraph edge probability p and is independent of the intra-subgraph architectures. This highlights that the inter-network connections in such graph models have an essential role.

In

Table 1. $EE\left(G_{m,n}\left(p\right)\right)$ for $m=2$, $n_{\left(1\right)}=n/2$ and $p=0.1$.

$\mathit{E}\mathit{E}\left({\mathit{G}}_{\mathit{m},\mathit{n}}\left(\mathit{p}\right)\right)$	Lower Bound	Upper Bound	Numerical Calculation
$n=2000$	$e^{100}(e^{O\left(\sqrt{2000}\right)}+o\left(1\right))$	$e^{1200}(e^{O\left(\sqrt{2000}\right)}+o\left(1\right))$	$e^{912}$
$n=4000$	$e^{200}(e^{O\left(\sqrt{4000}\right)}+o\left(1\right))$	$e^{2400}(e^{O\left(\sqrt{4000}\right)}+o\left(1\right))$	$e^{1843}$
$n=6000$	$e^{300}(e^{O\left(\sqrt{6000}\right)}+o\left(1\right))$	$e^{3600}(e^{O\left(\sqrt{6000}\right)}+o\left(1\right))$	$e^{2975}$
$n=8000$	$e^{400}(e^{O\left(\sqrt{8000}\right)}+o\left(1\right))$	$e^{4800}(e^{O\left(\sqrt{8000}\right)}+o\left(1\right))$	$e^{4036}$

we show the theoretical and experimental values for $EE\left(G_{m,n}\left(p\right)\right)$ with $m=2$ subgraphs, $n_{\left(1\right)}=n/2$ and $p=0.1$.

3. Laplacian Estrada Index of Random Interdependent Graphs

In this section, we consider Laplacian Estrada index $LEE\left(G_{m,n}\left(p\right)\right)$ for random interdependent graphs. The following result in regard to the Laplacian eigenvalues of random multipartite graphs has been essentially proved in [16]. We rephrase it as follows.

Lemma 2.

For $i=1,2,\cdots ,n-1$,

$$np-n_{\left(1\right)}p+o\left(n\right)\le {\lambda }_i\left(L\left(G\left(p\right)\right)\right)\le np+o\left(n\right)\hspace{1em}a.a.s.,$$

and

$$-n_{\left(1\right)}p+o\left(n\right)\le {\lambda }_n\left(L\left(G\left(p\right)\right)\right)\le o\left(n\right)\hspace{1em}a.a.s.$$

Theorem 2.

Let $G_{m,n}\left(p\right)$ be a random interdependent graph in $({\mathcal{S}}_n,{\mathcal{A}}_n,{\mathbb{P}}_n)$. We have

$$\begin{array}{cc}\hfill e^{np-n_{\left(1\right)}p}((n-1)e^{o\left(n\right)}+o\left(1\right))& \le LEE\left(G_{m,n}\left(p\right)\right)\hfill \\ \hfill & \le e^{np+n_{\left(1\right)}}(m{e}^{-n_{\left(1\right)}+o\left(n\right)}+(n-m)e^{o\left(n\right)})\hspace{1em}a.a.s.\hfill \end{array}$$

Proof. 

By the definition of random interdependent graphs, we observe that

$$\begin{array}{c}\hfill L\left(G_{m,n}\left(p\right)\right)=\overline{L}+L\left(G\left(p\right)\right),\end{array}$$

(6)

where $\overline{L}=\mathrm{diag}(L\left(G_1\right),L\left(G_2\right),\cdots ,L\left(G_m\right))$ is a block diagonal matrix with component blocks representing the Laplacian matrices of subgraphs.

To estimate $LEE\left(G_{m,n}\left(p\right)\right)$, we need to bound the eigenvalues of $L\left(G_{m,n}\left(p\right)\right)$. Note that the Laplacian eigenvalues vary monotonically with respect to edge addition or removal; see e.g., the interlacing theorem ([1] Theorem 7.1.5). Therefore, in view of Lemma 2, we have the following lower bounds

$$\begin{array}{c}\hfill np-n_{\left(1\right)}p+o\left(n\right)\le {\lambda }_i\left(L\left(G\left(p\right)\right)\right)\le {\lambda }_i\left(L\left(G_{m,n}\left(p\right)\right)\right)\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.\end{array}$$

(7)

for $i=1,2,\cdots ,n-1$, and

$$\begin{array}{c}\hfill -n_{\left(1\right)}p+o\left(n\right)\le {\lambda }_n\left(L\left(G\left(p\right)\right)\right)\le {\lambda }_n\left(L\left(G_{m,n}\left(p\right)\right)\right)\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.\end{array}$$

(8)

We obtained (7) and (8) by thinking all subgraphs $G_1,\cdots ,G_m$ as empty. On the other hand, if all these subgraphs are complete graphs, then ${\lambda }_1\left(\overline{L}\right),\cdots ,$${\lambda }_{n-m}\left(\overline{L}\right)\in \{(n_1-1)·n_1,(n_2-1)·n_2,\cdots ,(n_m-1)·n_m\}$, and ${\lambda }_{n-m+1}\left(\overline{L}\right)={\lambda }_{n-m+2}\left(\overline{L}\right)=\cdots ={\lambda }_n\left(\overline{L}\right)=0$. Here, we used the multiset representation $a_1·n_1$, meaning the multiplicity of element $n_1$ is $a_1$, etc. Therefore, by using Lemmas 1 and 2, and (6), we have

$$\begin{array}{c}\hfill {\lambda }_i\left(L\left(G_{m,n}\left(p\right)\right)\right)\le {\lambda }_i\left(\overline{L}\right)+{\lambda }_1\left(L\left(G\left(p\right)\right)\right)=np+o\left(n\right)\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.\end{array}$$

(9)

for $i=n-m+1,n-m+2,\cdots ,n$.

Recall that $n_{\left(1\right)}\ge n_{\left(2\right)}\ge \cdots \ge n_{\left(m\right)}$. We define $\langle n_{\left(1\right)}\rangle =\{1,2,\cdots ,n_{\left(1\right)}-1\}$, $\langle n_{\left(2\right)}\rangle =\{n_{\left(1\right)},n_{\left(1\right)}+1,\cdots ,n_{\left(1\right)}+n_{\left(2\right)}-2\}$, ⋯, $\langle n_{\left(m\right)}\rangle =\{n_{\left(1\right)}+n_{\left(2\right)}+\cdots +n_{(m-1)}-m+2,n_{\left(1\right)}+n_{\left(2\right)}+\cdots +n_{(m-1)}-m+3\cdots ,n_{\left(1\right)}+n_{\left(2\right)}+\cdots +n_{\left(m\right)}-m\}$. For $i\in \langle n_{\left(j\right)}\rangle$, $j=1,2,\cdots ,m$, it follows from Lemmas 1 and 2 that

$$\begin{array}{c}\hfill {\lambda }_i\left(L\left(G_{m,n}\left(p\right)\right)\right)\le {\lambda }_i\left(\overline{L}\right)+{\lambda }_1\left(L\left(G\left(p\right)\right)\right)\le n_{\left(j\right)}+np+o\left(n\right)\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.\end{array}$$

(10)

By the definition of Laplacian Estrada index and the estimates (7)–(10),

$$\begin{array}{cc}\hfill & (n-1)e^{np-n_{\left(1\right)}p+o\left(n\right)}+e^{-n_{\left(1\right)}p+o\left(n\right)}\hfill \\ \hfill \le & LEE\left(G_{m,n}\left(p\right)\right)\hfill \\ \hfill \le & m{e}^{np+o\left(n\right)}+(n_{\left(1\right)}-1)e^{n_{\left(1\right)}+np+o\left(n\right)}+\cdots +(n_{\left(m\right)}-1)e^{n_{\left(m\right)}+np+o\left(n\right)}\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.\hfill \end{array}$$

With a further look into the above lower and upper bounds, we have

$$\begin{array}{c}\hfill \frac{(n-1)e^{np-n_{\left(1\right)}p+o\left(n\right)}+e^{-n_{\left(1\right)}p+o\left(n\right)}}{e^{np-n_{\left(1\right)}p}}=(n-1)e^{o\left(n\right)}+o\left(1\right)\end{array}$$

and

$$\begin{array}{cc}\hfill & \frac{m{e}^{np+o\left(n\right)}+(n_{\left(1\right)}-1)e^{n_{\left(1\right)}+np+o\left(n\right)}+\cdots +(n_{\left(m\right)}-1)e^{n_{\left(m\right)}+np+o\left(n\right)}}{e^{np+n_{\left(1\right)}}}\hfill \\ \hfill \le & \frac{m{e}^{np+o\left(n\right)}+(n-m)e^{n_{\left(1\right)}+np+o\left(n\right)}}{e^{np+n_{\left(1\right)}}}\hfill \\ \hfill \le & m{e}^{-n_{\left(1\right)}+o\left(n\right)}+(n-m)e^{o\left(n\right)}.\hfill \end{array}$$

Hence,

$$\begin{array}{cc}\hfill e^{np-n_{\left(1\right)}p}((n-1)e^{o\left(n\right)}+o\left(1\right))& \le LEE\left(G_{m,n}\left(p\right)\right)\hfill \\ \hfill & \le e^{np+n_{\left(1\right)}}(m{e}^{-n_{\left(1\right)}+o\left(n\right)}+(n-m)e^{o\left(n\right)})\hspace{1em}\mathrm{a}.\mathrm{a}.\mathrm{s}.\hfill \end{array}$$

The proof is complete. ☐

Corollary 2.

Let $G_{m,n}\left(p\right)$ be a random interdependent graph in $({\mathcal{S}}_n,{\mathcal{A}}_n,{\mathbb{P}}_n)$. If $n_{\left(1\right)}=o\left(n\right)$, we have

$$\begin{array}{c}\hfill LEE\left(G_{m,n}\left(p\right)\right)=e^{np}(e^{o\left(n\right)}+o\left(1\right))\hspace{1em}a.a.s.\end{array}$$

Proof. 

If $n_{\left(1\right)}=o\left(n\right)$, by Theorem 2 we have

$$\begin{array}{cc}\hfill e^{np}((n-1+o\left(1\right))e^{o\left(n\right)}+o\left(1\right))& \le LEE\left(G_{m,n}\left(p\right)\right)\hfill \\ \hfill & \le e^{np+o\left(n\right)}(m{e}^{o\left(n\right)}+(n-m)e^{o\left(n\right)})\hfill \\ \hfill & \le e^{np+o\left(n\right)}(n{e}^{o\left(n\right)}+o\left(1\right)).\hfill \end{array}$$

Since $n=e^{o\left(n\right)}$, we obtain $LEE\left(G_{m,n}\left(p\right)\right)=e^{np}(e^{o\left(n\right)}+o\left(1\right))$ a.a.s. The proof is complete. ☐

Similar to Corollary 1, the Laplacian Estrada index $LEE\left(G_{m,n}\left(p\right)\right)$ only relies on the inter-subgraph edge probability p and is independent of the intra-subgraph topologies. If more information is available for the structure of subgraphs $G_1,\cdots ,G_m$, we may be able to derive sharper bounds for the (Laplacian) Estrada index of these interdependent graphs.

In

Table 2. $LEE\left(G_{m,n}\left(p\right)\right)$ for $m=2$, $n_{\left(1\right)}=n/2$ and $p=0.1$.

$\mathit{L}\mathit{E}\mathit{E}\left({\mathit{G}}_{\mathit{m},\mathit{n}}\left(\mathit{p}\right)\right)$	Lower Bound	Upper Bound	Numerical Calculation
$n=2000$	$e^{100}(1999e^{o\left(2000\right)}+o\left(1\right))$	$e^{1200}(1998e^{o\left(2000\right)}+o\left(1\right))$	$e^{935}$
$n=4000$	$e^{200}(3999e^{o\left(4000\right)}+o\left(1\right))$	$e^{2400}(3998e^{o\left(4000\right)}+o\left(1\right))$	$e^{1880}$
$n=6000$	$e^{300}(5999e^{o\left(6000\right)}+o\left(1\right))$	$e^{3600}(5998e^{o\left(6000\right)}+o\left(1\right))$	$e^{3127}$
$n=8000$	$e^{400}(7999e^{o\left(8000\right)}+o\left(1\right))$	$e^{4800}(7998e^{o\left(8000\right)}+o\left(1\right))$	$e^{4206}$

we show the theoretical and experimental values for $LEE\left(G_{m,n}\left(p\right)\right)$ with $m=2$ subgraphs, $n_{\left(1\right)}=n/2$ and $p=0.1$.

4. Conclusions

We have studied the Estrada index and Laplacian Estrada index of a class of random interdependent graphs $G_{m,n}\left(p\right)$. Some lower bounds and upper bounds are determined in the asymptotically almost sure limit as the number of vertices n tends to infinity. Our results have key connotations for such interdependent graphs: the interdependent edges between different subgraphs seem to play a key role in their spectral properties.

Although the random interdependent graph takes Erdös-Rényi random graph, random multipartite graph, and even stochastic block model as special cases, it does not cover the edge-independent random graph in general [25]. It seems interesting to further extend the results in this paper to random edge-independent interdependent graphs. It would also be desirable to explore more general spectral properties for such interdependent graphs. We leave these as future work.