Can one hear a matrix? Recovering a real symmetric matrix from its spectral data

Tomasz Maci ˛azek˙

E-mail: [email protected] School of Mathematics, University of Bristol, Fry Building, Woodland Road, Bristol BS8 1UG, UK

Uzy Smilansky Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 7610001, Israel

Abstract. The spectrum of a real and symmetric N × N matrix determines the matrix up to unitary equivalence. More spectral data is needed together with some sign indicators to remove the unitary ambiguities. In the first part of this work we specify the spectral and sign information required for a unique reconstruction of general matrices. More specifically, the spectral information consists of the spectra of the N nested main minors of the original matrix of the sizes 1, 2,...,N. However, due to the complicated nature of the required sign data, improvements are needed in order to make the reconstruction procedure feasible. With this in mind, the second part is restricted to banded matrices where the amount of spectral data exceeds the number of the unknown matrix entries. It is shown that one can take advantage of this redundancy to guarantee unique reconstruction of generic matrices, in other words, this subset of matrices is open, dense and of full measure in the set of real, symmetric and banded matrices. It is shown that one can optimize the ratio between redundancy and genericity by using the freedom of choice of the spectral information input. We demonstrate our constructions in detail for pentadiagonal matrices.

arXiv:2106.01819v3 [math-ph] 17 Aug 2021 1. Introduction

The answer to the question posed in the title is definitely: No! The spectrum determines the matrix up to unitary equivalence. The question is then, what additional spectral information is required for this purpose. This question accompanies Mathematical Physics already for a long time, and it is usually referred to as the spectral inversion or spectral reconstruction problems [1, 2]. The related reconstruction methods have a wide scope of applications in different areas of Physics, Materials Science and Engineering. They appear e.g., in studies of mechanical systems near equilibrium, where one tries to construct a quadratic Hamiltonian model [3], or in atomic, molecular and nuclear physics where spectra are measured, with the hope that they will provide information on the underlying interactions. In most cases the Hamiltonian operator is expressed as an ODE or PDE acting on an appropriate function Can one hear a matrix? 2 spaces. A prime example is the Sturm Liouville operator which appears in various guises [16]. In other cases, the systems are described in terms of discrete models, so that finite dimensional matrices are the operators under study. In the present paper we shall consider this case. Recent applications of inverse methods in experiments are numerous. We mention briefly a partial list. A model mechanical problem concerns a spring-mass system where the spring constants and masses can be deduced from the knowledge of system’s eigenfrequencies [1,4]. More complex problems include studying vibrating beams [5] or the composition of strings [6, 7]. Another paper [8] studies neutrino oscillations by reconstructing the so-called mixing matrix in matter from the spectra of its main (N − 1) × (N − 1)-minors. There, the relevant identity is the so-called eigenvector-eigenvalue identity [9]. Other types of eigenvector- eigenvalue identities have been used for parameter estimation in spin chain Hamiltonians [11] and Markov Models [12]. In the present paper the spectral information about the matrix comes from the spectra of its N nested main minors. Such a spectral information determines the matrix up to a finite (a prori large) number of possibilities. The reconstruction can be then made unique by providing additional discrete data encoded in the signs of certain expressions involving entries of the desired matrix. We devote much attention to the study of this sign problem. This may seem like a small subtlety, but the difficulty related to correctly choosing different types of signs in the study of the inverse matrix problems is often identified as one of the main factors that makes the practical use of the theoretical methods hard. G. M. L. Gladwell in his book [1] describes the subtle role played by different types of signs in the following way: “Signs, positive and negative, lie at the heart of any deep discussion of inverse problems.”. Arguably, the reconstruction methods have been most successful for symmetric tridiagonal matrices often referred to as Jacobi matrices [16]. For an N dimensional Jacobi matrix, the number of unknown entries is 2N − 1 and the spectral information necessary for the reconstruction is the spectrum with N eigenvalues, to which one adds the spectrum of the N − 1 dimensional minor obtained by removing the last column and row from the original matrix (also called a main minor of order N − 1). This data provides the diagonal elements (with their signs), but only the absolute value of the off diagonal entries can be recovered [13–15]. Thus, the related sign problem is particularly simple – it suffices to know the signs of the off-diagonal entries of the desired Jacobi matrix. When more general matrices which are less sparse are considered, the number of unknown entries is larger, and naturally, the number of input spectral data has to increase. This is accompanied by the sign problem becoming more acute. The sign-indicators needed in the general case are not necessarily the signs of the individual entries. However, they are absolutely needed to resolve the ambiguity due to spectral invariance of the main minors under similarity with diagonal matrices with entries arbitrarily chosen from {−1, +1}. Various methods were developed for recovering banded or full matrices [9,17–20,22,23,25]. Most of them provide only an exemplary matrix satisfying the spectral data and stop short of resolving the sign-related or other ambiguities. One of the main approaches in reconstructing banded symmetric matrices is to fix the spectra of a few largest main minors (specifically d + 1 main minors of orders N, N − 1, ..., N − d for a (2d + 1)-banded matrix) and using Can one hear a matrix? 3

the Lanczos algorithm, bring the reconstructed matrix to a (2d + 1)-banded form [17–20]. Importantly, the Lanczos algorithm changes the spectra of all remaining main minors and thus this method differs from the approach presented in this work. Similar problems have been considered for matrices with complex entries. A reconstruction algorithm using spectra of all 2N − 1 main minors has been developed [21]. The algorithm produces an example of a complex matrix whose main minors have the given spectra and, under some regularity assumptions, the result is unique up to similarity transformations by nonsingular diagonal matrices and the transposition operation [22]. Even more general settings concerning matrices with entries from an algebraically closed field have been considered. In particular, it has been shown that there exists a finite number of square N × N matrices with a prescribed spectrum and prescribed off-diagonal entries [23]. A different approach to the problem is presented in [24] where any Hermitian matrix can be reconstructed from its spectrum and the spectra of a suitable number of its perturbations. Recently, there has been some revived interest in the so-called eigenvector-eigenvalue identity [9] which allows one to find the amplitudes of eigenvector’s entries of a hermitian matrix using the spectra of its (N − 1) × (N − 1)-minors and the spectrum of the full matrix itself. As the authors of [9] point out, the eigenvector- eigenvalue identity appears in different guises in many places in the literature (see e.g. [10] – the authors of [9] also provide a thorough review of the relevant literature). In the present paper we derive (Appendix B) an identity which bears some similarity to the eigenvector- eigenvalue identity. This identity was proved previously [17–19] using a different method. The present paper consists of three sections. The first deals with the inverse problem of full, real and symmetric matrices of dimension N, where the number of unknown entries is 1 2 N(N + 1). The spectral data to be used is the union of the spectra of the first N nested main 1 minors of the matrix of the sizes 1, 2,...,N, and 2 N(N + 1) sign indicators needed for the complete reconstruction. The precise definition of the sign indicators will be given below. The actual construction is inductive: given a matrix A and a minor A(n) of dimension n, its next neighbor A(n+1) is obtained by computing the (n + 1)’th column from the given spectra and sign indicators. The uniqueness of the resulting solution is proved. Thus, the matrix unfolds like a telescope, hence its name - the telescopic construction. The fly in the ointment is that the computation of the sign indicators is rather cumbersome. The second section uses the telescopic method restricted to banded matrices with band width D = 2d+1 much smaller than N. The spectral input exceeds the number of unknowns, but this redundancy enables circumventing the need to compute the sign-indicators, and the only sign information required consists of the signs of the off-diagonal matrix elements in the diagonal which is d -steps away from the main diagonal, just as it was the case for the Jacobi matrices. The proposed method is proved to provide a unique solution only for generic D-diagonal matrices. Namely, this subset is open, dense and of full measure in the set of (N− 1 )(d+1) real, symmetric and banded matrices specified by entries in R 2 . An explicit criterion which distinguishes non-generic cases is proposed. Finally, in the last section it is shown that the large redundancy which exists in the telescopic approach can be reduced appreciably by studying a different construction: One considers only the d + 1 dimensional principal minors M (n) which are distinguished by the Can one hear a matrix? 4

position of their upper corner of M (n) along the diagonal. Their spectra and appropriate sign indicators are used to compute successive minors in a recursive process. This method will be referred to as the sliding-minor construction. Also in the sliding-minor construction we are able to introduce certain redundancy by considering (d + 2)-dimensional sliding minors. Then, we show that the required sign data typically reduces to the signs of the off-diagonal matrix elements in the first row. The application of the last two methods for banded matrices, requires a special treatment of the upper d × d main minor, as will be explained in detail for each case. To allow a smooth flow of the main ideas, some of the more technical proofs of lemmas and theorems which are stated in the text are deferred to the appendices. We believe that turning the focus to generic rather than to unconditionally applicable methods increases the domain of practical application and the scope of the study of inverse problems.

1.1. Notations and preliminaries

Before turning to the main body of the paper, the present subsection introduces notations and conventions that will be used throughout. Let A be a symmetric real N × N matrix. Denote by A(n) its n × n upper main diagonal (1) (N) minor, so that A = A1,1 and A = A. The upper suffix (n) will be used throughout for the dimension of the subspace under discussion. (n) (n) n (n) The Dirac notation for vectors is used. That is, given a list x = {xj }j=1, then x (n) denotes the column vector in dimension n with entries xj . Its transpose (row) will be denoted (n) (n) (n) by x so that the scalar product is x y . It is convenient also to introduce the unit (n) vectors e (j) , 1 ≤ j ≤ n whose entries all vanish but for the j’th which is unity. Thus, (n) (n) (n) (n) (n) (n) (n) clearly e (j) x = xj , and Ai,j := e (i) A e (j) . For every minor A(n), 1 ≤ n ≤ (N − 1) define the upper off-diagonal part of the next column  (n+1)  A1,n+1 n X (n+1) a(n) :=  .  = A e(n+1)(j) .  .  j,n+1 (n+1) j=1 An,n+1 (n) (n) (n) n (n) The spectra of A for 1 ≤ n ≤ N will be denoted by σ = {λj }j=1, with λk ≥ (n) λj ∀ k ≥ j, and the corresponding normalized eigenvectors (determined up to an overall  (n) n sign) are v (j) j=1. In general, the choice of eigenvectors is ambiguous. In the following, we will often remove this ambiguity by fixing the choice of eigenvectors of A(n). If the spectrum is non-degenerate, this is done, for instance by demanding the last entry of (n) (n) each v (j) to be positive. For a fixed choice of eigenvectors, the projections of a on (n) (n) (n) the eigenvectors of A will be denoted by sj ξj where ( n n (n) (n) (n)
o +1 for x ≥ 0 sj := Sign a v (j) , with Sign(x) = . (1) j=1 −1 for x < 0

(n) The sj will play the role of the sign-indicators in what follows. However, the arbitrary choice (n) of the overall signs of the vectors v (j) will not have any effect on the results, as long as Can one hear a matrix? 5

the chosen signs are not altered throughout the proof.

2. Inversion by the telescopic construction

Given the spectra σ(n) for 1 ≤ n ≤ N, one can use very simple arguments to derive the following information on the matrix A, without relaying on any sign information. They are known in different guises, see e.g. references [17–20]. We bring them here for completeness and as a background necessary for the ensuing developments. Theorem 1. The spectra σ(n) for 1 ≤ n ≤ N suffice to determine the diagonal elements of A and the norms of the vectors |a(n)i.

(n+1) (n) (n+1) Proof. The minor A is different from A by the added diagonal entry An+1,n+1 to be (n) denoted by h, and the off-diagonal column a . Clearly (n+1) (n) (n+1) tr A − tr A = An+1,n+1 = h , (2) Hence h is deduced directly from the spectral data. Also, n 2 X  (n+1) (n) 2 1 (n+1) 2
(n) 2
2
(n) A = a = tr (A ) − tr (A ) − h =: R . (3) k,n+1 2 k=1 2 Thus h, and the value of a(n) can be computed for all n. Remark 1. This result can be considered as the inverse of the Gerschgorin theorem [28], which for symmetric matrices states that given a symmetric matrix A, its spectrum lies in the PN union of real intervals centered at the points An,n and are of lengths 2 k=1 |An,k|. Theorem 1 states that given the spectra of the successive minors, the diagonal elements of the matrix are determined precisely, and the vectors of off-diagonal elements a(a) are restricted to a sphere of a radius determined by Equation (3). Remark 2. If A is a Jacobi (symmetric, tridiagonal) matrix, only the nth component of a(n) is different from zero, and therefore the off-diagonal entries are determined up to a sign. In other words, one has to provide their signs so that the Jacobi matrix could be heard. The following identities are similar in spirit to the former ones, and can be proven by similar ways. They will be used in the sequel.

n n (n) (n) (n) X (n) (n) (n) 2 X (n) (n) 2 a A a = λk |ha |v (k)i| = λk |ξk | = (4) k=1 k=1 1 = tr[(A(n+1))3] − tr[(A(n))3] − 3h|a(n)|2 − h3
. 3 and, n n (n+1) (n) (n)
−1 (n) X 1 (n) (n) 2 X 1 (n) 2 det A a A a = a v (k) = |ξ | = h − , (5) (n) (n) k det A(n) k=1 λk k=1 λk Can one hear a matrix? 6

which is valid if 0 is not in the spectrum of A(n). (If 0 ∈ σ(n) one can add to A a constant multiple of the identity which renders all the spectra non-vanishing, without affecting any of the proofs or the results derived in the sequel.) The identities (3,4,5) are quadratic in the n components of the vector |a(n)i. When A is assumed to be a pentadiagonal (D = 5), only the last two components of |a(n)i do not vanish, and the quadratic forms describe three concentric curves - a circle and two conic sections in R2. Their mutual intersections are candidates for the solution for the inversion problem. This is discussed in detail in section (3.2.2).

2.1. A general algorithm

Unlike the case of tridiagonal matrices where the signature indicators required are just the signs of the off diagonal entries (see Remark 2), for the general case one requires the sign n n (n)o indicators sj defined in (1) . j=1 In the following, it will often be convenient to assume certain regularity properties of the matrix A. Definition 2 (Regular matrix). A real N × N symmetric matrix is called regular if and only if the following conditions are satisfied. (i) The multiplicity of all eigenvalues of each A(n) is one. (ii) The successive spectra do not share a common value, σ(n) ∩ σ(n+1) = ∅ for n = 1,...,N − 1. Remark 3. The following theorem was kindly provided by Percy Deift [29], and it is quoted without proof. Theorem 3. (Deift) The subset of regular matrices is open, dense and of full measure in the set of real and symmetric matrices. Let A be an N × N real symmetric and regular matrix. We prove  (n) N  (n) N Theorem 4. Given the spectra σ n=1 and the sign indicators s n=1 as defined in (1.1) . Then, these data suffice to construct the original matrix A uniquely. ( Irregular matrices will be discussed in a subsequent subsection ). Consider two successive minors A(n) and A(n+1). Cauchy’s spectral interlacing theorem guarantees that (n+1) (n) (n+1) (n+1) (n) (n+1) λ1 ≤ λ1 ≤ λ2 ≤ · · · ≤ λn ≤ λn ≤ λn+1 . (6) (n) (n+1) Lemma 5. The vector a and the diagonal element h := An+1,n+1 , as well as the (n+1)  (n+1) n+1 (normalized) eigenvectors of the minor A , v (j) j=1 , are uniquely determined by the following data. (n) (n)  (n) n (i) The full spectral data of A , i.e. its spectrum σ and eigenvectors v (j) j=1. (ii) The spectrum σ(n+1) of A(n+1), (n) (iii) The sign indicators sj , j = 1, . . . , n. Can one hear a matrix? 7

Proof. The value of h is determined by (2). The next steps are based on the relation between the two spectra σ(n) and σ(n+1) which is used in the proof of (6). Define an auxiliary matrix ˜(n+1) (n) ˜(n+1) A with A being its main n × n minor, An+1,n+1 = h, and whose off diagonal entries of the (n + 1)’th column and row all vanish. Then, the two matrices of dimension (n + 1), A(n+1) and A˜(n+1) differ by a matrix of rank 2. (n+1) (n+1) (n+1) (n) (n) (n+1) A = A˜ + e (n + 1) a , 0 + a , 0 e (n + 1) (7) The spectrum of A˜(n+1) consists of σ(n) ∪ {h} . Correspondingly, the eigenvectors of A˜(n+1) are obtained by increasing the dimension of each of the eigenvectors of A(n) by adding a 0 in (n+1) the (n + 1)’th position and adding the eigenvector e (n + 1) . These eigenvectors will  (n+1) n+1 be denoted by v˜ (j) j=1 , where (n+1) (n) (n+1) (n+1) v˜ (j) = v (j), 0 for j = 1, . . . , n and v˜ (n + 1) = e (n + 1) . To express the spectrum and eigenvalues of A(n+1), expand any of its eigenvectors in the  (n+1) n+1 basis v˜ (j) j=1 , so that n+1 (n+1) X (n+1) (n+1) (n+1) (n+1) (n+1) v (k) = bk,r v˜ (r) and A v (k) = λk v (k) . (8) r=1 (n+1) n+1 Replacing A by its form (7) one easily finds that the expansion coefficients {bk,r}r=1 of the kth eigenvector satisfy

(n+1) (n) λk bk,r = (1 − δr,n+1)bk,rλr + δr,n+1h bk,n+1 (9) n X (n) (n+1) (n+1) (n) + δr,n+1 bk,j a , 0 v˜ (j) + (1 − δr,n+1) v˜ (r) a , 0 bk,n+1 . j=1 From these relations and the definitions above we conclude that: (i) for r = n + 1 n (n+1) X (n) (n) (λk − h)bk,n+1 = bk,j a v (j) , (10) j=1 (ii) for r ≤ n (n+1) (n) (n) (n) (λk − λr )bk,r = bk,n+1 v (r) a . (11) (n+1) (n) Since by the conditions of the theorem, λk 6= λr for all the pairs 1 ≤ r ≤ n, 1 ≤ k ≤ n + 1, one can write (n) (n) v (r) a bk,r = (n+1) (n) bk,n+1 (12) λk − λr

n (n) (n) 2 (n+1) X a v (r) λk − h = (n+1) (n) , ∀ 1 ≤ k ≤ n + 1 . (13) r=1 λk − λr Can one hear a matrix? 8

Moreover, all minors are assumed to be of full rank, and both σ(n+1) and σ(n) are given. Therefore, one can consider (13) as a set of linear equations from which the unknown n parameters (n) 2 (n) (n) 2 (ξr ) := a v (r) , r = 1, . . . , n could be computed and none vanishes as will be shown in (15) below. However, the number of equations exceeds by 1 the number of unknowns, and therefore there exists a solution only if the equations are consistent. To prove consistency, we start by obtaining an explicit solution of the first n equations in (13). Introducing the Cauchy matrix 1 Ck,r = (n+1) (n) , for, 1 ≤ k, r ≤ n , (14) λk − λr (n) 2 n the unknown parameters {(ξr ) }r=1 then read (see Appendix B): Qn+1  (n) (n+1) n λr − λ X (n+1) k=1 k (ξ(n))2 = C−1
(λ − h) = − . (15) r r,k k Qn  (n) (n) k=1 k=1,k6=r λr − λk Formula (15) can also be found in earlier works [17–19] and the recursive relation for the eigenvectors (12) has been also used in [18]. Similarly, one can combine Equation (12) with Equation (15) to obtain (after requiring the eigenvectors of A(n+1) to be normalized to one) 2 the following expression for the coefficients |bk,n+1| Qn  (n+1) (n) λ − λr 2 r=1 k |bk,n+1| = . (16) Qn  (n+1) (n+1) r=1,r6=k λk − λr Equation (16) is an instance of the eigenvector-eigenvalue identity from [9] which allows D E (n+1) (n+1) one to compute the amplitudes er v˜ (k) directly from the spectra of the n main minors of A(n+1) of the size n × n. Here, we are using only one n × n minor (namely, A(n)), thus we recover only the last squared entry of each eigenvector. By combining Equation (16) (n) with Equation (12) we can express the coefficients bk,r in terms of the spectra of A and A(n+1), but the resulting formulae will be different from the eigenvector-eigenvalue identity [9] as we use different spectral data and the entries refer to the basis of the eigenvectors of A(n) only. The consistency of the linear n+1 equations (10) and (11) could be proved by substituting (n) 2 the (ξr ) computed above in the last equation in (13), and showing that the resulting equation n n (n+1) X X 1 −1
(n+1) (λn+1 − h) = (n+1) (n) C r,k (λk − h), (17) r=1 k=1 λn+1 − λr −1 is satisfied identically. An explicit expression for (C )r,k is given in [33] : p(λ(n))q(λ(n+1)) C−1
= r k , r,k (n+1) (n) 0 (n+1) 0 (n) (λr − λk )p (λk )q (λr ) Can one hear a matrix? 9

n n Y (n+1) Y (n) where : p(x) = (x − λi ); q(x) = (x − λi ) . (18) i=1 i=1 The identities which are proved in Appendix A show that equation (17) is satisfied identically Pn+1 (n+1) Pn (n) provided that h = k=1 λk − r=1 λr as is the case in the present context. (n) 2 To recapitulate, the solution of equations (13) exists and it provides the (ξr ) in terms of the spectra in a unique way as shown in Equation (15). This information together with (n) (n) (n) (n) (n) the sign indicators determine ha |v (r)i = sr ξr , and the vector a is obtained by an orthogonal change of basis. Moreover, the possible 2n combinations of the sign indicators encapsulate all possible choices of the overall signs of the eigenvectors of A(n). Thus, an a priori choice of the overall phases of the eigenvectors does not result with any loss of generality or the loss of solutions.

Remark 4. A complementary result to the above was recently reported in [30]. Using another approach it is proved that given a real symmetric matrix A˜(n) whose spectrum σ˜(n) is simple (no multiplicities) and a list µ(n+1) of (n + 1) arbitrary numbers which interlace with σ˜(n), one can construct a column vector |a˜(n)i and a real h˜ which, when used to complete A˜(n) to A˜(n+1) yields A˜(n+1) with the spectrum µ(n+1). Moreover, |a˜(n)i can be constructed uniquely once an orthant (whose choice provides the sign indicators) is prescribed in the n dimensional space spanned by the eigenvectors of A˜(n). It is also proved that the assignment of |a˜(n)i and h˜ to µ(n+1) is a homeomorphism onto its image and a differomorphism. (n) The knowledge of the a v (r) with their signs enables the use of (12) to obtain (n+1) the coefficients of the eigenvectors of A up to the constant factor bk,n+1 which can be determined by requiring a proper normalization.

Proof. (of Theorem 4) The proof proceeds by a recursive construction: Start with a known (1) (2) (1) (1) A and with the known spectrum of A . Clearly, λ1 = A1,1, and v = 1. Furthermore,

(2) (2) (1) (1) (1) A1,2 h = λ1 + λ2 − λ1 = A2,2, a = A1,2, a v = A1,2, s1 = . |A1,2| Hence, Equation (15) gives

(1) 2 2  (1) (2)  (2) (1)  (2)  (2)  (ξ1 ) = (A1,2) = λ1 − λ1 λ2 − λ1 = A1,1 − λ1 λ2 − A1,1 .

(2) (1) Therefore, we get two possible solutions for A corresponding to choosing s1 = +1 or (1) s1 = −1, which written explicitly read r  (1)  (2)  (2)   A1,1 s1 A1,1 − λ1 λ2 − A1,1 (2) A =  r  . (19)  (1)  (2)  (2)  (2) (2)  s1 A1,1 − λ1 λ2 − A1,1 λ1 + λ2 − A1,1 The normalized eigenvectors can be computed from Equation (12). In the standard basis, their forms read v     u (2) (2) 1 α (2) 1 −1 (1)uλ2 − A1,1 v (1) = √ , v (2) = √ , α := −s t . (20) 2 1 2 α 1 (2) 1 + α 1 + α A1,1 − λ1 Can one hear a matrix? 10

Thus the theorem is valid for n = 1, 2. Note also that one can show that (13) is satisfied - which is not clear at first sight but can be easily checked. The lemma provides the final stage of the inductive proof.

The above inductive procedure is optimal for generic matrices whose all entries are typically nonzero. However, the requirement of providing signs of projections of the unknown column to the eigenvectors of A(n) seems awkward for practical use. However, when one wants to apply this result to banded matrices, it turns out that there is large redundancy which can be used to our advantage. In particular, as we argue in the following sections, the inductive reconstruction procedure for banded matrices can be improved so that the number of required sign data can be tremendously reduced.

2.1.1. Explicit reconstruction of A(3). It is instructive to work out the above inductive procedure explicitly for n = 3. The results will prove useful in the following sections. Consider the problem of reconstructing a regular 3 × 3 real symmetric matrix, A(3), given its (3) (3) (3) (2) spectrum λ1 < λ2 < λ3 and its 2 × 2 top-left minor A . It has a spectral decomposition (2) (2) (2) (2) (2) (2) (2) (2) (2) A = λ1 v (1) v (1) + λ2 v (2) v (2) , λ1 < λ2 , where r ! 1  2 λ(2) = A(2) + A(2) ± (A(2) − A(2))2 + 4 A(2) . 1,2 2 1,1 2,2 1,1 2,2 1,2 and v     u (2) (2) (2) 1 α (2) 1 −1 uλ2 − A1,1 v (1) = √ , v (2) = √ , α := −st (21) 2 1 2 α (2) (2) 1 + α 1 + α A1,1 − λ1

 (2) (3) with s = Sign A1,2 . The unknown entries of matrix A read

(3) (2) (2) (2)
(2) (2) (2)
(2) Ai,3 = s1 |ξ1 | v (1) i + s2 |ξ2 | v (2) i , si ∈ {+1, −1}, i ∈ {1, 2}, (22) (2) 2 where {ξr }r=1 are given by Equation (15). Thus, there are four possible solutions for the (3) last column of A corresponding to different choices of signs si. Next, let us show how to determine si directly from the sign of a certain expression involving only matrix elements of A(3). Using Equation (15) and Equation (21), we get that √ √ 2 (3) (2) (2) (2) 2 (3) (2) (2) (2) 1 + α A1,3 = s1|ξ1 |α − s2 |ξ2 |, 1 + α A2,3 = s1 |ξ1 | + αs2|ξ2 |. (2) (2) This allows us to express s1 and s2 as (3) (3) (3) (3) (2) αA1,3 + A2,3 (2) αA2,3 − A1,3 s1 = √ , s2 = √ . 2 (2) 2 (2) 1 + α |ξ1 | 1 + α |ξ2 | Because the denominators of the above expressions are positive numbers, we finally get (2)  (3) (3) (2)  (3) (3) s1 = Sign −sA1,3|α| + A2,3 , s2 = −Sign −sA2,3|α| + A1,3 . (23) The main advantage of these expressions is that they provide expressions for the signatures (2) (2) s1 , s2 in terms of data directly available from the matrix. Can one hear a matrix? 11

2.2. Non-regular matrices

In this subsection we show that non-regular matrices that have degeneracy in σ(n) or where the overlaps σ(n) ∩ σ(n+1) are nonempty can be effectively treated using the methods developed above for regular matrices. We consider the situation where there is a single block of (n) (n) (n) degeneracy in σ , i.e. for some 1 ≤ l ≤ n and m ≥ 0 we have λl = λl+1 = ... = (n) (n) λl+m := λ and other eigenvalues of A are non-degenerate. Let us denote the respective degeneracy indices by D(n)(l, m) := {l, l + 1, . . . , l + m}. The interlacing property (6) dictates the following possibilities for the respective degeneracy indices D(n+1)(l, m) (see Fig. 1): (n+1) DI (l, m) := {l, l + 1, . . . , l + m + 1}, (n+1) DII (l, m) := {l, l + 1, . . . , l + m}, (n+1) DIII (l, m) := {l + 1, l + 2, . . . , l + m + 1}, (n+1) (n+1) DIV (l, m) := {l + 1, l + 2, . . . , l + m} for m ≥ 1, DIV (l, 0) := ∅.

Figure 1. The four possibilities for the position of the degeneracy indices D(n+1)(l, m) relative to D(n)(l, m).

(n+1) (n+1) In what follows we assume that the eigenvalues λk with k∈ / D (l, m) are non- degenerate. Furthermore, the degenerate eigenspaces of A(n) and A(n+1) will be denoted by V (n)(λ) and V (n+1)(λ) respectively, and their respective orthogonal complements in Rn and Rn+1 by V (n)(λ)⊥ and V (n+1)(λ)⊥. We will also embed Rn into Rn+1 by appending a zero at the end of every vector. With a slight abuse of notation we will also denote by V (n)(λ) and V (n)(λ)⊥ the images of these spaces under the above embedding. Can one hear a matrix? 12

Theorem 6 below shows how to reduce each of the above four degenerate cases to a regular problem, where Theorem 4 can be applied. Interestingly, each case requires a different procedure. Theorem 6. Assume that the spectra of A(n) and A(n+1) as well as the eigenvectors  (n) n (n) v (r) r=1 of A are given and fixed. Assume that the spectra have single degeneracy blocks with the eigenvalue λ given by degeneracy indices D(n)(l, m) and D(n+1)(l, m) = (n+1) (n+1) Di (l, m) for some i ∈ {I,II,III,IV } respectively. The matrix A is reconstructed from the above spectral data as follows.

(n) (n+1) I. The vector a as well as the non-degenerate eigenvectors of A belong to the space V (n)(λ)⊥. They are constructed by applying Theorem 4 to the truncated regular

matrices A˜(n−m) := A(n+1) and A˜(n−m−1) := A(n) . The degenerate V (n)(λ)⊥ V (n)(λ)⊥ eigenvectors of A(n+1) n (n+1) (n+1) X (n) (n+1) v˜ (k) = bn+1,k e + bk,r v (r), 0 , k ∈ DI (l, m) r=1 have coefficients (n) (n) a v (r) (n) (n+1) bk,r = (n) bk,n+1, r∈ / D (l, m), k ∈ DI (l, m). λ − λr (n+1) (n) The remaining coefficients bk,n+1 and br,k with k ∈ DI (l, m) and r∈ / D (l, m) can be arbitrary so that the resulting eigenvectors form an orthonormal set. (n) (n+1) II. The vector a as well as the non-degenerate eigenvectors of A belong to the space V (n)(λ)⊥. They are constructed by applying Theorem 4 to the truncated regular

matrices A˜(n−m) := A(n+1) and A˜(n−m−1) := A(n) . The degenerate V (n)(λ)⊥ V (n)(λ)⊥ eigenvectors of A(n+1) form an orthonormal basis of V (n)(λ). III. The same as case II above. IV. The space V (n+1)(λ) is a subspace of V (n)(λ) of codimension one. The degenerate (n+1) (n+1) (n+1) eigenvectors of A , v˜ (k) , k ∈ D (l, m) can be freely chosen as any orthonormal subset of V (n)(λ). They determine the embedding φ : V (n+1)(λ) ⊂ (n) (n) (n+1) V (λ). The vector a , 0 as well as the remaining eigenvectors of A belong to the orthogonal complement of the image of the embedding φ in Rn+1, φ(V (n+1)(λ))⊥. They are constructed using Theorem 4 applied to the regular matrices A˜(n−m+1) :=

A(n+1) and A˜(n−m) := A(n) . V (n+1)(λ)⊥ φ(V (n+1)(λ))⊥ The proof of Theorem 6 is deferred to Appendix C.

3. Applications to banded matrices

In the present chapter we shall discuss the application of the above formalism for banded matrices. The fact that the spectral data exceeds the number of unknowns will be shown to provide stringent constraints on the sign distribution of the off diagonal elements, so that Can one hear a matrix? 13

generically their signs are determined up to an overall sign per column. If not stated otherwise, we will assume that the considered matrices are regular.

3.1. Tridiagonal matrices

In Remark 1 above we have shown that the present construction applies for the tridiagonal case. However it is certainly not efficient since the two spectra and the signs of the off- diagonal matrix elements suffice for the purpose - as is well known [27].

3.2. Pentadiagonal matrices

The next simple case which illustrate however the main ingredients of the general case are the D = 5 - banded matrices a.k.a. pentadioagonal matrices. Here,

(n) (n) (n) a1 = a2 = ··· = an−2 = 0 ,

hence the inductive step requires solving a number of equations in the two real variables (n) (n) an and an−1. We shall address the subject using two different approaches. In the first we shall apply the method based on Theorem 4. In the second we shall use the quadratic forms (4,5), which give an alternative point of view that applies exclusively for pentadiagonal matrices. Note that for both approaches, one should compute the upper 3 × 3 minor using the procedure outlined in subsection (2.1.1). Alternatively, it suffices to provide as input the spectrum and the eigenvectors of the upper 3 × 3 minor, and use this data as the starting data for the ensuing inductive steps.

3.2.1. Applying Theorem 4 to the pentadiagonal case In this subsection, we will focus on answering the following two questions (i) Given that the spectra σ(n), n = 1,...,N correspond to a pentadiagonal matrix, what additional information is needed to uniquely reconstruct the matrix? (ii) What are the necessary conditions for a matrix A(n) and the spectrum σ(n+1) of A(n+1) so that A(n+1) is a pentadiagonal matrix? Starting with the first question, the general answer is given by Theorem 4. However, we will make use of the redundant information to show that typically, one can do away with the computation of the sign indicators as prescribed in the theorem. It will be shown that the vectors |a(n)i are determined up to an overall sign, which is uniquely provided by the easily available signs of the entries in the upper diagonal. Consider the inductive step, where A(n) and the spectrum σ(n+1) are known. All vectors |a(n)i have the same norm determined by Equation (3). Hence,

 2 2 2 1 a(n) + a(n)
= R(n)
:= tr[(A(n+1))2] − tr[(A(n))2]
− h2). n−1 n 2 Furthermore, the solutions of Equation (15) imply that

(n) (n) (n) (n) (n) (n) (n) (n) (n) (n) either an−1v (r)n−1 + an v (r)n = ξr or an−1v (r)n−1 + an v (r)n = −ξr (24) Can one hear a matrix? 14

for every r = 1, . . . , n. Thus, every |a(n)i belongs to the intersection of one of the above lines (n) (n) (for every r) with the circle of radius R (recall also that ξr were defined to be positive). By the assumption of Theorem 4, σ(n) is non-degenerate and none of the elements of (n) (n+1) (n) σ belongs to σ . By Equation (15), this implies that ξr > 0 for all r = 1, . . . , n. Then, Equation (24) defines two distinct parallel lines for each r. Moreover, for a fixed r, the intersection of the two lines with the circle occurs at two distinct pairs of points, where a single pair consists of a point and its antipode (see Fig. 2).

Figure 2. Intersecting lines from Equation (24) with the circle. The blue line corresponds (n) (n) to the line with ξr as the constant while the green line corresponds to the choice of −ξr . Dashed lines connect the antipodal solutions.

Since the spectral data comes from a pentadiagonal matrix, each of the lines from Equation (24) must have an intersection point |a˜(n)i or − |a˜(n)i where |a˜(n)i is the nth column (with only two non zero entries) of the original matrix (see Fig. (3). However, there exist special situations when such intersections with the circle give more than just one solution up to a sign. These are highly degenerate cases where all lines are determined by a fixed pair of parallel lines (l+, l−). More precisely, for every r, the lines 0 0 from Equation (24) are either equal to (l+, l−) or to another fixed pair of lines (l+, l−) where 0 l+/− ⊥ l+/− (see Fig. 4). The precise conditions when that happens are as follows. Let us 0 1 denote the slope of l+/− by α. Then, the slope of l+/− is − α . Thus, for every r the slope of 1 the lines from Equation (24) is either α or − α , i.e. (n) (n) v (r)n−1 v (r)n−1 1 either (n) = α or (n) = − for all r = 1, . . . , n. (25) v (r)n v (r)n α Define the subsets I, J ⊂ {1, . . . , n}, I ∩ J = ∅ by  (n)   (n)  v (r)n−1 v (r)n−1 1 I := r : (n) = α , J := r : (n) = − . v (r)n v (r)n α Note that, we can never have I = ∅ or J = ∅, because this implies that the eigenvectors of Can one hear a matrix? 15

Figure 3. A typical situation showing the result of intersecting lines from Equation (24) with (n) the circle for r = 1, 2, 3. Blue lines correspond to the line with ξr as the constant while (n) green lines correspond to the choice of −ξr .

Figure 4. A typical situation where intersecting all lines from Equation (24) with the circle leads to four solutions for |a(n)i.

A(n) are linearly dependent, which is a contradiction. Furthermore, because the eigenvectors of A(n) are orthonormal, we necessarily have n n n X (n) (n) X (n)
2 X (n)
2 v (r)n−1v (r)n = 0, v (r)n−1 = 1, v (r)n = 1. r=1 r=1 r=1 (n) (n) Using α, we can express vr,n−1 by vr,n and obtain that 1 1 αS − S = 0, α2S + S = 1,S + S = 1, I α J I α2 J I J Can one hear a matrix? 16

where X (n)
2 X (n)
2 SI := v (i)n ,SJ := v (j)n . i∈I j∈J It is straightforward to see that S 1 α2 = J = − 1. (26) SI SI Finally, note that conditions (25) readily imply that the intersecting lines (24) with the circle results with only four points and it is not necessary to look at the free coefficients of the lines. This is because we assume that the spectral data comes from an existing pentadiagonal matrix and hence all lines must intersect in at least one common point. To sum up, we have the following theorem. Theorem 7. For a regular pentadiagonal matrix, the inductive step leads to at most four solutions for |a(n)i. Typically, by intersecting the lines from Equation (24) there exists only one solution for |a(n)i, up to an overall sign. However, there are two possible solutions for |a(n)i, up to an overall sign, if and only if the spectral data satisfies the following condition (which we call the α-condition): there exists a nonempty subset I ⊂ {1, . . . , n}, I= 6 {1, . . . , n}, (n) (n) v (r)n−1 v (r)n−1 1 such that for all r = 1, . . . , n we have (n) = α if r ∈ I and (n) = − if r∈ / I for v (r)n v (r)n α α defined by r 1 X (n)
2 α = ± − 1,SI := v (i)n . SI i∈I Remark 5. The set of matrices satisfying the α-condition from Theorem 7 is a subset of all pentadiagonal symmetric matrices of codimension at least one. This explains the results of simulations where sets of a few millions of randomly chosen pentadiagonal matrices of dimensions up to n = 10 where tested, and none satisfied the α-condition. This fact can be understood by noting that the matrices satisfying the α-condition satisfy certain polynomial equations, as explained in the following Section 3.2.2. Thus, the set of matrices which do not satisfy the α-condition is of a positive codimension and thus it is a full-measure set. Another outcome of the above analysis provides the answer to the second question stated at the beginning of the section: It provides a set of necessary conditions that ensures that the spectral data indeed belong to a pentadiagonal matrix. The conditions are that the distance of each line from Equation (24) from the origin must be no greater than R(n), i.e. (n) ξr (n) q ≤ R for all r = 1, . . . , n. (27) (n) 2 (n) 2 (v (r)n−1) + (v (r)n)

3.2.2. Using the conic sections for recovering pentadiagonal matrices Before proceeding to the general case, we shall show that using the three quadratic forms from Equation (3), Equation (4) and Equation (5) typically allow one to reconstruct a pentadiagonal matrix from its spectral data up to an overall sign of each of its columns. For a pentadiagonal matrix, (n) we have that a has only two non-zero entries and the quadratic forms from Equation (4) Can one hear a matrix? 17

and Equation (5) describe either ellipses or hyperbolae (a.k.a. conic sections). Their explicit forms read 2  (n)  (n)
2 (n)
2 an−1 + an = R , (28) 2 (n)  (n)  (n) (n) (n) (n) (n)
2 (n) α an−1 + 2β an−1an + γ an = ρ , (29) 2 (n)  (n)  ˜(n) (n) (n) (n) (n)
2 (n) α˜ an−1 + 2β an−1an +γ ˜ an =ρ ˜ , (30) where 2 1 R(n)
= tr (A(n+1))2
− tr (A(n))2
− h2
, 2 1 det A(n+1) ρ(n) = tr[(A(n+1))3] − tr[(A(n))3] − h3
, ρ˜(n) = h − , 3 det A(n) (n) (n) (n) (n) (n) (n) α = An−1,n−1 + h, β = An−1,n, γ = An,n + h, (n) (n)
−1 ˜(n) (n)
−1 (n) (n)
−1 α˜ = A n−1,n−1 , β = A n−1,n , γ˜ = A n,n . The three conic sections are centered at (0, 0), thus each of the curves from Equation (29) and Equation (30) typically intersects the circle given by Equation (28) at exactly four points (or none). There may also be only two intersection points in certain degenerate cases. If any of the intersections is empty or both intersections have zero overlap, then there is no pentadiagonal matrix A(n+1) whose top-left main minor is A(n). If such a pentadiagonal matrix A(n+1) exists, then the three conic sections intersect at at-least two points which yield a(n) up to an overall sign. However, in certain non-generic cases it may happen that the three conic curves intersect at four points (see Fig. 5). Then, one cannot decide which of the intersection points give a

Figure 5. Intersecting the conic sections (29)-(30) with the circle (28). a) A typical situation where the two conics intersect at two antipodal points which determine a(n) up to an overall sign. b) A degenerate situation when conditions (31) are satisfied and the conic sections intersect at four points.

correct solution. This happens if and only if the principal axes of the two conic sections are identical. This can be written as a condition for the eigenvectors of the quadratic forms Can one hear a matrix? 18

defining the conic sections. Namely, both sets of eigenvectors must be the same (up to a scalar factor) [26]. This in turn happens if and only if the matrices of both quadratic forms commute. This boils down to the following algebraic condition for the coefficients defining the conic intersections (29)-(30). α(n)β˜(n) − α˜(n)β(n) +γ ˜(n)β(n) − γ(n)β˜(n) = 0 (31) In the special case when n = 2, Equation (31) is automatically satisfied and all four points of intersection give correct solutions for a(n). For n > 2, the set of pentadiagonal n × n matrices satisfying condition (31) is a strict subset of the set of all n × n pentadiagonal matrices of codimension one. Thus, a typical (generic) pentadiagonal matrix can be completely reconstructed using the above conic curves’ method. However, considering non-generic cases brings in a few subtleties. In particular, satisfying the α-condition from Theorem 7 implies satisfying Equation (31), but the converse statement is not true. The set of n×n pentadiagonal matrices satisfying the α-condition is of codimension at least one, but it is a strict subset of the set of matrices satisfying Equation (31). Thus, in non-generic cases only Theorem 7 can give a definite answer about the possible solutions (see Fig. 6).

Figure 6. Conic curves (29)-(30) for a random 5 × 5 pentadiagonal matrix satisfying degeneracy condition (31). Green lines are lines from Equation (24) which determine the correct solutions for the column a(5). Despite Equation (31) being satisfied, only one pair of antipodal points gives solutions for a(5).

3.3. A finite algorithm for reconstructing a D-diagonal matrix.

Given that the spectra σ(n), n = 1,...,N correspond to a D-diagonal real symmetric matrix, we provide a finite algorithm which produces the set of possible D-diagonal matrices whose (n) N main minors have the given spectra {σ }n=1. By definition, the nth column of a D-diagonal (n) (n) 1 matrix satisfies a1 = ... = an−d = 0, where d := 2 (D − 1). (n) (n) The algorithm starts by computing |ξr| = a v (r) from Equation (15). By n (n) choosing the sign indicators s := (s1, . . . , sn), we obtain 2 potential solutions for |a i Can one hear a matrix? 19

which read n (n) X (n) a (s) := sr|ξr| v (r) . (32) r=1 For the inductive steps involving the main minors of size n = 1, . . . , d these are all valid solutions. However, if n > d, we have additional constraints for |a(n)i of the form (n) (n) e (k) a (s) = 0 for k = 1, . . . , n − d. The above equations are satisfied only by some choices of s which determine the valid solutions in the nth inductive step. Alternatively, one can verify the solutions by demanding (n) that vector a (s) belongs to the (d − 1)-dimensional sphere, i.e. n−1 X (n) (n) 2 (n) e (k) a (s) = R , (33) k=n−d where R(n) depends only on the spectral input data via Equation (3). Thus, we have the following algorithm for realizing the inductive step of constructing A(n+1) from A(n) (Algorithm 1). One can see that the main contribution to the numerical complexity of the

Algorithm 1 Algorithm for constructing A(n+1) from A(n) for a D-banded matrix. 1: Input: A(n) and its spectral decomposition, σ(n+1),  - numerical accuracy. 2: Output: The main minor A(n+1). 1 3: d = 2 (D − 1) ← size of the new column (n) n 4: {|ξr |}r=1 ← list of overlaps from Equation (15) 5: solutions ← list of solutions for the vector a(n) 6: for s ∈ {0, 1}n do (n) 7: a (s) ← vector from Equation (32) n 2 8: P e(n)(k) a(n)(s) − R(n) ≤  if k=n−d+1 then (n)
9: solutions.append a (s) 10: return solutions

above inductive step relies on the loop going through all 2n choices of the possible signs. n Thus, the complexity scales like polyn,d2 , where the polynomial factor polyn,d comes from the necessity of evaluating expressions from Equation (32) and Equation (33) in each step of the loop. What is more, heuristics shows that the required numerical accuracy  grows with d and n, as the potential solutions tend to bunch close to the surface of the sphere. Let us next argue that for generic input data, the result of the above Algorithm 1 for  sufficiently small are just two antipodal solutions for a(n). To this end, consider the hyperplanes

(n),± (n) (n) Hr : an−d+1v (r)n−d+1 + ... + anv (r)n = ±|ξr|. The procedure of intersecting lines with the circle from Subsection 3.2 for d > 2 generalizes (n),± here to a procedure involving the study of intersections of hyperplanes Hr with the Can one hear a matrix? 20

(hyper)sphere of radius R(n). As before, typically, as a result of such a procedure one obtains only a single pair of antipodal solutions for a(n). Non-typical situations arise when some hyperplanes are perpendicular to each other. This requires the existence of certain relations between the eigenvectors of A(n) analogously to the relations stated in Theorem 7. However, the precise form of these conditions seems to be much more complicated and we leave this problem open for future research. For an illustration of the intersection procedure for a typical heptadiagonal matrix (d = 3), see Fig. 7.

(n),± Figure 7. The result of intersecting eight hyperplanes Hr , r = 1, 2, 3, 4 for a typical random heptadiagonal matrix. Intersecting one of the hyperplanes with the sphere of radius R(n) results with a single circle. The common points where n circles meet (here n = 4) determine the solutions for the vector a(n) (the red dots). As one can see, typically there are just two such points which lie on the opposite sides of the sphere. The orange and blue colors were used to mark the relevant groups of intersecting circles.

Similarly to the pentadiagonal case, one can find necessary conditions for the spectral (n),± data to correspond to a D-diagonal matrix. Every hyperplane Hr must have a nontrivial intersection with the sphere of radius R(n), hence its distance from the origin must be no greater than R(n).

|ξr| (n) q ≤ R for all r = 1, . . . , n. Pd (n) 2 k=1 (v (r)n−d+k) Remark 6. Using Theorem 6 we can apply the above methodology mutatis mutandis to non- regular banded matrices. In particular, whenever the dimension of the smaller truncated matrix from Theorem 6 is greater than d, we can use the same redundancy effects to our advantage when dealing with typical matrices. Can one hear a matrix? 21

4. Improved strategies for D-diagonal matrices – the sliding minor construction

The number of independent non vanishing entries of a D-diagonal real symmetric matrix is 1 N = (2N − d)(d + 1), where D = 2d + 1. (34) D 2

The spectral data for all the main minors of a D-diagonal N × N matrix is NN−1, which exceeds ND by far. In this section, we provide two alternative methods for reconstructing D- diagonal matrices and which use much less spectral data. The first utilizes the minimum number of spectral parameters needed for the purpose, together with the necessary sign indicators. The second method makes use of N − d − 1 more spectral data. This introduces enough redundancy to render the method much more straight forward to use at the cost of being applicable only for generic matrices.

4.1. Inverse method with minimal spectral input

The number of non vanishing data needed to write a D-diagonal matrix, ND, can also be 1 written as ND = 2 d(d + 1) + (d + 1)(N − d). It can be interpreted as the sum of two terms: (1) (2) (d) 1 (i) The spectra of the minors A = A1,1, A , ..., A which give a total of 2 d(d + 1) numbers needed to reconstruct the first d × d minor. (ii) The spectra of (N − d) minors of size (d + 1) × (d + 1) with upper diagonal entry at (d+1) (d+1) successive positions along the diagonal of A which are denoted by M1 = A , (d+1) (d+1) M2 , ..., MN−d (see (35) for an example with N = 6 and d = 2). The "sliding N−d n (d+1)o minors" in the present construction are the Mk minors. k=2   A1,1 ? ? 0 0 0    ? A ? ? 0 0   2,2     ? ? A3,3 ? ? 0    (35)    0 ? ? A4,4 ? ?       0 0 ? ? A5,5 ?    0 0 0 ? ? A6,6

The reconstruction algorithm consists of the following steps. (1) Reconstruct A(d+1) from the spectral data of A(1) through A(d) and the sign indicators 1 ( 2 d(d + 1) signs in total) as described in the inductive procedure from the proof of Theorem 4. (d) (2) The matrix elements of the d × d minor M2 starting with the upper diagonal A2,2 are (d+1) (d+1) obtained from the known spectral decomposition of A = M1 . The eigenvalues (d) and eigenvectors of M2 are then computed. (d+1) (d) (3) Use Theorem 4 to reconstruct M2 from the spectral decomposition of M2 and appropriate sign indicator data. Can one hear a matrix? 22

(4) Repeat the previous two steps recursively to procedure the successive minors (d+1) (d+1) M3 , ··· ,MN−d . The main drawback of this method is the need to provide the sign indicators at each step. To overcome this, one can also introduce minimal redundancy to the above sliding-minor method so that for a typical matrix the required sign data reduces only to the overall signs of columns of A(N). This applies only for generic matrices, as explained in the previous section.

4.2. Inverse method with optimal spectral input

N−d−1 n (d+2)o Here, the sliding minors are the Mk minors (see (36)). The increase of their k=2 dimension amounts also to the vanishing of the (1, d + 2) and (d + 2, 1) entries. This enables turning the redundant information to impose constraints, which for generic matrices remove the need for computing the sign -indicators required in the previous construction.   A1,1 ? ? 0 0 0    ? A ? ? 0 0   2,2     ? ? A3,3 ? ? 0    (36)    0 ? ? A4,4 ? ?       0 0 ? ? A5,5 ?    0 0 0 ? ? A6,6

The inversion algorithm consists of the following steps. (1) Reconstruct A(d+1) from the spectral data of A(1) through A(d) and the sign indicators 1 ( 2 d(d + 1) signs in total) as described in the inductive procedure from the proof of Theorem 4. (d+2) (2) The matrix elements of the (d + 2) × (d + 2) minor M1 starting at the upper diagonal A1,1 are computed using the algorithm proposed in section (3.3). The only sign indicator needed is the sign of the entry A1,d+1. (d+1) (d+2) (3) The (d+1), (d+1) minor M2 is extracted from the computed M1 . Its eigenvalues (d+1) and eigenvectors are then computed. (M2 is the d + 1 × d + 1 minor with A2,2 as its upper diagonal entry) (d+2) (d+1) (4) Use Theorem 4 to reconstruct M3 from the spectral decomposition of M2 and the sign of the entry A1,d+1. (5) Repeat the previous two steps recursively to produce the successive minors (d+1) (d+2) M3 , ··· ,MN−d−1. In order to minimize further the set of exceptional matrices where this strategy fails, one could choose a larger dimension for the sliding minor. In choosing an "optimum" strategy, one has to weigh computational effort against minimizing the exceptional set. This is an individual decision left for the discretion of the practitioner. Can one hear a matrix? 23

4.3. Pentadiagonal matrices

The first inversion strategy for regular pentadiagonal matrices (d = 2) is a combination of two steps which we have already analyzed in detail in the previous sections. The first step of reconstructing A(2) has been described in the proof of Theorem 4 (see Equation (19) and Equation (20)). The second as well as all other steps rely on reconstructing 3 × 3 matrices N−2 n (3)o Mk from their top-left 2 × 2 minors. This has been done explicitly in subsection k=2 2.1.1. In total, the 2 × 2 → 3 × 3 steps require the following sign data.

(i) The sign of the matrix element A1,2. N−2 n (3)o (ii) A pair of signs from Equation (23) applied to each minor Mk . k=2 The second inversion strategy which uses the larger 4 × 4 sliding minor gives extra redundancy and, as explained above, it typically reduces the required sign data to the overall signs of columns of A(N). Here, the notion of typicality is particularly straightforward to state. Namely, all the sliding minors of a typical matrix do not satisfy the α-condition from Theorem 7 . The computation of the spectral decomposition of the (3 × 3) minors can be carried out analytically.

5. Summary

We revisited here the question of how much spectral data is needed to determine a real symmetric matrix. In principle, a matrix is determined by its spectrum only up to unitary equivalence. Hence, more spectral information is needed in order to reconstruct a non- diagonal matrix. Heuristically, the number of spectral data should match the number of unknowns. Hence, for a generic real symmetric N ×N matrix, we need N(N +1)/2 numbers coming from the spectral data. In our work, these numbers come from the spectra of the main minors. Given such spectra, we provide a finite algorithm that reconstructs the matrix up to a finite number of possibilities. The construction can be made unique by supplementing additional sign data as stated in Theorem 4 which comprises of N(N − 1)/2 signs in total. However, the signs are not easily deduced from the matrix itself, therefore the requirement of supplementing the signs is a notable limitation to our procedure. Thus, in further parts of this paper, we focus on ways of reducing the required sign data in the case of banded matrices. For (2d + 1)-banded matrices, there are less unknowns, hence one can find a smaller optimal set of (d + 1) × (d + 1) minors whose spectra allow one to reconstruct the matrix up to a finite number of possibilities. As before, some additional sign information is required to make the construction unique. By analysing the case of pentadiagonal matrices, we compare both methods of reconstruction. We find that the redundant information coming from the spectra of all main minors in a typical case (where the notion of typicality has been made precise in Subsection 3.2) allows one to reduce the number of required signs to N − 1 that are just the overall signs of columns of the upper-triangular part of A(N), while using the optimal amount of spectral information requires 2N −1 signs. In Section 4 we argue that this is a general fact, i.e. the redundancy coming from all main minors typically allows one to uniquely reproduce Can one hear a matrix? 24

a general banded matrix from the spectral data using only the overall signs of columns of the upper-triangular part of A(N).

Acknowledgment

We would like to thank Professor Raphael Loewy for suggestions and critical comments during the first stages of the present work. Professor Percy Deift accompanied the work since its beginning and offered critique and advise in many instances. We are obliged for his invaluable help and in particular for his permission to quote an unpublished theorem quoted in Remark 3. Professor Carlos Meite is acknowledged for bringing to our attention his work (Remark 4) before its publication. We also thank Professor Jon Keating for useful suggestions and Doctor Yotam Smilansky for helpful comments on the manuscript. The inspiration for the work came by listening to a seminar by Professor Christiane Tretter within the webinar series "Spectral Geometry in the Clouds". Thanks Christiane. The seminar is organized by Doctors Jean Lagacé and Alexandre Girouard whom we thank for initiating and running such successful seminar during the miserable Corona days.

Appendix A. Identities for Cauchy Matrices

n n For two sequences of real numbers, x := {xi}i=1 and y := {yi}i=1, define the Cauchy matrix C(x, y) as (C(x, y)) := 1 . The matrix C(x, y) has the followng symmetry i,j xi−yj (C(x, y))T = −C(y, x). (A.1) C will be used to denote C(x, y) by default. An explicit expression for C−1 reads [33] " n # " n # −1 Y yi − xk 1 Y xj − yk (C )i,j = (yi − xi) (xj − yj) (A.2) yi − yk yi − xj xj − xk k6=i k6=j The sum of the matrix elements of C−1 is known to be [35] n n X −1 X (C )i,j = (xk − yk) (A.3) i,j=1 k=1 We aim to prove the following identities Qn Qn X 1 yi X 1 xi C−1
= 1 − i=1 , C−1
= i=1 − 1, (A.4) i,j x Qn x y i,j Qn y i,j j i=1 i i,j i i=1 i Qn n Qn n X 1 yi X X 1 xi X y C−1
= i=1 (x − y ) , C−1
x = i=1 (x − y ) .(A.5) i i,j x Qn x j j y i,j j Qn y j j i,j j i=1 i j=1 i,j i i=1 i j=1 Note first, that the LHS identities imply the RHS identities by symmetry A.1. Hence, it is enough to prove the LHS identities. To this end, we write them down as the following matrix- vector products (n) (n) X −1
1 (n) −1 1 X −1
1 (n) −1 1 C = 1 C | i , yi C = y C | i , (A.6) i,j x x i,j x x i,j j i,j j Can one hear a matrix? 25

(n) 1 (n) (n) where vectors |1 i, | x i and |y i are defined as       1 1/x1 y1 (n)  1  1  1/x2   y2  |1(n)i :=   , | i :=   , |y(n)i :=   .  .   .   .   .  x  .   .  1 1/xn yn We will use the following lemma. Lemma 8 (Vandermonde Matrix Identity for Cauchy Matrix). C(x, y) has the following decomposition

−1 −1 C(x, y) = −PVx VyQ , where P and Q are diagonal matrices defined as

P = diag (p1(x1), . . . , pn(xn)) ,Q = diag (p(y1), . . . , p(yn)) with n n Y Y p(t) := (t − xi), pk(t) := (t − xi), 1 ≤ k ≤ n, i=1 i=1,i6=k and Vx, Vy are the Vandermonde matrices  1 1 ··· 1   1 1 ··· 1 

 x1 x2 ··· xn   y1 y2 ··· yn  V =   ,V =   . x  . . .. .  y  . . .. .   . . . .   . . . .  n−1 n−1 n−1 n−1 n−1 n−1 x1 x2 ··· xn y1 y2 ··· yn By Lemma 8, we have

−1 −1 −1 C = −QVY VxP .

Next, we define the vectors |φ1i, |φyi and |φxi as (n) (n) −1 (n) −1 −1 1 hφ1| := − 1 QV , hφy| := − y QV , |φxi := VxP | i . (A.7) Y Y x Calculation of A.6 boils down to computing overlaps X 1 X 1 C−1
= hφ | φ i , y C−1
= hφ | φ i . i,j x 1 x i i,j x y x i,j j i,j j

Let us start with finding |φxi. A straightforward calculation using Equation (A.7) shows  Pn 1  k=1 xkpk(xk)  Pn 1   k=1 pk(xk)  |φxi =  .  . (A.8)  .   n−2  Pn xk k=1 pk(xk) In order to compute the sums in Equation (A.8), we use the following Lemma [34]. Can one hear a matrix? 26

Lemma 9 (Summation of Powers over Product of Differences). Let (t1, t2, . . . , tn) be a sequence of real numbers. Then,  0 : 0 ≤ m < n − 1   n m !  1 : m = n − 1 X tk  Q = n (tk − ti) k=1 i=1,i6=k  X  tj : m = n  j=1 D E (n) By Lemma 9, we have ej φx = for j ≥ 2. The only nonzero entry is n n D (n) E X 1 (−1) e1 φx = = − , (A.9) xkpk(xk) x1· ... · xn k=1 which is evaluated in the following Proposition.

Proposition 10. Let (t1, t2, . . . , tn) be a sequence of real numbers. Then, n X 1 (−1)n Q = − . tk (tk − ti) t1· ... · tn k=1 i6=k

Proof. Consider the following complex integral over a circle of radius R > max{|ti| : 1 ≤ i ≤ n} 1 I dz IR(t1, . . . , tn) := . 2πi |z|=R z(z − t1)· ... · (z − tn)

Because the integrand is a regular function at z = ∞, we have IR(t1, . . . , tn) = 0. On the other hand, the integrand has simple poles at z = 0, z = t1, ..., z = tn. Hence, the residue theorem asserts that n n 1 X 1 IR(t1, . . . , tn) = (−1) + Q , t1· ... · tn tk (tk − ti) k=1 i6=k which finishes the proof.

Next, let us compute the vectors |φ1i and |φyi. Note first that because the first entry of |φxi is its only nonzero entry, we only need to find the first entries of |φ1i and |φyi. By Equation (A.7), the desired vectors are solutions to the following linear equations T T (n) T T (n) VY |φ1i = −Q |1 i ,VY |φyi = −Q |y i .     p(y1) y1p(y1)  p(y2)   y2p(y2)  QT |1(n)i =   ,QT |y(n)i =   .  .   .   .   .  p(y3) ynp(yn) We use now Crammer’s rule.  2 n−1  p(y1) y1 y1 ··· y1 2 n−1 D E 1  p(y2) y2 y ··· y  e(n) φ = −  2 2  1 1 T  ......  det Vy  . . . . .  2 n−1 p(yn) yn yn ··· yn Can one hear a matrix? 27

 2 n−1  y1p(y1) y1 y1 ··· y1 2 n−1 D E 1  y2p(y2) y2 y ··· y  e(n) φ = −  2 2  . 1 y T  ......  det Vy  . . . . .  2 n−1 ynp(yn) yn yn ··· yn Applying the Laplace expansion of both determinants with respect to their first columns. n D (n) E 1 X T e φ = − (−1)jp(y ) det (V (y ,..., yˆ , . . . , y )) , 1 1 det V T j 1 j n y j=1 n D (n) E 1 X T e φ = − (−1)jy p(y ) det (V (y ,..., yˆ , . . . , y )) , 1 y det V T j j 1 j n y j=1 where  2 n−1  y1 y1 ··· y1  . . .. .   . . . .   2 n−1  T  yj−1 yj−1 ··· yj−1  (V (y1,..., yˆj, . . . , yn)) :=    y y2 ··· yn−1   j+1 j+1 j+1   . . .. .   . . . .  2 n−1 yn yn ··· yn Using the expression for the determinant of a Vandermonde matrix, we have ! ! T Y Y det (V (y1,..., yˆj, . . . , yn)) = yk · (yl − yk) k6=j k

where en+1−l ({x1, . . . , xn}) is the (n+1−l)th Viete sum in variables {x1, . . . , xn}. Plugging the above expansion into Equation (A.10), we get n n−1 n n l−2 ! D (n) E n X yj X X yj e1 φ1 = (−1) y1· ... · yn + al . pj(yj) pj(yj) j=1 l=1 j=1

n n n n l−1 ! D (n) E n X yj X X yj e1 φy = (−1) y1· ... · yn + al . pj(yj) pj(yj) j=1 l=1 j=1 Can one hear a matrix? 28

By Lemma 9, we get that n ! D (n) E X 1 e φ = (−1)ny · ... · y 1 + a , 1 1 1 n 1 y p (y ) j=1 j j j n ! D E (n) n X e1 φy = (−1) y1· ... · yn yj + an . j=1 n Finally, after making use of Proposition 10 and expanding Viete sums a1 = (−1) x1· ... · xn and an = −(x1 + ... + xn) D E D E n (n) n (n) n X e1 φ1 = (−1) (y1· ... · yn − x1· ... · xn), e1 φy = (−1) y1· ... · yn (yj − xj). j=1 D E (n) The result follows now directly by multiplying the above expressions by e1 φx from Equation (A.9). The identity which is used in the discussion of (17) follows by combining Equation (A.4) and Equation (A.5) to give n X 1 −1 X s = (C )i,j(xj − s) , where s = (xk − yk) . (A.11) yi i,j k=1

(n) 2 Appendix B. Derivation of the formula for (ξr ) .

The aim is to prove the identity used in Equation (15) which reads n Qn+1 X (yr − xk) C−1
(x − h) = − k=1 , r,k k Qn (B.1) (yr − yk) k=1 k=1,k6=r Pn+1 Pn where C = C(x, y) is a Cauchy matrix and h := k=1 xk − k=1 yn. Note first that n n X −1
(n) −1 (n) X −1
(n) −1 (n) C r,k xk = er C |x i , C r,k = er C |1 i . (B.2) k=1 k=1 −1 −1 −1 Next, let us use Lemma 8 to find that C = −QVY VxP . By a reasoning relying on the use of Lemma 9 (which is analogous to the one used in Equation (A.9)), we obtain that   0   . 0  .  . −1 (n)   −1 (n)  .  VxP |x i =  0  ,VxP |1 i =   (B.3)    0   1  Pn 1 k=1 xk −1 The next step is the multiplication by Vy . To this end, we apply the following result concerning the inverse of a Vandermonde matrix [34]. (n) Lemma 11. Let Vt be the Vandermonde matrix of size n in variables t1, . . . , tn, i.e. j−1 Vi,j = ti , 1 ≤ i, j ≤ n. Then  −1 (−1)n−je ({t , . . . , t }/{t }) V (n) = n−j 1 n i , t Qk=n i,j k=1,k6=i(ti − tk) Can one hear a matrix? 29

where em denotes the mth Viete sum. −1 We only need the last two columns of (Vy) , because using Equation (B.3), we get n ! (n) −1 (n) −1 −1 X er C |x i = −p(yr) ((Vy) )r,n−1 + ((Vy) )r,n xk , (B.4) k=1 (n) −1 (n) −1 er C |1 i = −p(yr)((Vy) )r,n. (B.5) −1 By Lemma 11, the necessary matrix elements of Vy forms read Pn yk 1 ((V )−1) = − k=1,k6=r , ((V )−1) = . y r,n−1 Qk=n y r,n Qk=n k=1,k6=r(yr − yk) k=1,k6=r(yr − yk) Plugging the above results to the LHS of Equation (B.1) yields n n ! X p(yr) X C−1
(x − h) = h − y − (x − y ) , r,k k Qk=n r k k k=1 k=1,k6=r(yr − yk) k=1 which after substituting expressions for h and p(yr) yields the desired result.

Appendix C. Non-regular matrices: proof of Theorem 6

Let us start with the simplest situation where for some n we have σ(n) ∩ σ(n+1) = 1 and both spectra are non-degenerate. The other cases can be analyzed by essentially repeating and slightly adjusting the proof of the simple case. By the interlacing property (6), the scenario described above happens if and only if there exists a single l ∈ {1, . . . , n} such that either (n) (n+1) (n) (n+1) (n) λl = λl or λl = λl+1 . Equation (11) applied to k = r = l, gives bl,n+1ξl = 0, i.e. (n) i) bl,n+1 = 0 or ii) ξl = 0. The two cases are analyzed separately. (n+1) (n) (i) bl,n+1 = 0. Equation (11) applied to k = l and r 6= l gives (λl − λr )bl,r = 0, (n+1) (n) which in turn implies that bl,r for r 6= l (recall that λl 6= λr by the assumption (n) (n+1) (n+1) (n) that σ ∩ σ = 1). Thus, the eigenvector v˜ (l) = v (l), 0 . Equation (10) (n) (n) (n) applied for k = l implies v (l) a = 0, thus bl,n+1 = 0 leads to ξl = 0. (n) (ii) ξl = 0 and bl,n+1 6= 0. Here, we will show that this is in contradiction with the

assumption that σ(n) ∩ σ(n+1) = 1. To this end, we combine Equations (11) for r 6= l with Equations (10) to obtain the following reduced set of equations involving a Cauchy matrix 2  (n) n ξr (n+1) ˜ X λk − h = (n+1) (n) , ∀ 1 ≤ k ≤ n + 1, (C.1) r=1,r6=l λk − λr ˜ Pn+1 (n+1) Pn (n) where h = k=1,k6=l λk − k=1,k6=l λk . Note, that Equations (C.1) for k 6= l (n−1) (n) (n) can be viewed as a set equations for the regular spectra σ˜ := σ /{λl } and 2 (n) (n+1) (n+1)  (n) σ˜ := σ /{λl }. Thus, solutions for ξr , r 6= l read Qn+1  (n) (n+1) k=1,k6=l λr − λk (ξ(n))2 = − . (C.2) r Qn  (n) (n) k=1,k∈{ / r,l} λr − λk Can one hear a matrix? 30

2  (n) However, the above formula for ξr has to be consistent with Equation (C.1) for (n) (n+1) k = l. By denoting λl = λl := λ, this boils down to the following condition Qn+1  (n) (n+1) n λr − λ  ˜ X 1 k=1,k6=l k λ − h = − (n)  , (C.3) λ − λ Qn (n) (n) r=1,r6=l r k=1,k∈{ / r,l} λr − λk We view Equation (C.3) as a polynomial equation for λ. It turns out that Equation (C.3) Qn+1  (n+1) can be simplified to the equation k=1,k6=l λ − λk = 0. This is in contradiction (n+1) with the regularity assumption, because it implies that λ = λl+1 . (n) (n+1) Summing up the above analysis of the simplest case where λl = λl for a single l and (n) both spectra are non-degenerate, we have obtained that we necessarily have ξl = 0 and 2 (n+1) (n)  (n) v˜ (l) = v (l), 0 . The remaining solutions for ξr , r 6= l are given by Equation

(C.2). Moreover, Equation (11) applied for k 6= l yields bk,r = 0 whenever r = l and (n+1) (n+1) k 6= l. Thus, the expressions for the eigenvectors v˜ (r) , r 6= l, of the matrix A are formally identical with the expressions for the eigenvectors of a regular matrix A˜(n) obtained (n) (n) (n+1) by removing the eigenspaces corresponding to the eigenvalue λl in A and A . Let us next consider the situation where there is a single degeneracy block in σ(n) and σ(n+1). Note first that whenever k ∈ D(n+1)(l, m) and r ∈ D(n)(l, m), then Equation (11) (n) (n+1) implies that bk,n+1ξr = 0. Thus, we have two possibilities for each k ∈ D (l, m):

(i) bk,n+1 = 0 or (n) (n) (ii) bk,n+1 6= 0 and ξr = 0 for all r ∈ D (l, m). (n+1) (n+1) As the following Lemma states, case (ii) is only relevant when D (l, m) = DI (l, m). Lemma 12. Assume that the spectrum σ(n) has degeneracy indices D(n)(l, m) and denote the degenerate eigenvalue by λ. If the degeneracy indices in the spectrum σ(n+1) are of the (n+1) (n+1) (n+1) (n+1) form DII (l, m), DIII (l, m) or DIV (l, m), then for every k ∈ D (l, m) we have (n+1) (n+1) v˜ (k) e (n + 1) = 0, i.e. the coefficient bk,n+1 = 0.

(n+1) (n+1) Proof. (Ad absurdum.) Assume that bk,n+1 6= 0 for some k ∈ D (l, m) and D (l, m) (n+1) 0 (n+1) is of the form Di (l, m) with i ∈ {II,III,IV }. Denote K := {k ∈ D (l, m): 0 (n) bk,n+1 6= 0}. Then, Equation (11) applied for any k ∈ K and r ∈ D (l, m) (n) (n) 0 implies that ξr = 0 for all r ∈ D (l, m). On the other hand, for any k ∈ K ∪ {1, . . . , n + 1}/D(n+1)(l, m)
, Equation (11) gives

(n) ξr bk,r = (n+1) (n) bk,n+1. λk − λr This in turn combined with Equation (10) for the same k yields 2 2 P  (n)  (n) (n) ξr ξr  (n+1)  r∈D (l,m) X λk − h = (n+1) + (n+1) (n) . (C.4) λk − λ r∈D / (n)(l,m) λk − λr Can one hear a matrix? 31

In particular, we can use the Cauchy matrix method to solve the above equations for 2 2 P  (n)  (n) (n) (n) r∈D(n)(l,m) ξr or ξr with r∈ / D (l, m). However, each case with Di , i ∈ {II,III,IV } has to be treated separately.

(n+1) (n+1) (i) D (l, m) = DIV (l, m). Consider the set of Equations (C.4) for k ∈ {1, . . . , n + (n+1) 1}/DIV (l, m). They are Cauchy equations for the regular spectra of a reduced (n−m+1) (n+1) (n+1) (n−m) dimension σ˜ = {λk : k ∈ {1, . . . , n + 1}/DIV (l, m)} and σ˜ = 2 2 (n) (n) ˜(n−m) P  (n) {λ} ∪ {λk : k ∈ {1, . . . , n}/D4 (l, m)} with ξ1 := r∈D(n)(l,m) ξr . In analogy to Equation (15), we solve Equation (C.4) to obtain   Qn+1 λ − λ(n+1)  2 k=1,k∈D / (n+1)(l,m) k ξ˜(n−m) = − IV . (C.5) 1 Qn  (n) k=1,k∈D / (n)(l,m) λ − λk  2 0 ˜(n−m) As we mentioned earlier, K 6= ∅ implies ξ1 = 0, thus the above expression (n+1) (n+1) implies that λk = λ for some k∈ / DIV (l, m), which is in contradiction with the Lemma’s assumptions. (n+1) (n+1) (n+1) (n+1) (ii) D (l, m) = DII (l, m) or D (l, m) = DIII (l, m) Consider the set of (n+1) (n) Equations (C.4) for k ∈ {1, . . . , n + 1}/D (l, m) with ξr = 0 whenever r ∈ D(n)(l, m). They are Cauchy equations for the regular spectra of a reduced dimension (n−m) (n+1) (n+1) (n−m−1) (n) σ˜ = {λk : k ∈ {1, . . . , n + 1}/D (l, m)} and σ˜ = {λk : k ∈ {1, . . . , n}/D(n)(l, m)}. In analogy to Equation (15), we solve Equation (C.4) to obtain Qn+1  (n) (n+1) k=1,k∈D / (n+1)(l,m) λr − λk (ξ(n))2 = − , r ∈ {1, . . . , n}/D(n)(l, m). (C.6) r Qn  (n) (n) k=1,k∈D / (n)(l,m),k6=r λr − λk

(n) 2 The above expressions for (ξr ) have to be compatible with Equation (C.4) applied to k ∈ K0. This is in full analogy to Equation (C.3) which we treated as a polynomial equation for λ. This polynomial equation is equivalent to n+1 Y  (n+1) λ − λk = 0. k=1,k∈D / (n+1)(l,m)

(n+1) (n+1) This in turn implies that λk = λ for some k∈ / D (l, m), which is in contradiction with this Lemma’s assumptions.

(n+1) (n+1) Thus, whenever D (l, m) 6= DI (l, m), we have bk,n+1 = 0 for all k ∈ D(n+1)(l, m). Let us now focus on this particular case. By Equation (11) we get that for (n+1) (n) (n+1) every k ∈ D (l, m) and r∈ / D (l, m) bk,r = 0. Thus, vector v˜ (k) effectively belongs to the degenerate subspace of A(n), i.e. (n+1) X (n) (n+1) v˜ (k) = bk,r v (r), 0 , k ∈ D (l, m). (C.7) r∈D(n)(l,m) Can one hear a matrix? 32

Depending on the precise form of the set D(n+1)(l, m), the dimension of the degenerate subspace of A(n+1) can be: (i) equal to the dimension of the degenerate subspace of A(n); this happens when (n+1) (n+1) D (l, m) = Di (l, m) with i = II,III (ii) one dimension smaller than the degenerate subspace of A(n); this happens when (n+1) (n+1) D (l, m) = DIV (l, m). (n+1) (n+1) If D (l, m) = Di (l, m) with i = II,III, then Equation (C.7) effectively (n) describes just a change of basis. Thus, coefficients br,k with r ∈ D (l, m) and k ∈ (n+1) (n+1) (n+1) D (l, m) form an orthogonal matrix. If D (l, m) = DIV (l, m), then Equation (C.7) geometrically means that the degenerate eigenspace of A(n+1) can be any subspace of the degenerate eigenspace of A(n) of codimension one. Equation (10) yields the condition that (n+1) (n) (n+1) v˜ (k) a = 0 for all k ∈ D (l, m), i.e. X (n) (n+1) bk,rξr = 0, k ∈ D (l, m). (C.8) r∈D(n)(l,m)

(n+1) (n+1) If D (l, m) = Di (l, m) with i = II,III, the dimensionality of the degenerate (n) (n) eigenspaces together with Equation (C.8) forces ξr = 0 for all r ∈ D (l, m). This together with Equation (11) imply that any eigenvector of A(n+1) corresponding to a non- degenerate eigenvalue is orthogonal to the degenerate eigenspace of A(n). Thus, we have effectively reduced the dimension of the problem by m + 1 and we can apply Theorem 4 for the regular reduced matrices A˜(n−m−1) and A˜(n−m) with the degenerate eigenspaces removed. In particular, (n) ξr (n) (n+1) bk,r = (n+1) (n) bk,n+1, r∈ / D (l, m), k∈ / D (l, m), λk − λr Qn+1  (n) (n+1) (n+1) λr − λk (n) X (n) (n) (n)
2 k=1,k∈D / (l,m) a = ξ v (r), 0 , ξ = − r r Qn  (n) (n) r∈D / (n)(l,m) k=1,k∈D / (n)(l,m),k6=r λr − λk (n+1) and bk,r = 0 for k∈ / D (l, m) otherwise. (n+1) (n+1) Let us next move to the case when D (l, m) = DIV (l, m). The geometric (n) interpretation of Equation (C.8) is that a is orthogonal to the degenerate eigenspace of A(n+1) embedded in the degenerate eigenspace of A(n). Consider next Equation (C.4) (n+1) from the proof of Lemma 12 applied for k∈ / DIV (l, m). This is a set of n − m + 1 (n−m+1) (n+1) equations involving regular spectra of a reduced dimension σ˜ = {λk : k ∈ (n+1) (n−m) (n) (n) {1, . . . , n + 1}/DIV (l, m)} and σ˜ = {λ} ∪ {λk : k ∈ {1, . . . , n}/D (l, m)}. The reduced regular matrix A˜(n−m+1) is obtained from A(n+1) by removing the entire degenerate eigenspace with eigenvalue λ. The reduced regular matrix A˜(n−m) is determined by the embedding of the degenerate eigenspace of A(n+1) into the degenerate eigenspace of A(n) – the image of this embedding is removed. Next, we use the familiar method of Theorem 4 applied the above regular matrices. However, in contrast to the previous case the coefficients (n) (n) (n+1) ξr and bk,r for r ∈ D (l, m) and k∈ / D (l, m) may be nonzero and are determined by the choice of the embedding of eigenspaces. Can one hear a matrix? 33

(n+1) (n+1) (n+1) Finally, let us consider D (l, m) = DI (l, m). Firstly, note that |DI (l, m)| = |D(n)(l, m)| + 1 = m + 2, thus the degenerate eigenspace of A(n+1) cannot be spanned only (n) (n) by the vectors v (r), 0 with r ∈ D (l, m). Hence, we necessarily have bk,n+1 6= 0 for (n+1) (n) (n) some k ∈ DI (l, m). By Equation (11) this implies that ξr = 0 for all r ∈ D (l, m), (n) (n) i.e. a is orthogonal to the degenerate eigenspace of A . This in turn implies that for (n+1) (n) any k∈ / DI (l, m) we have bk,r = 0 whenever r ∈ D (l, m), i.e. the non-degenerate eigenvectors of A(n+1) are orthogonal to the degenerate eigenspace of A(n). Thus, we consider the regular matrices A˜(n−m) and A˜(n−m−1) which are the original matrices truncated to the space defined as the orthogonal complement of the degenerate eigenspace of A(n). They (n−m) (n+1) (n+1) have the regular spectra σ˜ = {λ} ∪ {λk : k ∈ {1, . . . , n + 1}/DI (l, m)} and (n−m−1) (n) (n) σ˜ = {λk : k ∈ {1, . . . , n}/D (l, m)}. Thus, we use Theorem 4 to find the (n) (n+1) vector a as well as the non-degenerate eigenvectors of A . What is left to find are the remaining m + 2 degenerate eigenvectors of A(n+1) corresponding to the eigenvalue λ. Their (n) coefficients bk,r with r∈ / D (l, m) are proportional to bk,n+1 as stated in Equation (11), i.e. (n) ξr (n) (n+1) bk,r = (n) bk,n+1, r∈ / D (l, m), k ∈ DI (l, m). λ − λr Note, that the equations n X (n) (n+1) (λ − h)bk,n+1 = bk,rξr , k ∈ DI (l, m). r=1,r∈D / (n)(l,m) are now automatically satisfied by the above expressions for bk,r. The unknown coefficients (n+1) (n) bk,n+1 and br,k with k ∈ DI (l, m) and r∈ / D (l, m) are completely free to choose so that  (n+1) the resulting vectors v˜ (k) (n+1) form an orthonormal basis of the degenerate k∈DI (l,m) eigenspace.

References

[1] G.M.L. Gladwell Inverse Problems in Vibration, Kluwer Academic Publishers, 2004 [2] V. Marchenko, Inverse Problems in the Theory of Small Oscillations, Translations of Mathematical Monographs vol 247, AMS, 2018 [3] F. R. Gantmacher and M. G. Krein, Oszillationsmatrizen, Oszillations-keme und kleine Schwingungen mechanischer Systeme, Akademie-Verlag, Berlin, 1960. [4] G. M. L. Gladwell, Antonino Morassi (Eds.), Dynamical Inverse Problems: Theory and Application, CISM International Centre for Mechanical Sciences, vol. 529, Springer-Verlag Wien, 2011 [5] G.M.L. Gladwell., The inverse problem for the vibrating beam, Proceedings of the Royal Society of London - Series A, 393:277–295, 1984 [6] Steven J. Cox, Mark Embree, Jeffrey M. Hokanson, One Can Hear the Composition of a String: Experiment with an Inverse Eigenvalue Problem SIAM REVIEW 54, 157–178, (2012 ) [7] V. Pivovarchik, N. Rozhenko and C. Tretter, Dirichlet–Neumann inverse spectral problem for a star graph of Stieltjes strings, Linear Algebra and its Applications 439, 2263–2292, (2013) [8] P. B. Denton, S. J. Parke, and X. Zhang, Neutrino oscillations in matter via eigenvalues, Phys. Rev. D 101, 093001 (2020) [9] P. B. Denton, S. J. Parke, T. Tao, Xining Zhang, Eigenvectors from eigenvalues: A survey of a basic identity in linear algebra, BULLETIN (New Series) OF THE AMERICAN MATHEMATICAL SOCIETY, 2021 Can one hear a matrix? 34

[10] P. Deift, L. Li, T. Nanda, and C. Tomei, The Toda flow on a generic orbit is integrable, Commun. Pure Appl. Math. 39, 183–232 (1986) [11] D. Burgarth and A. Ajoy, Evolution-Free Hamiltonian Parameter Estimation through Zeeman Markers, Phys. Rev. Lett. 119, 030402, 2017 [12] D. Burgarth, Identifying combinatorially symmetric Hidden Markov Models, Electronic Journal of Linear Algebra, Volume 34, pp. 393-398, 2018 [13] H. Hochstadt, On some inverse problems in matrix theory, Arch. Math. 18 (1967), 201-207. [14] H. Hochstadt, On the construction of a Jacobi matrix from spectral data , Linear Algebra Appl. 8 (1974), 435446. [15] Ole H. Hald, Inverse Eigenvalue Problems for Jacobi Matrices, LINEAR ALGEBRA AND ITS APPLICATIONS 14, 63-85 (1976) [16] Fritz Gesztesy, Inverse spectral theory as influenced by Barry Simon In Spectral Theory and Mathematical Physics: A Festschrift in Honor of Barry Simon’s 60th Birthday Fritz Gesztesy, Percy Deift, Cherie Galvez, Peter Perry, and Wilhelm Schlag, Editors ISBN-10: 0-8218-3783-4 (2010) Article electronically published on February 18, 2021 [17] Boley D., Golub G.H. Inverse eigenvalue problems for band matrices, In: Watson G.A. (eds) Numerical Analysis. Lecture Notes in Mathematics, vol 630. Springer, Berlin, Heidelberg, 1978 [18] F. W. Biegler-König, Construction of band matrices from spectral data, Linear Algebra and its Applications Vol. 40, Pages 79-87, 1981 [19] D. Boley and G. H. Golub, A survey of matrix inverse eigenvalue problems, Inverse Problems 3 595, 1987 [20] M. Chu and G. Golub, Inverse Eigenvalue Problems: Theory, Algorithms, and Applications, Oxford University Press, 2005 [21] Kent Griffin, Michael J. Tsatsomeros, Principal minors, Part II: The principal minor assignment problem, Linear Algebra and its Applications 419 125–171, 2006 [22] R. Loewy, Principal minors and diagonal similarity of matrices, Linear Algebra and its Applications 78 23–64, 1986 [23] Shmuel Friedland, Matrices with prescribed off-diagonal elements, lsrael J. Math. 11:184-189 (1972). [24] V. P. Kotlyarov, V. A. Marchenko, and V. V. Slavin, Reconstruction of the Hermitian Matrix by Its Spectrum and Spectra of Some Number of Its Perturbations, Doklady Mathematics, Vol. 94, No. 2, pp. 529–531, 2016 [25] R. S. Martin, J. H. Wilkinson, Symmetric Decomposition of Positive Definite Band Matrices, Numer. Math. 7, 355 - 361 (1965) [26] Ayoub, A. B., The central conic sections revisited, Mathematics Magazine, 66 (5): 322–325, 1993 [27] L. Anderson, On the effective determination of the wave operator from given spectral data in the case of a difference equation corresponding to a Sturm-Liouville differential equation Math. Anal. A 29 (1970), 467497. [28] Gerschgorin S Über die Abgrenzung der Eigenwerte einer Matrix Izv. Akad. Nauk. USSR Otd. Fiz.-Mat. Nauk 6 749–54 (1931) [29] Percy Deift, Private communication (2021). [30] Ricardo S. Leite and Carlos Tomei, The geometry of spectral interlacing, arXiv:2108.06285 (2021) [31] Daniel Perrin, Algebraic Geometry, an introduction. Springer-Verlag London Limited (2008) . [32] Q. I. Rahman and G. Schmeisser, Analytic theory of polynomials, London Mathematical Society Monographs. New Series, 26, The Clarendon Press Oxford University Press, Oxford, 2002, ISBN 0- 19-853493-0. MR 2004b:30015 [33] R. Gow Cauchy’s matrix, the Vandermonde matrix and polynomial interpolation Irish Math. Soc. Bull 28, 45-52 (1992) [34] D. E. Knuth, Art of Computer Programming, The: Volume 1: Fundamental Algorithms, Addison-Wesley; 3rd edition (1997) [35] See https://proofwiki.org/wiki/VandermondeMatrixIdentityforCauchyMatrix