Arxiv:1909.01291V3 [Math.SP]

Home , List of named matrices, Symmetric matrix

arXiv:1909.01291v3 [math.SP] 24 Jan 2020 ls fdul tcatcmtie;acrigt h Birk the to according of matrices; set stochastic the doubly inter of be matri class Permutation may a chains. matrices Markov stochastic symmetric finite-state Doubly of too). stochastic is [,5 ,7 0 1 1) erfrt ordsppr[]frag a for [9] ne paper of Mourad’s to plethora refer a We by 21]). problems. 11, reflected said 10, as 7, 6, impetus 5, new a ([2, Perfect gained by and has continuation 20] Sule˘ımanova subsequent problem [19, important of an works with the 15]) in origin stochast their doubly having or stochastic, entries, non-negative with arcso an of matrices doubl for problem eigenvalue inverse co the are consider called matrices th to Stochastic (hence is chains 1. Markov note to finite-state of up matrices adds transition row each and negative Sule˘ımanova spectrum. nes ievlepolm o lse fmtie uhas such matrices of classes for problems eigenvalue Inverse e od n phrases. and words Key 2010 Date qaemti ihra nre stermed is entries real with matrix square A NES RBESFRSMERCDUL STOCHASTIC DOUBLY SYMMETRIC FOR PROBLEMS INVERSE bistochastic aur 7 2020. 27, January : ahmtc ujc Classification. Subject Mathematics ulgnrlt.I spoe htwhenever that proved is m a It stochastic is doubly generality. this symmetric full for presented; problem is inverse matrix spectral stochastic doubly symmetric a Abstract. + 1 otecasclsffiin odtosdet Perfect–Mirsky to due Nader conditions to due sufficient refinements classical the to pcrmi precisely is spectrum results. λ n 2 × + ARCSWOESULE WHOSE MATRICES . . . n n eeetset. -element obysohsi arcsi h ovxhl ftesto pe of set the of hull convex the is matrices stochastic doubly + e ucetcniinfrals fra ubr ob h sp the be to numbers real of list a for condition sufficient new A nteltrtr;asohsi arxis matrix stochastic a literature; the in λ n > R ONE EO Y1/2 BY BELOW BOUNDED ARE 1 obysohsi arx itcatcmti,ivrepr inverse matrix, bistochastic matrix, stochastic doubly / nmmra:Mrsa ide (1926–2015) Fiedler Miroslav memoriam: In IHŁGAI N OAZKANIA TOMASZ AND GNACIK MICHAŁ (1 2 hnteeeit ymti,dul tcatcmti wh matrix stochastic doubly symmetric, a exists there then , λ , tal. et 2 λ , . . . , eas rvd oeeape n plctoso our of applications and examples some provide also We 1. n ) 51 piay,ad1A8 51 (secondary). 15A12 15A18, and (primary), 65F18 epitotta hscieini incomparable is criterion this that out point We . Introduction λ 1 2 λ , . . . , ˘ MNV SPECTRA IMANOVA stochastic n r o-oiiera ubr with numbers real non-positive are caewl-otdi h literature, the in well-rooted are ic obystochastic doubly tie hti tl pni its in open still is that atrices emnlg) h i fthis of aim The terminology). e ols n hi modern their and Soules, , hni a l nre non- entries all has it when , otiuint h classical the to contribution e r aaimeapeof example paradigm a are ces o—o emn theorem, Neumann hoff—von rtda rniinmatrices transition as preted dMrk 1] eety the Recently, [16]. Mirsky nd neednl,Pret([14, Perfect independently, , smerco o)matrices not) or (symmetric o vriwcnenn the concerning overview ood tcatcmtie (also matrices stochastic y vnetyitrrtdas interpreted nveniently ucetconditions sufficient w crmof ectrum fistranspose its if be,SDIEP, oblem, rmutation ose 2 MICHAŁGNACIKANDTOMASZKANIA

More explicitly, the symmetric doubly stochastic eigenvalue inverse problem (SDIEP) asks the following:

Let λ1,λ2,...,λn be a list of real numbers. In what circumstances does there exist a symmetric doubly stochastic matrix whose spectrum consists of these numbers? Quite naturally, the inverse problems for other classes of matrices are formulated totally analogously. In this short note, we focus only on providing a new sufficient condition to address SDIEP. Since any matrix A solving the above problem has non-negative entries, by the classical Frobenius–Perron theorem, A must have a non-negative eigenvalue λ1 (that is called the Perron eigenvalue of A) such that λ > λ for any other eigenvalue λ of A (and 1 | | the eigenvector associated to λ1 has non-negative entries, it is also known as the Perron eigenvector). It is to be noted that already for stochastic matrices, λ1 = 1 is the Perron eigenvalue to which corresponds the eigenvector comprising only 1s. Consequently, without loss of generality we shall restrict ourselves to λs from the interval [ 1, 1]. − Since the trace of (any power of) a square matrix with non-negative entries is non- negative, for a list (1,λ2,...,λn) of real numbers to form a spectrum of a solution to k k N SDIEP, it is necessary that 1+ λ2 + ... + λn > 0 for any k . Perfect and Mirsky [16] provided a fairly general sufficient condition for the possibility∈ of solving SDIEP for a list of real numbers 1= λ > λ > ... > 1. Namely, they proved that as long as 1 2 − 1 1 1 (1.1) + λ + ... + λ > 0, n n(n 1) 2 2 1 n − · there exists a symmetric doubly stochastic n n matrix whose spectrum is equal to × 1,λ2,...,λn . This condition was subsequently refined by Soules [18], who gave a finer condition{ depending} on the remainder of n when divided by 2. More specifically, let m be such that n =2m +1 in the case n is odd and n =2m +2 in the case where n is even. If n 1 n m 1 1 (1.2) + − − λ2 + λn 2k+2 > 0, n n(m + 1) (k + 1)k − k X=1 then there exists a symmetric doubly stochastic n n matrix whose spectrum coincides × with 1,λ ,...,λn . { 2 } Soules’ condition was refined further by Nader et al. who arrived at a more complicated condition that depends on the remainder n mod 4 ([11, Theorem 5]) that covers a large class of cases. Following the terminology introduced by Paparella [13], we call a list of real numbers σ = (λ1 = 1,λ2,...,λn) a Sule˘ımanova spectrum, whenever λj 6 0 for j = 2,...,n and 1+ λ2 + ... + λn > 0. Moreover, we say that Sule˘ımanova spectrum is normalised if λ > > λn. 2 ··· Fiedler ([4, Theorem 2.4]) showed that normalised Sule˘ımanova spectra are realisable by symmetric non-negative matrices. Soto and Ccapa ([17, Theorem 3.3]) proved that as INVERSE PROBLEMS FOR SYMMETRIC DOUBLY STOCHASTIC MATRICES 3

long as λj < 0 for j = 2,...,n every normalised Sule˘ımanova spectrum is realisable as a spectrum of a stochastic matrix (not necessarily symmetric). Johnson and Paparella ([6, Problem 6.2. and Theorem 6.3.]) proved, among other things, that every normalised Sule˘ımanova spectrum is realisable as a spectrum of a symmetric doubly stochastic matrix for all Hadamard orders (recall that the order of a Hadamard matrix must be 1, 2 or a multiple of 4). In particular, when n = 2k for some k, the resulting matrix A is trisymmetric ([6, Corollary 6.5]), that is, it satisﬁes any two of the following three conditions: symmetric, persymmetric (AK = KAT ), or centrosymmetric (AK = KA), where K is the exchange matrix (backward identity); satisfying any two of the above three conditions always implies that the remaining third condition holds.

It is to be noted that a normalised Sule˘ımanova spectrum with 1+ λ2 + ...+ λn =0 may fail to be realisable within the class of symmetric doubly stochastic matrices. Indeed, if n is an odd number, e.g., the list (1, 0, 0,..., 0, 1) cannot be a spectrum of any symmetric doubly stochastic matrix ([11, Corollary 1]). − For the sake of brevity, we introduce the following piece of terminology.

Deﬁnition 1. A list σ = (1,λ2,...,λn) is called a δ-Sule˘ımanova spectrum (δ > 0), whenever

(1.3) 1+ λ2 + ... + λn > δ

and λj 6 0 for j =2,...,n.

The main result of this note is to prove that 1/2-Sule˘ımanova spectra may be indeed realised as spectra of symmetric, doubly stochastic matrices. Note that we drop the as- sumption for the spectrum to be normalised.

Theorem. Let σ = (1,λ2,...,λn) be a 1/2-Sule˘ımanova spectrum. Then there exists a symmetric, doubly stochastic matrix whose spectrum is precisely σ. We discern that the idea for the proof has its origins in the theory of Markov chains. Roughly speaking, we consider a simple (symmetric) random walk on the discrete torus Z/nZ (n N) represented by the following graph that we denote by Sn: ∈ n 1 n 2 − −

1 1 ... 2 2 0 1 2 1 1 2 2 11 2 2

The crux of the proof is to use the eigenvectors of the transition probability matrix Pn of the above Markov chain as a generator for a class of new symmetric, doubly stochastic 4 MICHAŁGNACIKANDTOMASZKANIA

matrices. More specifically, we construct the desired matrices as Schur forms of the type T QΛQ , where the column vectors of Q are suitably chosen eigenvectors of Pn and Λ is a diagonal matrix containing the Sule˘ımanova spectrum. This very idea is not specific to Q, however it exploits the properties of Q, that is, symmetry and orthogonality. Such approach offers a room for improvement—we address this further in Remark 2.8. By appealing to this heuristics and backed with some numerical evidence, we raise the following problem. Problem 1. Let δ > 0. Can every δ-Sule˘ımanova spectrum be realised as the spectrum of a symmetric doubly stochastic matrix? The main result of the note asserts that the answer is positive for δ > 1/2.

2. Proof of the main result

We consider the symmetric random walk on the graph Sn described in the Introduction. The corresponding transition probability matrix Pn is then given by 1 1 0 2 0 0 0 0 2 1 1 ··· 2 0 2 0 0 0 0  . . . ···  Pn = ...... 00 0 0 1 0 1   2 2   1 0 0 0 ··· 0 1 0  2 2   ···  Note that Pn is a symmetric circulant matrix. Notation. In this section, we index the eigenvalues to start from 0 rather than 1 so that σ = (λ0,λ1,...,λn 1) rather than σ = (λ1,λ2,...,λn) as done in the previous section. This assures that the− formulae in this section are clearer and more compact.

A rote calculation shows that the eigenvalues λk and the corresponding (real) eigenvectors uk (k 0,...,n 1 ) of Pn are given by ∈{ − } 2πk (j) 2πkj (2.1) λk = cos , u = cos (k, j 0,...,n 1 ); n k n ∈{ − } we prove it, among other things, in Lemma 2.1. Due to symmetry, we have λk = λn k n − and uk = un k, so the eigenvectors uk (k 0,...,n 1 ) fail to span R ; in the next − ∈ { − } lemma we also ﬁnd another set of eigenvectors corresponding to Pn that do indeed form an orthonormal basis of Rn. For an elaborate discussion on eigenvalues of transition matrices we refer the reader to [8, chapter 12]. Lemma 2.1.

i) The eigenvalues λk and the corresponding eigenvectors uk (k 0,...,n 1 ) of ∈ { − } Pn are given by (2.1). ii) Set 2 2πkj π w(j) = sin + (k, j 0,...,n 1 ). k n n 4 ∈{ − } r INVERSE PROBLEMS FOR SYMMETRIC DOUBLY STOCHASTIC MATRICES 5

(j) Then wk = [wk ]06j6n 1 are eigenvectors of Pn corresponding to λ0,...,λn 1. Fur- thermore, these eigenvectors− form an orthonormal basis of Rn. − Proof.

i) Since Pn is a circulant matrix, it is diagonalisable by a Fourier matrix Fn, in par- 1 ticular, Pn = Fn∗diag(λ0,...,λn 1)Fn, where Fn∗ = [ϕ0 ϕ1 ... ϕn 1], − √n − − 2πik 2πikj 2πik(n 1) T ϕk = (1, e n ,...,e n ...,e n ) (k 0,...,n 1 ) ∈{ − } and λk are given by (2.1) (see, e.g., [3, Theorem 3.2.2]). Thus, each ϕk is a (com- plex) eigenvector of Pn corresponding to the eigenvalue λk. Since both the real and (j) 2πkj imaginary parts of ϕk are eigenvectors of Pn too, we conclude that uk = cos n and v(j) = sin 2πkj are the coordinates of the (real) eigenvectors u and v corre- k n k k sponding to λk for k, j 0,...,n 1 . ∈{ − } n 1 n 1 Rn ii) However, it is to be noted that neither of the systems (uk)k=0− , (vk)k=0− spans . As sums of eigenvectors corresponding to the same eigenvalue, if non-zero, are still eigenvectors, let us consider the eigenvector uk +vk corresponding to λk. Evidently, 2πkj π u(j) + v(j) = √2 sin + . k k n 4 We then deﬁne the vectors wk by 2 2πkj π w(j) = sin + (k, j 0,...,n 1 ). k n n 4 ∈{ − } r Conspicuously, wk is a (real) eigenvector of Pn corresponding to λk. n 1 In order to show that (wk)k=0− are pairwise orthogonal unit vectors, we invoke π the identities 2 sin(x) sin(y) = cos(x y) cos(x + y) and cos x + 2 = sin(x). Using these we arrive at − − − n 1 2 − 2πkj π 2πlj π wk,wl = sin + sin + h i n n 4 n 4 j=0 X n n 1 2πj(k l) 1 2πj(k + l) = 1+ cos − + sin , n − n n n j=0 ! j=0 X X where , denotes the standard inner product in Rn. h· ·i n 2πmj n 2πmj Let m be a positive integer. Then j=0 cos n and j=0 sin n are the real and imaginary parts, respectively, of the expression P P n 2iπ(m+ m ) 2πmj e n 1 i n e = m − =1 2iπ n j=0 e 1 X − respectively. Consequently, wk,wl = δk,l, where δk,l denotes the Kronecker delta, n 1 h i Rn and thus (wk)k=0− forms an orthonormal basis of . 6 MICHAŁGNACIKANDTOMASZKANIA

Before we proceed, a piece of notation is required. Let Q be the matrix whose columns are precisely the eigenvectors from the statement of Lemma 2.1, that is, Q = [w0 w1 ... wn 1]. − For a diagonal matrix Λ= diag(λ0 =1,λ1,...,λn 1) with λk [ 1, 1] (k 1,...,n 1 ) we consider the product matrix − ∈ − ∈{ − } n 1 − T T (2.2) P (Λ) := QΛQ = λjwjwj . j=0 X Remark 2.2. We note that the matrix P (Λ) is symmetric and its each row adds up to 1 since P (Λ)w0 = w0. Now we present the general form of the entries of QΛQT , with an arbitrary orthogonal Q 1 T n so that Q = [q0 q1 ... qn 1] with q0 = (1,..., 1) and qk R for each k =1,...,n 1. − √n ∈ − th th (k) Remark 2.3. Let us denote the entry of Q in the k row and l column by ql for (k) n 1 (k) j, k 0,...,n 1 . Note that Q = [q ] − and QΛ = [λlq ]k,l, where λ =1. Hence, ∈{ − } l k,l=0 l 0 n 1 n 1 − 1 − (2.3) (QΛQT ) = λ q(k)q(l) = q(k) + λ q(k)q(l). k,l j j j √n 0 j j j j=0 j=1 X X Now let Q be a matrix whose columns are precisely the eigenvectors of Pn from the statement of Lemma 2.1, namely, Q = [w0 w1 ... wn 1], then − 2 2πkj π (2.4) q(k) = w(k) = sin + (k, j 0,...,n 1 ), j j n n 4 ∈{ − } r (k) (k) 1 in particular, q0 = w0 = √n . To summarise let us record the form of the entries of the matrix P (Λ) in a lemma: n 1 Lemma 2.4. We have P (Λ) = [pkl]k,l−=0, where n 1 1 − 2πkj π 2πlj π pkl = 1+2 λj sin + sin + (k,l 0,...n 1 ). n n 4 n 4 ∈{ − } j=1 ! X Proof. Follows from (2.3) and (2.4). Proposition 2.5. Let P (Λ) be as in (2.2). The matrix P (Λ) is doubly stochastic if and only if n 1 − 2πkj π 2πlj π 1 (2.5) λj sin + sin + > n 4 n 4 −2 j=1 X for all k 0,...n 1 and l k,...,n 1 . ∈{ − } ∈{ − } Proof. The matrix P (Λ) is doubly-stochastic if and only if pkl > 0 for all k 0,...n 1 and l k,...,n 1 and this is indeed equivalent to (2.5) by Lemma 2.4.∈{ − } ∈{ − } INVERSE PROBLEMS FOR SYMMETRIC DOUBLY STOCHASTIC MATRICES 7

Finally, the main result follows directly from the next corollary to Proposition 2.5. Corollary 2.6. Let P (Λ) be as in (2.2).

(1) Suppose that λi 6 0 for all i 1,...,n 1 . Then P (Λ) is doubly stochastic as long as ∈ { − }

n 1 − 1 (2.6) λi > . −2 i=1 X (2) Suppose that λi > 0 for all i 1,...,n 1 . Then P (Λ) is doubly stochastic as long as ∈ { − }

n 1 − 1 (2.7) λ 6 . i 2 i=1 X 2πkj π Proof. First we show (1). Assume that (2.6) holds and set Sj(k) = sin n + 4 for any j 1,...,n 1 and k 0,...,n 1 . Clearly, Sj(k)Sj(l) 6 1 for all k,l 0,...n 1 ∈{ − } ∈{ − } ∈{ − } and since all λj 6 0, we have λjSj(k)Sj(l) > λj. Hence, we arrive at the estimate

n 1 n 1 − − λjSj(k)Sj(l) > λj j=1 j=1 X X 1 > . − 2

To show (2), assume (2.7) and use the fact that Sj(k)Sj(l) > 1. − Remark 2.7. Note that n 1 n 1 − − 1 1 0 6 δ := min λj : λjSj(k)Sj(l) > for all k,l 0,...n 1 6 . min −2 ∈{ − } 2 ( j=0 j=1 ) X X For the upper bound, let us note that Sj(k)=1 if and only if n is divisible by 8 and kj = n 1 + m for some m N 0 . Thus if n is not a multiple of 8 then δ < 1 . 8 ∈ ∪{ } min 2 For the lower bound we found that if n =5 then for σ = (1, 0.004, 0.002, 0.004, 0.51) 4 − − − − so that j=0 λj =0.48 we have a negative element in the third row, third column, namely P (Λ) 0.0005 < 0. These suggest that δ (0.48, 0.5). 22 ≈P − min ∈ We leave the task of finding the exact value of δmin an open problem. Remark 2.8. One may observe that our matrix Q is orthogonal and symmetric. Therefore, one may try to find different orthogonal symmetric matrix so that Problem 1 has a solution for δ (0, 1 ). This could be achieved by finding a constant M (1, 2) that is an upper ∈ 2 ∈ bound for the term nq(k)q(l) in equation (2.3) for each j, k, l 0, 1,...,n 1 . j j ∈{ − } 8 MICHAŁGNACIKANDTOMASZKANIA

A natural example of a matrix that is both orthogonal and symmetric, is a Householder matrix, i.e., vvT H(v)=I 2 − v 2 k k β β β ...β β β 1 α α . . . α α  − − . − −  β α 1 α .. α α =  − .− . −. −   ......   ......   .   β α α .. 1 α α   − − − −   β α α . . . α 1 α   − − − −  T 1 1  so that v = (1 √n, 1,..., 1) , α = √n(√n 1) and β = √n (see [22, Section 4.2]). However, − − given that λj 6 0 for each j = 1,...,n 1, the matrix H(v)ΛH(v), as n becomes large fails to have non-negative entries as −

n 1 − 1 2 2 (H(v)ΛH(v))k,k = + α λj + (1 α) λk. n − j=1 & j=k X 6 3. Examples and applications Examples. In this section we provide some examples of Sule˘ımanova spectra,

σn = (1,λ2,...,λn), 1 for which λ2,...,λn add up to 2 and, thus, yield symmetric doubly stochastic matrices obtained via our construction (2.2),− but do not satisfy any known suﬃcient conditions (e.g., (1.1), (1.2), etc.) to obtain symmetric doubly stochastic matrices. To wit, neither (1.1) nor (1.2) is satisﬁed for

σ5 = (1, 0.02, 0.03, 0.05, 0.4) (odd dimension); • σ = (1, −0.01, −0.02, −0.06, −0.08, 0.33) (even dimension), • 6 − − − − − respectively. Let σn = (1,λ2,...,λn), where λ2 > λ3 ... > λn. The improved Soules’ condition when n is even, [11, Theorem 3, Notation 1, Observation 1], that is,

n+2 1 n n+2 [ 4 ]− 1 1 2 4 λn 4k+4 + λ + − λ + − > 0 n n 2 n n+2 4 k(k + 1) 2 4 k=1 X is not satisﬁed as witnessed by σ = (1, 0.01, 0.01, 0.025, 0.03, 0.035, 0.04, 0.05, 0.08, 0.22) 10 − − − − − − − − − (the square brackets in the above formula denote the integral part of a real number). INVERSE PROBLEMS FOR SYMMETRIC DOUBLY STOCHASTIC MATRICES 9

Let n be odd, new condition 1 ([11, Theorem 4, Notation 1]; we adapt the naming conventions from the said paper), that is,

n+3 1 n+1 n+3 [ 4 ]− 1 n 1 2 4 λn 4k+4 + − λ + − λ + − > 0 n n(n + 1) 2 n+1 n+3 4 k(k + 1) 2 4 k=1 X is not satisﬁed as witnessed by the 1/2-Sule˘ımanova spectrum σ = (1, 0.03, 0.03, 0.04, 0.4). 5 − − − − Next we give examples of spectra that do not satisfy New condition 2 ([11, Theorem 5, Notation 2]) which is given by (3.1, 3.2, 3.3, 3.4) depending on the remainder n mod 4. Let m be an integer greater than 1. If (1) n =4m, then

n+4 1 n n+4 [ 8 ]− 1 1 2 4 8 λn 8k+8 (3.1) + λ + λ + − λ + − > 0 n n 2 n 4 n n+4 8 k(k + 1) 4 8 k=1 X is not satisﬁed by the 1/2-Sule˘ımanova spectrum σ = (1, 0.003, 0.003, 0.004, 0.007, 0.009, 0.02, 0.0209, 0.021, 16 − − − − − − − − 0.024, 0.026, 0.035, 0.042, 0.076, 0.0811, 0.128); − − − − − − − (2) n =4m +2, then

n+6 1 n+2 n+6 [ 8 ]− 1 1 2(n 2) 4 8 λn 8k+8 (3.2) + λ + − λ + − λ + − > 0 n n 2 n(n + 2) 4 n+2 n+6 8 k(k + 1) 4 8 k=1 X is not satisﬁed by the 1/2-Sule˘ımanova spectrum σ = (1, 0.01, 0.01, 0.01, 0.02, 0.02, 0.04, 0.07, 0.1, 0.22); 10 − − − − − − − − − (3) n =4m +3, then

n+5 1 n+1 n+5 [ 8 ]− 1 n 1 2 4 8 λn 8k+8 (3.3) + − λ + λ + − λ + − > 0 n n(n + 1) 2 n +1 4 n+1 n+5 8 k(k + 1) 4 8 k=1 X is not satisﬁed by the 1/2-Sule˘ımanova spectrum σ = (1, 0.001, 0.004, 0.01, 0.01, 0.012, 0.013, 0.05, 0.09, 0.11, 0.2); 11 − − − − − − − − − − (4) n =4m +1, then

n+7 1 n+3 n+7 [ 8 ]− 1 n 1 2(n 1) 4 8 λn 8k+8 (3.4) + − λ + − λ + − λ + − > 0 n n(n + 1) 2 (n + 1)(n + 3) 4 n+3 n+7 8 k(k + 1) 4 8 k=1 X is not satisﬁed by the 1/2-Sule˘ımanova spectrum σ = (1, 0.006, 0.018, 0.02, 0.028, 0.028, 0.053, 0.105, 0.242). 9 − − − − − − − − 10 MICHAŁGNACIKANDTOMASZKANIA

New condition 3 ([11, Conjecture 1, Example 1]) that for n = 26 takes the form 1 1 6 3 1 1 + λ + λ + λ + λ + λ > 0 26 26 2 13 7 4 28 8 4 16 2 26 · is not satisﬁed by the 1/2-Sule˘ımanova spectrum σ = (1, 0.004, 0.005, 0.006, 0.007, 0.01, 0.01, 0.011, 0.011, 0.011, 0.012, 26 − − − − − − − − − − λ2 λ4 λ8

0.012, |0.015{z ,} 0.015, |0.016{z ,} 0.017, 0.019, 0.02, | 0{z.022}, 0.022, 0.025, 0.028, − − − − − − − − − − − λ16

0.028| ,{z0.}032, 0.069, 0.073). − − − − λ26

Applications to random generation of doubly stochastic| {z } matrices. Let n N and α 1 , 1 . Our construction provides a simple way to randomly generate symmetric∈ ∈ − 2 2 doubly stochastic matrices via their spectra. Namely, let X1,...,Xn 1 be independent random variables having probability distributions supported on [0, 1] and− let us consider Sn := X1 + ... + Xn 1. Set − Xi λi = α (i 1,...,n 1 ). Sn ∈{ − }

Then σ = (1,λ1,...,λn 1) is a spectrum of a symmetric doubly stochastic matrix and the corresponding matrix− may be obtained via (2.2). More algorithms to generate doubly stochastic matrices (not necessarily symmetric) can be found in [1]. For an elaborate discussion on spectral properties of random doubly stochastic matrices we refer the reader to [12].

Supplementary material. We supplemnt the material with a Python code organised in a JupyterNotebook available at https://github.com/Nty24/DoublyStochasticMatricesGenerator that generates further examples and counterexamples in the spirit of Section 3.

Acknowledgements. The ﬁrst-named author would like to thank his colleagues, James Burridge, for the idea to investigate the eigenvectors of symmetric random walks, and Ittay Weiss, for great discussions that led to Remark 2.8. The authors would like to thank Pietro Paparella for clarifying the history of the results on realising normalised Sule˘ımanova spactra. M.G.’s visit to Prague in August 2019, during which some part of the project was completed, was supported by funding received from GAČR project 19-07129Y; RVO 67985840, which is acknowledged with thanks. The work was supported by the Iuventus Plus grant IP2015 035174. INVERSE PROBLEMS FOR SYMMETRIC DOUBLY STOCHASTIC MATRICES 11

References

[1] V. Cappellini, H.J. Sommers, W. Bruzda, and K. Życzkowski, Random bistochastic matrices, J. Phys. A: Math. Theor. 42 (2009) 365209, 23 pp. [2] J. Ccapa, R.L. Soto, On spectra perturbation and elementary divisors of positive matrices, Electron. J. Linear Algebra 18 (2009) 462-–481. [3] P. J. Davis, Circulant Matrices, Wiley, New York, 1979. [4] M. Fiedler, Eigenvalues of nonnegative symmetric matrices, Linear Algebra Appl. 9 (1974) 119–142. [5] S.G. Hwang, S.S. Pyo, The inverse eigenvalue problem for symmetric doubly stochastic matrices, Linear Algebra Appl. 379 (2004) 77–83. [6] C. R. Johnson, P. Paparella, Perron spectratopes and the real nonnegative inverse eigenvalue problem, Linear Algebra Appl. 493 (2016) 281–300. [7] Y.-J. Lei, W.-R. Xu, Y. Lu, Y.-R. Niu, and X.-M. Gu, On the symmetric doubly stochastic inverse eigenvalue problem, Linear Algebra Appl. 445 (2014), 181–205. [8] D. A. Levin, Y. Peres, and E. L. Wilmer, Markov chains and mixing times, American Mathematical Society, Providence, RI, 2009 [9] B. Mourad, An inverse problem for symmetric doubly stochastic matrices, Inverse Problems 19 (2003), 821–831 [10] B. Mourad, H. Abbas, A. Mourad, A. Ghaddar, and I. Kaddoura, An algorithm for constructing doubly stochastic matrices for the inverse eigenvalue problem, Linear Algebra Appl, 439 (2013), 1382–1400. [11] R. Nader, B. Mourad, A. Bretto, and H. Abbas, A note on the real inverse spectral problem for doubly stochastic matrices, Linear Algebra Appl, 569 (2019), 206–240. [12] H.H Nguyen, Random doubly stochastic matrices: the circular law, Ann. Probab. 42, 3 (2014), 1161– 1196. [13] P. Paparella, Realizing Sule˘ımanova spectra via permutative matrices, Electron. J. Linear Algebra, 31 (2016), 306–312. [14] H. Perfect, On positive stochastic matrices with real characteristic roots, Proc. Cambridge Philos. Soc., 48 (1952), 271–276. [15] H. Perfect, Methods of constructing certain stochastic matrices I, Duke Math. J., 20 (1953), 395–404. [16] H. Perfect, L. Mirsky, Spectral properties of doubly-stochastic matrices, Monatsh. Math. 69 (1965) 35–57. [17] R.L. Soto and J. Ccapa, Nonnegative matrices with prescribed elementary divisors. Electron. J. Linear Algebra, 17 (2008), 287—303. [18] G.W. Soules, Constructing symmetric nonnegative matrices, Linear Multilinear Algebra 13 (1983), 241–251. [19] H. R. Sule˘ımanova, Stochastic matrices with real characteristic numbers. (Russian), Doklady Akad. Nauk SSSR (N.S.) 66, (1949). 343–345 [20] H. R. Sule˘ımanova, The question of a necessary and suﬃcient condition for the existence of a stochastic matrix with prescribed characteristic numbers (in Russian) Trudy Vsesojuz. Zaočn. Ènerget. Inst. Vyp. 28 (1965), 33–49. [21] W.-R. Xu, Y.-J. Lei, X.-M. Gu, Y. Lu, and Y.-R. Niu, Comment on "A note on the inverse eigenvalue problem for symmetric doubly stochastic matrices", Linear Algebra Appl. 439 (2013), 2256–2262. [22] S. Zhu, T. Gu, X. Liu, Solving inverse eigenvalue problems via Householder and rank-one matrices, Linear Algebra Appl, 430 (2009), 318–334. 12 MICHAŁGNACIKANDTOMASZKANIA

School of Mathematics and Physics, Lion Gate Building, Lion Terrace, University of Portsmouth, Portsmouth, PO1 3HF, United Kingdom E-mail address: [email protected] Institute of Mathematics, Czech Academy of Sciences, Žitná 25, 115 67 Prague 1, Czech Republic and Institute of Mathematics, Jagiellonian University, Łojasiewicza 6, 30-348 Kraków, Poland E-mail address: [email protected], [email protected]