Open access

Nonnegative Inverse Eigenvalue Problem

Written By

Ricardo L. Soto

Submitted: 17 November 2011 Published: 11 July 2012

DOI: 10.5772/48279

From the Edited Volume

Linear Algebra - Theorems and Applications

Edited by Hassan Abid Yasser

Chapter metrics overview

3,373 Chapter Downloads

View Full Metrics

1. Introduction

Nonnegative matrices have long been a sorce of interesting and challenging mathematical problems. They are real matrices with all their entries being nonnegative and arise in a number of important application areas: communications systems, biological systems, economics, ecology, computer sciences, machine learning, and many other engineering systems. Inverse eigenvalue problems constitute an important subclass of inverse problems that arise in the context of mathematical modeling and parameter identification. A simple application of such problems is the construction of Leontief models in economics [1]-[2].

The nonnegative inverse eigenvalue problem (NIEP) is the problem of characterizing those lists Λ={λ 1 ,λ 2 ,...,λ n } of complex numbers which can be the spectra of n×n entrywise nonnegative matrices. If there exists a nonnegative matrix A with spectrum Λ we say that Λ is realized by A and that A is the realizing matrix. A set Kof conditions is said to be a realizability criterion if any list Λ={λ 1 ,λ 2 ,...,λ n }, real or complex, satisfying conditions K is realizable. The NIEP is an open problem. A full solution is unlikely in the near future. The problem has only been solved for n=3 by Loewy and London ([3], 1978) and for n=4 by Meehan ([4], 1998) and Torre-Mayo et al.([5], 2007). The case n=5 has been solved for matrices of trace zero in ([6], 1999). Other results, mostly in terms of sufficient conditions for the problem to have a solution (in the case of a complex list Λ), have been obtained, in chronological order, in [7]-[8].

Two main subproblems of the NIEP are of great interest: the real nonnegative inverse eigenvalue problem (RNIEP), in which Λ is a list of real numbers, and the symmetric nonnegative inverse eigenvalue problem (SNIEP), in which the realizing matrix must be symmetric. Both problems, RNIEP and SNIEP are equivalent for n4 (see [9]), but they are different otherwise (see [10]). Moreover, both problems remains unsolved for n5. The NIEP is also of interest for nonnegative matrices with a particular structure, like stochastic and doubly stochastic, circulant, persymmetric, centrosymmetric, Hermitian, Toeplitz, etc.

The first sufficient conditions for the existence of a nonnegative matrix with a given real spectrum (RNIEP) were obtained by Suleimanova ([11], 1949) and Perfect ([12], [13], 1953 and 1955). Other sufficient conditions have also been obtained, in chronological order in [14]-[15], (see also [16], [17], and references therein for a comprehensive survey).

The first sufficient conditions for the SNIEP were obtained by Fiedler ([18], 1974). Other results for symmetric realizability have been obtained in [19], [7] and [20]-[21]. Recently, new sufficient conditions for the SNIEP have been given in [22]-[23].

1.1. Necessary conditions

Let A be a nonnegative matrix with spectrum Λ={λ 1 ,λ 2 ,...,λ n }. Then, from the Perron Frobenius theory we have the following basic necessary conditions

(1)Λ ¯={λ 1 ¯,...,λ n ¯}=Λ(2)max j {λ j }Λ(3)s m (Λ)= j=1 n λ j m 0,m=1,2,...,uid2

where Λ ¯=Λ means that Λ is closed under complex comjugation.

Moreover, we have

(4)(s k (Λ)) m n m-1 s km (Λ),k,m=1,2,...(5)(s 2 (Λ)) 2 (n-1)s 4 (Λ),nodd,tr(A)=0.uid3

Necessary condition (4) is due to Loewy and London [3]. Necessary condition (5), which is a refinement of (4), is due to Laffey and Meehan [24]. The list Λ={5,4,-3,-3,-3} for instance, satisfies all above necessary conditions, except condition (5). Therefore Λ is not a realizable list. In [25] it was obtained a new necessary condition, which is independent of the previous ones. This result is based on the Newton's inequalities associated to the normalized coefficients of the characteristic polynomial of an M-matrix or an inverse M-matrix.

The chapter is organized as follows: In section 2 we introduce two important matrix results, due to Brauer and Rado, which have allowed to obtain many of the most general sufficient conditions for the RNIEP, the SNIEP and the complex case. In section 3 we consider the real case and we introduce, without proof (we indicate where the the proofs can be found), two sufficient conditions with illustrative examples. We consider, in section 4, the symmetric case. Here we introduce a symmetric version of the Rado result, Theorem , and we set, without proof (see the appropriate references), three sufficient conditions, which are, as far as we know, the most general sufficient conditions for the SNIEP. In section 5, we discuss the complex (non real) case. Here we present several results with illustrative examples. Section 6 is devoted to discuss some Fiedler results and Guo results, which are very related with the problem and have been employed with success to derive sufficient conditions. Finally, in section 7, we introduce some open questions.

Advertisement

2. Brauer and Rado Theorems

A real matrix A=(a ij ) i=1 n is said to have constant row sums if all its rows sum up to the same constant, say, α, that is, j=1 n a ij =α, i=1,...,n. The set of all real matrices with constant row sums equal to α is denoted by 𝒞𝒮 α . It is clear that any matrix in 𝒞𝒮 α has eigenvector 𝐞=(1,1,...,1) T corresponding to the eigenvalue α. Denote by 𝐞 k the n-dimensional vector with one in the k-th position and zeros elsewhere.

It is well known that the problem of finding a nonnegative matrix with spectrum Λ={λ 1 ,...,λ n } is equivalent to the problem of finding a nonnegative matrix in 𝒞𝒮 λ 1 with spectrum Λ (see [26]). This will allow us to exploit the advantages of two important theorems, Brauer Theorem and Rado Theorem, which will be introduced in this section.

The spectra of circulant nonnegative matrices have been characterized in [27], while in [28], a simple complex generalization of Suleimanova result has been proved, and efficient and general sufficient conditions for the realizability of partitioned spectra, with the partition allowing some of its pieces to be nonrealizable, provided there are other pieces, which are realizable and, in certain way, compensate the nonnrealizability of the former, have been obtained. This is the procedure which we call negativity compensation. This strategy, based in the use of the following two perturbation results, together with the properties of real matrices with constant row sums, has proved to be successful.

Theorem 1 Brauer [29] Let A be an n×n arbitrary matrix with eigenvalues λ 1 ,...,λ n . Let 𝐯=(v 1 ,...,v n ) T an eigenvector of A associated with the eigenvalue λ k and let 𝐪=(q 1 ,...,q n ) T be any n-dimensional vector. Then the matrix A+𝐯𝐪 T has eigenvalues λ 1 ,...,λ k-1 ,λ k +v T q,λ k+1 ,...,λ n .

Let U be an n×n nonsingular matrix such that

U -1 AU=λ 1 **λ 2 *λ n

is an upper triangular matrix, where we choose the first column of U as 𝐯 (U there exists from a well known result of Schur). Then,

U -1 (A+𝐯𝐪 T )U=U -1 AU+q 1 q 2 q n U=λ 1 +𝐪 T 𝐯**λ 2 *λ n .

and the result follows.

This proof is due to Reams [30].

Theorem 2 Rado [13] Let A be an n×n arbitrary matrix with eigenvalues λ 1 ,...,λ n and let Ω=diag{λ 1 ,...,λ r } for some rn. Let X be an n×r matrix with rank r such that its columns x 1 ,x 2 ,...,x r satisfy Ax i =λ i x i , i=1,...,r. Let C be an r×n arbitrary matrix. Then the matrix A+XC has eigenvalues μ 1 ,...,μ r ,λ r+1 ,...,λ n , where μ 1 ,...,μ r are eigenvalues of the matrix Ω+CX.

Let S=XY a nonsingular matrix with S -1 =[]UV. Then UX=I r , VY=I n-r and VX=0, UY=0. Let C=C 1 C 2 , X=[]X 1 X 2 , Y=[]Y 1 Y 2 . Then, since AX=XΩ,

S -1 AS=[]UVXΩAY=ΩUAY0VAY

and

S -1 XCS=[]I r 0C 1 C 2 S=C 1 C 2 00X 1 Y 1 X 2 Y 2 =CXCY00.

Thus,

S -1 (A+XC)S=S -1 AS+S -1 XCS=Ω+CXUAY+CY0VAY,

and we have σ(A+XC)=σ(Ω+CX)+σ(A)-σ(Ω).

Advertisement

3. Real nonnegative inverse eigenvalue problem.

Regarding the RNIEP, by applying Brauer Theorem and Rado Theorem, efficient and general sufficient conditions have been obtained in [13], [31], [32], [33].

Theorem 3 [32] Let Λ={λ 1 ,λ 2 ,...,λ n } be a given list of real numbers. Suppose that:

i) There exists a partition Λ=Λ 1 ...Λ t , where

Λ k ={λ k1 ,λ k2 ,...λ kp k },λ 11 =λ 1 ,λ k1 λ kp k ,λ k1 0,

k=1,...,t, such that for each sublist Λ k we associate a corresponding list

Γ k ={ω k ,λ k2 ,...,λ kp k },0ω k λ 1 ,

which is realizable by a nonnegative matrix A k CS ω k of order p k .

ii) There exists a nonnegative matrix BCS λ 1 with eigenvalues λ 1 ,λ 21 ,...,λ t1 (the first elements of the lists Λ k ) and diagonal entries ω 1 ,ω 2 ,...,ω t (the first elements of the lists Γ k ).

Then Λ is realizable by a nonnegative matrix ACS λ 1 .

Perfect [13] gave conditions under which λ 1 ,λ 2 ,...,λ t and ω 1 ,ω 2 ,...,ω t are the eigenvalues and the diagonal entries, respectively, of a t×t nonnegative matrix BCS λ 1 . For t=2 it is necessary and sufficient that λ 1 +λ 2 =ω 1 +ω 2 , with 0ω i λ 1 . For t=3 Perfect gave the following result:

Theorem 4 [13] The real numbers λ 1 ,λ 2 ,λ 3 and ω 1 ,ω 2 ,ω 3 are the eigenvalues and the diagonal entries, respectively, of a 3×3 nonnegative matrix BCS λ 1 , if and only if:

i)0ω i λ 1 ,i=1,2,3ii)λ 1 +λ 2 +λ 3 =ω 1 +ω 2 +ω 3 iii)λ 1 λ 2 +λ 1 λ 3 +λ 2 λ 3 ω 1 ω 2 +ω 1 ω 3 +ω 2 ω 3 iv)max k ω k λ 2 uid8

Then, an appropriate 3×3 nonnegative matrix B is

B=ω 1 0λ 1 -ω 1 λ 1 -ω 2 -pω 2 p0λ 1 -ω 3 ω 3 ,uid9

where

p=1 λ 1 -ω 3 (ω 1 ω 2 +ω 1 ω 3 +ω 2 ω 3 -λ 1 λ 2 +λ 1 λ 3 +λ 2 λ 3 ).

For t4, we only have a sufficient condition:

i)0ω k λ 1 ,k=1,2,...t,ii)ω 1 +ω 2 +ω t =λ 1 +λ 2 +λ t ,iii)ω k λ k ,ω 1 λ k ,k=2,3,...,t,uid10

with the following matrix BCS λ 1 having eigenvalues and diagonal entries λ 1 ,λ 2 ,...,λ t and ω 1 ,ω 2 ,...,ω t , respectively:

B=ω 1 ω 2 -λ 2 ω r -λ t ω 1 -λ 2 ω 2 ω r -λ t ω 1 -λ t ω 2 -λ 2 ω t .uid11

Example 1 Let us consider the list Λ={6,1,1,-4,-4} with the partition

Λ 1 ={6,-4},Λ 2 ={1,-4},Λ 3 ={1}

and the realizable associated lists

Γ 1 ={4,-4},Γ 2 ={4,-4},Γ 3 ={0}.

From () we compute the 3×3 nonnegative matrix

B=4023 241 2060

with eigenvalues 6,1,1, and diagonal entries 4,4,0. Then

A=0400040000000400040000000+100100010010001000023 20001 200600=04002400023 20041 23 20401 200600

is nonnegative with spectrum Λ.

A map of sufficient conditions for the RNIEP it was constructed in [17], There, the sufficient conditions were compared to establish inclusion or independence relations between them. It is also shown in [17] that the criterion given by Theorem contains all realizability criteria for lists of real numbers studied therein. In [33], from a new special partition, Theorem is extended. Now, the first element λ k1 of the sublist Λ k need not to be nonnegative and the realizable auxiliar list Γ k ={ω k ,λ k1 ,...,λ kp k } contains one more element. Moreover, the number of lists of the partition depend on the number of elements of the first list Λ 1 , and some lists Λ k can be empty.

Theorem 5 [33] Let Λ={λ 1 ,λ 2 ,...,λ n } be a list of real numbers and let the partition Λ=Λ 1 Λ p 1 +1 be such that

Λ k ={λ k1 ,λ k2 ,...λ kp k },λ 11 =λ 1 ,λ k1 λ k2 λ kp k ,

k=1,...,p 1 +1, where p 1 is the number of elements of the list Λ 1 and some of the lists Λ k can be empty. Let ω 2 ,...,ω p 1 +1 be real numbers satisfying 0ω k λ 1 , k=2,...,p 1 +1. Suppose that the following conditions hold:

i) For each k=2,...,p 1 +1, there exists a nonnegative matrix A k CS ω k with spectrum Γ k ={ω k ,λ k1 ,...,λ kp k },

ii) There exists a p 1 ×p 1 nonnegative matrix BCS λ 1 , with spectrum Λ 1 and with diagonal entries ω 2 ,...,ω p 1 +1 .

Then Λ is realizable by a nonnegative matrix ACS λ 1 .

Example 2 With this extension, the authors show for instance, that the list

{5,4,0,-3,-3,-3}

is realizable, which can not be done from the criterion given by Theorem . In fact, let the partition

Λ 1 ={5,4,0,-3},Λ 2 ={-3},Λ 3 ={-3}withΓ 2 ={3,-3},Γ 3 ={3,-3},Γ 4 =Γ 5 ={0}.

The nonnegative matrix

B=3020030230020320

has spectrum Λ 1 and diagonal entries 3,3,0,0. It is clear that

A 2 =A 3 =0330realizesΓ 2 =Γ 3 .

Then

A=A 2 A 3 00+100010000100010000100001000020000002300002003020=030020300020000302003002300002003020

has the desired spectrum {5,4,0,-3,-3,-3}.

Advertisement

4. Symmetric nonnegative inverse eigenvalue problem

Several realizability criteria which were first obtained for the RNIEP have later been shown to be symmetric realizability criteria as well. For example, Kellogg criterion [14] was showed by Fiedler [18] to imply symmetric realizability. It was proved by Radwan [7] that Borobia's criterion [34] is also a symmetric realizability criterion, and it was proved in [21] that Soto's criterion for the RNIEP is also a criterion for the SNIEP. In this section we shall consider the most general and efficient symmetric realizability criteria for the SNIEP (as far as we know they are). We start by introducing a symmetric version of the Rado Theorem:.

Theorem 6 [22] Let A be an n×n symmetric matrix with spectrum Λ={λ 1 ,λ 2 ,...,λ n } and for some rn, let {𝐱 1 ,𝐱 2 ,...,𝐱 r } be an orthonormal set of eigenvectors of A spanning the invariant subspace associated with λ 1 ,λ 2 ,...,λ r . Let X be the n×r matrix with i-th column 𝐱 i , let Ω=diag{λ 1 ,...,λ r }, and let C be any r×r symmetric matrix. Then the symmetric matrix A+XCX T has eigenvalues μ 1 ,μ 2 ,...,μ r ,λ r+1 ,...,λ n , where μ 1 ,μ 2 ,...,μ r are the eigenvalues of the matrix Ω+C.

Since the columns of X are an orthonormal set, we may complete X to an orthogonal matrix W=[X Y], that is, X T X=I r , Y T Y=I n-r , X T Y=0, Y T X=0. Then

W -1 AW=[]X T Y T AXY=ΩX T AY0Y T AYW -1 (XCX T )W=[]I r 0CI r 0=C000.

Therefore,

W -1 (A+XCX T )W=Ω+CX T AY0Y T AY

and A+XCX T is symmetric with eigenvalues μ 1 ,...,μ r ,λ r+1 ,...,λ n .

By using Theorem , the following sufficient condition was proved in [22]:

Theorem 7 [22] Let Λ={λ 1 ,λ 2 ,...,λ n } be a list of real numbers with λ 1 λ 2 λ n and, for some tn, let ω 1 ,...,ω t be real numbers satisfying 0ω k λ 1 , k=1,...,t. Suppose there exists:

i) a partition Λ=Λ 1 Λ t with

Λ k ={λ k1 ,λ k2 ,...λ kp k },λ 11 =λ 1 ,λ k1 0,λ k1 λ k2 λ kp k ,

such that for each k=1,...,t, the list Γ k ={ω k ,λ k2 ,...,λ kp k } is realizable by a symmetric nonnegative matrix A k of order p k , and

ii) a t×t symmetric nonnegative matrix B with eigenvalues λ 11 ,λ 21 ,...,λ t1 } and with diagonal entries ω 1 ,ω 2 ,...,ω t .

Then Λ is realizable by a symmetric nonnegative matrix.

Since A k is a p k ×p k symmetric nonnegative matrix realizing Γ k , then A=diag{A 1 ,A 2 ,...,A t } is symmetric nonnegative with spectrum Γ 1 Γ 2 Γ t . Let {𝐱 1 ,...,𝐱 t } be an orthonormal set of eigenvectors of A associated with ω 1 ,...,ω t , respectively. Then the n×t matrix X with i-th column 𝐱 i satisfies AX=XΩ for Ω=dig{ω 1 ,...,ω t }. Moreover, X is entrywise nonnegative, since each 𝐱 i contains the Perron eigenvector of A i and zeros. Now, if we set C=B-Ω, the matrix C is symmetric nonnegative and Ω+C has eigenvalues λ 1 ,...,λ t . Therefore, by Theorem the symmetric matrix A+XCX T has spectrum Λ. Besides, it is nonnegative since all the entries of A,X, and C are nonnegative.

Theorem not only ensures the existence of a realizing matrix, but it also allows to construct the realizing matrix. Of course, the key is to know under which conditions does there exists a t×t symmetrix nonnegative matrix B with eigenvalues λ 1 ,...,λ t and diagonal entries ω 1 ,...,ω t .

The following conditions for the existence of a real symmetric matrix, not necessarily nonnegative, with prescribed eigenvalues and diagonal entries are due to Horn [35]: There exists a real symmetric matrix with eigenvalues λ 1 λ 2 λ t and diagonal entries ω 1 ω 2 ω t if and only if

i=1 k λ i i=1 k ω i ,k=1,...,t-1 i=1 t λ i = i=1 t ω i uid17

For t=2, the conditions () become

λ 1 ω 1 ,λ 1 +λ 2 =ω 1 +ω 2 ,

and they are also sufficient for the existence of a 2×2 symmetric nonnegative matrix B with eigenvalues λ 1 λ 2 and diagonal entries ω 1 ω 2 0, namely,

B=ω 1 (λ 1 -ω 1 )(λ 1 -ω 2 )(λ 1 -ω 1 )(λ 1 -ω 2 )ω 2 .

For t=3, we have the following conditions:

Lemma 1 [18] The conditions

λ 1 ω 1 λ 1 +λ 2 ω 1 +ω 2 λ 1 +λ 2 +λ 3 =ω 1 +ω 2 +ω 3 ω 1 λ 2 uid19

are necessary and sufficient for the existence of a 3×3 symmetric nonnegative matrix B with eigenvalues λ 1 λ 2 λ 3 and diagonal entries ω 1 ω 2 ω 3 0.

In [22], the following symmetric nonnegative matrix B, satisfying conditions (), it was constructed:

B=ω 1 μ-ω 3 2μ-ω 2 -ω 3 sμ-ω 2 2μ-ω 2 -ω 3 sμ-ω 3 2μ-ω 2 -ω 3 sω 2 (μ-ω 2 )(μ-ω 3 )μ-ω 2 2μ-ω 2 -ω 3 s(μ-ω 2 )(μ-ω 3 )ω 3 ,uid20

where μ=λ 1 +λ 2 -ω 1 ; s=(λ 1 -μ)(λ 1 -ω 1 ).

For t4 we have only a sufficient condition:

Theorem 8 Fiedler [18] If λ 1 λ t and ω 1 ω t satisfy

i) i=1 s λ i i=1 s ω i ,s=1,...,t-1ii) i=1 t λ i = i=1 t ω i iii)ω k-1 λ k ,k=2,...,t-1,uid22

then there exists a t×t symmetric nonnegative matrix with eigenvalues λ 1 ,...,λ t and diagonal entries ω 1 ,...,ω t .

Observe that

B=521 21 2251 21 21 21 2521 21 225

has eigenvalues 8,6,3,3, but λ 2 =6>5=ω 1 .

Example 3 Let us consider the list Λ={7,5,1,-3,-4,-6} with the partition

Λ 1 ={7,-6},Λ 2 ={5,-4},Λ 3 ={1,-3}withΓ 1 ={6,-6},Γ 2 ={4,-4},Γ 3 ={3,-3}.

We look for a symmetric nonnegative matrix B with eigenvalues 7,5,1 and diagonal entries 6,4,3. Then conditions () are satisfied and from () we compute

B=63 52 53 5462 563andC=B-Ω,

where Ω=diag{6,4,3}. The symmetric matrices

A 1 =0660,A 2 =0440,A 3 =0330

realize Γ 1 ,Γ 2 ,Γ 3 . Then

A=A 1 A 2 A 3 +XCX T ,whereX=2 2002 20002 2002 20002 2002 2,

is symmetric nonnegative with spectrum Λ.

In the same way as Theorem was extended to Theorem (in the real case), Theorem was also extended to the following result:

Theorem 9 [33] Let Λ={λ 1 ,λ 2 ,...,λ n } be a list of real numbers and let the partition Λ=Λ 1 Λ p 1 +1 be such that

Λ k ={λ k1 ,λ k2 ,...λ kp k },λ 11 =λ 1 ,λ k1 λ k2 λ kp k ,

k=1,...,p 1 +1, where Λ 1 is symmetrically realizable, p 1 is the number of elements of Λ 1 and some lists Λ k can be empty. Let ω 2 ,...,ω p 1 +1 be real numbers satisfying 0ω k λ 1 , k=2,...,p 1 +1. Suppose that the following conditions hold:

i) For each k=2,...,p 1 +1, there exists a symmetric nonnegative matrix A k with spectrum Γ k ={ω k ,λ k1 ,...,λ kp k },

ii) There exists a p 1 ×p 1 symmetric nonnegative matrix B with spectrum Λ 1 and with diagonal entries ω 2 ,...,ω p 1 +1 .

Then Λ is symmetrically realizable.

Example 4 Now, from Theorem , we can see that there exists a symmetric nonnegative matrix with spectrum Λ={5,4,0,-3,-3,-3}, which can not be seen from Theorem . Moreover, we can compute a realizing matrix. In fact, let the partition

Λ 1 ={5,4,0,-3},Λ 2 ={-3},Λ 3 ={-3}withΓ 2 ={3,-3},Γ 3 ={3,-3},Γ 4 =Γ 5 ={0}.

The symmetric nonnegative matrix

B=3060030660020620

has spectrum Λ 1 and diagonal entries 3,3,0,0. Let Ω=diag{3,3,0,0} and

X=2 20002 200002 20002 20000100001,A 2 =A 3 =0330,C=B-Ω.

Then, from Theorem we obtain

A=A 2 A 3 00+XCX T ,

which is symmetric nonnegative with spectrum Λ.

The following result, although is not written in the fashion of a sufficient condition, is indeed a very general and efficient sufficient condition for the SNIEP.

Theorem 10 [36] Let A be an n×n irreducible symmetric nonnegative matrix with spectrum Λ={λ 1 ,λ 2 ,...,λ n }, Perron eigenvalue λ 1 and a diagonal element c. Let B be an m×m symmetric nonnegative matrix with spectrum Γ={μ 1 ,μ 2 ,...,μ m } and Perron eigenvalue μ 1 .

i) If μ 1 c, then there exists a symmetric nonnegative matrix C, of order (n+m-1), with spectrum {λ 1 ,...,λ n ,μ 2 ,...,μ m }.

ii) If μ 1 c, then there exists a symmetric nonnegative matrix C, of order (n+m-1), with spectrum {λ 1 +μ 1 -c,λ 2 ,...,λ n ,μ 2 ,...,μ m }.

Example 5 The following example, given in [36], shows that

{7,5,0,-4,-4,-4} with the partition

Λ={7,5,0,-4},Γ={4,-4},

satisfies conditions of Theorem , where

A=40b0040db0060d60withb 2 +d 2 =23,bd=46,

is symmetric nonnegative with spectrum Λ. Then there exists a symmetric nonnegative matrix C with spectrum {7,5,0,-4,-4} and a diagonal element 4. By applying again Theorem to the lists {7,5,0,-4,-4} and {4,-4}, we obtain the desired symmetric nonnegative matrix.

It is not hard to show that both results, Theorem and Theorem , are equivalent (see [37]). Thus, the list in the Example is also realizable from Theorem , while the list in the example is also realizable from Theorem .

Advertisement

5. List of complex numbers

In this section we consider lists of complex nonreal numbers. We start with a complex generalization of a well known result of Suleimanova, usually considered as one of the important results in the RNIEP (see [11]): The list λ 1 >0>λ 2 λ n is the spectrum of a nonnegative matrix if and only if λ 1 +λ 2 ++λ n 0.

Theorem 11 [28] Let Λ={λ 0 ,λ 1 ,...,λ n } be a list of complex numbers closed under complex conjugation, with

Λ ' ={λ 1 ,...,λ n }{z:Rez0;RezImz}.

Then Λ is realizable if and only if i=0 n λ i 0.

Suppose that the elements of Λ ' are ordered in such a way that λ 2p+1 ,...,λ n are real and λ 1 ,...,λ 2p are complex nonreal, with

x k =Reλ 2k-1 =Reλ 2k andy k =Imλ 2k-1 =Imλ 2k

for k=1,...,p. Consider the matrix

B=000.-x 1 +y 1 x 1 -y 1 .-x 1 -y 1 y 1 x 1 .-x p +y p 00.x p -y p -x p -y p 00.y p x p -λ 2p+1 00.λ 2p+1 .-λ n 0.λ n .

It is clear that B𝒞𝒮 0 with spectrum {0,λ 1 ,...,λ n } and all the entries on its first column are nonnegative. Define 𝐪=(q 0 ,q 1 ,...,q n ) T with q 0 =λ 0 + i=1 n λ i and

q k =-Reλ k fork=1,...,2pandq k =-λ k fork=2p+1,...,n.

Then, from the Brauer Theorem A=B+𝐞𝐪 T is nonnegative with spectrum Λ.

In the case when all numbers in the given list, except one (the Perron eigenvalue), have real parts smaller than or equal to zero, remarkably simple necessary and sufficient conditions were obtained in [38].

Theorem 12 [38] Let λ 2 ,λ 3 ,...,λ n be nonzero complex numbers with real parts less than or equal to zero and let λ 1 be a positive real number. Then the list Λ={λ 1 ,λ 2 ,...,λ n } is the nonzero spectrum of a nonnegative matrix if the following conditions are satisfied:

i)Λ ¯=Λii)s 1 = i=1 n λ i 0iii)s 2 = i=1 n λ i 2 0uid30

The minimal number of zeros that need to be added to Λ to make it realizable is the smallest nonnegative integer N for which the following inequality is satisfied:

s 1 2 (n+N)s 2 .

Furthermore, the list {λ 1 ,λ 2 ,...,λ n ,0,...,0} can be realized by C+αI, where C is a nonnegative companion matrix with trace zero, α is a nonnegative scalar and I is the n×n identity matrix.

Corollary 1 [38] Let λ 2 ,λ 3 ,...,λ n be complex numbers with real parts less than or equal to zero and let λ 1 be a positive real number. Then the list Λ={λ 1 ,λ 2 ,...,λ n } is the spectrum of a nonnegative matrix if and only if the following conditions are satisfied:

i)Λ ¯=Λii)s 1 = i=1 n λ i 0iii)s 2 = i=1 n λ i 2 0iv)s 1 2 ns 2 uid32

Example 6 The list Λ={8,-1+3i,-1-3i,-2+5i,-2-5i} satisfies conditions (). Then Λ is the spectrum of the nonnegative companion matrix

C=01000001000001000001232049427812.

Observe that Theorem gives no information about the realizability of Λ.

The list {19,-1+11i,-1-11i,-3+8i,-3-8i} was given in [38]. It does not satisfy conditions (): s 1 =11, s 2 =11 and s 1 2 ns 2 . The inequality 11 2 (5+N)11 is satisfied for N6. Then we need to add 6 zeros to the list to make it realizable.

Theorem (in section 3), can also be extended to the complex case:

Theorem 13 [8] Let Λ={λ 2 ,λ 3 ,...,λ n } be a list of complex numbers such that Λ ¯=Λ, λ 1 max i λ i , i=2,...,n, and i=1 n λ i 0. Suppose that:

i) there exists a partition Λ=Λ 1 Λ t with

Λ k ={λ k1 ,λ k2 ,...λ kp k },λ 11 =λ 1 ,

k=1,...,t, such that Γ k ={ω k ,λ k2 ,...,λ kp k } is realizable by a nonnegative matrix A k 𝒞𝒮 ω k , and

ii) there exists a t×t nonnegative matrix B𝒞𝒮 λ 1 , with eigenvalues

λ 1 ,λ 21 ,...,λ t1 (the first elements of the lists Λ k ) and with diagonal entries ω 1 ,ω 2 ,...,ω t (the Perron eigenvalues of the lists Γ k ).

Then Λ is realizable.

Example 7 Let Λ={7,1,-2,-2,-2+4i,-2-4i}. Consider the partition

Λ 1 ={7,1,-2,-2},Λ 2 ={-2+4i},Λ 3 ={-2-4i}withΓ 1 ={3,1,-2,-2},Γ 2 ={0},Γ 3 ={0}.

We look for a nonnegative matrix B𝒞𝒮 7 with eigenvalues 7,-2+4i,-2-4i and diagonal entries 3,0,0, and a nonnegative matrix A 1 realizing Γ 1 . They are

B=30441 708 7070andA 1 =0201200101020120.

Then

A=A 1 00+10010010010001000100000441 700008 7000070=02010420010401020401200441 700008 7000070

has the spectrum Λ.

Advertisement

6. Fiedler and Guo results

One of the most important works about the SNIEP is due to Fiedler [18]. In [18] Fiedler showed, as it was said before, that Kellogg sufficient conditions for the RNIEP are also sufficient for the SNIEP. Three important and very useful results of Fiedler are:

Lemma 2 [18] Let A be a symmetric m×m matrix with eigenvalues α 1 ,...,α m , A𝐮=α 1 𝐮, 𝐮=1. Let B be a symmetric n×n matrix with eigenvalues β 1 ,...,β n , B𝐯=β 1 𝐯, 𝐯=1. Then for any ρ, the matrix

C=Aρ𝐮𝐯 T ρ𝐯𝐮 T B

has eigenvalues α 2 ,...,α m ,β 2 ,...,β n ,γ 1 ,γ 2 , where γ 1 ,γ 2 are eigenvalues of the matrix

C ˜=α 1 ρρβ 1 .

Lemma 3 [18] If {α 1 ,...,α m } and {β 1 ,...,β n } are lists symmetrically realizable and α 1 β 1 , then for any t0, the list

{α 1 +t,β 1 -t,α 2 ,...,α m ,β 2 ,...,β n }

is also symmetrically realizable.

Lemma 4 [18] If Λ={λ 1 ,λ 2 ,...,λ n } is symmetrically realizable by a nonnegative matrix and if t>0, then

Λ t ={λ 1 +t,λ 2 ,...,λ n }

is symmetrically realizable by a positive matrix.

Remark 1 It is not hard to see that Lemma can be obtained from Theorem . In fact, it is enough to consider

C=AB+𝐮00𝐯0ρρ0𝐮 T 0 T 0 T 𝐯 T =Aρ𝐮𝐯 T ρ𝐯𝐮 T B,

which is symmetric with eigenvalues γ 1 ,γ 2 ,α 2 ,...,α m ,β 2 ,...,β n , where γ 1 ,γ 2 are eigenvalues of

B=α 1 ρρβ 1 .

Now we consider a relevant result due to Guo [39]:

Theorem 14 [39] If the list of complex numbers Λ={λ 1 ,λ 2 ,...,λ n } is realizable, where λ 1 is the Perron eigenvalue and λ 2 , then for any t0 the list Λ t ={λ 1 +t,λ 2 ±t,λ 3 ,...,λ n } is also realizable.

Corollary 2 [39] If the list of real numbers Λ={λ 1 ,λ 2 ,...,λ n } is realizable and t 1 = i=2 n t i with t i , i=2,...,n, then the list Λ t ={λ 1 +t 1 ,λ 2 +t 2 ,...,λ n +t n } is also realizable.

Example 8 Let Λ={8,6,3,3,-5,-5,-5,-5} be a given list. Since the lists Λ 1 =Λ 2 ={7,3,-5,-5} are both realizable (see [31] to apply a simple criterion, which shows the realizability of Λ 1 =Λ 2 ), then

Λ 1 Λ 2 ={7,7,3,3,-5,-5,-5,-5}

is also realizable. Now, from Theorem , with t=1, Λ is realizable.

Guo also sets the following two questions:

Question 1: Do complex eigenvalues of nonnegative matrices have a property similar to Theorem ?

Question 2: If the list Λ={λ 1 ,λ 2 ,...,λ n } is symmetrically realizable, and t>0, is the list Λ t ={λ 1 +t,λ 2 ±t,λ 3 ,...,λ n } symmetrically realizable?.

It was shown in [40] and also in [41] that Question 1 has an affirmative answer.

Theorem 15 [40] Let Λ={λ 1 ,a+bi,a-bi,λ 4 ,...,λ n } be a realizable list of complex numbers. Then for all t0, the perturbed list

Λ t ={λ 1 +2t,a-t+bi,a-t-bi,λ 4 ,...,λ n }

is also realizable.

Question 2, however, remains open. An affirmative answer to Question 2, in the case that the symmetric realizing matrix is a nonnegative circulant matrix or it is a nonnegative left circulant matrix, it was given in [42]. The use of circulant matrices has been shown to be very useful for the NIEP [27], [32]. In [32] it was given a necessary and sufficient condition for a list of 5 real numbers, which corresponds to a even-conjugate vector, to be the spectrum of 5×5 symmetric nonnegative circulant matrix:

Lemma 5 [32] Let λ=(λ 1 ,λ 2 ,λ 3 ,λ 3 ,λ 2 ) T be a vector of real numbers (even-conjugate) such that

λ 1 λ j ,j=2,3λ 1 λ 2 λ 3 λ 1 +2λ 2 +2λ 3 0uid45

A necessary and sufficient condition for {λ 1 ,λ 2 ,λ 3 ,λ 3 ,λ 2 } to be the spectrum of a symmetric nonnegative circulant matrix is

λ 1 +(λ 3 -λ 2 )5-1 2-λ 2 0.uid46

Example 9 From Lemma we may know, for instance, that the list

{6,1,1,-4,-4} is the spectrum of a symmetric nonnegative circulant matrix.

Advertisement

7. Some open questions

We finish this chapter by setting two open questions:

Question 1: If the list of real numbers Λ={λ 1 ,λ 2 ,...,λ n } is symmetrically realizable, and t>0, is the list Λ t ={λ 1 +t,λ 2 ±t,λ 3 ,...,λ n } also symmetrically realizable?

Some progress has been done about this question. In [42], it was given an affirmative answer to Question 1, in the case that the realizing matrix is symmetric nonnegative circulant matrix or it is nonnegative left circulant matrix. In [43] it was shown that if 1>λ 2 λ n 0, then Theorem holds for positive stochastic, positive doubly stochastic and positive symmetric matrices.

Question 2: How adding one or more zeros to a list can lead to its symmetric realizability by different symmetric patterned matrices?

The famous Boyle-Handelman Theorem [44] gives a nonconstructive proof of the fact that if s k =λ 1 k +λ 2 k ++λ n k >0, for k=1,2,..., then there exists a nonnegative number N for which the list {λ 1 ,...,λ n ,0,...,0}, with N zeros added, is realizable. In [38] Laffey and Šmigoc completely solve the NIEP for lists of complex numbers Λ={λ 1 ,...,λ n }, closed under conjugation, with λ 2 ,...,λ n having real parts smaller than or equal to zero. They show the existence of N0 for which Λ with N zeros added is realizable and show how to compute the least such N. The situation for symmetrically realizable spectra is different and even less is known.

Advertisement

8. Conclusion

The nonnegative inverse eigenvalue problem is an open and difficult problem. A full solution is unlikely in the near future. A number of partial results are known in the literature about the problem, most of them in terms of sufficient conditions. Some matrix results, like Brauer Theorem (Theorem ), Rado Theorem (Theorem ), and its symmetric version (Theorem ) have been shown to be very useful to derive good sufficient conditions. This way, however, seems to be quite narrow and may be other techniques should be explored and applied.

References

  1. 1. A. Berman, R. J. Plemmons (1994) Nonnegative Matrices in the Mathematical Sciences. In: Classics in Applied Mathematics 9, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA.
  2. 2. M. T. Chu, G. H. Golub (2005) Inverse eigenvalue problems: theory, algorithms and applications, Oxford University Press, New York.
  3. 3. H. Minc (1988) Nonnegative Matrices, John Wiley & Sons, New York.
  4. 4. R. Loewy, D. London (1978) A note on an inverse problem for nonnegative matrices. In: Linear and Multilinear Algebra 6 83-90.
  5. 5. M.E. Meehan (1998) Some results on matrix spectra, Ph. D. Thesis, National University of Ireland, Dublin.
  6. 6. J. Torre-Mayo, M.R. Abril-Raymundo, E. Alarcia-Estévez, C. Marijuán, M. Pisonero (2007) The nonnegative inverse eigenvalue problem from the coefficients of the characteristic polynomial. EBL digraphs In: Linear Algebra Appl. 426 729-773.
  7. 7. T. J. Laffey, E. Meehan (1999) A characterization of trace zero nonnegative 5x5 matrices. In: Linear Algebra Appl. 302-303 295-302.
  8. 8. N. Radwan (1996) An inverse eigenvalue problem for symmetric and normal matrices. In: Linear Algebra Appl. 248 101-109.
  9. 9. O. Rojo, R. L. Soto (2003) Existence and construction of nonnegative matrices with complex spectrum. In: Linear Algebra Appl. 368 53-69
  10. 10. A. Borobia, J. Moro, R. L. Soto (2004) Negativity compensation in the nonnegative inverse eigenvalue problem. In: Linear Algebra Appl. 393 73-89.
  11. 11. T. J. Laffey, H. Šmigoc (2006) Nonnegative realization of spectra having negative real parts. In: Linear Algebra Appl. 416 148-159.
  12. 12. T. J. Laffey (2005) Perturbing non-real eigenvalues of nonnegative real matrices. In: Electronic Journal of Linear Algebra 12 73-76.
  13. 13. R.L. Soto, M. Salas, C. Manzaneda (2010) Nonnegative realization of complex spectra. In: Electronic Journal of Linear Algebra 20 595-609.
  14. 14. W. Guo (1996) An inverse eigenvalue problem for nonnegative matrices. In: Linear Algebra Appl. 249 67-78.
  15. 15. C. R. Johnson, T. J. Laffey, R. Loewy (1996) The real and the symmetric nonnegative inverse eigenvalue problems are different. In: Proc. AMS 124 3647-3651.
  16. 16. H. R. Suleimanova (1949) Stochastic matrices with real characteristic values. In: Dokl. Akad. Nauk SSSR 66 343-345.
  17. 17. H. Perfect (1953) Methods of constructing certain stochastic matrices. In: Duke Math. J. 20 395-404.
  18. 18. H. Perfect (1955) Methods of constructing certain stochastic matrices II. In: Duke Math. J. 22 305-311.
  19. 19. R. Kellogg (1971) Matrices similar to a positive or essentially positive matrix. In: Linear Algebra Appl. 4 191-204.
  20. 20. F. Salzmann (1972) A note on eigenvalues of nonnegative matrices. In: Linear Algebra Appl. 5.329-338.
  21. 21. A. Borobia (1995) On the Nonnegative Eigenvalue Problem. In: Linear Algebra Appl. 223-224 131-140.
  22. 22. R. L. Soto (2003) Existence and construction of nonnegative matrices with prescribed spectrum. In: Linear Algebra Appl. 369 169-184.
  23. 23. H. Šmigoc (2005) Construction of nonnegative matrices and the inverse eigenvalue problem. In: Linear and Multilinear Algebra 53 88-96.
  24. 24. R. L. Soto, O. Rojo (2006) Applications of a Brauer Theorem in the nonnegative inverse eigenvalue problem. In: Linear Algebra Appl. 416 (2-3) (2006) 844-856.
  25. 25. A. Borobia, J. Moro, R. L. Soto (2008) A unified view on compensation criteria the real nonnegative inverse eigenvalue problem. In: Linear Algebra Appl. 428 2574-2584.
  26. 26. R.L. Soto (2011) A family of realizability criteria for the real and symmetric nonnegative inverse eigenvalue problem. In: Numerical Linear Algebra with Appl. (2011). DOI: 10.1002/nla.835.
  27. 27. P. Egleston, T. Lenker, S. Narayan (2004) The nonnegative inverse eigenvalue problem. In: Linear Algebra Appl. 379 475-490.
  28. 28. C. Marijuán, M. Pisonero, R. L. Soto (2007) A map of sufficient conditions for the real nonnegative inverse eigenvalue problem. In: Linear Algebra Appl. 426 (2007) 690-705.
  29. 29. M. Fiedler (1974) Eigenvalues of nonnegative symmetric matrices. In: Linear Algebra Appl. 9 119-142.
  30. 30. G. Soules (1983) Constructing symmetric nonnegative matrices. In: Linear and Multilinear Algebra 13 241-251.
  31. 31. R. Reams (2002) Constructions of trace zero symmetric stochastic matrices for the inverse eigenvalue problem. In: Electronic Journal of Linear Algebra 9 270-275.
  32. 32. R. Loewy, J. J. McDonald (2004) The symmetric nonnegative inverse eigenvalue problem for 5×5 matrices. In: Linear Algebra Appl. 393 275-298.
  33. 33. R. L. Soto (2006) Realizability criterion for the symmetric nonnegative inverse eigenvalue problem. In: Linear Algebra Appl. 416 (2-3) 783-794.
  34. 34. R. L. Soto, O. Rojo, J. Moro, A. Borobia (2007) Symmetric nonnegative realization of spectra. In: Electronic Journal of Linear Algebra 16 1-18.
  35. 35. T. J. Laffey, H. Šmigoc (2007) Construction of nonnegative symmetric matrices with given spectrum. In: Linear Algebra Appl. 421 97-109.
  36. 36. R.L. Soto, O. Rojo, C.B. Manzaneda (2011) On nonnegative realization of partitioned spectra. In: Electronic Journal of Linear Algebra 22 557-572.
  37. 37. O. Spector (2011) A characterization of trace zero symmetric nonnegative 5×5 matrices. In: Linear Algebra Appl. 434 1000-1017.
  38. 38. T. J. Laffey, E. Meehan (1998) A refinament of an inequality of Johnson, Loewy and London on nonnegative matrices and some applications. In: Electronic Journal of Linear Algebra 3 119-128.
  39. 39. O. Holtz (2004) M-matrices satisfy Newton's inequalities. In: Proceedings of the AMS 133 (3) (2004) 711-717.
  40. 40. C. R. Johnson (1981) Row stochastic matrices similar to doubly stochastic matrices. In: Linear and Multilinear Algebra 10 113-130.
  41. 41. A. Brauer (1952) Limits for the characteristic roots of a matrix. IV: Aplications to stochastic matrices. In: Duke Math. J., 19 75-91.
  42. 42. R. Reams (1996) An inequality for nonnegative matrices and the inverse eigenvalue problem. In: Linear and Multilinear Algebra 41 367-375.
  43. 43. R. A. Horn, C. R. Johnson (1991) Matrix Analysis, Cambridge University Press, Cambridge.
  44. 44. R.L. Soto, A.I. Julio (2011) A note on the symmetric nonnegative inverse eigenvalue problem. In: International Mathematical Forum 6 N 50, 2447-2460.
  45. 45. W. Guo (1997) Eigenvalues of nonnegative matrices. In: Linear Algebra Appl. 266 261-270.
  46. 46. S. Guo, W. Guo (2007) Perturbing non-real eigenvalues of nonnegative real matrices. In: Linear Algebra Appl. 426 199-203.
  47. 47. O. Rojo, R. L. Soto (2009) Guo perturbations for symmetric nonnegative circulant matrices. In: Linear Algebra Appl. 431 594-607.
  48. 48. J. Ccapa, R.L. Soto (2009) On spectra perturbation and elementary divisors of positive matrices. In: Electron. J. of Linear Algebra 18 462-481.
  49. 49. M. Boyle and D. Handelman, The spectra of nonnegative matrices via symbolic dynamics, Ann.of Math. 133: 249-316 (1991).

Written By

Ricardo L. Soto

Submitted: 17 November 2011 Published: 11 July 2012