arxiv: v1 [math.co] 15 Jan 2019

Transcription

1 arxiv: v1 [math.co] 15 Jan 2019 DEFORMED FOURIER MATRICES IN BISTOCHASTIC FORM TEO BANICA Abstract. Many interesting examples of complex Hadamard matrices H M N (C) can be put, up to the standard equivalence relation for such matrices, in bistochastic form. We discuss here this phenomenon, with a number of computations for some of the most basic such matrices, namely the deformations of the Fourier matrix F N M N (C). Introduction A complex Hadamard matrix is a square matrix H M N (C) whose entries are on the unit circle, H ij = 1, and whose rows are orthogonal with respect to the scalar product of C N. The basic example is the Fourier matrix, F N = (w ij ) ij with w = e 2πi/N : w w 2... w N 1 F N = w N 1 w 2(N 1)... w (N 1)2 More generally, the Fourier matrix F G =< i,j > ij with i G,j Ĝ of any finite abelian group G is a complex Hadamard matrix. As basic examples here, for the cyclic group G = Z N we obtain the above standard Fourier matrix F N, while for group G = Z n 2, having order N = 2 n, we obtain the n-th Walsh matrix W N = F2 n. In general, the complex Hadamard matrices are best thought of as being generalized Fourier matrices, of somewhat exotic type. These matrices can successfully replace the Fourier matrices, in a number of key situations, and due to this fact, they appear in a wide array of questions, in mathematics and physics. See [19]. Our starting observation here regards F 2. By multiplying the first row by i, and the second column by i, we can put this matrix in a better-looking form, as follows: ( ) ( ) 1 1 i 1 F 2 = = F i Mathematics Subject Classification. 15B34 (05B20). Key words and phrases. Complex Hadamard matrix, Fourier matrix. 1

2 2 TEO BANICA Indeed, the matrix on the right F 2 is Hadamard as well, and looks more balanced than F 2. To be more precise, this matrix has the extra property of being bistochastic, in the sense that the sums on all rows and all columns are equal (to 1+i). Many other interesting Hadamard matrices can be put in bistochastic form. This is true for instance for the Walsh matrices, simply by taking tensor powers: ( ) n ( ) n 1 1 i 1 W N = = W i N In fact, any Fourier matrix F N can be put in bistochastic form, F N F N, and the same remains true for the generalized Fourier matrices F G. Inaddition, several interesting deformations of the Fourier matrices F G can be put in bistochastic form as well. Our aim here is to discuss this phenomenon. We have three motivations, as follows: (1) Classical Hadamard matrices. Certain real Hadamard matrices H M N (±1) can be put in bistochastic form, over the field of reals R. However, this is quite restrictive, because in this case the order must be a square, N = M 2. The point now is that, as shown by the above equivalence F 2 F 2, when allowing complex scalars the situation drastically changes. So, it would be very interesting to know which real matrices H M N (±1) can be put in bistochastic form, over the field of complex numbers C. This would shed some new light on these matrices. (2) Almost Hadamard conjecture. This conjecture, recently formulated in [6], following some previous work from [7], states that the rescaled versions U = H/ N of the complex Hadamard matrices H M N (C) are precisely the local maximizers of the 1-norm U 1 = ij U ij on the unitary group U N. This might seem a bit technical, but as explained in [6], [7], having such a result would be something very useful. And the point is that the few available verifications of the conjecture, from [6], are in the bistochastic case, and this suggests a relation with [13]. (3) Quantum algebra questions. Associated to any Hadamardmatrix H M N (C) is a quantum permutation group G S + N, and the computation of the representation theory invariants of G is a key question [4], [5]. These invariants are believed to be related to the geometry of the Hadamard manifold around H, and to the notion of defect [2], [15], [20], and also to the notion of glow from [3], inspired from the Gale-Berlekamp game [11], [16]. And the point is that the support of the glow is maximal precisely when H can be put in bistochastic form. Summarizing, understanding which complex Hadamard matrices H M N (C) can be put in bistochastic form appears to be an interesting question, related to many things. We will discuss here this problem, with some results concerning the most basic examples of such matrices, namely the Fourier matrices and their deformations.

3 FOURIER MATRICES 3 The paper is organized as follows: 1 contains some preliminaries, and a straightforward studyofthecomplexhadamardmatricesofordern 5,usingelementarymethods, in2-3 we use more advanced technology, namely Björck s cyclic roots, in order to investigate the general Fourier matrices and some of their deformations, 4 contains some supplementary preliminary material, and in 5-6 we state and prove our main results, regarding the general deformed tensor products, and the tensor powers of Fourier matrices. 1. Hadamard matrices One of the most interesting objects in combinatorics are the matrices H M N (±1) whose rows are pairwise orthogonal. These are called Hadamard matrices [17], [18]. We are interested in what follows in the complex generalization of such matrices: Definition 1.1. A complex Hadamard matrix is a square matrix H M N (C) whose entries are on the unit circle, H ij T, and whose rows are pairwise orthogonal. Here, and in what follows, the scalar product is the usual one on C N. Observe that orthogonality condition tells us that we must have U N U N. Since the transpose of a unitary matrix is unitary too, the columns of H must be orthogonal as well. The basic example is provided by the Fourier coupling G Ĝ T of a finite abelian group G, regarded as a square matrix, F G M G Ĝ (T) M N(T), with N = G. In the cyclic group case G = Z N this matrix is denoted F N, and is given by F N = (w ij ) ij with w = e 2πi/N. In general, with G = Z N1... Z Ns, we have F G = F N1... F Ns. There are many other interesting examples of such matrices. See [19]. We are interested in such matrices, up to the following equivalence relation: Definition 1.2. Two complex Hadamard matrices H,K M N (C) are called equivalent, and we write in this case H K, if one can pass from H to K by permuting the rows and columns, and by multiplying the rows and columns by numbers in T. As a basic example of such an equivalence, which will be of key importance in what follows, consider the Fourier matrix F 2. By multiplying the first row by i, and the second column by i, we can put this matrix in a better-looking form, as follows: ( ) ( i 1 1 i In general, most of the questions regarding the complex Hadamard matrices, including their classification, can be considered up to the above equivalence relation. It was shown by Haagerup in [12] that the only Hadamard matrices at N = 2,3,4,5 are, up to equivalence, the matrices F 2,F 3,F s 4,F 5, where F s 4 with s T is a certain one-parameter deformation of F 4. The classification problem is open at N = 6. See [12], [19]. We will be mostly interested here in the Hadamard matrices which are bistochastic, or which are equivalent to a bistochastic matrix, in the following sense: )

4 4 TEO BANICA Definition 1.3. A matrix H M N (C) is called bistochastic when the N sums on its rows and the N sums on its columns are all equal, to a certain number λ C. In other words, H is bistochastic when Hξ = H t ξ = λξ, where ξ C N is the all-one vector. Equivalently, we must have HJ = JH = λj, where J M N (C) is the all-one matrix. We should mention that these matrices go under different names, including also doubly stochastic, in Sinkhorn normal form, and regular. See [13], [14], [18]. Let us first discuss some examples. As a first, elementary result, we have: Proposition 1.4. The following are bistochastic complex Hadamard matrices, ( ) i 1 F 2 =, F 1 i 3 = w ν 1 ν 1 w 1, F 4 = ν 1 ν 1 1 ν 1 ν 1 1 w ν 1 ν 1 where w = e 2πi/3,ν = e πi/4, equivalent to the Fourier matrices F 2,F 3,F 4. Proof. The orthogonality between rows being clear, we have here complex Hadamard matrices. These matrices being circulant, they are bistochastic as well. The fact that we have an equivalence F 2 F 2 is something that we already know. At N = 3 now, the equivalence F 3 F 3 can be constructed as follows: w w w 1 w 1 w w 1 1 w 2 w w w As for the case N = 4, here the equivalence F 4 F 4 can be constructed as follows, where we use the logarithmic notation [k] s = e 2πki/s, with respect to s = 8: Finally, observe that the matrices in the statement are all circulant and symmetric, which is stronger than our bistochasticity claim. We will be back later to this. Also, we will explain as well later the reasons for denoting the last matrix F 4, instead of F 4. In order to further build on these examples, let us introduce: Definition 1.5. A matrix H M N (C), with matrix indices i,j {0,1,...,N 1}, taken modulo N, is called circulant when we have H ij = γ j i, for some γ C N. 8 8

5 FOURIER MATRICES 5 In the complex Hadamard matrix case, we already have as examples the above matrices F 2,F 3,F 4. Here is another example, which is fundamental as well: K 4 = Observe that this matrix is, once again, both circulant and symmetric. In fact, according to the Circulant Hadamard Conjecture (CHC) this matrix is, up to equivalence, the only circulant Hadamard matrix H M N (±1), regardless of the value of N N. In relation with our Fourier matrix considerations, K 4 is equivalent to the Fourier matrix F G = F 2 F 2 of the Klein group G = Z 2 Z 2, as shown by: More generally, we have the following construction of circulant and symmetric Hadamard matrices at N = 4, which involves an extra parameter q T: Proposition 1.6. The following circulant and symmetric matrix is Hadamard, 1 q 1 q K q 4 = q 1 q 1 1 q 1 q q 1 q 1 for any q T. At q = 1,e πi/4 recover respectively the matrices K 4,F 4. Proof. The rows of the above matrix are pairwise orthogonal for any q C, and so at q T we obtain a complex Hadamard matrix. The last assertion is clear. We can now formulate our main result is this section, as follows: Theorem 1.7. Any complex Hadamard matrix of order N = 2,3,4,5 can be put, up to equivalence, in circulant and symmetric form. Proof. As already mentioned, by a fundamental result of Haagerup [12], the only complex Hadamard matrices at N = 2,3,4,5 are, up to equivalence, the matrices F 2,F 3,F p 4,F 5, where F p 4 with p T is a certain one-parameter deformation of F 4. At N = 2,3 the problem is solved by Proposition 1.4 above.

6 6 TEO BANICA AtN = 4now, ourclaimisthatwehavek q 4 Fs 4, withs = q 2. Indeed, bymultiplying the rows of K q 4, and then the columns, by suitable scalars, we have: 1 q 1 q 1 q 1 q K q 4 = q 1 q 1 1 q 1 q 1 q 1 q 1 q 1 q 1 s 1 s q 1 q 1 1 q 1 q 1 s 1 s On the other hand, by permuting the second and third rows of F4, s we obtain: F4 s = s 1 s 1 s 1 s s 1 s 1 s 1 s Thus these matrices are equivalent, and the result follows from Proposition 1.6. At N = 5, the matrix that we are looking for is as follows, with w = e 2πi/5 : w 2 1 w 4 w 4 1 F 5 = 1 w 2 1 w 4 w 4 w 4 1 w 2 1 w 4 w 4 w 4 1 w w 4 w 4 1 w 2 It is indeed clear that this matrix is circulant, symmetric, and complex Hadamard, and the fact that we have F 5 F 5 follows either directly, or by using [12]. 2. Cyclic roots We have seen in the previous section that the Fourier matrices of order N 5 can be put in symmetric and circulant form. We discuss here the general case, by putting any Fourier matrix F N in circulant and symmetric form. In addition, we will present a unification of this result with Proposition 1.6, by using a construction from [1]. We use Björck s cyclic root formalism [8], which is as follows: Proposition 2.1. Assume that H M N (T) is circulant, H ij = γ j i. Then H is Hadamard if and only if the vector (z 0,z 1,...,z N 1 ) given by z i = γ i /γ i 1 satisfies: z 0 +z z N 1 = 0 z 0 z 1 +z 1 z z N 1 z 0 = 0... z 0 z 1...z N z N 1 z 0...z N 3 = 0 z 0 z 1...z N 1 = 1 If so is the case, we say that z = (z 0,...,z N 1 ) is a cyclic N-root.

7 FOURIER MATRICES 7 Proof. This follows indeed from a direct computation, and we refer here to [8]. Observe that, up to a global multiplication by a scalar w T, the first row vector γ = (γ 0,...,γ N 1 ) of the matrix H M N (T) constructed above is as follows: γ = (z 0,z 0 z 1,z 0 z 1 z 2,...,z 0 z 1...z N 1 ) Now back to the Fourier matrices, we have the following result: Proposition 2.2. Given N N, set ν = e πi/n and q = ν N 1,w = ν 2. Then (q,qw,qw 2,...,qw N 1 ) is a cyclic N-root, and the corresponding complex Hadamard matrix F N is circulant and symmetric, and equivalent to the Fourier matrix F N. Proof. Given q,w T, let us find out when (q,qw,qw 2,...,qw N 1 ) is a cyclic root: (1) In order for the = 0 equations in Proposition 2.1 to be satisfied, the value of q is irrelevant, and w must be a primitive N-root of unity. (2) As for the = 1 equation in Proposition 2.1, this states in our case that we must have q N w N(N 1) 2 = 1, and so that we must have q N = ( 1) N 1. We conclude that with the values of q,w T in the statement, we have indeed a cyclic N-root. Now construct H ij = γ j i as in Proposition 2.1. We have: γ k = γ k, k q k+1 w k(k+1) 2 = q k+1 w k(k 1) 2, k q 2k w k = 1, k q 2 = w 1 But this latter condition holds indeed, because we have q 2 = ν 2N 2 = ν 2 = w 1. We conclude that our circulant matrix H is symmetric as well, as claimed. It remains to construct an equivalence H F N. In order to do this, observe that, due to our conventions q = ν N 1,w = ν 2, the first row vector of H is given by: γ k = q k+1 w k(k+1) 2 = ν (N 1)(k+1) ν k(k+1) = ν (N+k 1)(k+1) Thus, the entries of H are given by the following formula: H i,j = H 0,i+j = ν (N+i+j 1)(i+j+1) = ν i2 +j 2 +2ij+Ni+Nj+N 1 = ν N 1 ν i2 +Ni ν j2 +Nj ν 2ij With this formula in hand, we can now finish. Indeed, the matrix H = (H ij ) is equivalent to the matrix H = (H i,j ). Now regarding H, observe that in the above

8 8 TEO BANICA formula, thefactorsν N 1, ν i2 +Ni, ν j2 +Nj correspondrespectively toaglobalmultiplication byascalar, andtorowandcolumnmultiplications byscalars. ThusH isequivalent tothe matrix H obtained by deleting these factors. But this latter matrix, given by H ij = ν2ij with ν = e πi/n, is precisely the Fourier matrix F N, and we are done. As an illustration, let us work out the cases N = 2,3,4,5. We have here: Proposition 2.3. The matrices F N are as follows: (1) At N = 2,3 we obtain the old matrices F 2,F 3, from section 1. (2) At N = 4 we obtain the following matrix, with ν = e πi/4 : ν 3 1 ν 7 1 F 4 = 1 ν 3 1 ν 7 ν 7 1 ν ν 7 1 ν 3 (3) At N = 5 we obtain the old matrix F 5, from section 1. Proof. With notations from Proposition 2.2, the proof goes as follows: (1) At N = 2 we have ν = i,q = i,w = 1, so the cyclic root is (i, i), the first row vector is (i,1), and we obtain indeed the old matrix F 2, from section 1. At N = 3 we have ν = e πi/3 and q = w = ν 2 = e 2πi/3, the cyclic root is (w,w 2,1), the first row vector is (w,1,1), and we obtain indeed the old matrix F 3, from section 1. (2) At N = 4 we have ν = e πi/4 and q = ν 3,w = ν 2, the cyclic root is (ν 3,ν 5,ν 7,ν), the first row vector is (ν 3,1,ν 7,1), and we obtain the matrix in the statement. (3) At N = 5 we have ν = e πi/5 and q = ν 4 = w 2, with w = ν 2 = e 2πi/5, the cyclic root is therefore (w 2,w 3,w 4,1,w), the first row vector is (w 2,1,w 4,w 4,1), and we obtain in this way the old matrix F 5 from section 1, as claimed. Regarding the above matrix F 4, observe that this is equivalent to the matrix F 4 from section 1, with the equivalence F 4 F 4 being obtained by multiplying everything by ν = e πi/4. Whileboththesematricesarecirculantandsymmetric, andofcourseequivalent tof 4, oneofthem, namelyf 4, is better thantheother, becausethecorresponding cyclic root comes from a progression. This is the reason for our notations F 4,F 4. We recall now that the tensor product of the usual matrices is given by (H K) ia,jb = H ij K ab. In matrix form, by using the lexicographic order on the double indices: H 11 K... H 1M K H K =. H M1 K... H MM K We have the following elementary result, regarding such products:

9 FOURIER MATRICES 9 Proposition 2.4. If H M M (C) and K M N (C) are both complex Hadamard matrices, then so is their tensor product, constructed as follows: H K M MN (C), (H K) ia,jb = H ij K ab In addition, if both H,K are symmetric or bistochastic, then so is H K. Proof. All these assertions are well-known, and follow from definitions. We can now discuss the case of the generalized Fourier matrices F G, with G being a finite abelian group. These decompose as tensor products, and we have: Theorem 2.5. The generalized Fourier matrices F G, associated to the finite abelian groups G, can be put in symmetric and bistochastic form. Proof. We recall that F G is by definition the Fourier coupling G Ĝ T of the abelian group G, regarded as a square matrix, F G M (T) M G Ĝ N(T), with N = G. It is standard to check that with G = H K we have F G = F H F K, and also that in the cyclic group case, G = Z N, the Fourier matrix is F G = F N. Thus, when decomposing our group as G = Z N1... Z Ns, we obtain an identification as follows: F G = F N1... F Ns With this formula in hand, the result follows from the observations in Proposition 2.4 above, with input coming from Proposition 2.2 above. Observe that the property of F N of being circulant does not extend via tensor products. In fact, being circulant means something with respect to the cyclic group G = Z N, and in the case of the generalized Fourier matrices F G, the underlying group is G itself. 3. Backelin matrices We have the following alternative generalization of Proposition 2.2, coming this time from Backelin s work in [1], and remaining in the circulant and symmetric setting: Theorem 3.1. Let M N, and set w = e 2πi/N. We have a cyclic root as follows, (q 1,...,q }{{ M,q } 1 w,...,q M w,...,q }{{} 1 w N 1,...,q M w N 1 ) }{{} M M M provided that q 1,...,q M T satisfy (q 1...q M ) N = ( 1) M(N 1). Moreover, assuming q 1 q 2 = 1, q 3 q M = q 4 q M 1 =... = w which imply (q 1...q M ) N = ( 1) M(N 1), the Hadamard matrix is symmetric.

10 10 TEO BANICA Proof. Let us first check the = 0 equations for a cyclic root. Given arbitrary numbers q 1,...,q M T, if we denote by (z i ) the vector in the statement, we have: ( ) q1...q z i+1...z i+k = K +q 2...q K q M K+1...q M +q M K+2...q M q 1 w +...+q M q 1...q K 1 w K 1 i (1+w K +w 2K +...+w (N 1)K ) Now since the sum on the right vanishes, the = 0 conditions are satisfied. Regarding now the = 1 condition, the total product of the numbers z i is given by: z i = (q 1...q M ) N (1 w w 2...w N 1 ) M = (q 1...q M ) N w MN(N 1) 2 i By using w = e 2πi/N we obtain that the coefficient on the right is: w MN(N 1) 2 = e 2πi N MN(N 1) 2 = e πim(n 1) = ( 1) M(N 1) Thus, if (q 1...q M ) N = ( 1) M(N 1), we obtain a cyclic root, as stated. See [1], [10]. The corresponding first row vector can be written as follows: V = q w M 1 w 2 1,q 1 q 2,...,q 1...q }{{ M,...,,...,, w,1 } q 2...q M q M 1 q M q M M }{{} M Thus, the corresponding circulant complex Hadamard matrix is as follows: q 1 q 1 q 2 q 1 q 2 q 3 q 1 q 2 q 3 q 4 q 1 q 2 q 3 q 4 q q 1 q 1 q 2 q 1 q 2 q 3 q 1 q 2 q 3 q 4... w q M 1 q 1 q 1 q 2 q 1 q 2 q 3... H = w 2 w q M 1 q M q M 1 q 1 q 1 q 2... w 3 w 2 w q M 2 q M 1 q M q M 1 q M q M 1 q We are therefore led to the symmetry conditions in the statement, and we are done. Observe that the symmetry conditions q 1 q 2 = 1,q 3 q M = q 4 q M 1 =... = w in the above statement imply in fact the cyclic root condition, namely: (q 1...q M ) N = ( 1) M(N 1) It is quite unclear, however, on how to further simplify the formulation of these symmetry conditions, as to reach to an overall improvement of the statement of Theorem 3.1. We believe that this should be possible, but we have no further results here. Observe also that Theorem 3.1 generalizes both Proposition 2.2, and the construction in Proposition 1.6. Thus, we have here a full generalization of Theorem 1.7. The question

11 FOURIER MATRICES 11 of further generalizing this result, as to cover for instance the Fourier matrices F G = F N1... F Ns of the arbitrary abelian groups G = Z N1... Z Ns, remains open. We will be back to the Fourier matrices and their deformations, with some further results on the subject, regarding more general deformations, in sections 5-6 below. 4. Bistochastic matrices We have seen so far that the Fourier matrices and some of their deformations can be put in bistochastic form, with the bistochasticity property actually coming from some much stronger properties, such as the fact of being circulant, and symmetric too. In general, we cannot expect such stronger properties to appear. So, let us develop us now some general theory for the general bistochastic Hadamard matrices. The question that we are mostly interested in is that of deciding, given an Hadamard matrix H M N (C), whether H can be put or not in bistochastic form. That is, we would like to know if we have an equivalence of type H H, with H being bistochastic. Let us begin with the following standard result: Proposition 4.1. For an Hadamard matrix H M N (C), the following are equivalent: (1) If H is bistochastic, with sums λ. (2) If H is row-stochastic, with sums λ, and λ 2 = N. Proof. Both the implications are elementary, as follows: (1) = (2) If we denote by H 1,...,H N T N the rows of H, we have indeed: N = < H 1,H 1 >= < H 1,H i >= H 1j Hij i i j = H 1j H ij = H 1j λ = λ λ = λ 2 j i j (2) = (1) Consider the all-one vector ξ = (1) i C N. The fact that H is rowstochastic with sums λ, respectively column-stochastic with sums λ, reads: H ij = λ, i H ij ξ j = λξ i, i Hξ = λξ j j H ij = λ, j H ij ξ i = λξ j, j H t ξ = λξ j i We must prove that the first condition implies the second one, provided that the row sum λ satisfies λ 2 = N. But this follows from the following computation: Hξ = λξ = H Hξ = λh ξ = N 2 ξ = λh ξ = N 2 ξ = λh t ξ = H t ξ = λξ Thus, we have proved both the implications, and we are done.

12 12 TEO BANICA Now back to our questions, we first have: Proposition 4.2. For a complex Hadamard matrix H M N (C), and a number λ C satisfying λ 2 = N, the following are equivalent: (1) We have H H, with H being bistochastic, with sums λ. (2) K ij = a i b j H ij is bistochastic with sums λ, for some a,b T N. (3) The equation Hb = λā has solutions a,b T N. Proof. Once again, this is an elementary result, the proof being as follows: (1) (2) Since the permutations of the rows and columns preserve the bistochasticity condition, the equivalence H H that we are looking for can be assumed to come only from multiplying the rows and columns by numbers in T. Thus, we are looking for scalars a i,b j T such that K ij = a i b j H ij is bistochastic with sums λ, as claimed. (2) (3) The row sums of the matrix K ij = a i b j H ij are given by: K ij = a i b j H ij = a i (Hb) i j j Thus K is row-stochastic with sums λ precisely when Hb = λā, and by using the equivalence in Proposition 4.1 above, we obtain the result. With these observations in hand, we can now formulate a first conceptual result, regarding the bistochastic complex Hadamard matrices, as follows: Proposition 4.3. For a complex Hadamard matrix H M N (C), the excess, E(H) = ij H ij satisfies E(H) N N, with equality if and only if H is bistochastic. Proof. In terms of the all-one vector ξ = (1) i C N, we have: E(H) = H ij = H ij ξ j ξi = (Hξ) i ξi =< Hξ,ξ > ij ij i Now by using the Cauchy-Schwarz inequality, along with the fact that U = H/ N is unitary, and hence of norm 1, we obtain, as claimed: E(H) Hξ ξ H ξ 2 = N N Regarding now the equality case, this requires the vectors Hξ, ξ to be proportional, and so our matrix H to be row-stochastic. Moreover, if we assume Hξ = λξ, the above computation gives λ 2 = N, and by Proposition 4.1, we obtain the result. In what follows we will not exactly be interested in the complex Hadamard matrices which are bistochastic, but rather in those which can be put in bistochastic form. We can summarize the above results, in terms of these latter matrices, as follows:

13 FOURIER MATRICES 13 Proposition 4.4. An Hadamard matrix H M N (C) can be put in bistochastic form when one of the following conditions is satisfied: (1) The equations Ha i = N, with i = 1,...,N, have solutions a T N. (2) The quantity E attains its maximum N N over the equivalence class of H. Proof. This follows indeed from Proposition 4.2 and Proposition 4.3 above. Observe that system of N equations in (1) has as unknowns N real variables, namely the phases of a 1,...,a N. However, this system is highly non-linear, and at this level of generality, we have no results. As for the method in (2), here a few more things can be said at the general level, and explaining this will be our purpose, in what follows. Let us discuss indeed the relation with the glow. We recall from [3] that we have: Definition 4.5. The glow of a complex Hadamard matrix H M N (C) is the complex probability measure µ given by the formula ( ϕ(x)dµ(x) = a i b j H ij )d(a,b) C T N T N ϕ valid for any continuous function ϕ : C C. In other words, µ is the law of the excess E(K) = ij K ij, over the Hadamard equivalence class of H. Here the fact that µ is the law of E comes from the fact that E is invariant under permutations of rows and columns. Indeed, the law of E over the equivalence class of H is therefore the same as the law of the variable f : T N T N C given by: ] f(a,b) = E [(a i b j H ij ) ij = a i b j H ij ij As a first observation regarding the glow, we first have: Proposition 4.6. The glow µ P(C) has the following properties: (1) µ is invariant under rotations around the origin, and we have µ = ε µ +, where µ + is the law of E [0, ), and ε is the uniform measure on T. (2) supp(µ) D 0 (N N) and supp(µ + ) [0,N N] and in both cases, the covering radius R = N N is optimal precisely when H can be put in bistochastic form. Proof. These results are all elementary, as follows: (1) Here the first assertion follows by using the transformations H zh, with z T, which leave invariant µ, and the second assertion follows from it. (2) Here the assertion regarding µ follows from Proposition 4.3 above, and the assertion regarding µ + follows from it, and from (1). As a basic example, at N = 2 we have the following result: ij

14 14 TEO BANICA Proposition 4.7. For the Fourier matrix F 2 we have E 2 4+2Re(α β) for certain variables α,β T which are uniform, and independent. Proof. Thematrixthatweinterestedin, namelythefouriermatrixf 2 alteredbyavertical switching vector (a, b) and an horizontal switching vector (c, d), is: ( ) ac ad F 2 = bc bd With this notation, we have the following formula: E 2 = ac+ad+bc bd 2 = 4+ ad bc + bc ad bd ac ac bd For proving that α = ad bd and β = are independent, we use the moment method: bc ac ( ) p ( ) q ad bd = a p q b q p c p q d p+q = δ p,q,0 T bc ac 4 T T T T Thus α,β are indeed independent, and we are done. As a consequence, the excess over the equivalence class of F 2 can descrease continuously from 2 2 to 0, the endpoints being realised for instance by the following matrices: ( ) ( ) 1 i 1 i F 2 =, F i 1 2 = i 1 The glow was studied in [3] by using the moment method: Proposition 4.8. For H M N (T) the even moments of E are given by H i1 j a i b j H ij 2pd(a,b) = 1...H ipj p H k1 l 1...H kpl p T N T N ij [i]=[k],[j]=[l] where the sets between brackets are by definition sets with repetition. Proof. We have indeed the following computation: E 2p ( 2p = H ij a i b j = T N T N T N T N ij T N T N ijkl = H i1 j 1...H ipj p a i1...a ip H k1 l ijkl 1...H kpl p T a N k1...a kp H ij H kl aib j a k b l b j1...b jp T b N l1...b lp Now since the integrals at right equal respectively the Kronecker symbols δ [i],[k] and δ [j],[l], we are led indeed to the formula in the statement. As explained in [3], we can deduce from this an asymptotic result, as follows: ) p

15 FOURIER MATRICES 15 Theorem 4.9. The glow of a complex Hadamard matrix H M N (C) is given by ( 1 ij a ) ib j H ij 2p ( ) p d(a,b) = 1 N 1 +O(N 2 ) p! N 2 T N T N so in particular, E/N becomes complex Gaussian in the N limit. Proof. This follows indeed from Proposition 4.8 above, by letting N, and by using some standard combinatorial estimates. See [3]. In relation with our present considerations, which concern fixed matrices, this kind of result cannot really be used, in order to derive some explicit consequences. However, this provides evidence for the fact that many complex Hadamard matrices H M N (C) can be probably put in bistochastic form, at least in some kind of N limit sense. 5. Deformed tensor products We will be interested in what follows in the following general construction of complex Hadamard matrices, due to Diţă [9], which produces many interesting examples: Proposition 5.1. If H M M (C) and K M N (C) are Hadamard, then so is H Q K = (Q ib H ij K ab ) ia,jb for any choice of a parameter matrix Q M M N (T). Proof. This is indeed elementary, and known since [9]. As a basic example, consider the matrix F 2 F 2. With binary indices, and by using the lexicographic ordering {00, 01, 10, 11} on the double indices, we have: ( ) ( ) F 2 F 2 = = By deforming now this matrix as in Proposition 5.1 above, we obtain: ( ) ( ) p q p q 1 1 F 2 Q F 2 = 1 1 p q 1 1 = p q p q 1 1 r s r s r s r s r s Observe that when Q is dephased, p = q = r = 1, we obtain the matrix F4 s from section 1. In general, by dephasing, we obtain a matrix which is equivalent to F ps/qr 4. In order to formulate our results, we will use the following notion: Definition 5.2. We say that a complex Hadamard matrix H M N (C) is in almost bistochastic form when all the row sums belong to N T.

16 16 TEO BANICA Observe that, assuming that this condition holds, the matrix H can be of course put in bistochastic form, just by multiplying its rows by suitable numbers. We will be particularly interested here in the special situation where the affine deformations H q M N (C) of a given complex Hadamard matrix H M N (C) can be put in almost bistochastic form, independently of the value of the parameter q. For the simplest deformations, namely those of F 2 F 2, this is indeed the case: Proposition 5.3. The deformations of F 2 F 2, with parameter matrix Q = ( p r q s), p q p q F 2 Q F 2 = p q p q r s r s r s r s can be put in almost bistochastic form, independently of the value of Q. Proof. By multiplying the columns of the matrix in the statement with 1,1, 1,1 respectively, we obtain the following matrix: p q p q F 2 Q F 2 = p q p q r s r s r s r s The row sums of this matrix being 2q, 2q,2r,2r 2T, we are done. We will see later on that the above matrix F 2 Q F 2 is equivalent to a certain matrix F 2 F 2, which looks a bit more complicated, but is part of a series F N F N. Now back to the general case, we have the following result: Theorem 5.4. A deformed tensor product H Q K can be put in bistochastic form when there exist numbers x i a T such that with we have (H G) ib = MN, for any i,b. G ib = (K x i ) b Q ib Proof. The deformed tensor product L = H Q K is given by L ia,jb = Q ib H ij K ab. By multiplying the columns by scalars R jb T, this matrix becomes: L ia,jb = R jbq ib H ij K ab

17 The row sums of this matrix are given by: FOURIER MATRICES 17 S ia = jb = b = b R jb Q ib H ij K ab K ab Q ib H ij R jb K ab Q ib (HR) ib j In terms of the variables C i b = Q ib(hr) ib, these rows sums are given by: S ia = b K ab C i b = (KC i ) a Thus H Q K can be put in bistochastic form when we can find scalars R jb T and x i a T such that, with Ci b = Q ib(hr) ib, the following condition is satisfied: (KC i ) a = MNx i a, i,a But this condition is equivalent to KC i = MNx i for any i, and by multiplying to the left by the adjoint matrix K, we are led to the following condition: NC i = MK x i, i Now by recalling that Cb i = Q ib(hr) ib, this condition is equivalent to: NQib (HR) ib = M(K x i ) b, i,b With G ib = (K x i ) b /Q ib as in the statement, this condition reads: N(HR)ib = MG ib, i,b But this condition is equivalent to NHR = MG, and by multiplying to the left by the adjoint matrix H, we are led to the following condition: MNR = H G Thus, we have obtained the condition in the statement. As an illustration for this result, assume that H,K can be put in bistochastic form, by using vectors y T M,z T N. If we set x i a = y iz a, with Q = 1 we have: G ib = (K x i ) b = [K (y i z)] b = y i (K z) b

18 18 TEO BANICA We therefore obtain the following formula: (H G) ib = j (H ) ij G jb = j (H ) ij y j (K z) b = (H y) i (K z) b Thus the usual tensor product H K can be put in bistochastic form as well. 6. Fourier matrices Generally speaking, going beyond Theorem 5.4 seems impossible, at this level of generality. However, in the case H = F M the equations simplify, and we have: Proposition 6.1. A deformed tensor product F M Q K can be put in bistochastic form when there exist numbers x i a T such that with G ib = (K x i ) b Q ib we have the following formulae, with l being taken modulo M: G jb Ḡ j+l,b = MNδ l,0, l,b j Moreover, the M N matrix G jb 2 is row-stochastic with sums N 2, and the l = 0 equations state that this matrix must be column-stochastic, with sums MN. Proof. With notations from Theorem 5.4, and with w = e 2πi/M, we have: (H G) ib = j w ij G jb The absolute value of this number can be computed as follows: (H G) ib 2 = jk w i(k j) G jb Ḡ kb = jl w il G jb Ḡ j+l,b = l w il j G jb Ḡ j+l,b If we denote by v b l the sum on the right, we obtain: (H G) ib 2 = l w il v b l = (F M v b ) i

19 FOURIER MATRICES 19 Now if we denote by ξ the all-one vector in C M, the condition (H G) ib = MN for any i,b found in Theorem 5.4 above reformulates as follows: F M v b = MNξ, b By multiplying to the left by FM /M, this condition is equivalent to: MN v b = NFM ξ = 0., b 0 Let us examine the first equation, v0 b = MN. By definition of vb l, we have: v b 0 = j G jb Ḡ jb = j G jb 2 Now recall from Theorem 5.4 that we have G jb = (K x j ) b /Q jb, for certain numbers x j b T. Since we have Q jb T and K / N U N, we obtain: G jb 2 = (K x j ) b 2 = K x j 2 2 = N xj 2 2 = N2 b b Thus the M N matrix G jb 2 is row-stochastic, with sums N 2, and our equations v0 b = MN for any b state that this matrix must be column-stochastic, with sums MN. Regardingnow theother equations thatwe found, namelyvl b = 0forl 0, by definition of vl b and of the variables G jb, these state that we must have: G jb Ḡ j+l,b = 0, l 0, b Thus, we are led to the conditions in the statement. j As an illustration for this result, let us go back to the Q = 1 situation, explained at the end of section 5. By using the formula G ib = y i (K z) b there, we have: G jb Ḡ j+l,b = y j (K z) b y j+l (K z) b j j y j = (K z) b 2 j y j+l = M Nδ l,0 Thus, if K can be put in bistochastic form, then so can be put F M K. As a second illustration, let us go back to the matrices F 2 Q F 2 from the proof of Proposition 5.3 above. The vector of the row sums being S = (2q, 2q,2r,2r), we have

20 20 TEO BANICA x = (q, q,r,r), and so we obtain the following formulae for the entries of G: [( )( )] ( 1 1 q q 2q) b b G 0b = = Q 0b Q 0b [( )] ( ) 1 1 r 2r 1 1)( r 0 b b G 1b = = Q 1b Q 1b Thus, in this case the matrix G is as follows, independently of Q: ( ) 0 2 G = 2 0 In particular, we see that the conditions in Proposition 6.1 are satisfied. As a main application now, we have the following result: Theorem 6.2. The Diţă deformations F N Q F N can be put in almost bistochastic form, independently of the value of the parameter matrix Q M N (T). Proof. We use Proposition 6.1 above, with M = N, and with K = F N. Let w = e 2πi/N, and consider the vectors x i T N given by: x i = (w (i 1)a ) a Since K K = N1 N, and x i are the column vectors of K, shifted by 1, we have: 0 N 0 K x 0 = 0. 0, K x 1 = 0. 0,..., K x N 1 = 0. N N 0 0 We conclude that we have (K x i ) b = Nδ i 1,b, and so the matrix G is given by: G ib = Nδ i 1,b Q ib With this formula in hand, the sums in Proposition 6.1 are given by: G jb Ḡ j+l,b = j j Nδ j 1,b Q jb Nδ j+l 1,b Q j+l,b In the case l 0 we clearly get 0, because the products of Kronecker symbols are 0. In the case l = 0 the denominators are Q jb 2 = 1, and we obtain: G jb Ḡ jb = N 2 δ j 1,b = N 2 j j Thus, the conditions in Proposition 6.1 are satisfied, and we obtain the result.

21 Here is an equivalent formulation of the above result: FOURIER MATRICES 21 Proposition 6.3. The matrix F N Q F N, with Q M N (T), defined by (F N Q F N) ia,jb = wij+ab w bj+j Q ib Q b+1,b where w = e 2πi/N is almost bistochastic, and equivalent to F N Q F N. Proof. OurclaimisthatthisisthematrixconstructedintheproofofTheorem6.2. Indeed, let us first go back to the proof of Theorem 5.4. In the case M = N and H = K = F N, the Diţă deformation L = H Q K studied there is given by: L ia,jb = Q ib H ij K ab = w ij+ab Q ib As explained in the proof of Theorem 5.4, if the conditions in the statement there are satisfied, then the matrix L ia,jb = R jbl ia,jb is almost bistochastic, where: MN R = H G In our case now, M = N and H = K = F N, we know from the proof of Proposition 6.1 that the choice of G which makes work Theorem 5.4 is as follows: G ib = Nδ i 1,b Q ib With this formula in hand, we can compute the matrix R, as follows: R jb = 1 N (H G) jb = 1 w ij G ib = w ij δi 1,b = w (b+1)j N Q i i ib Q b+1,b Thus, the modified version of F N Q F N which is almost bistochastic is given by: L ia,jb = R jb L ia,jb = w (b+1)j Q b+1,b w ij+ab Q ib = wij+ab w bj+j Q ib Q b+1,b Thus we have obtained the formula in the statement, and we are done. As an illustration, let us work out the case N = 2. Here we have w = 1, and with Q = ( p r q s), and then with u = p,v = s, we obtain the following matrix: r q p q p q r q r q u 1 u 1 p F 2 Q F 2 = q p q r q r q r s r s = u 1 u 1 1 v 1 v r q r q 1 v 1 v r r s q r r s q

22 22 TEO BANICA Observe that this matrix is indeed almost bistochastic, with row sums 2, 2,2,2. We conjecture that the converse of Theorem 6.2 holds, in the sense that if the matrices H Q K can be put in almost bistochastic form, independently of the value of the parameter matrix Q M N (T), then we must have H = K = F N, for some N N. References [1] J. Backelin, Square multiples n give infinitely many cyclic n-roots, preprint [2] T. Banica, The defect of generalized Fourier matrices, Linear Algebra Appl. 438 (2013), [3] T. Banica, The glow of Fourier matrices: universality and fluctuations, Oper. Matrices 9 (2015), [4] T. Banica, Deformed Fourier models with formal parameters, Studia Math. 239 (2017), [5] T. Banica and J. Bichon, Random walk questions for linear quantum groups, Int. Math. Res. Not. 24 (2015), [6] T. Banica and I. Nechita, Almost Hadamard matrices with complex entries, Adv. Oper. Theory 3 (2018), [7] T. Banica, I. Nechita and K. Życzkowski, Almost Hadamard matrices: general theory and examples, Open Syst. Inf. Dyn. 19 (2012), [8] G. Björck, Functions of modulus 1 on Z n whose Fourier transforms have constant modulus, and cyclic n-roots, NATO Adv. Sci. Inst. Ser. C Math. Phys. Sci. 315 (1990), [9] P. Diţă, Some results on the parametrization of complex Hadamard matrices, J. Phys. A 37 (2004), [10] J.-C. Faugère, Finding all the solutions of Cyclic 9 using Gröbner basis techniques, Lecture Notes Ser. Comput. 9 (2001), [11] P.C. Fishburn and N.J.A. Sloane, The solution to Berlekamp s switching game, Discrete Math. 74 (1989), [12] U. Haagerup, Orthogonal maximal abelian -subalgebras of the n n matrices and cyclic n-roots, in Operator algebras and quantum field theory, International Press (1997), [13] M. Idel and M.M. Wolf, Sinkhorn normal form for unitary matrices, Linear Algebra Appl. 471(2015), [14] H. Kharaghani and J. Seberry, The excess of complex Hadamard matrices, Graphs Combin. 9 (1993), [15] R. Nicoara and J. White, Analytic deformations of group commuting squares and complex Hadamard matrices, J. Funct. Anal. 272 (2017), [16] R. Roth and K. Viswanathan, On the hardness of decoding the Gale-Berlekamp code, IEEE Trans. Inform. Theory 54 (2008), [17] J. Seberry and M. Yamada, Hadamard matrices, sequences, and block designs, Wiley (1992). [18] D.R. Stinson, Combinatorial designs: constructions and analysis, Springer-Verlag (2006). [19] W. Tadej and K. Życzkowski, A concise guide to complex Hadamard matrices, Open Syst. Inf. Dyn. 13 (2006), [20] W. Tadej and K. Życzkowski, Defect of a unitary matrix, Linear Algebra Appl. 429 (2008), T.B.: Department of Mathematics, University of Cergy-Pontoise, F Cergy- Pontoise, France. teo.banica@gmail.com