arxiv: v3 [math.pr] 18 Oct 2017

Transcription

1 MEIXNER MATRIX ENSEMBLES arxiv: v3 [math.pr] 18 Oct 017 W LODZIMIERZ BRYC AND GÉRARD LETAC Abstract. We construct a family of matrix ensembles that fits Anshelevich s regression postulates for Meixner laws on matrices, namely the distribution with an invariance property of X when E(X X+Y) = a(x+y) +b(x+y)+ci n where X and Y are i.i.d. on symmetric matrices of order n. We show that the Laplace transform of a general n n Meixner matrix ensemble satisfies a system of partial differential equations which is explicitly solvable for n =. We rely on these solutions to identify the six types of Meixner matrix ensembles. 1. Introduction The classical Meixner laws are the one dimensional binomial, Poisson, negative binomial, gamma, normal, and hyperbolic Meixner laws that include the hyperbolic secant law. These laws make their appearance as the orthogonality measures of a certain family of orthogonal polynomials [35], as the laws characterized by quadratic regression property [7], and as the exponential families corresponding to the quadratic variance functions [36], see also [18, 1, 8, 40, 41]. The Laha-Lukacs [5, 7] characterization is the description of the possible distributions of a square-integrable random variable X such that there exist real numbers a,b,c satisfying E(X S) = as + bs + c where S = X +Y, where Y X and where Y is independent of X. Another convenient form of this condition is to postulate the existence of real numbers A,B,C such that (1.1) E((X Y) S) = AS +BS +C Date: Created: June 006. Revised: Sept 010, and March 011. Printed: September 13, 018 file: BandB-0 cln.tex. 000 Mathematics Subject Classification. 60B0. Key words and phrases. Meixner laws; random projections; quadratic conditional moments; matrix ensembles; systems of PDEs; Jack polynomials. Department of Mathematical Sciences, University of Cincinnati, Cincinnati, OH , USA. Laboratoire de Statistique et Probabilités, Université Paul Sabatier, 3106, Toulouse, France. 1

2 W LODZIMIERZ BRYC AND GÉRARD LETAC (use(x Y) = X +Y S to see that(a,b,c) = (4a 1,4b,4c)). In006M.Anshelevich[4, pagesand5]observesthat(1.1)makes sense when the real random variable X is replaced by a random square matrix X. He then raises the question of defining analogues of Meixner distributions on matrices. In other words, he asks for the matrix version of Laha-Lukacs [7] result. This question on random matrices is so general that it is suitable to restrict this problem to random variables X valued in the space H n,1 of all symmetric matrices and such that UXU X for all matrices U in the orthogonal group O(n) = K n,1. On the other hand, it is natural to extend this simpler framework to the space H n, of Hermitian complex matrices for random variables X such that UXU X for all matrices U in the unitary group U(n) = K n, and even to the space H n,4 of Hermitian-quaternionic matrices, again for random variables X such that UXU X for matrices U in the symplectic group Sp(n) = K n,4. There are also connections with free probability (see [10] and Section 6.3 below), and connections with Jordan algebras (see Section 6.4 below). A probability law µ on H n,β (where β = 1, or 4), or a random variable X with law µ, is said to be rotation invariant if the law of UXU is also µ for all U in K n,β. In such case, following a physicists tradition we shall call µ an ensemble. Definition 1.1. We will say that the probability law µ on H n,β is a Meixner ensemble with parameters A,B,C R if µ is rotation invariant, the second moments exist, and µ has the following property: if X,Y are independent with the same law µ and S = X+Y, then (1.) E((X Y) S) = AS +BS+CI n. This is a different concept of a Meixner matrix ensemble than the one introduced in [16] under the name multivariate Meixner classes of invariant distributions of random matrices. In particular, we note that the oldest known ensemble, namely the Wishart distribution with mean proportional to I n is not a Meixner ensemble in the sense of Definition 1.1, since if X, Y are i.i.d. with Wishart distribution, then there are constants A,B such that (1.3) E((X Y) S) = AS +BStrS, which is quadratic in S but not of the desired form (1.). For a proof of (1.3), see [30, page 58]. To complete this introduction, let us make a number of simple remarks about the Meixner ensembles. It is easy to check that if the law of X is a Meixner ensemble with parameters A,B,C, then an affine

3 MEIXNER MATRIX ENSEMBLES 3 transformation X = αx+t with real α,t is also Meixner, with parameters à = A, B = αb 4At, C = α C +4At Bαt. In particular, since passing to X changes only the sign of B, without loss of generality we will consider only B 0. Suppose now that E(X) and E(X ) are given; in this case the constants A,B,C satisfy an equation with coefficients that depend on E(X), E(X ). The equation is obtained from taking the expected values of both sides of (1.). Finally, rotation invariance implies that both E(X) and E(X ) are constant multiples of I n, so in the non-degenerate case E(X ) 0 we can normalize the matrices so that E( X) = 0 and E( X ) = I n, and express relation (1.) in terms of two real constants a,b as (1.4) E ( ( X Ỹ) ) S) = (I n +b S+a S. 1+a In this normalized setting one can classify the Meixner ensembles into (i) the elliptic case b > 4a, which in Laha-Lukacs theorem on R corresponds to negative binomial (a > 0), Poisson (a = 0), or binomial (a < 0) laws. (ii) the parabolic case b = 4a, which corresponds to gamma law (a > 0) or Gaussian law (a = 0) (iii) the hyperbolic case b < 4a, which corresponds to the family of hyperbolic Meixner laws that include hyperbolic secant law. Of course, for non-degenerate ensembles (1.) and (1.4) are equivalent. If (1.) holds with constants A 1,B,C, then the standardization X of X with mean µi n and variance Var(X) := E(X ) (E(X)) = σ I n satisfies (1.4) with A (1.5) a = (1 A), b = B +4Aµ σ(1 A). What we call trivial Meixner ensembles are random multiples of the identity, X = ξi n where ξ is R-valued with one of the six Meixner laws. The aim of the paper is the study of non-trivial Meixner ensembles for alln. Wewillbeabletofindallofthemforn =, wewillgiveexamples of Meixner ensembles for n 3 and we will write a system of n linear partial differential equations (PDEs) of the second order satisfied by the Laplace transform of a Meixner ensemble of given parameters (a, b) as defined by (1.4). Let us emphasize the fact that several different ensembles correspond to one set of parameters a,b in (1.4). The zoo of the known Meixner ensembles includes obviously the Gaussian ensembles which are described by their cumulant transform c 1 trθ+ c (trθ) + c 3 tr(θ ) with the familiar GOE, GUE and GSE (see

4 4 W LODZIMIERZ BRYC AND GÉRARD LETAC [] or [34]) corresponding to c 3 = 1 and c 1 = c = 0 among them. The Gaussian ensembles yield a family of Meixner ensembles corresponding to a = b = 0 in (1.4), but the existence of non Gaussian ensembles for a = b = 0 and n 3 has not been proved or disproved. Another important example is a Bernoulli ensemble: if µ m is the distribution of the projection on a uniformly random m dimensional space of the Euclidean space of dimension n, a Bernoulli ensemble is any mixing of µ 0,...,µ n ; for n = 1 we get back the ordinary Bernoulli distribution. The Bernoulli ensembles can be seen as the only distributions of random projections P such that UPU P for all U K n,β. The convolution of N identical Bernoulli ensembles provides the binomial ensemble, and randomizing N by Poisson or negative binomial distribution yields the Poisson and the negative binomial ensembles. The construction of the full family of Meixner ensembles with b = 4a or with b < 4a is available only for n = and their extension to n 3 is a challenge. The paper is organized as follows. In Section we describe Laplace transforms of Meixner ensembles. In Section 3 we give examples of Meixner ensembles. In Section 4 we derive the system of PDEs that determine the Laplace transforms of the general n n Meixner ensembles. The general solutionof this system is elusive for n 3. InSection 5 we use the system of PDEs to show that under a natural integrability assumption the examples of Meixner ensembles from Section 3 exhaust all Meixner ensembles. In Section 6 we collected additional material: Proposition 6.1 answers in negative a question raised in [10, Remark 5.8] for random matrices with finite Laplace transform. Proposition 6.3 gives a series expansion for a Laplace transform of a random projection which can be used to derive (some) solutions to the system of PDEs. Section 6.4 makes some comments about the Jordan algebra context of the problem.. Laplace transforms In this section we rewrite the Meixner property (1.) in terms of the Laplace transform. We will work simultaneously with symmetric, unitary and symplectic ensembles, thus with probability laws on H n,β which are invariant under the action X UXU where U is in the group K n,β. We call β = 1,,4 the Peirce constant (see [17, page 97], where β is denoted as d and n is denoted by r). The space H n,β is a real vector space equipped with the Euclidean structure defined by the real inner product (x,y) x y := R(tr(xy))

5 MEIXNER MATRIX ENSEMBLES 5 (the real part of the trace is needed for the symplectic case β = 4 only). For an H n,β -valued random variable X, the Laplace transform is a mapping from H n,β into (0, ], defined by Let L(θ) = E(exp( θ X )). (.1) Θ X = int{θ : L(θ) < }. (We will also use notation Θ µ when µ is the law of random variable X.) Throughout this paper we consider only random matrices X with values in H n,β with non-empty Θ X. We need a result of linear algebra that we are not going to prove. Consider the space L s (H n,β ) of symmetric endomorphisms of the Euclidean space H n,β. To each y H n,β we associate the elements y y and P y of L s (H n,β ) defined respectively by h (y y)(h) = ytr(yh), h P y (h) = yhy. TheyprovideimportantexamplesofelementsofL s (H n,β ). Now,L s (H n,β ) is itself a linear space, and the result that we are going to admit as a black box is the following (see [15, Lemma 6.3] and [30, Proposition 3.1] for a proof): Proposition A. There exists a unique endomorphism Ψ of L s (H n,β ) such that for all y H n,β one has Ψ(y y) = P y. Furthermore (.) Ψ(P y ) = β ( y y + 1 β ) P y. With this notation we have the simple but crucial result. (Note that in this result we could have dispensed with the hypothesis of invariance under rotations.) Proposition.1. Let µ be an ensemble on H n,β such that its Laplace transform L is finite on an open non empty set. The following are equivalent: (i) The ensemble µ is Meixner with real parameters A,B,C, i.e. (1.) holds for X µ. (ii) With k(θ) = lnl(θ), we have (.3) (1 A)Ψ(k (θ))(i n ) = 4A(k (θ)) +Bk (θ)+ci n for all θ Θ X. Proof. To prove (i) (ii), we first observe that (.4) E ( (X Y) (X Y)e θ S ) = L (θ)l(θ) L (θ) L (θ).

6 6 W LODZIMIERZ BRYC AND GÉRARD LETAC To see this, note that since X, Y are i.i.d, the left hand side of (.4) is L(θ)E ( X Xe θ X ) E ( Xe θ X ) E ( Ye θ Y ). The same reasoning gives (.5) E ( S Se θ S ) = L (θ)l(θ)+l (θ) L (θ). Since E ( Se θ S ) = L(θ)L (θ), applying Ψ to (.4) and (.5), from (1.) we get (.6) (1 A)Ψ(L (θ))(i n )L(θ) = (1+A)(L (θ)) +BL(θ)L (θ)+c(l(θ)) I n. Since L = k L and L = (k +k k )L, this implies (.3). Implication (ii) (i) also follows from the above identities. Applying again Ψ to (.4) and (.5), from (.6) we infer that E ( (X Y) e θ S ) = E ( (AS +BS+CI)e θ S ) for all θ Θ X. It is known that this is equivalent to (1.), see [5, Section 1.1.3]. Remark.1. In the above proposition, the use of the endomorphism Ψ may seem surprising. For a random variable X valued in a linear space E, the link between E(X X) and the Laplace transform L of X is easy through the second differential of L. However, if E is an algebra, like the space of square real matrices, relating E(X ) to L is difficult. Here our reasoning was based on P X (A) = XAX and thus P X (I n ) = X. Some other facts about Ψ are known; it is a symmetric operator on L s (H n,β ) with only two eigenspaces: the one generated by all the β y y +P y for the eigenvalue 1 and the one generated by all the y y P y for the eigenvalue β/. Details are in [3, Section 5] where the dimensions of these two subspaces are also computed. Remark.. Let µ be a Meixner ensemble with cumulant function k that satisfies (.3). If θ 0 Θ µ, then k θ0 (θ) = k(θ + θ 0 ) k(θ 0 ) is the cumulant function of the probability measure µ θ0 (dx) = e θ 0 x k(θ 0 ) µ(dx). Inotherterms, µ θ0 isamember ofthenaturalexponential familygenerated by µ. Clearly, k θ0 also satisfies (.3). However, µ θ0 is an ensemble i.e. it is rotation invariant, if and only if θ 0 = ti n for some real t. In this case θ 0 x = ttrx. Remark.3. Note that if A = 1 then X is degenerate. Indeed, from (.3) we see that A = 1 implies k (θ) = const.

7 MEIXNER MATRIX ENSEMBLES 7 We now show that Meixner ensembles are preserved under convolution power. Definition.1. Let µ be a probability measure on a finite dimensional linear space such that its Laplace transform is finite on a non empty open set. The Jørgensen set Λ(µ) of µ is the set of α > 0 such that there exists a probability measure µ α with L µα = L α µ and Θ µ α = Θ µ, see e.g. [15, page 767]. For instance, Λ(µ) = (0, ) if µ is infinitely divisible, and Λ(µ) is the set of positive integers if µ is the Bernoulli distribution on {0,1}. Proposition.. Suppose Θ µ. If µ is a Meixner ensemble with parameters A 1,B,C and if α Λ(µ), then µ α is a Meixner ensemble with parameters A α = A/(A+α(1 A)), B α = αb/(a+α(1 A)), C α = α C/(A+α(1 A)). Proof. It is clear that µ α is invariant under rotations. Since k α (θ) := k µα (θ) = αk µ (θ), multiplying (.3) by α(1 A α )/(1 A) we get (1 A α )Ψ(k α(θ))(i n ) = 4 A(1 A α) α(1 A) (k α(θ)) + B(1 A α) 1 A k α(θ)+ Cα(1 A α) 1 A. Solving the equation A α = A(1 A α) α(1 A), we get the formulas for A α,b α,c α such that (.3) holds. It will be convenient to have a version of Proposition.1 for standardized ensembles. Corollary.3. Suppose E( X) = 0, E( X ) = I n. Denote k(θ) = ln E(e θ X ). If (1.4) holds, then (.7) Ψ(k (θ))(i n ) = I n +bk (θ)+4a(k (θ)). Conversely, if the logarithm of the Laplace transform of X satisfies (.7) (together with the initial conditions) then (1.4) holds, and E( X) = 0, E( X ) = I n. Proof. We apply (.3) with A = a/(1 + a), B = b/(1 + a), C = /(1+a). 3. Examples of Meixner Ensembles In this section we give examples of Meixner ensembles. In Theorem 5.1 we will show that the examples exhaust the family of all Meixner ensembles on H,β with Θ µ.

8 8 W LODZIMIERZ BRYC AND GÉRARD LETAC 3.1. Bernoulli and binomial ensembles. The Bernoulli ensembles (which are different from the Bernoulli random matrices with independent two-valued entries) will be our basic building blocks for more complicated ensembles. Suppose P m is the orthogonal projection onto the random and uniformly distributed m-dimensional subspace, i.e. P m = UI m,n U where I m,n is the n n matrix with m ones on the diagonal followed by n m zeroes, and U is uniformly distributed on the group K n,β. Definition 3.1. Denote by µ m the distribution of P m. A Bernoulli ensemble with parameters q 1,...,q n 0 such that q 1 + +q n 1, is the law µ = q 0 δ 0 +q 1 µ 1 + +q n 1 µ n 1 +q n δ In, where q 0 = 1 (q 1 + +q n ) and δ I denotes a point-mass at I. In terms of random variables, Bernoulli ensemble can be realized by the random variable 0 with probability q 0 = 1 (q 1 + +q n ), P 1 with probability q 1, (3.1) X =. I n with probability q n. Proposition 3.1. A Bernoulli ensemble is a Meixner ensemble with parameters A = 1, B =, C = 0. Proof. We note that for any pair of projections (3.) (P Q) = (P +Q) (P +Q). Applying this algebraic identity to random projections X and Y, where X is given by (3.1) and Y is its independent copy, we see that (1.) holds. We note that for the Bernoulli ensemble formula (.6) takes a particularly simple form: its Laplace transform L B (θ) = E(exp( θ X )) satisfies (3.3) Ψ(L B (θ))(i n) = L B (θ). The corresponding version of (.3) is (3.4) Ψ(k B (θ))(i n) = k B (θ)(i n k B (θ)). The fact that all Bernoulli ensembles share the same A,B,C will be convenient for calculations. But it will be easier to classify the resulting Meixner laws into the familiar families from the standardized form that appears in formula (1.4). To do so, we find the moments of X. We first

9 MEIXNER MATRIX ENSEMBLES 9 note that by rotation invariance, E(P k ) = ci n, so Etr(P k ) = nc. But tr(p k ) = tr(i k,n ) = k, so (3.5) E(P k ) = E(P k ) = k n I n. Therefore, with q = 1(q n 1+q + +nq n ), we have E(X ) = E(X) = qi n. If q 0,1, the standardized matrix is X 1 = q(1 q) (X qi n), and we get (3.6) E(( X Ỹ) S) = 4 ( I n + (1 q) S q(1 q) S ) 4 So (1.4) holds with parameters a = 1/4 and b = (1/ q)/ q(1 q). An intrinsic characterization of Bernoulli ensembles with an elementary proof is as follows. Proposition 3.. If X with values in H n,β has a law invariant under rotations and X = X then its law is a Bernoulli ensemble. Proof. Since X is symmetric, it is an orthogonal projection. Denote by d(x) the dimension of its image. Invariance under rotations means that for any bounded measurable function f, we have E(f(UXU ) f(x)) = 0 for all U K n,β. If d = 0,1,...,n and if f is zero outside of matrices of trace d, then E(f(UXU ) f(x) d(x) = d) = E(f(UXU ) f(x)) = 0. This implies that the conditional law L(X d(x) = d) is invariant under rotations, so it is a law of P d. Then X is a Bernoulli ensemble with parameters q m = Pr(d(X) = m). This implies the following converse of Proposition 3.1. Proposition 3.3. Suppose that a law µ belongs to a Meixner ensemble with parameters A = 1, B =, C = 0, and that the first four moments are finite. Then µ is a Bernoulli ensemble. Proof. Let X,Y be independent random variables with law µ, and let S = X+Y. Since (X Y) = X +Y S, and since the law of (X,S) is the same as (Y,S) from (1.) we get E(X S) = 1 S..

10 10 W LODZIMIERZ BRYC AND GÉRARD LETAC From E(X S) = 1 S we see that (3.7) E(X S) = E(X S). From this, we deduce that E(X) = E(X ), and that which implies that E(X 3 ) = E(X). Using (3.7) again, we have E(X S) = E(XS), (3.8) tr E(X S ) = tr E(XS ). Since E tr = tr E, from (3.8) using independence and cyclic invariance of the trace we get tr E(X 4 )+tr(ex 3 EX)+tr((EX ) ) = tr E(X 3 )+tr(ex ) +tr(ex EX). Since we already proved that E(X 3 ) = E(X ) = E(X), we see that tr E(X 4 ) = tr E(X 3 ). Let Λ 1,Λ,...Λ n be the (random) eigenvalues of X. From the equality of the first four moments of X, we see that n E(Λ j (1 Λ j )) = tr E(X (I n X) ) = tr E(X )+tr E(X 4 ) tr E(X 3 ) = 0. j=1 Therefore, all eigenvalues of X are either 0 or 1, and X = X. So by Proposition 3., the law µ of X is a Bernoulli ensemble. The Jørgensen set of a Bernoulli ensemble is of interest due to connections with free probability which allows continuous values for the analog of parameter N. For n =, we will show that the Jørgensen set is N. The following sufficient condition shows that this is a generic case. Proposition 3.4. If α Λ(µ) and µ is a Bernoulli ensemble with parameters q 1,...,q n on H n,β, then (q 0 +q 1 z+ +q n z n ) α is a polynomial in variable z. In particular, if q n > 0 then nα is an integer. Proof. For θ s = diag(s,s,...,s), the Laplace transform is L(θ s ) = n r=0 q re rs = j (es z j ) m j, where z j C are the distinct roots of the polynomial q 0 + q 1 z + + q n z n taken with their multiplicities m j. Since the Laplace transform must be an analytic function on its domain, we see that α must be rational, and that if α = p/q with relatively prime p, q N, then q must be a common divisor of all multiplicities m 1,m,... of the roots. So q divides also their sum m 1 +m + = n.

11 MEIXNER MATRIX ENSEMBLES Binomial ensemble. The real, complex or quaternionic binomial ensembles are measures on H n,β which are parametrized by integer N and a discrete probability law on {0,1,...,n}. We will follow the tradition of not listing the probability of 0 among the parameters. FixintegerN andnon-negativenumbersq 1,...,q n withq 1 + +q n 1. Let X 1,...,X N be independent random matrices with the same Bernoulli distribution (3.1). Definition 3.. The binomial ensemble Bin(N,q 1,...,q n ) with parameters N = 1,,... and q 1,...,q n, is the law of X = N j=1 X j. Proposition 3.5. A binomial ensemble with parameter N is a Meixner ensemble with parameters A = 1/(N 1), B = N/(N 1), C = 0. Proof. This is a special case of Proposition., applied to the law µ of X 1 with parameters described in Proposition 3.1, and to α = N. For standardization, we will need to know that (3.9) E(X) = N q and Var(X) = N q(1 q), where q is defined just below (3.5). From (1.5) we then get that (1.4) holds with a = 1/(4N) and b = (1/ q)/ N q(1 q). Remark 3.1. Applying Remark. to the Bernoulli and binomial ensemble with parameters q 1,...,q n gives the new binomial ensemble with parameters q 1 (t),q (t),...,q n (t), where q j (t) = q j r j /P(r), r = e t and P(r) = q 0 +q 1 r+ +q n r n. To see this, we observe that if µ(dx) = q 0 δ 0 + q 1 µ 1 + +q n δ In then e ttrx µ(dx) = q 0 +q 1 rµ 1 + +q n r n δ In since trx = j on the support of µ j. Remark 3.. Let us mention here a geometric interpretation of the binomial distribution for n =. For simplicity we explain this interpretation for β = 1; its extension to β =,4 is fairly straightforward, compare (3.0). All real symmetric matrices are parametrized by (a,b+ic) R C as [ ] a+b c M = v c a b M = (a,b+ic). In particular, matrix M corresponding to v M = (a,z) is semipositive definite if andonlyif a z. We have trm = aanddetm = a z. If P is a projection matrix of rank 1 then v P = 1 (1,eiT ) and P has distribution µ 1 from Definition 3. if and only if e it is uniformly

12 1 W LODZIMIERZ BRYC AND GÉRARD LETAC distributed on the unit circle. If X 1,...,X N are i.i.d. with distribution µ 1 then S N = X X N is Bin(N,1,0) = µ N 1. So we have v SN = 1 (N,eiT 1 +e it + +e it N ), where v Xj = 1 (1,eiT j ). Therefore the eigenvalues of S N are λ ± = 1 (N ± eit 1 +e it + +e it N ). The distribution of e it 1 +e it + +e it N was studied by Kluyver and Rayleigh, see [44, pages ] and it is known that P( e it 1 +e it + +e it N r) = r 0 J 1 (rt)(j 0 (t)) N dt for 0 r N and N. (Here J 0, J 1 are Bessel functions.) In particular, P(λ + λ < 1) = 1, see [4, page 104]. N+1 The general binomial ensemble with parameters N and q 1,q is N! ν 0!ν 1!ν! qν 0 0 q ν 1 1 q ν µ ν 1 1 δ ν I, ν 0,ν 1,ν 0,ν 0 +ν 1 +ν =N andthetermµ ν 1 1 δ ν I isgivenbythedistributionofv = ( ν 1 +ν,e it 1 + e it + +e it N ). 3.. Poissonensemble. OurPoissonensembleshaveparametersλ 1,λ,...,λ n 0 with λ = λ 1 + +λ n > 0. Definition 3.3. Let N be a Poisson real random variable, Pr(N = j) = e λ λ j /j!, j = 0,1,..., and let q 1 = λ 1 /λ,...,q n = λ n /λ 0. Let X 1,X,... be independent Bernoulli matrices with the same parameters q 1,...,q n, and X 0 = 0. The Poisson ensemble is the law of X = N k=0 X k. We will use notation Poiss(λ 1,...,λ n ). Proposition 3.6. The Poisson ensemble is a Meixner ensemble, as it satisfies (1.) with parameters A = 0, B = 1, C = 0. Proof. LetL B (θ)bethelaplacetransformofthecorrespondingbernoulli ensemble. Then the log of the Laplace transform for the Poisson ensemble is k(θ) = λ(l B (θ) 1). From (3.3) we get Ψ(k (θ))(i n ) = k (θ), which is (.3) with A = 0, B = 1, C = 0. Therefore (1.) holds by Proposition.1.

13 MEIXNER MATRIX ENSEMBLES 13 For standardization, we will need to know that, with λ = (λ 1 +λ + +nλ n )/n, (3.10) E(X) = λi n and Var(X) = λi n. Indeed, E(X) = E(E(X N)) and Var(X) = E(Var(X N))+Var(E(X N)) = E(N q(1 q))i n +Var(N q)i n. One can then check that (1.4) holds with a = 0 and b = λ. 1 It is also easy to see that if the law of X is Poisson with parameters λ 1,...,λ n > 0, then X = n m=1 X m is the sum of independent random variables X m from Poisson ensembles with parameters (0,...,0,λ m,0...,0). Furthermore, real random variable tr(x) has the compound Poisson law with the law of summands given as n λ k k=0 λ δ k. We remark that the Poisson model with exactly one parameter λ 1 0 appears explicitly in [13, page 638], as part of the construction of matrix models for all free-infinitely divisible laws; see also[8]. However, regression properties of the model were not analyzed Negative binomial ensemble. We now use Bernoulli ensembles to construct the negative binomial ensemble. Let N be a negative binomial random variable with parameters r > 0, 0 < p < 1, i.e. (3.11) P(N = j) = Γ(r+j) p r q j, j = 0,1,,..., q = 1 p. Γ(r)j! Definition 3.4. Thenegativebinomialensemble NB(r,q 1,...,q n )with parameters r > 0 and q 1,...,q n 0 such that q 1 + +q n < 1, is the law of the random sum N X = X k, k=0 where N has distribution (3.11) with p = 1 (q 1 + +q n ), X 1,X,..., areindependentbernoulliensembleswithparametersq 1 /(1 p),...,q n /(1 p), and X 0 = 0. Proposition 3.7. The negative binomial ensemble is a Meixner ensemble, as it satisfies (1.) with parameters A = 1, B = r, r+1 r+1 C = 0. Proof. LetL B (θ)bethelaplacetransformofthecorrespondingbernoulli ensemble. The generating function E(z N ) = p r /(1 qz) r gives (3.1) L(θ) = E(e θ X ) = p r (1 ql B (θ)) r, so 1 ql B (θ) = p(l(θ)) 1/r.

14 14 W LODZIMIERZ BRYC AND GÉRARD LETAC Differentiating, we get and ql B(θ) = p r L (θ)(l(θ)) (1+r)/r ql B (θ) = p r L (θ)(l(θ)) (1+r)/r p(1+r) L (θ) L (θ)(l(θ)) (1+r)/r. r From (3.3), we get 1 r Ψ(L (θ)(i n ))L(θ) = 1+r (L (θ)) + 1 r r L (θ)l(θ). We re-write this as r r +1 Ψ(L (θ)(i n ))L(θ) = which is (.6) with A = 1, B = r r+1 by Proposition.1. (r +1) r +1 (L (θ)) + r r +1 L (θ)l(θ), r+1, C = 0. Therefore (1.) holds For standardization, we will need to know that with q = (q 1 +q + +nq n )/n, (3.13) E(X) = r q p I n and Var(X) = r q p (p+ q)i n. Indeed, E(X) = E(E(X N)) = q/(1 p)e(n)i n = qr/pi n and Var(X) = E(Var(X N))+Var(E(X N)) = E(N q(1 q))i n +Var(N q)i n = q(1 q 1 p )r p I n + q r (1 p)p I n. Then from (1.5) we see that (1.4) holds with a = 1 4r and b = p+ q r q(p+ q) Gaussian ensemble. Since X Y and S are independent for any Gaussian independent identically distributed pair X, Y, formula (1.) holds with A = B = 0 and C =. The requirement of rotational invariance reduces the choices of the Gaussian law to the three-parameter family: (3.14) k(θ) = c 1 trθ+c (trθ) /+c 3 tr(θ )/ with c 3 0 and nc + c 3 0. To see this, without loss of generality we take c 1 = 0. Let θ 1,...,θ n be the eigenvalues of θ, and consider the matrix of the quadratic form c (θ 1 + +θ n ) +c 3 (θ 1+ +θ n) whose characteristic polynomial is (z (nc +c 3 ))(z c 3 ) n 1. This shows that c 3 0 and nc +c 3 0. The explicit construction is to start with auxiliary GUE/GOE/GSE matrix Z and independent real standard normal ζ.

15 MEIXNER MATRIX ENSEMBLES 15 Definition 3.5. For c 1 R, c 3 0 and c R such that nc +c 3 0, the Gaussian Meixner ensemble is the law of X = c 3 (Z 1 ) n tr(z)i n + c + c 3 n ζi n +c 1 I n. Proposition 3.8. The logarithm of the Laplace transform of the Gaussian Meixner ensemble is given by (3.14). Proof. Since θ Z atrz = θ atrθ Z and tr((θ 1 n trθ) ) = trθ 1 n (trθ), we see that the answer follows from the well known GUE/GOE/GSE formula Eexp(tr(θZ)) = exp(trθ /) Gamma ensemble. The remaining types of Meixner ensembles will be constructed only for matrices of size n =. One difficulty we encounter is lack of continuity, which we now explain. It is natural to expect that gamma ensemble arises as lim p 0 px p of a sequence of negative binomial ensembles with varying parameter p while r and the ratios q j /(1 p) are kept fixed. However from (3.1) we see that lim p 0 Ee θ pxp = limp r (1 ql B (θ)) r = ( 1 tr(θl p 0 B(0)) r = (1 qtrθ) ) r. SothisisatrivialensembleoftheformξI n withrealgamma-distributed ξ. We now show that there are non-trivial gamma ensembles for n =. This construction is based on a more detailed analysis of the system of PDEs that arises from (.7). Definition 3.6. The Gamma ensemble on H,β with parameters p > β/, c > 1 is defined by its Laplace transform ( (3.15) E(e θ X ) = 1 4c p, 1+βtrθ+βtr θ +4detθ) defined on Θ X = {θ H,β : 1 4c 1+βtrθ+βtr θ+4detθ > 0}. Of course, this definition requires a proof that the required law on H,β exists. The Gamma ensemble on H,β is constructed by choosing theappropriatelawon R β+, andarrangingthe corresponding real random variables into the random matrix. The construction is based on the laws analyzed by Letac and Weso lowski [31, Theorem 3.1]. According to this result, for p > β/ and c > 1 there is a probability measure β+1 ν p,c (dx) on the open Lorentz cone Ω β = {x R β+ : x 0 > r=1 x r }

16 16 W LODZIMIERZ BRYC AND GÉRARD LETAC with the Laplace transform ( (3.16) e β+1 p s rx r ν p,c (dx) = 1 cs 0 +s 0 sr). Ω β Proposition 3.9. Consider (ξ 0,...,ξ β+1 ) with joint distribution ν p,c. The gamma ensemble with parameters p, c is: (i) for β = 1, (3.17) X = (ii) for β =, (3.18) X = r=1 [ ] 1+βξ0 +ξ 1 ξ. ξ 1+βξ0 ξ 1 [ ] 1+βξ0 +ξ 1 ξ +iξ 3, ξ iξ 3 1+βξ0 ξ 1 where i = 1. (iii) for β = 4, [ ] 1+βξ0 +ξ (3.19) X = 1 ξ +iξ 3 +jξ 4 +kξ 5, ξ iξ 3 jξ 4 kξ 5 1+βξ0 ξ 1 where i,j,k are the standard quaternion basis. Proof. We give the proof for β = 4. Writing [ ] s (3.0) θ = 0 +s 1 s +is 3 +js 4 +ks 5 s is 3 js 4 ks 5 s 0 s 1 we have θ X = 5s 0 ξ r=1 s rξ r, so (3.15) follows from (3.16). Proposition A gamma ensemble on H,β with parameters p > β/, c > 1 is a Meixner ensemble with parameters A = 1, B = 0, C = 0. 1+p The mean and the variance are E(X) = pc 1+βI, Var(X) = 4pc (1+β)I. For the standardized version, (1.4) holds with a = 1/(4p) and b = 1/ p. The proof relies on the system of PDEs derived in Theorem 4.1; it appears after Proposition 5.4. Remark 3.3. Somewhat more generally, one can define gamma ensembles on H n,β as the ensembles with Laplace transform (3.1) L(θ) = (1 ctrθ +(β(n 1)+)(trθ) trθ ) p.

17 MEIXNER MATRIX ENSEMBLES 17 Such an ensemble is a Meixner ensemble, as it follows from Theorem 4.1 that (1.4) holds with a = 1/(4p), b = 1/ p. Since this is only a one-parameter subset of the possible solutions, we do not pursue this construction further Exceptional hyperbolic ensemble. It turns out that there is just one exceptional class of non-trivial hyperbolic ensembles that, unlike in all previous cases, cannot be continuously deformed into the trivial ensemble. Recall that the Bessel functions and the modified Bessel functions are [44, 3.1 (8), 3.7 ()] (3.) J ν (z) = n=0 ( 1) n (z/) ν+n n!γ(n+ν +1), I ν(z) = n=0 (z/) ν+n n!γ(n+ν +1). Following [6], we work with differently normalized Bessel functions (3.3) (3.4) J ν (x) := Γ(ν +1)(/x) ν J ν (x), I ν (x) := Γ(ν +1)(/x) ν I ν (x). which are entire functions of z C for ν > 1. Lemma For every α > 0, ν 0 and m N, there exist a unique probability measure µ α on R m with the Laplace transform (3.5) e t,y 1 µ α (dy) = R (J m ν ( t )) α for all t R m small enough. Proof. According to [44, 15.7], Bessel function J ν has simple real zeros that come in opposite pairs. Arranging the positive zeros in increasing order, 0 < j 1 < j <..., from [44, (3)] we get J ν (z) = (1 z /jk ), so for 0 < r < j 1, k=1 1 (J ν ( r)) α = k=1 1 (1 r/j k )α is a Laplace transform (of the series of independent Gamma random variables). Therefore, by Schoenberg s theorem [39, Theorem ], for every m there exists a measure µ α on Borel sets of R m such that for t < j 1, 1 (J ν ( t )) α = exp t,y µ α (dy). R m

18 18 W LODZIMIERZ BRYC AND GÉRARD LETAC The hyperbolic Meixner ensembles have two disconnected components: the two-parameter family of trivial ensembles ξi, and the following exceptional one-parameter non-trivial ensemble. Definition 3.7. For α > 0, consider (ξ 1,...,ξ β+1 ) with joint distribution µ α from Lemma The exceptional hyperbolic ensemble on H,β with parameter α is: (i) for β = 1, (ii) for β =, X = [ ] ξ1 ξ, ξ ξ 1 [ ] ξ X = 1 ξ +iξ 3, ξ iξ 3 ξ 1 where i = 1, (iii) for β = 4, [ ] ξ (3.6) X = 1 ξ +iξ 3 +jξ 4 +kξ 5, ξ iξ 3 jξ 4 kξ 5 ξ 1 where i,j,k is the standard basis of quaternions. We remark that the Laplace transform of the exceptional hyperbolic ensemble is ( ( α, (3.7) E(e θ X ) = J (β 1)/ tr θ 4detθ)) defined for all θ with eigenvalues θ 1 θ < l β, the first positive zero of J (β 1)/. To see this, we use (3.0) (or its appropriate modifications for β < 4) with θ X = s 0 ξ 0 + β+1 r=1 s rξ r. From (3.5) with s = s 0 and t = (s 1,s,...,s β+1 ) we get (3.8) E(e β+1 j=0 s jξ j ) = ( J (β 1)/ ( t ) ) α, whichgives(3.7). (Wealsousedtheelementary observation θ 1 θ = tr θ 4detθ.) Proposition 3.1. The exceptional hyperbolic ensemble with parameter α > 0 is a Meixner ensemble with parameters 1 4α A =, B = 0, C = 1+α α+1. The mean and the variance are E(X) = 0, Var(X) = αi. For the standardized version, (1.4) holds with a = 1/(4α) and b = 0.

19 MEIXNER MATRIX ENSEMBLES 19 Our proof uses the system of PDEs derived in Theorem 4.1, so it appears after Proposition The system of PDEs Inthis section, we derive asystem ofpdes forthelaplacetransform of general Meixner ensembles. As previously, we consider simultaneously the real, complex, and quaternionic cases. We introduce some notation. For θ H n,β and j = 0,1,..., we consider σ j (θ) defined by n (4.1) det(i n +zθ) = σ j (θ)z j. Recall that for β = 4, the determinants of quaternionic matrices are defined only for hermitian matrices, see [17, page 9] and [14]. Note that σ 0 = 1, σ 1 (θ) = trθ, σ (θ) = 1 (trθ) 1 tr(θ ), σ n (θ) = detθ, and that σ j = 0 for j > n. We also adopt the convention that σ j = 0 for j < 0. In general, σ j (θ) = e j (θ 1,θ,...) is the j-th elementary symmetric function of the eigenvalues θ 1,θ,... of θ H n,β ; this notation and most of our calculations do not depend on the dimension n. Denote by U the open subset of R n such that if (σ 1,...,σ n ) U then the polynomial in x defined by n i=0 ( 1)n i σ i x n i has only distinct real roots. Themainresult ofthissectionisthesystem ofpdesthatdetermines the Laplace transform of a Meixner ensemble on a nonempty open set Θ X H n,β. This system is written in terms of an auxiliary function g(σ 1,...,σ n ) which is defined on an appropriate non-empty open set U X U. The link between L and this auxiliary function g varies according to the fact that a or b are zero or not. To define the set U X, we need Θ X to be closed under conjugation (θ Θ X implies UθU Θ X for all U K nβ ), and to have 0 in its closure. Consider the n- dimensional real subspace D n H n,β consisting of diagonal matrices, and let Dn 0 D n denote the set of matrices with distinct eigenvalues. Since Θ X is invariant under rotations, Θ X D n is a nonempty open set with 0 in its closure. The set Θ X Dn 0 is obtained from Θ X D n by removing a finite number of hyperplanes, so it is also non-empty, and has 0 in its closure. Denoteby σ : H n,β R n the mapping θ (σ 1 (θ),...,σ n (θ)). If θ Dn 0 then by the Vandermonde s determinant, B θ = {I n,θ,θ,...,θ n 1 } is a basis of D n. From formula (4.10) below, we see that in that basis, j=0

20 0 W LODZIMIERZ BRYC AND GÉRARD LETAC the matrix representation of the derivative of σ restricted to D n is triangular at θ Dn 0, with σ 0 = 1 on the diagonal, σ σ 1 (θ) σ [σ (θ)] Bθ = σ (θ) σ 1 (θ) σ ( 1) n 1 σ n 1 (θ) ( 1) n σ n (θ) σ 0 Therefore, the Jacobian of σ is 1 and σ : D 0 n Rn is an open mapping. Thus, the set (4.) U X = σ(θ X D 0 n ) is open, non-empty, and since σ(0) = 0, it has 0 R n in its closure. To shorten the formulas, for i,j = 1,,...,n we write g i = g σ i and g ij = g σ i σ j. For g : U X R consider the vector n g 11 σ r+s g rs (n 1) β g r,s= n σ 1 g +g 1 σ r+s 3 g rs (n ) β g 3 r,s=3. j n (4.3) D(g) = σ r+s 1 j g rs σ r+s 1 j g rs (n j) β g j+1 r,s=1 r,s=j+1. n 1 σ r+s n g rs σ n g nn β g n r,s=1 n σ r+s n 1 g rs r,s=1 Recall the convention that σ j = 0 for j < 0 or j > n. Theorem 4.1. Let X be a random variable with values in H n,β with the law invariant under rotations. Let L(θ) = e k(θ) = E(e θ X ) be its Laplace transform, and assume that Θ X 0. Suppose that E(X) = 0, E(X ) = I n and that X satisfies (1.4) with some a and b. Let U X be the associated open subset of R n, see (4.).

21 MEIXNER MATRIX ENSEMBLES 1 (i) If a 0, defineg : U X Rbyg(σ 1 (θ),...,σ n (θ)) = exp( 4ak(θ) btrθ). Then g satisfies the system of PDEs (b 4a)g D(g) = 0. 0 Furthermore, (4.4) g(σ 1,...,σ n ) 1, g(σ 1,...,σ n ) σ 1 b as (σ 1,...,σ n ) 0 over U X. (ii) If a = 0,b 0, define g : U X R by g(σ 1 (θ),...,σ n (θ)) = k(θ)+ 1 trθ. Then g satisfies the system of PDEs b Furthermore, (4.5) g(σ 1,...,σ n ) 0, g 1 g D(g) = b. g n g(σ 1,...,σ n ) σ 1 1 b as (σ 1,...,σ n ) 0 over U X. (iii) If a = 0,b = 0, define g : U X R by g(σ 1 (θ),...,σ n (θ)) = k(θ). Then g satisfies the system of PDEs 1 D(g) = 0. 0 Furthermore, (4.6) g(σ 1,...,σ n ) 0, as (σ 1,...,σ n ) 0 over U X. g(σ 1,...,σ n ) σ 1 0 Conversely, let µ be an ensemble on H n,β such that Θ µ. If one can find constants a,b that lead to one of the equations in (i), (ii) or (iii), then the corresponding version of equation (.7) holds for all θ Θ µ, and hence µ is a Meixner ensemble with parameters a,b.

22 W LODZIMIERZ BRYC AND GÉRARD LETAC The left hand side D(g) of our systems of the PDE s resembles the systemofpde sforthebesselfunctionofmatrixargumentgivenin[3, Eqtn (18)] when β = 1, see also [37, page 75 Eqtn (31)]. However, as we will see in Proposition 5.4, in contrast to [3] our system has multiple analytic solutions Proof of Theorem 4.1. The Cayley Hamilton theorem (for the case β = 4, see [17, Proposition II..1]) implies (4.7) θ n = n ( 1) i 1 σ i (θ)θ n i. i=1 We use also without proof the Newton formula [7, X.4] m 1 (4.8) mσ m (θ) = ( 1) m 1 i σ i (θ)tr(θ m i ); i=0 we remark that under the convention that σ m = 0 for m > n, formula (4.8) holds true for all m 1, i.e. also for m > n. Since tri n = n, it can be rewritten as m (4.9) (m n)σ m (θ) = ( 1) m 1 i σ i (θ)tr(θ m i ). i=0 We also need a differentiation formula. Proposition 4.. Fix n 1. Then for all m 1, m 1 (4.10) σ m(θ) = ( 1) m 1 i σ i (θ)θ m 1 i. i=0 Proof. We prove (4.10) by induction with respect to m. Since σ 1 = I n = σ 0 θ 0, the formula holds true for m = 1. Suppose (4.10) holds for some m 1. With t j = tr(θ j ), we differentiate Newton s formula (4.8) written for σ m+1 and use the induction assumption. We get (4.11) (m+1)σ m+1 = ( 1)m (m+1)θ m m ( j 1 ) + ( 1) m j (m+1 j)σ j θ m j +t m+1 j ( 1) j 1 s σ s θ j 1 s j=1 = ( 1) m (m+1)θ m +(m+1) s=0 m ( 1) m j σ j θ m j +R, j=1

23 MEIXNER MATRIX ENSEMBLES 3 where R = = j 1 m m ( 1) m j jσ j θ m j + ( 1) m j t m+1 j ( 1) j 1 s σ s θ j 1 s j=1 m ( 1) m j jσ j θ m j + j=1 j=1 s=0 m s 1 ( 1) m s θ m s ( 1) s 1 j σ j t s j. s=1 UsingNewton sformulaagain,weseethat s 1 j=0 ( 1)s 1 j σ j t s j = sσ s, so R = 0 and the formula follows from (4.11). We now consider separately the following cases, in which we re-write (.7) by the indicated substitutions: j=0 (i) If a 0, then f(θ) = exp( 4ak(θ) btr(θ)) solves (4.1) Ψ(f (θ))(i n ) = (b 4a)f(θ)I n. The limits as θ 0 in Θ X are: (4.13) f(θ) 1, f (θ) bi n. (ii) If a = 0,b 0, then f(θ) = k(θ)+ 1 tr(θ) solves b (4.14) Ψ(f (θ))(i n ) = bf (θ)i n. The limits as θ 0 in Θ X are: (4.15) f(θ) 0, f (θ) 1 b I n. (iii) If a = 0,b = 0, then f(θ) = k(θ) solves (4.16) Ψ(f (θ))(i n ) = I n. The limits as θ 0 in Θ X are: (4.17) f(θ) 0, f (θ) 0. The remaining part of the proof consists of finding a suitable expression for f. Now the fact that the law of X is invariant under the rotations is equivalent to k being invariant under the substitution θ UθU for all U K n,β, and is equivalent to saying that in each of the three cases mentioned above, there exists a real analytic function g on U X such that f(θ) = g(σ 1 (θ),...,σ n (θ))

24 4 W LODZIMIERZ BRYC AND GÉRARD LETAC for all θ Θ X with distinct eigenvalues. (This is Weyl s formula, see [34, Lemma.6.1].) Using (4.10) we get (4.18) f (θ) = n n g m σ m (θ) = m=1 m=1 ( 1) m 1 i σ i (θ)θ m 1 i. g m m 1 The second derivative of f is trickier: we write f = A+B with n n A = g mi σ m (θ) σ i (θ), B = m=1 i=1 n g m σ m (θ). m=1 i=0 The calculation of Ψ(A)(I n ) is quick and gives n n Ψ(A)(I n ) = g mi σ mσ i. m=1 i=1 For calculating Ψ(B)(I n ) we prove the surprising formula: Proposition 4.3. Ψ(σ m)(i n ) = β (m 1 n)σ m 1. Proof. We split σ m (θ) into σ m (θ) = C m +D m where (4.19) (4.0) C m = = D m = m 1 i 1 ( 1) m 1 i i=0 ( 1) i 1 j σ j θ i 1 j θ m 1 i j=0 m m j ( 1) m j σ j θ s θ m j s, j=0 s=0 m 1 ( 1) m 1 j σ j [θ m 1 j ]. j=0 Since z y+y x = (x+y) (x+y) x x y y, pairing up the appropriate powers from Proposition A we get ( m j (4.1) Ψ θ s θ )(I m j s n ) = (m j 1)θ m j, s=0 which leads to the calculation of Ψ(C m )(I n ) : m (4.) Ψ(C m )(I n ) = ( 1) m j σ j (m j 1)θ m j. j=0

25 MEIXNER MATRIX ENSEMBLES 5 Having in mind the calculation of Ψ(D m )(I n ) we write s 1 [θ s ] (h) = θ i hθ s 1 i = 1 s 1 [P θ i +θ s 1 i P θ i P θ s 1 i](h). i=0 i=0 From Proposition A, we have Ψ(P x ) = ( 1 ) β Px + β x x which leads to Ψ(P x )(I n ) = (1 β )x + β xtr(x) and to Ψ(P x+y P x P y )(I n ) = ( 1 β ) (xy +yx)+β xtr(y)+ytr(x). Thus we can compute Ψ([θ s ] )(I n ) = 1 = ( 1 β ) s 1 [θ i θ s i 1 +θ s i 1 θ i ]+ β 4 i=0 ( 1 β ) sθ s 1 + β s 1 θ i tr(θ s 1 i ). i=0 s 1 [θ i tr(θ s 1 i )+θ s 1 i tr(θ i )] i=0 We apply this to s = m 1 j and compute Ψ(D m )(I n ) : (4.3) Ψ(D m )(I n ) = ( 1 β ) m ( 1) m 1 j σ j (m 1 j)θ m j j=0 + β m m j ( 1) m 1 j σ j θ i tr(θ m j i ). j=0 i=0

26 6 W LODZIMIERZ BRYC AND GÉRARD LETAC Note that we replaced m 1 by m in the first sum of the right hand side. We obtain Ψ(σ m )(I n) = (4.4) (4.5) (4.6) = = m β ( 1) m j σ j (m 1 j)θ m j j=0 + β m m i ( 1) i θ i ( 1) m 1 j i σ j tr(θ m i j ) i=0 j=0 m β ( 1) m j σ j (m 1 j)θ m j j=0 + β m ( 1) i θ i (m i n)σ m i i=0 m β ( 1) m j σ j (m 1 j)θ m j j=0 + β m ( 1) m j θ m j (j n)σ j j=0 = β m 1 n m ( 1) m j σ j θ m j j=0 = β m 1 n σ m 1. In this chain of equalities, the first is obtained by combining (4.) and (4.3) and reversing the summation in i and j in the last sum; (4.4) is obtained by applying (4.9) (replacing m by m i); (4.5) is obtained by changing the variable j = m i. The last equality (4.6) uses the equality (4.10). Thus the Proposition 4.3 is proved. Gathering all terms, using the change of variable m = j 1 in Ψ(B) and the fact that σ 0 = 0, the expression Ψ(f )(I n ) = Ψ(A + B)(I n ) becomes, Ψ(f )(I n ) = n m=1 i=1 n g mi σ mσ i + β n 1 (j n)g j+1 σ j. j=1

27 MEIXNER MATRIX ENSEMBLES 7 Since σ 1 (θ) = I n, the three equations (4.1),(4.14) and (4.16) respectively become n n (4.7) g mi σ m σ i + β n 1 (j n)g j+1 σ j = (b 4a)gσ 1 (θ), m=1 i=1 j=1 (4.8) n n g mi σ m σ i + β n 1 n (j n)g j+1 σ j = b g m σ m (θ), m=1 i=1 j=1 m=1 (4.9) n m=1 i=1 n g mi σ m σ i + β n 1 (j n)g j+1 σ j = σ 1 (θ). Noting that σ 1,σ,...,σ n is a basis of the algebra generated by θ with distinct eigenvalues, the corresponding systems of PDEs arise by comparing the coefficients at σ j in ( ), respectively. To end the proof, we only need to find the coefficients P j (m,i) in the basis expansion n (4.30) σ m σ i = P j (m,i)σ j. j=1 j=1 This is accomplished by the following formula. Proposition 4.4. For all r,s 1, r 1 (4.31) σ r (θ)σ s (θ) = [σ j,σ r+s 1 j ], where j=0 [f,g] := f(θ)g (θ) f (θ)g(θ). Proof. From (4.10) we get σ j+1 = σ jθ 0 θσ j. Since the algebra generatedby θ is commutative, σ j+1 σ n σ j σ n+1 = (σ jθ 0 θσ j )σ n σ j (σ nθ 0 θσ n) simplifies to (4.3) σ j+1σ n σ jσ n+1 = [σ j,σ n ]. We now prove the proposition by induction. We assume that (4.31) holds for some fixed r and all s 1. This is trivially true for r = 1 since σ 1 σ s = σ s = [σ 0,σ s ]. If (4.31) holds true for some r 1 and all s 1, then by (4.3) and the induction assumption, r 1 σ r+1σ s = σ rσ s+1 +[σ r,σ s ] = [σ j,σ r+s j ]+[σ r,σ s ], j=0 which shows that (4.31) holds for r +1 and all s 1.

28 8 W LODZIMIERZ BRYC AND GÉRARD LETAC Proof of Theorem 4.1. Cases(i)-(iii) correspond to equations( ). The left hand sides of the resulting three equations ( ) are then re-written using (4.30). From (4.31) we see that coefficients P j (m,i) are linear with respect to σ 1,σ,...,σ n. Rewriting (4.31) as j r+s 1 P j (r,s)σ j = j=s r 1 σ r+s 1 j σ j σ r+s 1 j σ j, we get the explicit formula (4.33) σ r+s 1 j if max{r,s} j r +s 1 n, P j (r,s) = σ r+s 1 j if 1 j < min{r,s} and j r +s 1 n, 0 otherwise. With this notation equations( ) become the equations asserted in Theorem 4.1. To derive the corresponding boundary limits, we use (4.18) and the respective limits (4.13), (4.15), or (4.17). It is also clear that these calculations can be reversed. From the equations in listed in (i) (iii), we see that (.7) holds for all θ Θ X with distinct eigenvalues. By uniqueness of analytic extension, (.7) extends to the whole open set Θ X. Remark 4.1. Note that to have a solution of the system in each of the three cases is not a guarantee of existence of an ensemble: the function L(θ) still have to be a Laplace transform of a probability measure. For example, for n = there are solutions of the system for all a < 0, but according to Theorem 5.1(i) they are the Laplace transforms of a matrix ensemble only if 1/(4a) N. Remark 4.. Particular Meixner ensemble can be found from the system of PDEs by restricting g to be a function of σ 1,σ only. In other worlds, we lookforanensemble such thatthelaplacetransforml(θ) is a function of trθ, trθ. Such strategy leads us to the gamma ensemble on H n,β mentioned in Remark 3.3 and to Definition 3.5 of the Gaussian ensemble. This strategy leads only to trivial ensembles for b 4a. j=1 5. Meixner ensembles on matrices The goal of this section is to use Theorem 4.1 to prove that the constructions from Section 3 exhaust all Meixner ensembles on H,β. Theorem 5.1. Up to affine transformations, the only Meixner ensembles on H,β with the Laplace transform defined on a non-empty open subset of H,β are of the six types, depending on the values of a,b from (1.4):

29 MEIXNER MATRIX ENSEMBLES 9 (i) If µis a Meixner ensembleon H,β with a < 0 andθ µ, then 1/(4a) = N N and µ is a binomial ensemble Bin(N,q 1,q ) from Definition 3. for some q 1,q. (ii) If µ is a Meixner ensemble on H,β with a = 0, b > 0 and Θ µ, then µ is a Poisson ensemble from Definition 3.3 for some λ > 0. (iii) If µ is a Meixner ensemble on H,β with b > 4a > 0 and Θ µ, then µ is one of the negative binomial ensembles from Definition 3.4. (iv) If µ is a Meixner ensemble on H,β with a = b = 0 and Θ µ, then µ is one of the Gaussian ensembles from Definition 3.5. (v) If µ is a Meixner ensemble on H,β with b = 4a > 0 and Θ µ, then µ is one of the gamma ensembles from Definition 3.6. (vi) If µ is a Meixner ensemble on H,β with b < 4a and Θ µ, then either µ is the trivial ensemble of the form ξi or b = 0 and µ is the exceptional hyperbolic Meixner ensemble from Definition 3.7. (The seventh case of non-random ensembles arises from degenerate affine transformations.) 5.1. Laplace transforms of binomial, negative binomial and Poisson ensembles. Lemma 5.. If P is a random projection in Bin(1,1,0), then with σ j = σ j (θ), we have ) (5.1) E(e θ P ) = e σ 1/ I (β 1)/ ( σ 1 4 σ where I ν is a Bessel function (3.4). [ X1 ] Proof. To construct P, we choose a direction at random by taking X two independent standard normal random variables as the components ofthevector. NotethatthismeansthatX 1 = Z 1 +iz inthecaseβ = or X 1 = Z 1 +iz +jz +kz 4 in the case β = 4, where Z 1,Z,Z 3,Z 4 are the standard real-valued independent normal random variables. The matrix representation of P is [ ] 1 X1 P = X X 1 X 1 + X X 1 X X.,

30 30 W LODZIMIERZ BRYC AND GÉRARD LETAC Due to invariance [ under ] rotations, without loss of generality we take θ1 0 diagonal θ =, so that 0 θ θ P = tr(θp) = 1 X 1 + X ( θ1 X 1 +θ X ) = θ 1 U+θ (1 U). Since X 1, X are independent gamma distributed with shape parameterβ/andscaleparameter1/,thedistributionofu = X 1 /( X 1 + X ) is Beta(β/,β/). This gives E(e θ P ) = Γ(β) Γ (β/) 1 0 e uθ 1+(1 u)θ u β/ 1 (1 u) β/ 1 du. Substituting u = (1+v)/ so that u(1 u) = (1 v )/4 we get (5.) L P (θ) = e θ 1/+θ / Γ(β) β 1 Γ (β/) 1 1 e (θ 1 θ )v/ (1 v ) β/ 1 dv. The integral is expressed in terms of the Bessel functions in [44, 3.71 (9)]. Remark 5.1. Denoting by θ 1,θ R the eigenvalues of θ, I 0 (θ 1 θ ) if β = 1, E(e θ R sinh(θ 1 θ ) ) = θ 1 θ if β =, 3 ( ) cosh(θ1 θ ) (θ 1 θ sinh(θ 1 θ ) ) (θ 1 θ ) 3 if β = 4. Combining this with the Laplace transforms derived in the proofs of Propositions , and 3.7, we get. Corollary 5.3. (i) If X has Bin(N,q 1,q ) law, then its Laplace transform is well defined for all θ H β, and is given by (5.3) ( ) E(e θ X ) = e Nσ 1/ q 0 e σ 1/ +q 1 I (β 1)/ ( σ 1 4 σ +q e σ 1/) N, where q 0 = 1 q 1 q. (ii) If X has Poiss(λ 1,λ ) law, then its Laplace transform is well defined for all θ H β, and is given by ( ) ) (5.4) E(e θ X ) = exp where λ = λ 1 +λ. λ 1 e σ 1/ I (β 1)/ ( σ 1 4 σ +λ e σ 1 λ,

31 MEIXNER MATRIX ENSEMBLES 31 (iii) If X has NB(r,q 1,q ) law, then its Laplace transform is well defined for θ in a neighborhood of 0 in H,β, and is given by (5.5) E(e θ X ) = ( 1 (1 (q 1 +q )) ( (q r ) 1 e σ 1/ σ1 I (β 1)/ 4 σ +q e σ 1 )) r. Remark 5.. A calculation shows that for a given b > 4a > 0, parameter p = 1 q 1 q in (5.5) can range over the set {p (0,1) : p/( p) κ/b p}, while q = (bp κ)/κ [0,1 p] and q 1 = 1 p q are uniquely determined by p. (Here κ = b 4a, b > 0.) 5.. Solutions of the system of PDEs for n =. We consider in more detail the case n = with (5.6) U = {(σ 1,σ ) : 4σ < σ 1 }. We require that the solutions extend by continuity to σ 1 = σ = Case a 0 (Theorem 4.1(i)). The system of PDEs simplifies to: (5.7) (5.8) g 11 σ g = β g +(b 4a)g, g 1 +σ 1 g = 0. We will use as the initial conditions (5.9) g(0,0) = 1, Denote κ = b 4a. g(σ 1,0) σ 1 σ1 =0 = b. Proposition 5.4. Consider the system (5.7), (5.8) with initial conditions (5.9). (i) For b = 4a > 0 all solutions are g(σ 1,σ ) = 1 bσ 1 +C(βσ 1 +4σ ), where C is an arbitrary constant. Accordingly, (5.10) L(θ) = exp( σ 1 (θ)/b) (1 bσ 1 (θ)+c(βσ 1 (θ)+4σ (θ))) 1/b. (ii) For b > 4a all solutions are (5.11) g(σ 1,σ ) = C 1 e κσ 1 +C e κσ 1 +C 3 I (β 1)/ (κ σ1/4 σ ),