Loading [MathJax]/jax/output/SVG/fonts/TeX/fontdata.js
Review

A review on epileptic foci localization using resting-state functional magnetic resonance imaging

  • Epilepsy is a brain syndrome caused by synchronous abnormal discharge of brain neurons. As an effective treatment for epilepsy, successful surgical resection requires accurate localization of epileptic foci to avoid iatrogenic disability. Previous studies have demonstrated the potential of restingstate functional magnetic resonance imaging (rs-fMRI) technique to localize epileptic foci though clinical applications of rs-fMRI are still at an early stage of development. fMRI data analysis approaches seek pre-defined regressors modeling contributions to the voxel time series, including the BOLD response following neuronal activation. In present study, localization strategies of epileptic foci in rs-fMRI technology were classified and summarized. To begin with, data-driven approaches attempting to determine the intrinsic structure of the data were discussed in detail. Then, as novel fMRI data analysis methods, deconvolution algorithms such as total activation (TA) and blind deconvolution were discussed, which were applied to explore the underlying activity-inducing signal of the BOLD signal. Lastly, effective connectivity approaches such as autocorrelation function method and Pearson correlation coefficient have also been proposed to identify the brain regions driving the generation of seizures within the epileptic network. In the future, fMRI technology can be used as a supplement of intraoperative subdural electrode method or combined with traditional epileptic focus localization technologies, which is one of the most attractive aspect in clinic. It may also play an important role in providing diagnostic information for epilepsy patients.

    Citation: Yue Shi, Xin Zhang, Chunlan Yang, Jiechuan Ren, Zhimei Li, Qun Wang. A review on epileptic foci localization using resting-state functional magnetic resonance imaging[J]. Mathematical Biosciences and Engineering, 2020, 17(3): 2496-2515. doi: 10.3934/mbe.2020137

    Related Papers:

    [1] Zhen Lin . On the sum of the largest Aα-eigenvalues of graphs. AIMS Mathematics, 2022, 7(8): 15064-15074. doi: 10.3934/math.2022825
    [2] Bilal A. Rather, Fawad Ali, Nasim Ullah, Al-Sharef Mohammad, Anwarud Din, Sehra . Aα matrix of commuting graphs of non-abelian groups. AIMS Mathematics, 2022, 7(8): 15436-15452. doi: 10.3934/math.2022845
    [3] Alfredo Sotelo-Pejerrey . Traces of certain integral operators related to the Riemann hypothesis. AIMS Mathematics, 2023, 8(10): 24971-24983. doi: 10.3934/math.20231274
    [4] Shuo Wang, Juan Liu, Xindong Zhang . Properties of solutions for fractional-order linear system with differential equations. AIMS Mathematics, 2022, 7(8): 15704-15713. doi: 10.3934/math.2022860
    [5] Efruz Özlem Mersin . Sturm's Theorem for Min matrices. AIMS Mathematics, 2023, 8(7): 17229-17245. doi: 10.3934/math.2023880
    [6] Samer Al-Ghour, Hanan Al-Saadi . Soft weakly connected sets and soft weakly connected components. AIMS Mathematics, 2024, 9(1): 1562-1575. doi: 10.3934/math.2024077
    [7] Meraj Ali Khan, Ali H. Alkhaldi, Mohd. Aquib . Estimation of eigenvalues for the α-Laplace operator on pseudo-slant submanifolds of generalized Sasakian space forms. AIMS Mathematics, 2022, 7(9): 16054-16066. doi: 10.3934/math.2022879
    [8] Jun He, Yan-Min Liu, Jun-Kang Tian, Ze-Rong Ren . A note on the inclusion sets for singular values. AIMS Mathematics, 2017, 2(2): 315-321. doi: 10.3934/Math.2017.2.315
    [9] Xiaofei Cao, Yuyue Huang, Xue Hua, Tingyu Zhao, Sanzhang Xu . Matrix inverses along the core parts of three matrix decompositions. AIMS Mathematics, 2023, 8(12): 30194-30208. doi: 10.3934/math.20231543
    [10] Xiaomei Zhang, Xiang Chen . Uniqueness of difference polynomials. AIMS Mathematics, 2021, 6(10): 10485-10494. doi: 10.3934/math.2021608
  • Epilepsy is a brain syndrome caused by synchronous abnormal discharge of brain neurons. As an effective treatment for epilepsy, successful surgical resection requires accurate localization of epileptic foci to avoid iatrogenic disability. Previous studies have demonstrated the potential of restingstate functional magnetic resonance imaging (rs-fMRI) technique to localize epileptic foci though clinical applications of rs-fMRI are still at an early stage of development. fMRI data analysis approaches seek pre-defined regressors modeling contributions to the voxel time series, including the BOLD response following neuronal activation. In present study, localization strategies of epileptic foci in rs-fMRI technology were classified and summarized. To begin with, data-driven approaches attempting to determine the intrinsic structure of the data were discussed in detail. Then, as novel fMRI data analysis methods, deconvolution algorithms such as total activation (TA) and blind deconvolution were discussed, which were applied to explore the underlying activity-inducing signal of the BOLD signal. Lastly, effective connectivity approaches such as autocorrelation function method and Pearson correlation coefficient have also been proposed to identify the brain regions driving the generation of seizures within the epileptic network. In the future, fMRI technology can be used as a supplement of intraoperative subdural electrode method or combined with traditional epileptic focus localization technologies, which is one of the most attractive aspect in clinic. It may also play an important role in providing diagnostic information for epilepsy patients.


    Let

    G=(V(G),E(G))

    be a simple undirected graph on n vertices such that V(G) and E(G) denote the vertex set and edge set of G, respectively. Let V(G)={v1,v2,,vn}, for viV(G), and dG(vi)=d(vi)=di denotes the degree of vertex vi and the set of vertices adjacent to vi, denoted by N(vi), refers to the neighborhood of vi. The maximum and the minimum degree of G are denoted by Δ=Δ(G) and δ=δ(G), respectively. λ(M), λk(M) and λmin(M) denote the largest eigenvalue, the k-th largest eigenvalue, and the smallest eigenvalue of a matrix M, respectively. The set of all eigenvalues for a matrix M together with their multiplicities is called the M-spectrum. A real symmetric matrix M is called positive semidefinite if λmin(M)0. Denote by Kn, Pn, Cn and Ks,ns the complete graph, path, cycle, and complete bipartite graph with n vertices, respectively. Let A(G) and D(G) be the adjacency matrix and the diagonal matrix of the degrees of G, respectively. The signless Laplacian Q(G) of G is defined as

    Q(G)=D(G)+A(G).

    The Laplacian L(G) of G is defined as

    L(G)=D(G)A(G).

    We denote the eigenvalues of Q(G) and L(G) by ρk(G) and λk(G), respectively, where 1kn (we drop G when it is clear from the context). In particular, the largest eigenvalue of Q(G) is called the signless Laplacian spectral radius of G, denoted by ρ(G), and the largest eigenvalue of L(G) is called the Laplacian spectral radius of G, denoted by λ(G).

    In [1], Nikiforov introduced the Aα-matrix of G, which is the convex linear combination of A(G) and D(G) defined by

    Aα(G)=αD(G)+(1α)A(G),0α1.

    Clearly,

    A0(G)=A(G),  A1(G)=D(G)and2A12(G)=Q(G).

    Thus, the matrix Aα(G) is a generalization of the adjacency matrix and the signless Laplacian matrix. The largest eigenvalue of Aα(G) is called the Aα-spectral radius of G, denoted by ρα(G). Nikiforov [1] had studied the matrix Aα(G), and he has investigated many properties on Aα(G), including bounds on the k-th largest (especially, the largest, the second largest, and the smallest) eigenvalue of Aα(G), the positive semidefiniteness of Aα(G), etc. From then on, the study of Aα-spectra and the Aα-spectral radius for the graph has attracted the attention of many researchers. In [2], Nikiforov et al. gave an upper bound of the spectral radius of Aα(TΔ), where TΔ is the tree of maximal degree Δ. Also, he obtained several bounds on the spectral radius of Aα of general graphs. The graphs with

    λk(Aα(G))=αn1,2kn

    had been characterized by Lin et al. in [3]. In [4], Guo and Zhou gave upper bounds for the Aα-spectral radius for some graphs under some conditions. Guo and Zhang [5]; obtained a sharp upper bound on the Aα-spectral radius for α[12,1), and proved that for two connected graphs G1 and G2, ρα(G1)>ρα(G2) under certain conditions. For further related studies, one may see [6,7,8].

    The research on L(G) has shown that it is a remarkable matrix in many respects. Since; L(G) is just the difference between D(G) and A(G), then to understand to what extent each of the summands A(G) and D(G) determines the properties of L(G), and motivated by the above works, we introduce the Aα-matrix defined by

    Aα(G)=αD(G)+(α1)A(G),0α1. (1.1)

    Through our study, several facts indicate that the study of the Aα(G) family is of a unique importance. First of all, obviously

    A(G)=A0(G),D(G)=A1(G)andL(G)=2A12(G).

    Since A12(G) is essentially equivalent to L(G), then the matrix Aα(G) can be regarded as a generalization of the Laplacian matrix L(G). Since

    A1(G)=A1(G)=D(G);

    unless otherwise specified, we only consider the case of 0α<1 in this paper.

    In this paper, we study some properties on Aα(G) and obtain some bounds of the largest eigenvalue of Aα(G). For the relation between the Q-spectrum and L-spectrum of a connected graph G, it is well known that G is bipartite if and only if Q(G) and L(G) has the same spectrum ([9, Proposition 1.3.10]), and we extend this relation on the Aα- and Aα-spectrum. In [10], Akbari et al. proved that, for a graph G of order n and any real number β, if 0<β1 or 2β3, then

    nk=1ρβknk=1λβk;

    and if 1β2, then

    nk=1ρβknk=1λβk.

    In this work we present the relation between the sum of powers of Aα and Aα eigenvalues, which is a generalization of the relation between the sum of powers of Q(G) and L(G) eigenvalues in [10].

    The rest of this paper is organized as follows. In Section 2, we give some preliminaries and lemmas used later. In Section 3, we derive new basic properties of Aα(G). In Section 4, we present the Aα-spectra of the complete graphs and the complete split graphs. Section 5 gives some bounds of the largest eigenvalue of Aα(G). In Section 6, we deduce the relation between the Aα- and Aα-spectral radius. Finally, in Section 7, we prove the relation between the sum of powers of Aα and Aα eigenvalues.

    Although Weyl's inequalities have been known for almost a century (see, e.g., [11]), their equality case was first established by So in [12], and for convenience we state Weyl and So's complete theorem as follows:

    Theorem 2.1. Let A and B be Hermitian matrices of order n, and let 1in and 1jn. Then,

    λi(A)+λj(B)λi+jn(A+B),ifi+jn+1, (2.1)
    λi(A)+λj(B)λi+j1(A+B),ifi+jn+1. (2.2)

    In either of these inequalities, the equality holds if and only if there exists a nonzero n-vector that is an eigenvector to each of the three eigenvalues involved.

    A simplified version of (2.1) and (2.2) gives

    λk(A)+λmin(B)λk(A+B)λk(A)+λ(B). (2.3)

    We shall need the following lemmas for our new results.

    Lemma 2.2. [3] Let G be a graph of order n. If eE(G) and α12, then,

    ρi(Aα(G))ρi(Aα(Ge)),for1in.

    Lemma 2.3. [13] Let B be a real symmetric nonnegative irreducible matrix and λ be the largest eigenvalue of B. ZRn. If ZtBZ=λ and Z=1, then BZ=λZ.

    Lemma 2.4. [5] Let G be a connected graph with n4 vertices and m edges. If α[12,1), then,

    ρα(G)max{αΔ(G),(1α)(mn12)}+2α.

    Equality holds if and only if α=n1n+1 and G=Kn.

    For a graph G of order n, suppose that λ is an eigenvalue of Aα(G) and x is an eigenvector of Aα(G) with respect to λ. We use x(v) to denote the entry of x corresponding to the vertex vV(G). It is clear that the system of eigenequations for the matrix Aα(G) is

    λx(v)=αdG(v)x(v)+(α1)uN(v)x(u). (3.1)

    If G is a graph of order n with Aα(G)=Aα, and x is a real vector, the quadratic form Aαx,x can be represented in several equivalent ways:

    Aαx,x=uvE(G)(αx(u)2+2(α1)x(u)x(v)+αx(v)2), (3.2)
    Aαx,x=(2α1)uV(G)x(u)2dG(u)+(1α)uvE(G)(x(u)x(v))2, (3.3)
    Aαx,x=αuV(G)x(u)2dG(u)+2(α1)uvE(G)x(u)x(v). (3.4)

    Each of these representations can be useful in proofs.

    Now, we give some of the spectral properties of the Aα-matrix. Let us call the largest eigenvalue of Aα(G) the Aα-spectral radius of G, and denote it as λα(G). Let us also denote the smallest eigenvalue of Aα(G) as μα(G). Since Aα(G) is a real symmetric matrix, and by using Rayleigh's principle, the following result holds:

    Proposition 3.1. If α[0,1) and G is a graph of order n, then

    λα(G)=maxx2=1Aα(G)x,xandμα(G)=minx2=1Aα(G)x,x. (3.5)

    Moreover, if x is a unit n-vector, then,

    λα(G)=Aα(G)x,x

    if and only if x is an eigenvector to λα(G), and

    μα(G)=Aα(G)x,x

    if and only if x is an eigenvector to μα(G).

    By using these relations, the following result is evident:

    Proposition 3.2. If α[0,1) and G is a graph, then

    λα(G)=max{λα(H):HisacomponentofG},
    μα(G)=min{μα(H):HisacomponentofG}.

    It is clear that if G is a d-regular graph of order n, then

    Aα(G)=αdIn+(α1)A(G),

    and so there is a linear correspondence between the spectra of Aα(G) and of A(G),

    λk(Aα(G))=αd+(α1)λnk+1(A(G)), 1kn. (3.6)

    In particular, if k=n, then

    μα(G)=(2α1)d

    for any α[0,1]. As a consequence of Weyl's inequality (2.3), the following result is immediate:

    Proposition 3.3. If α[0,1] and G is a graph with

    A(G)=AandAα(G)=Aα,

    then,

    αδ+(α1)λnk+1(A)λk(Aα)αΔ+(α1)λnk+1(A).

    An important property of the Laplacian L(G) is that it is positive semidefinite. This is certainly not true for Aα(G) if α<12 and G is a regular graph, but if α12, then Aα(G) is like L(G). We give this result in the following:

    Proposition 3.4. If α12, and G is a graph, then Aα(G) is positive semidefinite, and if α>12 and G has no isolated vertices, then Aα(G) is positive definite. Moreover, if α<12 and G is a regular graph, then Aα(G) is not positive semidefinite.

    Proof. Let x be a nonzero vector. If α12, then for any edge uvE(G), we see that

    Aα(G)x,x(2α1)(x(u)2+x(v)2)+(1α)(x(u)x(v))20. (3.7)

    Hence, Aα(G) is positive semidefinite. Now, assume that G has no isolated vertices. Choose a vertex u with x(u)0 and let uvE(G). Then, (3.7) becomes a strict inequality, and so Aα(G) is positive definite. Finally, suppose that G is a dregular graph, α<12, and let A be its adjacency matrix, then λ(A)=d. Thus

    μα(G)=(2α1)d<0.

    By the above proposition we get the following lemma which gives a relation between the Aα eigenvalues of G and the Aα eigenvalues of spanning subgraphs of G.

    Lemma 3.5. Let G be a graph of order n and let α[12,1). If G=Ge, where eE(G), then,

    λi(G)λi(G)for1in.

    Proof. Let e=uv such that u,vV(G). It is easy to see that

    Aα(G)=Aα(G)+M,

    where M is the matrix of order n indexed by the vertices of G having (u,v)th and (v,u)th entries both equal to α1, and the (u,u)th and (v,v)th entries both equal to α and all other entries equal to zero, hence M is an Aα-matrix of a graph containing only one edge. Since α[12,1), then Aα(G), Aα(G) and M are positive semidefinite and Weyl's inequalities (2.3) imply that λi(G)λi(G).

    Equality (3.6) and the fact the eigenvalues of A(Kn) are {n1,1,,1} give the spectrum of Aα(Kn) as follows:

    Proposition 4.1. The eigenvalues of Aα(Kn) are

    μα(Kn)=(2α1)(n1)andλk(Aα(Kn))=α(n2)+1for1kn1.

    If SV(G), then we use G[S] to denote the subgraph of G induced by S. Recall that G[S] is an independent set if no two vertices of S are adjacent and G[S] is a clique if it is a complete subgraph of G. The graph Kr(nr)K1 is called a complete split graph, denoted by CSr,nr. The work in the following proposition is motivated by the proof of [3, Proposition 2.4].

    Proposition 4.2. Let G be a graph with Aα(G), and α[0,1). Let SV(G) and |S|=k. Suppose that dG(u)=d for each vertex uS, and N(v){w}=N(w){v} for any two vertices v,wS. Then, we have the following statements:

    (i) If G[S] is a clique, then α(d1)+1 is an eigenvalue of Aα(G) with multiplicity at least k1.

    (ii) If G[S] is an independent set, then αd is an eigenvalue of Aα(G) with multiplicity at least k1.

    Proof. Let S={v1,v2,,vk}. Clearly, d1==dk=d. Let z1,z2,,zk1 be vectors such that

    {zi(v1)=1,zi(vi+1)=1,zi(v)=0,if vV(G){v1,vi+1},

    for i=1,,k1. Assume that G[S] is a clique. It is easy to obtain that

    Aα(G)zi=(α(d11)+1,0,,0,α(1di)1,0,,0)=(α(d1)+1)zi,

    for i=1,,k1. Hence, α(d1)+1 is an eigenvalue of Aα(G) and z1,z2,,zk1 are eigenvectors of Aα(G) corresponding to α(d1)+1. In addition, since z1,z2,,zk1 are linearly independent, the multiplicity of α(d1)+1 is at least k1. Now, supposing that G[S] is an independent set, we have

    Aα(G)zi=(αd1,0,,0,αdi,0,,0)=αdzi,

    for i=1,,k1. Since z1,z2,,zk1 are linearly independent, it follows that αd is an eigenvalue of Aα(G) with multiplicity at least k1. Thus, the proof is completed.

    Consider an n×n real symmetric matrix

    S=(S11S12S1tS21S22S2tSt1St2Stt),

    whose rows and columns are partitioned according to a partitioning P1, P2, , Pt of {1,2,,n}. The quotient matrix B of the matrix S is the t×t matrix whose entries are the average row sums of the blocks Sij of S. The partition is equitable if each block Sij of S has constant row sum.

    Lemma 4.3. [14] Let B be an equitable quotient matrix of a symmetric real matrix S. If λ is an eigenvalue of B, then λ is also an eigenvalue of S.

    Now, we can determine all Aα-eigenvalues of CSr,nr as follows:

    Proposition 4.4. The Aα-spectrum of CSr,nr contains α(n2)+1 with multiplicity r1, αr with multiplicity nr1, and the remaining two Aα-eigenvalues are

    α(n+2(r1))+1±(α(n+2(r1))+1)2+4r(12α)(nr+α(r1))2. (4.1)

    Proof. We can write Aα(CSr,nr) as

    Aα(CSr,nr)=((1α)Jr,r+(αn1)Ir(1α)Jr,nr(1α)Jnr,rαrInr,nr).

    Then, the quotient matrix of Aα(CSr,nr) is equitable and it can be written in the form

    B(Aα(CSr,nr))=((nr)α+r1(1α)(nr)(1α)rαr).

    Thus, by Lemma 4.3, the eigenvalues of B(Aα(CSr,nr)) are eigenvalues of Aα(CSr,nr), and according to Proposition 4.2, we get the result.

    In this section we give a few general bounds on λα(G).

    Proposition 5.1. Let G be a graph, with

    Δ(G)=Δ,  A(G)=AandD(G)=D.

    The following inequalities hold for λα(G):

    λα(G)αΔ+(α1)λ(A), (5.1)
    λα(G)(2α1)λ(A). (5.2)

    Proof. Inequality (5.1) follows by Weyl's inequalities (2.3) because

    αΔ+(α1)λ(A)=αλ(D)+(α1)λ(A)=λ(αD)+λmin((α1)A)λ(αD+(α1)A)=λα(G).

    To prove the inequality (5.2), let H be a component of G such that λ(A)=λ(A(H)). Let x be a positive unit vector to λ(A(H)). For every edge uvE(H), the AM-GM inequality implies that

    2x(u)x(v)=2αx(u)x(v)+2(1α)x(u)x(v)α(x(u)+x(v))22+2(1α)x(u)x(v)=12(αx(u)2+2(α1)x(u)x(v)+αx(v)2)+(32α)x(u)x(v).

    Summing this inequality over all edges uvE(H), and using (3.2), we get

    λ(A)=λ(A(H))=A(H)x,x12Aα(H)x,x+12(32α)λ(A);

    and then we get

    2λ(A)(32α)λ(A)Aα(H)x,xλα(G);

    hence

    (2α1)λ(A)λα(G).

    Having inequality (5.2) in hand, if α12, then every lower bound of λ(A) gives a lower bound on λα(G), but if α<12, then every upper bound of λ(A) gives a lower bound on λα(G). We mention just two such bounds.

    Corollary 5.2. Let G be a graph such that α12. If G is of order n and has m edges, then

    λα(G)(2α1)2mn.

    Corollary 5.3. Let G be a connected graph such that α<12. If G is of order n and has m edges, then

    λα(G)(2α1)2mn+1.

    Proposition 5.4. Let G be a graph of order n. If α(12,1], then,

    λα(G)α(n2)+1.

    Moreover,

    λα(G)=α(n2)+1

    with multiplicity k1 if G has k vertices of degree n1.

    Proof. Applying Lemma 3.5 leads to

    λα(G)λα(Kn)=α(n2)+1.

    If G has k vertices of degree n1, then it follows from Proposition 4.2 that α(n2)+1 is an eigenvalue of Aα(G) with multiplicity at least k1, and since

    λα(G)α(n2)+1,

    we get

    λα(G)=α(n2)+1

    with multiplicity k1.

    Let G be a connected graph. Merris [15] pointed out λ(L(G))ρ(Q(G)), and the equality holds if G is a bipartite graph. In the next result, we generalize this result to the Aα- and Aα-spectral radius of a connected graph.

    Theorem 6.1. Let G be a graph of order n, α(0,1), λα(G)=λα and ρα(G)=ρα. We have λαρα. Moreover, if G is connected, then the equality holds if and only if G is bipartite.

    Proof. Let

    V(G)={v1,v2,,vn},  x=(x1,x2,,xn)tRn

    be an arbitrary vector such that x=1. Let

    y=(y1,y2,,yn)tRn

    be a unit eigenvector of Aα(G) belonging to λα and

    y=(y1,y2,,yn)tRn

    be a unit eigenvector of Aα(G) belonging to ρα. Let

    |y|=(|y1|,|y2|,,|yn|)t.

    First, we prove that λαρα.

    λα=maxxtAαx=maxxt(αD+(α1)A)x=max(αni=1x2idi+2(α1)vivjE(G)xixj)=ytAαy=αni=1y2idi+2(α1)vivjE(G)yiyj

    and

    ρα=maxxtAαx=maxxt(αD+(1α)A)x=max(αni=1x2idi+2(1α)vivjE(G)xixj)=ytAαy=αni=1yi2di+2(1α)vivjE(G)yiyj.

    Thus,

    λα=αni=1y2idi+2(α1)vivjE(G)yiyjαni=1y2idi+2(1α)vivjE(G)|yiyj|=|yt|Aα|y|maxxtAαx=ρα. (6.1)

    Now, if G is bipartite, then the matrix Aα and the matrix Aα are similar by a diagonal matrix D with diagonal entries ±1, that is, Aα=DAαD1. Therefore, Aα and Aα have the same spectrum, and thus we get λα=ρα. Finally, when G is connected and λα=ρα, all inequalities (6.1) must be equalities. By Lemma 2.3 and the equality

    |yt|Aα|y|=ρα;

    we know that |y| is an eigenvector of Aα belonging to ρα. So, |y|=±y. Using the Perron-Frobenius' theorem for Aα(G), we have y>0, |y|=y, and |yi|>0, i=1,2,,n.

    Since

    αni=1y2idi+2(α1)vivjE(G)yiyj=αni=1y2idi+2(1α)vivjE(G)|yiyj|;

    we get

    vivjE(G)yiyj=vivjE(G)|yiyj|,

    hence, |yiyj|=yiyj when vivjE(G). Therefore, yiyj<0 if vivjE(G).

    Let

    U={vi:yi>0}andW={vj:yj<0}.

    For each edge e=vivj, we have yiyj<0, one of the vertices of edge e is in U, and the other is in W. So, G is a bipartite graph.

    Remark 6.2. If G is not bipartite and is not connected, and λα is the largest eigenvalue of Aα for a bipartite connected component of G, then the equality in Theorem 6.1 still holds.

    Example 6.3. Take

    G=C3C4.

    Then, G is not bipartite and not connected. Now, A14(G) has a spectrum 2[1],1.25[2],0.5[2],1[2] (where λ[i] means the eigenvalue λ is repeated i times in the spectrum). On the other hand, A14(G) has a spectrum 2[2],0.5[2],0.25[2],1[1]. Then, we have

    λ14(G)=ρ14(G)=2.

    Note that 2 is the largest eigenvalue of A14(C4) as well.

    Now, we introduce the relationship between the Aα- and Aα-spectra of bipartite graphs, which is a generalization [9, Proposition 1.3.10], and it follows from the proof of the above theorem.

    Corollary 6.4. Let G be a connected graph. Then, G is bipartite if and only if the Aα-spectrum and Aα-spectrum are equal.

    Remark 6.5. In fact, if G is bipartite and is not connected, the Aα-spectrum still equals the Aα-spectrum.

    Example 6.6. Take G=P3C4. Then G is bipartite and not connected. We have that A34(G) has a spectrum 2[1],1.64039[1],1.5[2],1[1],0.75[1],0.609612[1]. In contrast, A34(G) has the same spectrum as A34(G), although G is not connected.

    According to Corollary 6.4 and [1, Propositions 38 and 39], we get the following two results:

    Corollary 6.7. The Aα-spectrum and the Aα-spectrum of the complete bipartite graph Ka,b are equal, that is, if ab1 and α(0,1), the eigenvalues of Aα(Ka,b) are

    λα(Ka,b)=12(α(a+b)+α2(a+b)2+4ab(12α)),
    μα(Ka,b)=12(α(a+b)α2(a+b)2+4ab(12α)),
    λk(Aα(Ka,b))=αafor1<kb,
    λk(Aα(Ka,b))=αbforb<k<a+b.

    Corollary 6.8. The Aα-spectrum and the Aα-spectrum of the star K1,n1 are equal, that is, the eigenvalues of Aα(K1,n1) are

    λα(K1,n1)=12(αn+α2n2+4(n1)(12α)),
    μα(K1,n1)=12(αnα2n2+4(n1)(12α)),
    λk(Aα(K1,n1))=αfor1<k<n.

    Indeed, many practical results of Theorem 6.1 can be deduced, and here are some of them.

    Proposition 6.9. Let G be a connected graph with n4 vertices and m edges. If α[12,1), then

    λα(G)<max{αΔ(G),(1α)(mn12)}+2α. (6.2)

    Proof. Let

    λα(G)=λαandρα(G)=ρα.

    Then, Theorem 6.1 and Lemma 2.4 lead to

    λαmax{αΔ(G),(1α)(mn12)}+2α. (6.3)

    Suppose that the equality in (6.3) holds, thus ρα=λα, and so G is bipartite and

    ρα=max{αΔ(G),(1α)(mn12)}+2α.

    Therefore, by Lemma 2.4, G=Kn, and thus G is bipartite if and only if n=2, but n4, and hence the inequality is strict.

    Theorem 6.1 and [2, Theorem 2] lead directly to the next result.

    Proposition 6.10. If T is a tree of order n and α[0,1], then

    λα(T)nα+n2α2+4(n1)(12α)2.

    The equality holds if and only if T is the star K1,n1.

    By Theorem 6.1 and [1, Proposition 20], we get the next result.

    Proposition 6.11. If G is a graph with no isolated vertices, then

    λα(G)maxuV(G){αd(u)+1αd(u)uvE(G)d(v)}.

    If α(12,1) and G is connected, the equality holds if and only if G is a regular bipartite graph.

    In this part we give two explicit expressions for the sums and the sum of squares of the eigenvalues of Aα and Aα, considering [1, Propositions 34 and 35].

    Proposition 6.12. If G is a graph of order n and has m edges, then

    ni=1λi(Aα(G))=ni=1λi(Aα(G))=trAα(G)=trAα(G);

    where

    ni=1λi(Aα(G))=αuV(G)d(u)=2αm.

    A similar formula for the sum of the squares of the Aα and Aα-eigenvalues is given as follows:

    Proposition 6.13. If G is a graph of order n and has m edges, then

    ni=1λ2i(Aα(G))=ni=1λ2i(Aα(G))=trA2α(G)=trA2α(G);

    where

    ni=1λ2i(Aα(G))=α2uV(G)d2(u)+2(1α)2m.

    Proof. Let

    Aα=Aα(G),  A=A(G)andD=D(G).

    Calculating the square A2α and taking its trace, we find that

    trA2α=tr(α2D2+(1α)2A2+α(α1)DA+α(α1)AD)=α2trD2+(1α)2trA2α(1α)trDA+α(α1)trAD=α2uV(G)d2(u)+2(1α)2m.

    Pirzada et al. [16] introduced the sum of the βth powers of the Aα eigenvalues of G as

    Sαβ(G)=ni=1ρβi.

    Now, we have the notation

    Sαβ(G)=ni=1λβi

    for the sum of the βth powers of the Aα eigenvalues of G. The following theorem is a generalization [10, Theorem 2].

    Theorem 7.1. Let G be a graph of order n>1, α[12,1), and let β be a real number.

    (i) If 0<β1 or 2β3, then

    Sαβ(G)Sαβ(G).

    (ii) If 1β2, then

    Sαβ(G)Sαβ(G).

    For β(0,1)(2,3), the equality holds in (i) if and only if G is a bipartite graph. Moreover, for β(1,2), the equality holds in (ii) if and only if G is a bipartite graph.

    Proof. We recall that, for any real number r, the binomial series

    k=0(rk)xk

    converges to (1+x)r if |x|<1. This also remains true for x=1 if r>0 (see, e.g., [17, p. 419]). By Lemma 3.5, we find that,

    λα(G)λα(Kn)=α(n2)+1andρα(G)ρα(Kn)=n1.

    Since Aα is positive semidefinite when α[12,1), then if ρi>0, we have

    |ρin1||n1n1|=|1n|<1

    and if ρi=0, we get

    ρin1=1.

    Therefore,

    Sαβ(G)nβ=(ρ1n)β++(ρnn)β=k=0(βk)(ρ1n1)k++k=0(βk)(ρnn1)k=k=0(βk)tr(1n(αD+(1α)A)I)k.

    Also, since Aα is positive semidefinite when α[12,1), then if λi>0,

    |λin1||α(n2)+1n1|=|(α1)+12αn|<1

    and if λi=0, we have

    λin1=1.

    Thus, in a similar manner as above, we obtain that

    Sαβ(G)nβ=k=0(βk)tr(1n(αD+(α1)A)I)k.

    We claim that

    ifkis even, tr(αD+(1α)AnI)ktr(αD+(α1)AnI)k;
    ifkis odd, tr(αD+(1α)AnI)ktr(αD+(α1)AnI)k.

    When ((αDnI)+(1α)A)k and ((αDnI)+(α1)A)k are expanded in terms of the powers of αDnI and (1α)A, respectively, the terms appearing in both expansions, regardless of their signs, are the same. To prove this claim, we identify the sign of each term in both expansions. Consider the terms in the expansion of ((αDnI)+(1α)A)k, where there are exactly j factors equal to αDnI, for some j=0,1,,k. The sign of the trace for each of these terms is (1)j or 0 because all entries of αDnI and (1α)A are non-positive and non-negative, respectively. On the other hand, in each term in the expansion of ((αDnI)+(α1)A)k all factors are matrices with non-positive entries, hence the sign of the trace of each term is (1)k or 0. Therefore, the claim has been proven.

    Now, note that if 0<β<1 or 2<β<3, then the sign of (βk) is (1)k1, except that (β2)>0, for 2<β<3. According to this, for 0<β<1 and every k,

    (βk)tr(αD+(1α)AnI)k(βk)tr(αD+(α1)AnI)k.

    This inequality remains true for 2β3 as

    tr(αD+(1α)AnI)2=tr(αD+(α1)AnI)2,

    since trA2α=trA2α. Thus, Part (i) is proved. For 1<β<2, the sign of (βk) is (1)k1, except that (β1)>0. Since trAα=trAα, we have

    tr(αD+(1α)AnI)=tr(αD+(α1)AnI),

    and so part (ii) is similarly proved.

    Now, we examine the equality case. Since Aα and Aα are similar if G is bipartite, it follows that the equality holds in both (i) and (ii). Since for any positive integer i, trAi equals the total number of closed walks of length i in G, then if G is not bipartite, there exists an odd integer r such that trAr>0 (see [18, Lemma 2.5]). Hence,

    tr(αD+(1α)AnI)r>tr(αD+(α1)AnI)r;

    and so the inequalities in both (i) and (ii) are strict.

    We know that the Aα-spectra and Aα-spectra of Kn are {(n1)[1],(αn1)[n1]} and {((2α1)(n1))[1],(α(n2)+1)[n1]}, respectively. Therefore,

    Sαβ(Kn)Sαβ(Kn)=(n1)β+(n1)(αn1)β((2α1)(n1))β(n1)(α(n2)+1)β.

    By Theorem 7.1 and based on our numerical experiments, we propose the following conjecture:

    Conjecture 7.2. For every α(12,1) and each integer n3, we have

    Sαβ(Kn)Sαβ(Kn)0

    for any β[0,1][2,) and

    Sαβ(Kn)Sαβ(Kn)0

    for any β(,0)(1,2).

    In this paper, we have introduced the Aα-matrix of a graph G, which is a generalization of the Laplacian matrix L(G), and we have studied the basic properties of Aα; and derived some bounds for its spectral radius. Furthermore, we have determined the Aα-spectra for the complete graph and the complete split graph. Building upon previous results, we have extended findings related to the spectral radius of L(G) and Q(G) matrices to the Aα- and Aα-spectral radius. Specifically, in Theorem 6.1, we have generalized Merris' [15] observation that λ(L(G))ρ(Q(G)) with equality holding for bipartite graphs. Additionally, we have extended the known relation that G is bipartite if and only if Q(G) and L(G) share the same spectrum, as demonstrated in Corollary 6.4. Finally, in Theorem 7.1; we have generalized a relation established by S. Akbari et al. in [10]; which relates the sum of powers of the eigenvalues of Q(G) and L(G). In conclusion, these findings have implications for various applications involving graph analysis.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The authors are grateful to the referees for their valuable suggestions and comments.

    The authors declare no conflicts of interest.



    [1] N. Jetté, J. W. Sander, M. R. Keezer, Surgical treatment for epilepsy: The potential gap between evidence and practice, Lancet Neurol., 15 (2016), 982-994.
    [2] J. S. Ebersole, S. M. Ebersole, Combining MEG and EEG source modeling in epilepsy evaluations, J. Clin. Neurophysiol., 27 (2010), 360-371.
    [3] S. Baxendale, The role of functional MRI in the presurgical investigation of temporal lobe epilepsy patients: a clinical perspective and review, J. Clin. Exp. Neuropsychol., 24 (2002), 664-676.
    [4] K. Lee, S. Tak, J. C. Ye, A data-driven sparse GLM for fMRI analysis using sparse dictionary learning with MDL criterion, IEEE Trans. Med. Imaging, 30 (2010), 1076-1089.
    [5] F. Grouiller, L. Vercueil, A. Krainik, C. Segebarth, P. Kahane, O. David, Characterization of the hemodynamic modes associated with interictal epileptic activity using a deformable model‐based analysis of combined EEG and functional MRI recordings, Hum. Brain Mapp., 31 (2010), 1157-1173.
    [6] V. L. Morgan, R. R. Price, A. Arain, P. Modur, B. Abou-Khalil, Resting functional MRI with temporal clustering analysis for localization of epileptic activity without EEG, Neuroimage, 21 (2004), 473-481.
    [7] P. P. Mitra, S. Ogawa, X. Hu, K Uǧurbil, The nature of spatiotemporal changes in cerebral hemodynamics as manifested in functional magnetic resonance imaging, Magn. Reson. Med., 37 (1997), 511-518.
    [8] M. J. McKeown, S. Makeig, G. G. Brown, T. Jung, S. S. Kindermann, A. J. Bell, et al., Analysis of fMRI data by blind separation into independent spatial components, Hum. Brain Mapp., 6 (1998), 160-188.
    [9] J. H. Gao, S. H. Yee, Iterative temporal clustering analysis for the detection of multiple response peaks in fMRI, Magn. Reson. Imaging, 21 (2003), 51-53.
    [10] S. H. Yee, J. H. Gao, Improved detection of time windows of brain responses in fMRI using modified temporal clustering analysis, Magn. Reson. Imaging, 20 (2002), 17-26.
    [11] K. Hamandi, A. S. Haddadi, A. Liston, H. Laufs, D. R. Fish, L. Lemieux, fMRI temporal clustering analysis in patients with frequent interictal epileptiform discharges: comparison with EEG-driven analysis, Neuroimage, 26 (2005), 309-316.
    [12] V. L. Morgan, J. C. Gore, B. Abou-Khalil, Cluster analysis detection of functional MRI activity in temporal lobe epilepsy, Epilepsy Res., 76 (2007), 22-33.
    [13] Q. Song, H. Chen, D. Yao, G Lu, Z Zhang, Detecting Epileptic Activities from resting fMRI Time-course by using PCA algorithm, 2005 Int. Conf. Neural Networks Brain, (2005), 1552-1555.
    [14] W. Backfrieder, R. Baumgartner, M. Samal, E. Moser, H. Bergmann, Quantification of intensity variations in functional MR images using rotated principal components, Phys. Med. Biol., 41 (1996), 1425.
    [15] R. Rodionov, F. De Martino, H. Laufs, D. W. Carmichael, E. Formisano, M. Walker, et al., Independent component analysis of interictal fMRI in focal epilepsy: Comparison with general linear model-based EEG-correlated fMRI, Neuroimage, 38 (2007), 488-500.
    [16] C. H. Zhang, Y. Lu, B. Brinkmann, K Welker, G. Worrell, B. He, Using functional MRI alone for localization in focal epilepsy, 2014 36th Annu. Int. Conf. IEEE Eng. Med. Biol. Soc., (2014), 730-733.
    [17] V. D. Calhoun, T. Eichele, T. Adalı, E. A. Allen, Decomposing the brain: Components and modes, networks and nodes, Trends Cogn. Sci.., 16 (2012), 255-256.
    [18] B. Hunyadi, S. Tousseyn, P. Dupont, S. V. Huffel, M. D. Vos, W. V. Paesschen, A prospective fMRI-based technique for localising the epileptogenic zone in presurgical evaluation of epilepsy, Neuroimage, 113 (2015), 329-339.
    [19] M. J. McKeown, T. P. Jung, S. Makeig, G. Brown, S. S. Kindermann, T. Lee, et al., Spatially independent activity patterns in functional MRI data during the Stroop color-naming task, Proc. Natl. Acad. Sci., 95 (1998), 803-810.
    [20] C. H. Zhang, Y. Lu, B. Brinkmann, K. Welker, G. Worrell, B. He, Lateralization and localization of epilepsy related hemodynamic foci using presurgical fMRI, Clin. Neurophysiol., 126 (2015), 27-38.
    [21] F. D. Martino, F. Gentile, F. Esposito, M. Balsi, F. D. Salle, R. Goebel, et al., Classification of fMRI independent components using IC-fingerprints and support vector machine classifiers, Neuroimage, 34 (2007), 177-194.
    [22] D. Cordes, V. M. Haughton, K. Arfanakis, J. D. Carew, P. A. Turski, C. H. Moritz, et al., Frequencies contributing to functional connectivity in the cerebral cortex in "resting-state" data, Am. J. Neuroradiol., 22 (2001), 1326-1333.
    [23] M. Negishi, R. Martuzzi, E. J. Novotny, D. D. Spencer, R. T. Constable, Functional MRI connectivity as a predictor of the surgical outcome of epilepsy, Epilepsia, 52 (2011), 1733-1740.
    [24] H. Chen, D. Yao, G. Lu, Z. Zhang, Q. Hu, Localization of latent epileptic activities using spatio-temporal independent component analysis of fMRI data, Brain Topogr., 19 (2006), 21-28.
    [25] X. J. Chai, A. N. Castañón, D. Öngür, S. Whitfield-Gabrieli, Anticorrelations in resting state networks without global signal regression, Neuroimage, 59 (2012), 1420-1428.
    [26] J. S. Archer, D. F. Abbott, A. B. Waites, G. D. Jackson, fMRI "deactivation" of the posterior cingulate during generalized spike and wave, Neuroimage, 20 (2003), 1915-1922.
    [27] E. Kobayashi, A. P. Bagshaw, C. Grova, F. Dubeau, J. Gotman, Negative BOLD responses to epileptic spikes, Hum. Brain Mapp., 27 (2010), 488-497.
    [28] S. Wang, Z. Zhang, G. Lu, L. Luo, Localization of brain activity by temporal anti-correlation with the posterior cingulate cortex, 2007 29th Annu. Int. Conf. IEEE Eng. Med. Biol. Soc., (2007), 5227-5230.
    [29] J. A. Maldjian, P. J. Laurienti, J. H. Burdette, Precentral gyrus discrepancy in electronic versions of the Talairach atlas, Neuroimage, 21 (2004), 450-455.
    [30] J. A. Maldjian, P. J. Laurienti, R. A. Kraft, J. H. Burdetteet, An automated method for neuroanatomic and cytoarchitectonic atlas-based interrogation of fMRI data sets, Neuroimage, 19 (2003), 1233-1239.
    [31] D. M. Cole, S. M. Smith, C. F. Beckmann, Advances and pitfalls in the analysis and interpretation of resting-state fMRI data, Front. Syst. Neurosci., 4 (2010), 8.
    [32] F. I. Karahanoglu, C. Caballero-Gaudes, F. Lazeyras, D. V. D. Ville, Total activation: fMRI deconvolution through spatio-temporal regularization, Neuroimage, 73 (2013), 121-134.
    [33] F. I. Karahanoğlu, F. Grouiller, C. C. Gaudes, M. Seeck, S. Vulliemoz, D. V. D. Ville, Spatial mapping of interictal epileptic discharges in fMRI with total activation, 2013 IEEE 10th Int. Symp. Biomed. Imaging, (2013), 1500-1503.
    [34] C. G. Bénar, D. W. Gross, Y. Wang, V. Petre, B. Pike, F. Dubeau, et al., The BOLD response to interictal epileptiform discharges, Neuroimage, 17 (2002), 1182-1192.
    [35] R. Lopes, J. M. Lina, F. Fahoum, J. Cotman, Detection of epileptic activity in fMRI without recording the EEG, Neuroimage, 60 (2012), 1867-1879.
    [36] I. Khalidov, J. Fadili, F. Lazeyras, D. V. D. Ville, M. Unser, Activelets: Wavelets for sparse representation of hemodynamic responses, Signal Process., 91 (2011), 2810-2821.
    [37] G. Bettus, E. Guedj, F. Joyeux, S. Confort-Gouny, E. Soulier, V. Laguitton, et al., Decreased basal fMRI functional connectivity in epileptogenic networks and contralateral compensatory mechanisms, Hum. Brain Mapp., 30 (2009), 1580-1591.
    [38] G. Bettus, F. Bartolomei, S. Confort-Gouny, E. Guedj, P. Chauvel, P. J. Cozzone, et al., Role of resting state functional connectivity MRI in presurgical investigation of mesial temporal lobe epilepsy, J. Neurol., Neurosurg., 81 (2010), 1147-1154.
    [39] F. Pittau, C. Grova, F. Moeller, F. Dubeau, Jean Gotman, Patterns of altered functional connectivity in mesial temporal lobe epilepsy, Epilepsia, 53 (2012), 1013-1023.
    [40] F. R. Pereira, A. Alessio, M. S. Sercheli, T. Pedro, E. Bilevicius, J. M. Rondina, et al., Asymmetrical hippocampal connectivity in mesial temporal lobe epilepsy: evidence from resting state fMRI, BMC Neurosci., 11 (2010), 66.
    [41] W. Liao, Z. Zhang, Z. Pan, D. Mantini, J. Ding, X. Duan, et al., Altered functional connectivity and small-world in mesial temporal lobe epilepsy, PLoS One, 5 (2010), e8525.
    [42] V. L. Morgan, B. P. Rogers, H. H. Sonmezturk, J. C. Gore, B Abou-Khalil, Cross hippocampal influence in mesial temporal lobe epilepsy measured with high temporal resolution functional magnetic resonance imaging, Epilepsia, 52 (2011), 1741-1749.
    [43] L. Maccotta, B. J. He, A. Z. Snyder, L. N. Eisenman, T. L. Benzinger, B. M. Ances, et al., Impaired and facilitated functional networks in temporal lobe epilepsy, Neuro Image Clin., 2 (2013), 862-872.
    [44] S. M. Stufflebeam, H. Liu, J. Sepulcre, N. Tanaka, R. L. Buckner, J. R. Madsen, Localization of focal epileptic discharges using functional connectivity magnetic resonance imaging, J. Neurosurg., 114 (2011), 1693-1697.
    [45] S. Nedic, S. M. Stufflebeam, C. Rondinoni, T. R. Velasco, A. C. Santos, J. P. Leite, et al., Using network dynamic fMRI for detection of epileptogenic foci, BMC Neurol., 15 (2015), 262.
    [46] D. Hartman, J. Hlinka, M. Paluš, D. Mantini, M. Corbetta, The role of nonlinearity in computing graph-theoretical properties of resting-state functional magnetic resonance imaging brain networks, Chaos: Int. J. Nonlinear Sci., 21 (2011), 013119.
    [47] M. Ke, X. Duan, F. Zhang, X. Yang, The research of the generalized tonic-clonic seizure epilepsy by resting state fMRI, First Int. Conf. Inf. Sci., Mach., Mater. Energy, (2015).
    [48] R. L. Buckner, J. Sepulcre, T. Talukdar, F. M. Krienen, H. Liu, T. Hedden, et al., Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer's disease, J. Neurosci., 29 (2009), 1860-1873. doi: 10.1523/JNEUROSCI.5062-08.2009
    [49] J. Sepulcre, H. S. Liu, T. Talukdar, I. Martincorena,B. T. T. Yeo, R. L. Buckner, The organization of local and distant functional connectivity in the human brain, PLoS Comput. Biol., 6 (2010), e1000808.
    [50] R. P. Carne, T. J. O'Brien, C. J. Kilpatrick, L. R. MacGregor, R. J. Hicks, M. A. Murphy, et al., MRI-negative PET-positive temporal lobe epilepsy: a distinct surgically remediable syndrome, Brain, 127 (2004), 2276-2285. doi: 10.1093/brain/awh257
    [51] K. S. Hong, S. K. Lee, J. Y. Kim, D. S. Lee, C. K. Chung, Pre-surgical evaluation and surgical outcome of 41 patients with non-lesional neocortical epilepsy, Seizure,11 (2002), 184-192.
    [52] R. L. Kutsy, Focal extratemporal epilepsy: clinical features, EEG patterns, and surgical approach, J. Neurol. Sci., 166 (1999), 1-15.
    [53] J. T. Kim, S. J. Bai, K. O. Choi, Y. J. Lee, H. J. Park, D, S. Kim, et al., Comparison of various imaging modalities in localization of epileptogenic lesion using epilepsy surgery outcome in pediatric patients, Seizure, 18 (2009), 504-510.
    [54] I. Y. Capraz, G. Kurt, O. Akdemir, T. Hirfanoglu, Y. Oner, T. Sengezer, et al., Surgical outcome in patients with MRI-negative, PET-positive temporal lobe epilepsy, Seizure, 29 (2015), 63-68.
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5507) PDF downloads(691) Cited by(5)

Figures and Tables

Tables(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog