Symmetries of the Feinberg-Zee Random Hopping Matrix Raffael Hagger November 6, 7 arxiv:4937v [mathsp] 3 Sep 5 Abstract We study the symmetries of the spectrum of the Feinberg-Zee Random Hopping Matrix introduced in [6] and studied in various papers therafter (eg [], [], [3], [7], [9]) In [3], Chandler-Wilde and Davies proved that the spectrum of the Feinberg-Zee Random Hopping Matrix is invariant under taking square roots, which implied that the unit disk is contained in the spectrum (a result already obtained slightly earlier in []) In a similar approach we show that there is an infinite sequence of symmetries at least in the periodic part of the spectrum (which is conjectured to be dense) Using these symmetries and the result of [3], we can exploit a considerably larger part of the spectrum than the unit disk As a further consequence we find an infinite sequence of Julia sets contained in the spectrum These facts may serve as a part of an explanation of the seemingly fractal-like behaviour of the boundary Mathematics Subject Classification: Primary 47B8; Secondary 47A, 47B36 Keywords: random hopping, random operator, spectrum, tridiagonal, periodic, symmetry Introduction In recent years some progress has been made in the study of non-self-adjoint random operators (see [], [], [3], [7], [9] and references therein) However, still a lot of questions remain open even in the tridiagonal case In particular, the spectrum of most tridiagonal random operators is yet unknown Although the techniques used here are possibly also suited for more general types of operators, we focus on the Feinberg-Zee Random Hopping Matrix [6] here It is defined as follows: A b := b b L(l (Z)) for a random sequence b {±} Z, ie a sequence with randomly (iid) distributed entries Similarly we define A c for arbitrary sequences c {±} Z The spectrum sp(a b ) := { λ C : A b λi is not invertible } of A b is independent of the sequence b in the following sense: raffaelhagger@tuhhde
Proposition (eg [, Lemma 3, Theorem 5]) For a sequence of iid random variables (b j ) j Z taking values in {±} with non-zero probability each, the spectrum of A b is given by sp(a b ) = Σ := sp(a c ) () c {±} Z almost surely In [] Chandler-Wilde, Chonchaiya and Lindner discovered a beautiful connection to the Sierpinski triangle This connection was then used to show that the unit disk is contained in Σ, which disproved earlier conjectures that Σ may be of fractal dimension In [] the same authors also gave an upper bound using the numerical range This upper bound was further improved in [7] An image of Σ (or more precisely: what it is conjectured to be) and the upper bounds computed in [] and [7] are provided in Figure Figure : The conjectured shape of Σ (black), the upper bound computed in [] (red) and the upper bound computed in [7] (blue) Another way of deriving lower bounds is using equation () directly The spectra of many operators on the right hand side can be computed explicitly This is in particular the case for periodic operators We call an operator A L(l (Z)) (not necessarily tridiagonal) m-periodic if A i,j = A i+m,j+m for all i, j Z In the tridiagonal case, we will also use the notation A k per for
k {±} m, ie A k per = k m k L(l (Z)) k m k Note that k is not unique because (k,, k m ) {±} m and (k,, k m, k,, k m ) {±} m define the same operator We will make use of this fact later on The good thing about periodic operators is that they can be diagonalized in some sense Proposition (eg [4, Theorem 449]) Let m N, let A L(l (Z)) be an m-periodic operator and define A rm+, A rm+,m a r := A rm+m, A rm+m,m for all r Z If r Z a r <, then A is unitarily equivalent to the (generalized) multiplication operator M a L(L ([, π), C m )) The function a L ([, π), C m m ) is called the symbol of A and given by a(ϕ) = r Z a r e irϕ (ϕ [, π)) Moreover, sp(a) = {λ C : det(a(ϕ) λi m ) = for some ϕ [, π)} = ϕ [,π) sp(a(ϕ)) In the following, we will always have a finite number of non-vanishing coefficients a r Thus we will not have to worry about the convergence of r Z a r Proposition enables us to compute the spectra of periodic operators explicitly (numerically) We call the set π := m N k {±} m sp(a k per) the periodic part of Σ Another related set is the set of eigenvalues of finite matrices of this kind For n N and k {±} n we define A k fin := k C(n+) (n+) k n 3
Then we call the set σ := n N k {±} n sp(a k fin) the finite part of Σ It is clear by equation () that π is a subset of Σ Furthermore, it was shown in [] that σ π holds It was then conjectured in [] that clos(σ ) = clos(π ) = Σ and that Σ is a simply connected set which is the closure of its interior and which has a fractal boundary That the first equality holds was recently shown in [8] whereas the second equality and the other assertions remain open Some light was shed on this question in [3], where it was proved that D clos(π ) Combined with the result of [8] and equation (), this implies D clos(σ ) = clos(π ) Σ () The result of Chandler-Wilde and Davies in [3] is based on the observation that π (and also Σ) is invariant under taking square roots, ie λ π = λ π In this paper we extend this result to an infinite number of symmetries We prove that for p S, where S is an infinite set of polynomials made precise below, the following holds: p(λ) π = λ π In other words, π (and hence also clos(π )) is invariant under taking roots of polynomials p S This implies that () can be extended to for all p S Thus p (D) p (clos(σ )) = p (clos(π )) clos(π ) Σ p (D) clos(σ ) = clos(π ) Σ (3) p S This improvement in comparison with () is significant as Figure shows Clearly, this construction can also be iterated, ie and hence p n (D) p n (clos(σ )) = p n (clos(π )) clos(π ) Σ n N p n (D) n N p n (clos(σ )) = n N p n (clos(π )) clos(π ) Σ This implies that Σ contains an infinite sequence of (presumably filled) Julia sets (see Remark below), eg the set indicated in Figure 3 These two approaches can also be combined as follows Let T be the closure of S with respect to composition, ie T = {q : C C : q = p p n for p,, p n S, n N} Then q (D) clos(σ ) = clos(π ) Σ q T 4
5 5 5 5 Figure : The left-hand side of (3) (polynomials in S up to order 5 to be precise), the upper bound to Σ computed in [7] and the unit circle as a reference In this way one can construct even more Julia sets that are contained in Σ This richness of symmetries might be a part of an explanation of the seemingly fractal boundary of Σ Surely, this observation needs further investigation Note that in [3] it was also shown that λ Σ = ± λ Σ Here we only have the (possibly weaker) statement λ Σ = p ({λ}) Σ 6= for p S In addition to the polynomial symmetries mentioned above, there are also the following symmetries (see [, Lemma 34]): λ π = iλ, λ π and λ Σ = iλ, λ Σ We start with some well-known preparatory results and end up with the two main theorems of this paper A short list of polynomials in S and selected pictures of new subsets of Σ are provided at the end 5
Figure 3: The filled Julia set corresponding to p(λ) = λ 3 λ Symmetries Let m N and k {±} m Then we denote the corresponding m-periodic operator k m k L(l (Z)) k m k by A k per By Proposition, we can use the symbol a k to compute the spectrum of A k per In our case, due to tridiagonality, the formula can be simplified as follows Lemma 3 Let m N, k {±} m and let a k denote the symbol of A k per Then the only ϕ- dependent term in the characteristic polynomial of a k (ϕ) is the term of order zero More precisely, 6
there exists a polynomial p k : C C of degree m such that det(a k (ϕ) λi m ) = ( ) m p k (λ) e iϕ m j= k j e iϕ (4) for all ϕ [, π) The polynomial p k is monic and given by λ λ p k (λ) = ( ) m det k k m det k λ λ k m Furthermore, p k is an even (odd) function if m is even (odd) Proof The symbol a k is given by k m e iϕ a k (ϕ) = k e iϕ k m Using Laplace s formula, we get λ det(a k (ϕ) λi m ) = λ det k k m λ λ k det k 3 + ( )m k m e iϕ k m λ λ + ( ) m e iϕ + ( )m k m e iϕ det λ λ = det k k m det λ k m m + ( ) m e iϕ k j + ( ) m e iϕ j= k m m j= k j k λ k m k λ k m 7
That p k is an even (odd) function if m is even (odd) follows easily by induction over m The most important part of Lemma 3 is that there are no mixed terms of λ and ϕ in equation (4) This leads to the fact (see Corollary 5 below) that the spectrum of every periodic operator A k per can be written as the preimage of the interval [, ] under some polynomial p k In this way the various parts of π are connected We will make great use of this fact in Theorem But first observe the following The term in (4) involving ϕ can be simplified to cos(ϕ) or i sin(ϕ) depending on the product m k j To avoid unnecessary paperwork, we give the following definition j= Definition 4 Let m N and k {±} m Then we call k even if m k j = and odd if m k j = Note that we can always assume that a periodic operator A k per has an even period k This is because we can always double the period, ie take (k,, k m, k,, k m ) {±} m instead of (k,, k m ) {±} m as mentioned in the introduction Lemma 3 has the following important corollary Corollary 5 Let m N and k {±} m Then sp(a k per) = p k ([, ]) if k is even and sp(ak per) = p (i[, ]) if k is odd k Proof Let k be even first By Proposition and Lemma 3 we have j= sp(a k per) = {λ C : det(a(ϕ) λi m ) = for some ϕ [, π)} = {λ C : p k (λ) cos(ϕ) = for some ϕ [, π)} = {λ C : p k (λ) [, ]} = p k ([, ]) If k is odd, just replace cos(ϕ) by i sin(ϕ) The next proposition not only shows that sp(p k (A k per)) = [, ] (i[, ]) for all even (odd) k {±} m but also that p k (A k per) has a very simple form Proposition 6 Let m N, k {±} m and let p k be the corresponding polynomial given by Lemma 3 If k is even, then p k (a k (ϕ)) = cos(ϕ)i m and p k (A k per) is the Laurent operator with on its m-th sub- and superdiagonal (and everywhere else): p k (A k per) = (5) j= 8
If k is odd, then p k (a k (ϕ)) = i sin(ϕ)i m and p k (A k per) is the Laurent operator with on its m-th subdiagonal and on its m-th superdiagonal (and everywhere else): p k (A k per) = (6) Proof In both cases the first part follows immediately from Lemma 3 and the theorem of Cayley- Hamilton For the second part observe that cos(ϕ)i m = (e iϕ + e iϕ )I m is the symbol of (5) and i sin(ϕ)i m = ( e iϕ + e iϕ )I m is the symbol of (6) The assertion thus follows by Proposition This simple observation now enables us to prove the first of our main results Theorem 7 Let m, k := (k,, k m,, ) {±} m and ˆk := (k,, k m,, ) {±} m Furthermore, let b {±} Z and A b := b, b where the box indicates A b, If the corresponding polynomials p k and pˆk are equal, then there exist c {±} Z and B L(l (Z)) such that p k (A c ) = B A b, where we consider the following decomposition of the Hilbert space l (Z): l (Z) = l (Z \ mz) l (mz) = l (Z) l (Z) In particular, sp(a b ) sp(p k (A c )) Proof Let c {±} Z be the sequence defined by: c =, c =, c rm+j = k j for j {,, m }, r Z, 9
m j= c rm j = b r for r Z, c rm+ = c rm for r Z Note that A c is very similar to A k per and Aˆk per The difference is that, depending on the sequence b, the entries c rm and c rm+ are swapped for some r Z For m = this is exactly the same construction as in [3, Lemma 5] First we will prove the following claim by induction: Claim : ((A c ) s ) i,i s+j only depends on the coefficients c i s+j+,, c i+j for s N, i Z and j {,, s } Furthermore, ((A c ) s ) i,i+s = for all s N, i Z and all other entries are For s = we have (A c ) i,i +j = c i+j for j =, i Z, and (A c ) i,i+ = for i Z All other entries are So assume that the claim holds for s Then ((A c ) s ) i,i s+j = (A c ) i,i+ ((A c ) s ) i+,i s+j + (A c ) i,i ((A c ) s ) i,i s+j = ((A c ) s ) i+,i s+j + c i ((A c ) s ) i,i s+j (7) = f(c i s+j+,, c i+j ) + c i g(c i s+j+,, c i+j ), where f and g are some polynomials Observe that c i is contained in {c i s+j+,, c i+j } for all i Z, j {,, s } Thus ((A c ) s ) i,i s+j only depends on the coefficients c i s+j+,, c i+j for all i Z, j {,, s } Plugging j = s into (7) yields ((A c ) s ) i,i+s = for all i Z by the same induction argument Plugging in j Z + and j Z \ {,, s} shows that all other entries are This finishes the proof of the claim Using Claim with s = m and i = rm for r Z, we get that ((A c ) m ) rm,(r )m+j only depends on the coefficients c (r )m+j+,, c rm+j for all j {,, m }, ((A c ) m ) rm,rm+m = and all other entries in row rm are Similarly, ((A c ) m ) (r )m+j,rm only depends on the coeffcients c (r )m+j+,, c rm+j for all j {,, m }, ((A c ) m ) rm+m,rm = and all other entries in column rm are Moreover, Claim also implies that the same is true for p k (A c ) because p k is an even/odd monic polynomial of degree m by Lemma 3 Claim : ((A c ) m ) i,i m = m j= This again follows easily by induction: c i j for all i Z ((A c ) m ) i,i m = (A c ) i,i ((A c ) m ) i,i m + (A c ) i,i+ ((A c ) m ) i+,i m = c i ((A c ) m ) i,i m + = c i c i (m ) Since ((A c ) s ) i,i m = for all s < m, also (p k (A c )) i,i m = m c, it thus follows (p k (A c )) rm,(r )m = m j= j= c rm j = b r c i j for all i Z By definition of
for all r Z Furthermore, (p k (A c )) rm,(r )m+j only depends on the coefficients c (r )m+j+,, c rm+j for j {,, m } In particular, these numbers all depend on c rm and c rm+ but not on c (r )m, c (r )m+, c (r+)m or c (r+)m+ This implies or (using p k = pˆk) (p k (A c )) rm,(r )m+j = (p k (A k per)) rm,(r )m+j = (if c rm =, c rm+ = ) (p k (A c )) rm,(r )m+j = (pˆk(aˆk per )) rm,(r )m+j = (if c rm =, c rm+ = ) for j {,, m } by Proposition 6 In other words, the entries (p k (A c )) rm,(r )m+j (j {,, m }) can not know whether we swapped some of the entries c lm and c lm+ (l Z) or not Thus they have to remain zero Similarly, the entries (p k (A c )) (r )m+j,rm (j {,, m }) remain Therefore p k (A c ) looks like this (where means some unimportant entries ): p k (A c ) = b r b r+ }{{} m columns Decomposing our Hilbert space l (Z) = l (Z \ mz) l (mz) = l (Z) l (Z), we get the following decomposition of p k (A c ): p k (A c ) = B A b for some B L(l (Z)) In particular, sp(a b ) sp(p k (A c )) By construction of the sequence c, we also have the following important corollary for periodic operators Corollary 8 Under the same assumptions as in Theorem 7, we have that if b {±} Z is an n-periodic sequence with even period (b,, b n ), then c, as defined in the proof of Theorem 7, is an nm-periodic sequence with even period (c,, c nm ) B L(l (Z)) then is a periodic operator, too Furthermore, if we denote the symbols of A b, A c and B by a b, a c and a B, then also p k (a c ) can be decomposed as p k (a c ) = a B a b In particular, sp(a b (ϕ)) sp(p k (a c (ϕ))) for every ϕ [, π)
Proof The first part follows by contruction of c It remains to prove that the symbol p k (a c ) can be decomposed in a similar way By Proposition, A c is unitarily equivalent to the multiplication operator M a c L(L ([, π), C nm )) Let us denote this equivalence by F nm, ie A c = F nmm a cf nm It follows p k (A c ) = p k (F nmm a cf nm ) = F nmp k (M a c)f nm = F nmm pk (a c )F nm Furthermore, A b is unitarily equivalent to the multiplication operator M a b L(L ([, π), C m )) and B is unitarily equivalent to the multiplication operator M a B L(L ([, π), C (n )m )) Let us denote these equivalences by F m and F (n )m Furthermore, let us denote the decomposition l (Z) = l (Z\mZ) l (mz) by U and the decomposition L ([, π), C nm ) = L ([, π), C (n )m ) L ([, π), C m ) (in the obvious way) by V It is not hard to see that holds Thus M pk (a c ) = F nm p k (A c )F nm (F (n )m F m )UF nm = V = F nm U (B A b )UF nm = F nm U (F (n )m M a BF (n )m F mm a bf m )UF nm = V (M a B M a b)v This implies p k (a c ) = a B a b (as functions of ϕ) Definition 9 We define S := {p k : p k is a polynomial such that the assumptions of Theorem 7 are satisfied} as our set of symmetries of π In the case of periodic operators we can prove the following stronger version of Theorem 7 that justifies the definition of S Theorem Let n N, let b {±} Z be an n-periodic sequence and let k {±} m be such that p := p k S Moreover, let c {±} Z be the sequence as constructed in the proof of Theorem 7 Then the following assertion holds: p(λ) sp(a b ) λ sp(a c ) Proof Wlog we can assume that (b,, b n ) is even By Corollary 8, (c,, c nm ) is even, too Let q and r be the polynomials given by Lemma 3 corresponding to b and c Also denote by a b and a c the symbols of A b and A c Fix some ϕ [, π) and let µ C be an eigenvalue of a b (ϕ) Again by Corollary 8, there exists some λ sp(a c (ϕ)) such that p(λ) = µ Since q(µ) = cos(ϕ) by Proposition 6, we have q(p(λ)) = cos(ϕ) On the other hand, also r(λ) = cos(ϕ) by Proposition 6 Thus (q p)(λ) = cos(ϕ) = r(λ) (cf Figure 4) Since both q p and r are polynomials and the above argument is valid for every ϕ [, π), we conclude that q p and r are equal It follows sp(a c ) = r ([, ]) = p (q ([, ])) = p (sp(a b )) by Corollary 5
µ sp(a b ) p q sp(a c ) λ r - cos ϕ π Figure 4: Schematic picture of the maps involved in the proof of Theorem Theorem combined with the result D clos(π ) from [3] implies the following corollary Corollary Let p S Then p (π ) π and hence p (D) p (clos(π )) clos(π ) Recall that k = (k,, k m ) {±} m generates a polynomial p = p k S if k m k m, so that ˆk = (k,, k m, k m, k m ) k but still p k = pˆk Thus it is immediate that all k of the form (,,,, ) and (,,,, ) generate a polynomial p k S Indeed, if k = (,,,, ) and ˆk = (,,,, ), then A k per and Aˆk per are unitarily equivalent by a simple shift This implies that S contains a countable number of polynomials However, there are also a lot more than these trivial examples as the following table shows We conjecture that there are approximately m polynomials of degree m in S 3
No k p k (λ) (, ) λ 3 (,, ) λ 3 λ 3 (,, ) λ 3 + λ 4 (,,, ) λ 4 λ 4 (,,, ) λ 4 + λ 5 (,,,, ) λ 5 3λ 3 + λ 5 (,,,, ) λ 5 λ 3 + λ 53 (,,,, ) λ 5 + λ 3 + λ 54 (,,,, ) λ 5 + 3λ 3 + λ 6 (,,,,, ) λ 6 4λ 4 + 3λ 6 (,,,,, ) λ 6 λ 63 (,,,,, ) λ 6 + 4λ 4 + 3λ 7 (,,,,,, ) λ 7 5λ 5 + 6λ 3 λ 7 (,,,,,, ) λ 7 3λ 5 + λ 3 + λ 73 (,,,,,, ) λ 7 λ 5 + λ 3 λ 74 (,,,,,, ) λ 7 + λ 5 λ 3 + λ 75 (,,,,,, ) λ 7 λ 5 λ 3 λ 76 (,,,,,, ) λ 7 + λ 5 + λ 3 + λ 77 (,,,,,, ) λ 7 + 3λ 5 + λ 3 λ 78 (,,,,,, ) λ 7 + 5λ 5 + 6λ 3 + λ Table : Short list of elements in S 4
Remark As mentioned in the introduction, we can iterate Corollary to get U := n N p n (D) n N p n (clos(π )) clos(π ) for every p S In other words, z U if and only if p n (z) for some n N Thus there is clearly a connection to the filled Julia set corresponding to p which is given by J f (p) := {z C : (p n (z)) n N is bounded} (see [, Lemma 7]) Indeed, the boundary J(p) := J f (p) (which is usually just called the Julia set corresponding to p) is contained in the closure of p n (z) for every z C except for at most one point (see [5, Corollary 48(a)]) Hence J(p) clos(u) clos(π ) Considering the pictures (viii) and (ix) below, it seems natural to conjecture that even the filled Julia set J f (p) is contained in clos(u) We conclude with some pictures of subsets of Σ The red unit circle serves as a reference n N 5 5 5 5 5 5 5 5 (i) p k (D) for k = (,, ) (ii) p (D) for k = (,,,, ) k 5 5 5 5 5 5 5 5 (iii) p k (D) for k = (,,,,,, ) (iv) p k (D) for k = (,,,,,,, ) 5
5 5 5 5 5 5 5 5 (v) S (vi) p (D) S p (D) p S deg(p)= p S deg(p) 7 5 5 5 5 (vii) S p S deg(p) 5 6 p (D)
5 5 5 5 5 5 5 5 (viii) p 9 k (D) for k = (,, ) (ix) p 6 k (D) for k = (,,,, ) (x) p 5 (D) p S deg(p) 8 7
References [] SN Chandler-Wilde, R Chonchaiya and M Lindner: Eigenvalue Problem meets Sierpinski Triangle: Computing the Spectrum of a Non-Self-Adjoint Random Operator, Operators and Matrices, 5 (), 633-648 [] SN Chandler-Wilde, R Chonchaiya and M Lindner: On the Spectra and Pseudospectra of a Class of non-self-adjoint Random Matrices and Operators, Operators and Matrices, 7 (3), 739-775 [3] SN Chandler-Wilde and EB Davies: Spectrum of a Feinberg-Zee Random Hopping Matrix, Journal of Spectral Theory, (), 47-79 [4] EB Davies: Linear Operators and their Spectra, Cambridge Studies in Advanced Mathematics 6, Cambridge University Press, Cambridge, 7 [5] K Falconer: Fractal Goemetry: Mathematical Foundations and Applications, nd Ed, John Wiley, Chichester, 3 [6] J Feinberg and A Zee: Spectral Curves of Non-Hermitean Hamiltonians, Nucl Phys B 55 (999), 599-63 [7] R Hagger: On the Spectrum and Numerical Range of Tridiagonal Random Operators, Preprint at arxiv: 475486 [8] R Hagger: The Eigenvalues of Tridiagonal Sign Matrices are Dense in the Spectra of Periodic Tridiagonal Sign Operators, Preprint at arxiv: 474 [9] DE Holz, H Orland and A Zee: On the Remarkable Spectrum of a Non-Hermitian Random Matrix Model, J Phys A Math Gen 36 (3), 3385-34 [] J Milnor: Dynamics in one Complex Variable Introductory Lectures, Vieweg & Sohn, Braunschweig, 999 Author s address: Raffael Hagger Institute of Mathematics Hamburg University of Technology Schwarzenbergstr 95 E D-73 Hamburg GERMANY raffaelhagger@tuhhde 8