arXiv:1705.03510v1 [math.PR] 9 May 2017

The middle-scale asymptotics of
Wishart matrices
Didier Chételat and Martin T. Wells
Department of Applied Mathematics
and Industrial Engineering
École Polytechnique
Montréal, Québec H3T 1J4
Canada
e-mail: didier.chetelat@polymtl.ca
Department of Statistical Science
Cornell University
1198 Comstock Hall
Ithaca, NY 14853-3801
USA
e-mail: mtw1@cornell.edu
Abstract: We study the behavior of a real p-dimensional Wishart random
matrix with n degrees of freedom when n, p Ñ 8 but p{n Ñ 0. We establish
the existence of phase transitions when p grows at the order npK`1q{pK`3q
for every k P N, and derive expressions for approximating densities between
every two phase transitions. To do this, we make use of a novel tool we call
the G-transform of a distribution, which is closely related to the characteristic function. We also derive an extension of the t-distribution to the
real symmetric matrices, which naturally appears as the conjugate distribution to the Wishart under a G-transformation, and show its empirical
spectral distribution obeys a semicircle law when p{n Ñ 0. Finally, we discuss how the phase transitions of the Wishart distribution might originate
from changes in rates of convergence of symmetric t statistics.
MSC 2010 subject classifications: Primary 60B20, 60B10; secondary
60E10.

1. Introduction
The roots of random matrix theory lies in statistics, with the work of Wishart
[1928] and Bartlett [1933], and in numerical analysis, with the work of Von Neumann and Goldstine [1947]. In this early period, many well-known matrix distributions were introduced. This includes the real Gaussian matrix ensemble
Gpp, qq, a p ˆ q matrix with independent standard Gaussian entries, the Gaussian orthogonal
ensemble GOEppq, the distribution of a symmetric matrix pX `
?
X t q{ 2 with X „ Gpp, pq, and the Wishart (also known as Laguerre) distribution Wp pn, Ip {nq, the distribution of a symmetric matrix XX t {n with
X „ Gpp, nq. During that time, the main concern was to derive properties
of these distributions for a fixed dimension. Some asymptotics of the Wishart
distribution were considered, but only as n Ñ 8 for fixed p.
1

Chételat and Wells/Mid-scale Wishart asymptotics

2

Starting with the pioneering work of Wigner [1951, 1955, 1957], Porter and
Rosenzweig [1960], Gaudin [1961] and Mehta [1960a,b], researchers began investigating the asymptotics of Gaussian ensembles as their dimension grew to
infinity. As a result of decades of work, the behavior of a GOEppq matrix is now
well understood both in the classical setting where p is fixed, and in the setting
where p Ñ 8.
However, the situation asymptotics of the Wishart distribution is more complicated, as it depends on two parameters, n and p, and initial progress was slow.
The work of Marchenko and Pastur [1967] clearly established that the analogue
of a Gaussian orthogonal ensemble matrix whose dimension p grows to infinity
is a Wishart matrix whose degrees of freedom n and dimension p jointly grow
to infinity in such a way that p{n Ñ c P p0, 1q. Since then, we gained a very
good understanding of the behavior of Wishart matrices in this regime.
But this body of work left open the question as to what happens to a Wishart
matrix when n, p Ñ 8 with p{n Ñ 0. Since such asymptotics are middle-scale
between the classical regime where p is fixed as n Ñ 8 and the high-dimensional
regime where p{n Ñ c P p0, 1q, we might refer to them as middle-scale regimes.
Hence, we might ask: what is the asymptotic behavior of a Wishart matrix
Wp pn, Ip {nq in the middle-scale regimes? This question is addressed this article.
To gain some intuition, it is instructive to look at the eigenvalues λ1 ą ¨ ¨ ¨ ą
λp ą 0 of a Wp pn, Ip {nq Wishart matrix. In the classical regime where p is
fixed as n Ñ 8, the eigenvalues must all almost surely tend to 1 by the strong
law of large numbers. In constrast, in the high-dimensional regime where both
n, p Ñ 8 with p{n Ñ c P p0, 1q, the Marchenko-Pastur law states that for any
bounded, continuous f ,
a
ż c`
n
pc` ´ lqpc´ ´ lq
1 ÿ
dl
a.s.,
f pλi q Ñ
f plq
n i“1
2πcl
c´
?
where c˘ “ p1 ˘ cq2 . Thus the eigenvalues do not all tend to 1, but rather
distribute themselves in the shape of a Marchenko-Pastur law with parameter
c.
What happens between these two extremes? When c Ñ 0, the MarchenkoPastur law converges weakly to a Dirac measure with mass at 1. This suggests
that whenever n, p Ñ 8 with p{n Ñ 0
n
1 ÿ
f pλi q Ñ f p1q
n i“1

a.s.,

or in other words that the eigenvalues converge almost surely to 1, as in the
classical case.
This motivates a binary view of Wishart asymptotics. It appears that the
behavior of a Wishart matrix in the middle-scale regimes is the same as in
the classical regime, and therefore that there really are only two regimes: lowdimensional where p{n Ñ 0, and high-dimensional where p{n Ñ c P p0, 1q.
This binary view has very concrete repercussions. For example, in statistics,
many covariance matrix estimators have been developed that leverage high-

Chételat and Wells/Mid-scale Wishart asymptotics

3

dimensional Wishart asymptotics (see Pourahmadi [2013] for a review). When
faced with a problem where p is large with respect to n, it has been argued
that the high-dimensional asymptotics, rather than the classical, constitute the
correct model. The binary view provides a useful rule of thumb: small p’s call
for classical covariance estimators, while large p’s call for high-dimensional covariance estimators.
Unfortunately, recent results establish that this binary view is incorrect. In
the classical regime where p is fixed, the central limit theorem implies that
ı
? ”
n Wp pn, Ip {nq ´ Ip ñ GOEppq,
as n Ñ 8, where the arrow stands for weak convergence. In fact, something
better is known: recent work has extended this result to the case where p
tends to infinity. Recall that the total variation distance between two absolutely continuous distributions şF1 and F2 with densities f and g is given by
dTV pF1 , F2 q “ dTV pf1 , f2 q “ |f1 pxq ´ f2 pxq|dx. With different approaches,
Jiang and Li [2015] and Bubeck et al. [2016] independently established that
˙
ˆ ”
ı
?
n Wp pn, Ip {nq ´ Ip , GOEppq Ñ 0
dTV
whenever p3 {n Ñ 0. Thus, when p3 {n Ñ 0, the same asymptotics hold as in the
p fixed case, and we might regard these regimes as rightfully belonging to the
classical setting.
The surprising part is that the converse is true! When p3 {n Û 0, results of
Bubeck et al. [2016] and Rácz and Richey [2016] show that
ˆ ”
˙
ı
?
dTV
n Wp pn, Ip {nq ´ Ip , GOEppq Û 0.
?
Thus a phase transition occurs when p is of the order 3 n. This begs the?question:
if a normal approximation fails to hold when p grows faster than 3 n, what
asymptotics hold? Is there a uniform asymptotic behavior that holds whenever
p{n Ñ 0 with p3 {n Û 0, or are there further phase transitions as the growth
rate of p is increased?
The results of this paper offers a mostly complete answer to this question.
Namely, we establish that when p3 {n Û 0 but p2 {n Ñ 0,
ˆ ”
˙
ı
?
dTV
n Wp pn, Ip {nq ´ Ip , F1 Ñ 0,
where F1 is a continuous distribution on the space of real symmetric matrices
whose density is given when n ě 3p ´ 3 by
?
ˇ „
"
ˇ
2i
i
i
1
3
4
f1 pXq 9 ˇˇE exp ? tr XZ ´ ? tr Z `
tr Z ` 3{2 tr Z 5
4n
5n
8
3 2n
?
*ˇ2
ˇ
1
ipp ` 1q
p`1
4 2ipp ` 1q
6
2
3
ˇ ,
´ 2 tr Z ` ?
tr Z ´
tr Z ´
tr Z
(1.1)
ˇ
3{2
3n
4n
3n
2 2n

Chételat and Wells/Mid-scale Wishart asymptotics

4

Fig 1. Correct picture of Wishart asymptotics. This contrasts with the binary view, where no
phase transitions occur between p held constant and p growing like n.

?
for a Z „ GOEppq. When p grows like n, another phase transition occurs.
Namely, we establish that when p2 {n Û 0 but p5{3 {n Ñ 0,
˙
ˆ ”
ı
?
n Wp pn, Ip {nq ´ Ip , F2 Ñ 0,
dTV
where F2 is a continuous distribution on the space of real symmetric matrices
whose density is given when n ě 3p ´ 3 by
?
ˇ „
"
ˇ
i
i
1
2i
3
4
ˇ
f2 pXq 9 ˇE exp ? tr XZ ´ ? tr Z `
tr Z ` 3{2 tr Z 5
4n
5n
8
3 2n
?
?
ipp`1q
1
2 2i
p`1
4 2ipp`1q
7
6
2
tr Z 3
´ 2 tr Z ´ 5{2 tr Z ` ?
tr Z ´
tr Z ´
3n
4n
7n
3n3{2
2 2n
*ˇ2
ˇ
p`1
1024pp`1q
512pp`1q
5
4
6
ˇ ,
`
tr
Z
´
tr
Z
`
i
tr
Z
(1.2)
ˇ
4n2
3n3
5n5{2
again for a Z „ GOEppq.
In general, for every K P N we find a continuous distribution FK on the space
of real symmetric matrices, with density given when n ě 3p ´ 3 by
2K`3`

1rK
odds
ˇ «
#
ÿ
ˇ
´ 2 ¯ k2 tr Z k
i
n
ˇ
fK pXq 9 ˇE exp ? trpXZq `
ik
ˇ
4
n
k
8
k“3
2K`2´
1rK odds

p `1
`
4

ÿ

k

i

+ffˇ2
ˇ
ˇ
ˇ
ˇ
k

´ 2 ¯ k2 tr Z k
n

(1.3)

k“1

for a Z „ GOEppq, which approximates the normalized Wishart distribution in
some (but not all) middle-scale regimes. Namely, we prove the following, which
can be regarded as the main result of this paper.
Theorem 1.For any K P N,?the total variation distance between the the normalized Wishart distribution nrWp pn, Ip {nq ´ Ip s and the K th degree density
fK satisfies
ˆ ”
˙
ı
?
dTV
n Wp pn, Ip {nq ´ Ip , FK Ñ 0

Chételat and Wells/Mid-scale Wishart asymptotics

5

as n Ñ 8 with pK`3 {nK`1 Ñ 0.
The definition of fK and proof of Theorem 1 are found in Section 6, and
follow from definitions and results from Sections 3, 4 and 5 that constitute the
bulk of this paper.
The main consequence of this theorem is the existence of an infinite countable
number of phase transitions, occurring when p grows like npK`1q{pK`3q for K P
N. A diagram is provided at Figure 1. This naturally groups“ the middle-scale
˘
log p
K
K`1
regimes satisfying lim log
n ă 1 by which semi-open interval K`2 , K`3 their
nÑ8

log p
nÑ8 log n

limit lim

belongs to. We might refer to this grouping as the degree of the

log p
regime. In other words, we will say an middle-scale regime satisfying lim log
n ă
nÑ8
“ K K`1 ˘
log p
1 has degree K when lim log n P K`2 , K`3 .
nÑ8
The main result of this paper, Theorem 1, may then be summarized as saying
that the normalized Wishart distribution can be approximated by the distribution with density fK in every middle-scale regime of degree K or less. The
0th degree case corresponds to the classical setting, while the higher degrees
correspond to previously unknown behavior. In fact, we show that our 0th degree approximation F0 is asymptotically equivalent to the Gaussian orthogonal
ensemble. The results of this paper can therefore be regarded as a wide generalization of the Wishart asymptotics results of Jiang and Li [2015], Bubeck et al.
[2016], Bubeck and Ganguly [2016] and Rácz and Richey [2016].
Our approach relies on a novel technical tool we call the G-transform. It
turns out that to understand middle-scale regime behavior of Wishart matrices,
densities are less clear than characteristic functions (that is, Fourier transforms
of densities). Unfortunately, characteristic functions are difficult to relate to
metrics like the total variation distance. To remedy this problem, we develop
the G-transform and some associated theory in Section 3. An interesting aspect
of G-transform theory is that to every distribution we can associate a closely
related distribution called its G-conjugate. In fact, the G-conjugate of a Wishart
matrix is essentially a generalization of the t distribution to real symmetric
matrices. In Section 4, we define and derive several results concerning this new
distribution, including a semicircle law. From these results, we derive in Section
5 approximations to the Wishart distribution for middle-scale regimes of every
degree. Since these approximations are given using the language of G-transforms,
we derive in Section 6 density approximations, from which Theorem 1 follows.
We briefly discuss what concrete effects the phase transitions might have on
Wishart asymptotics in Section 7. Finally, we compile auxiliary results in Section
8, while we discuss in Section 9 open questions that arise from these results.
Although the results of this paper explain a large part of the behavior of
Wishart matrices when p{n Ñ 0, there exists regimes for which p{n Ñ 0 yet
log p
p R OpnpK`1q{pK`3q q for all K P N, or in other words for which lim log
n “ 1.
nÑ8

?

An example is when p grows at the order n1´1{ log n . Although the results of
our paper characterize almost all middle-scale regimes in the sense that among
log p
log p
those regimes satisfying lim log
n ď 1, those such that lim log n “ 1 represent
nÑ8

nÑ8

Chételat and Wells/Mid-scale Wishart asymptotics

6

a negligible set, they nonetheless exist. One might regard regimes such as those
as having infinite degree. Beyond this, however, it is difficult to say anything
about the behavior of Wishart matrices in these regimes. More work in that
direction is clearly needed.
2. Notation and definitions
The transpose of a matrix is denoted t , and the identity matrix of dimension
p is Ip . As is standard, we take the trace operator to have lower priority than
the power operator: thus for a matrix X, tr X k means the trace of X k . We will
write trk X when we mean the k th power of the trace of X. The Kronecker delta
is the symbol δkl “ 1rk “ ls.
The space of all real-valued symmetric matrices is denoted Sp pRq “ tX P
Mp pRq|X “ X t u. For a symmetric matrix X, we define the symmetric differentiation operator Bs {Bs Xkl by
#
1 B
for k ­“ l
1 ` δkl B
Bs
kl
“
“ 2 BX
B
Bs Xkl
2 BXkl
for k “ l.
BXkk
This operator has the elegant property that BSBXS kl trpXY q “ Ykl for any two
symmetric matrices X, Y .
The space of symmetric matrices Sp pRq can be assimilated to Rppp`1q{2 by
mapping a symmetric matrix to its upper triangle. By integration over Sp pRq,
we mean integration with respect to the pullback Lebesgue measure under this
isomorphism, that is,

ż

ż
f pXq dX “

Sp pRq

p
` ˘ź
f X
dXij .

Rppp`1q{2 iďj

We say a real symmetric matrix follows the Gaussian orthogonal ensemble
GOEppq distribution if Xkl , k ď l are all independent, with diagonal elements
Xkk „ Np0, 2q and off-diagonal elements Xkl „ Np0, 1q.
Let X be a n ˆ p matrix of i.i.d. Np0, 1q random variables, and let Σ be a
p ˆ p positive-definite matrix. The Wishart distribution Wp pn, Σq is the distri1
1
bution of the random matrix Σ 2 X t XΣ 2 . This is a special case of the matrix
gamma distribution. Following Gupta and Nagar [1999, Section 3.6], we say
a positive-definite matrix X has a matrix gamma distribution Gp pα, Σq with
shape parameter α ą pp ´ 1q{2 and scale parameter Σ if it has density over
Sp pRq given by
!
)
ˇ ˇα´ p`1
1
2
ˇX ˇ
exp ´ trpΣ´1 Xq 1rX ą 0s ,
f pXq “
α
|Σ| Γp pαq
where Γp is the multivariate gamma function. With this definition, the Wishart
distribution Wp pn, Σq is a matrix gamma with shape n2 and scale 2Σ.
While studying the Wishart distribution, the expression n ´ p ´ 1 comes up
so often that it makes sense to give it its own symbol. We will therefore write
m “ n ´ p ´ 1.

Chételat and Wells/Mid-scale Wishart asymptotics

7

The Hellinger distance is metric between absolutely continuous probability measures. For two distributions F1 and F2 with densities f1 and f2 , their
Hellinger distance is defined as
ˇ2 ¯ 12
´ż ˇ
ˇ
ˇ 1{2
1{2
HpF1 , F2 q “ Hpf1 , f2 q “
ˇf1 pxq ´ f2 pxqˇ dx .
The Hellinger distance is closely related to the total variation distance by the
inequalities
1
1{2
dTV pf1 , f2 q ď Hpf1 , f2 q ď dTV pf1 , f2 q.
(2.1)
2
In particular, Hpf1 , f2 q Ñ 0 if and only if dTV pf1 , f2 q Ñ 0. Thus they can
be seen as inducing the same topology on absolutely continuous probability
measures, called the strong topology, in contrast to the topology induced by
weak convergence of measures called the weak topology. One can show that if a
sequence of measures converges in the strong sense (i.e. in the dTV or H metrics),
then it converges weakly.
3. G-transforms
Our analysis of Wishart matrices relies heavily on a tool we call the G-transform
of a probability measure. To do so, we first need to define the Fourier transform
over symmetric matrices.
In Section 2, we clarified what we meant by integration over Sp pRq. For a
function f : Sp pRq Ñ C in L1 pSp pRqq, we define its Fourier transform to be

ż
Ftf upT q “

1
p

22 π

ppp`1q
4

e´i trpT Xq f pXq dX.

(3.1)

Sp pRq

It is more common to define
the Fourier
transform on symmetric(matrices with
(
ř
the integrand exp ´ i kďl Tkl Xkl , but choosing exp i trpT Xq considerably
simplifies our computations.
We extend this definition to f P Lr pSp pRqq, 1 ă r ď 2 in the usual manner. Because of the specific normalization chosen, this definition obeys a simple
version of Plancherel’s theorem, namely

ż

ż
f pXqs
g pXq dX “

Sp pRq

Ę q dT.
Ftf upT qFtgupT
Sp pRq

We now define the G-transform. In itself, the definition has nothing to do
with symmetric matrices and could have been perfectly well defined on any
other space endowed with a Fourier transform.
Definition 1 (G-transform of a density). Let f be an integrable function
Sp pRq Ñ C. Its G-transform is the complex-valued function Gtf u : Sp pRq Ñ C
defined by
Gtf u “ Ftf 1{2 u2 ,
where z 1{2 stands for the principal branch of the complex logarithm.

(3.2)

Chételat and Wells/Mid-scale Wishart asymptotics

8

In the same way that the Fourier transform maps L2 pSp pRqq to itself, the
G-transform maps L1 pSp pRqq to itself.
By extension, for an absolutely continuous distribution on Sp pRq with density
f , we will define its G-transform to be the G-transform of its density. (This
usage mirrors other transforms, such as the Stietjes transform.) We will usually
denote the G-transform of f by ψ. Since a density is integrable, this is always
well-defined. Moreover, f “ F ´1 tψ 1{2 u2 , so the density can be recovered from
the G-transform, and therefore to understand a distribution it is equivalent to
study its density or its G-transform.
Two comments are in order. First, for many densities, f 1{2 P L1 pSp pRqq. In
this case, the G-transform can be written explicitly as
ż
ˆ
˙2
1
´i trpT Xq 1{2
(3.3)
e
f pXq dX .
ψpT q “ Gtf upT q “
ppp`1q
2p π 2
Sp pRq
Second, throughout this article we will often talk about “the” square root of a
G-transform. To be clear, by ψ 1{2 we will always mean Ftf 1{2 u.
Now, in many ways, the G-transform behaves similarly to the characteristic
function (Fourier transform of a density), but it has unique features. First,
Plancherel’s theorem yields that

ż

ż
|ψ 1{2 pT q|2 dT

|ψpT q| dT “
Sp pRq

żS pRq
p

ż

|f 1{2 pXq|2 dX “

“

Sp pRq

|f pXq| dX “ 1.

(3.4)

Sp pRq

Thus |ψ| is itself a density, which we will call the G-conjugate of f . (In particular, ψ 1{2 is much like a quantum-mechanical wavefunction.) We will also
use an asterisk notation, so that the G-conjugate of a Np0, 1q distribution
will be denoted Np0, 1q˚ . For example, straightforward computations yield that
2
2
?1
Np0, 1q˚ “ Np0, 1{8q, χ2˚
n “ 8n tn{2 (where χν and tν are the univariate χ and
˚
´1 ˚
t distributions with ν degrees of freedom, respectively) and paF ` bq “ a F
for any distribution F and scalars a ­“ 0, b P R. Studying the G-conjugate of the
Wishart distribution will play a key part in deriving results about the Wishart
distribution itself. We should note that, in general, the double G-conjugate F ˚˚
is not the same as F . For example, χ2˚˚
is a density involving modified Bessel
n
functions of the first kind, not a χ2n .
A second feature that distinguishes G-transforms from characteristic functions is that they are easy to relate to the Hellinger distance between probability
measures. Consider two densities f1 , f2 with G-transforms ψ1 , ψ2 . By analogy,
we could define the “total variation” and “Hellinger” distances of ψ1 and ψ2
by

ż
dTV pψ1 , ψ2 q “

|ψ1 pT q ´ ψ2 pT q| dT,
Sp pRq

(3.5)

Chételat and Wells/Mid-scale Wishart asymptotics

ż
´
and

1{2

1{2

|ψ1 pT q ´ ψ2 pT q|2 dT

Hpψ1 , ψ2 q “

¯ 12

.

9

(3.6)

Sp pRq

Since the modulus of the G-transforms integrate to one, their total variation and
Hellinger distances are related to each other in the same way as in Equation 2.1
for densities, namely
1
1{2
dTV pψ1 , ψ2 q ď Hpψ1 , ψ2 q ď dTV pψ1 , ψ2 q.
(3.7)
2
Thus dTV pψ1 , ψ2 q Ñ 0 if and only if Hpψ1 , ψ2 q Ñ 0. But the Hellinger distance
between G-transforms is much more useful. Indeed, by the Plancherel theorem,
for any two densities f1 , f2 with G-transforms ψ1 , ψ2 , their Hellinger distance
satisfies
ż
H2 pf1 , f2 q “

1{2

1{2

|f1 pXq ´ f2 pXq|2 dX

żS pRq
p

“

1{2

1{2

|ψ1 pT q ´ ψ2 pT q|2 dT “ H2 pψ1 , ψ2 q.

(3.8)

Sp pRq

Thus to compute the Hellinger distance H2 pf1 , f2 q between two densities, we
can instead compute the Hellinger distance H2 pψ1 , ψ2 q of their G-transforms.
In contrast, there is no explicit way to express the Hellinger distance in terms
of characteristic functions. And no such connection exists between the total
variation distances of densities and G-transforms.
The G-transform does have some disadvantages compared to the Fourier
transform. It is a non-linear transformation (and therefore not a true transform),
and it does not behave well with respect to convolution. For our purposes,
however, the advantages listed above outweigh these problems.
In practice, it is not aways easy to control the Hellinger distance directly, and
one often focuses on the Kullback-Leibler divergence instead. The two quantities
are related through the well known inequality

„
f1 pXq
for X „ f1 .
H2 pf1 , f2 q ď E log
f2 pXq
For G-transforms, the following analog holds, which clarifies our interest in Gconjugates:
ˇ2
„

„ˇ
ˇ
ψ1 pT q
ψ1 pT q ˇˇ
H2 pψ1 , ψ2 q ď E < Log
` 2E ˇˇ= Log
for T „ F1˚ ,
ψ2 pT q
ψ2 pT q ˇ
where Log stands for the principal branch of the complex logarithm. In fact, in
this article we will need a further generalization, where ψ2 does not need to be
a G-transform of a density.
Proposition 1 (Kullback-Leibler inequality for G-transforms). Let ψ1 be the
G-transform of an absolutely continuous distribution F1 on Sp pRq, and let ψ2

Chételat and Wells/Mid-scale Wishart asymptotics

10

be an integrable function Sp pRq Ñ C. Then
ż

„

„
ψ1 pT q
H2 pψ1 , ψ2 q ď
|ψ2 |pT q dT ´ 1 ` E < Log
ψ2 pT q
Sp pRq
ż
ˇ 1
„ˇ
ˇ
1
ψ1 pT q ˇˇ 2
ˇ
2
`2
|ψ2 |pT q dT ¨ E ˇ= Log
ψ2 pT q ˇ
Sp pRq

for T „ F1˚ , where Log stands for the principal branch of the complex logarithm.
Proof. We can write

ż
1{2

H2 pψ1 , ψ2 q “

Sp pRq

„

1{2 1{2
´ ψ1 ψs2 dT

ż

ż

«


|ψ2 |pT q| ´ 1 ` 2 ´

“

1{2

|ψ1 |pT q ` |ψ2 |pT q| ´ ψs1 ψ2

1{2

ψ1 pT q

Sp pRq

1{2

1{2

ψ2 pT q

`

ψ̄2 pT q
1{2

ψ̄1 pT q

ff
|ψ1 |pT q dT

Sp pRq

ż
„

ż
„
|ψ2 |pT q| ´ 1 ` 2 1 ´

#



“

żS pRq

żS pRq

p

„

<
p

1{2

ψ2 pT q
1{2

ψ1 pT q

+


|ψ1 |pT q dT


„
|ψ2 |pT q| ´ 1 ` 2 1 ´

"
*
1
ψ1 pT q
exp ´ < Log
2
ψ2 pT q

Sp pRq

Sp pRq

“

˙

ψ1 pT q
1
= Log
|ψ1 |pT q dT .
2
ψ2 pT q
a
Now using the inequality ´ cospxq ď ´1 ` 2|x| that holds for any x P R. The
last quantity is bounded as
ˆ

¨ cos

ż
„

ż

„
|ψ2 |pT q| ´ 1 ` 2 1 ´

"
*

1
ψ1 pT q
exp ´ < Log
|ψ1 |pT q dT
2
ψ2 pT q

Sp pRq

Sp pRq

ď

ż
`2

ˇ
"
*dˇ
ˇ
1
ψ1 pT q
ψ1 pT q ˇˇ
ˇ
exp ´ < Log
ˇ= Log ψ2 pT q ˇ|ψ1 |pT q dT.
2
ψ2 pT q

Sp pRq

In the second term, use 1 ´ x ď ´ logpxq for x ě 0, while in the third term, use
the Cauchy-Schwarz inequality to obtain

ż

ż

„
ď


|ψ2 |pT q| ´ 1 ´ 2 log

"
*
ψ1 pT q
1
|ψ1 |pT q dT
exp ´ < Log
2
ψ2 pT q

Sp pRq

Sp pRq

Chételat and Wells/Mid-scale Wishart asymptotics

ż
`2

11

ż
"
*
1
ψ1 pT q
exp ´< Log
|ψ1 |pT q dT 2
ψ2 pT q

ˇ
1
ψ1 pT q ˇˇ
ˇ
ˇ= Log
ˇ|ψ1 |pT q dT 2 .
ψ2 pT q

Sp pRq

Sp pRq

Now use Jensen’s inequality in the second term and the algebraic identity
expt´< Log ψ1 pT q{ψ2 pT qu “ |ψ2 |pT q{|ψ1 |pT q in the third term to obtain

ż

ż

„
ď


|ψ2 |pT q| ´ 1 `

< Log

Sp pRq

Sp pRq

ψ1 pT q
|ψ1 |pT q dT
ψ2 pT q

ż

ż
|ψ2 |pT q dT 2

ˇ
1
ψ1 pT q ˇˇ
ˇ
ˇ= Log
ˇ|ψ1 |pT q dT 2 ,
ψ2 pT q

Sp pRq

Sp pRq

1

`2
as desired.

Let us now compute the G-transform of the Gaussian Orthogonal Ensemble
and the normalized Wishart distribution, which will be needed in our proofs.
The density of a GOEppq matrix over Sp pRq is
! 1
)
1
fGOE pXq “ ppp`3q ppp`1q exp ´ tr X 2 .
(3.9)
4
2 4 π 4
To compute its G-transform, we will make use of the fact that the elements
of a GOEppq matrix are independent to reduce the expression to a product of
characteristic functions.
Proposition 2. The G-transform of the Gaussian Orthogonal Ensemble density
on Sp pRq is
pp3p`1q
!
)
2 4
ψGOE pT q “ ppp`1q exp ´ 4 tr T 2 .
π 4
?
1{2
Proof. From Equation (3.9), fGOE is proportional to the density of the 2 GOEppq
distribution, so it is integrable. Therefore, we can apply Equation (3.3) to find
that

ż
1

1{2

ψGOE pT q “

p
2

2 π

ppp`1q
4

!
)
1{2
exp ´ i trpT Xq fGOE pXq dX

żS pRq
p

1

“
2

ppp`7q
8

π

3ppp`1q
8

!
)
1
exp ´ i trpT Xq ´ tr X 2 dX
8

żS pRq
p

1

“
2

ppp`7q
8

π

3ppp`1q
8

p
p
!
ÿ
1ÿ 2
exp ´ 2i
Tkl Xkl ´
X
4 kăl kl
kăl

Rppp`1q{2

´i

p
ÿ
k“1

Tkk Xkk ´

p
p
1 ÿ 2 )ź
Xkk
dXkl
8 k“1
kďl

Chételat and Wells/Mid-scale Wishart asymptotics

12

ż
“

2

pp3p`1q
8

π

ppp`1q
8

p
ź
kăl

) exp ´ 1 X 2 (
? 4 kl dXkl
exp ´ 2iTkl Xkl
4π
!

R

ż
p
ź

) exp ´ 1 X 2 (
? 8 kk dXkk
exp ´ iTkk Xkk
8π
!

¨
k“1

“

2

pp3p`1q
8

π

ppp`1q
8

R

p
p
”
!
)ı
”
! ?
)ı ź
ź
E exp ´ 2iTkk Z
E exp ´ 8iTkl Z
kăl

k“1

for Z „ Np0, 1q. The characteristic function of a Np0, 1q is expp´t2 {2q, so
pp3p`1q
p
p
!
)ź
!
)
ź
2 8
2
2
“ ppp`1q
exp ´ 4Tkl
exp ´ 2Tkk
π 8 kăl
k“1
pp3p`1q
!
)
2 8
“ ppp`1q exp ´ 2 tr T 2 .
π 8
Squaring this result yields the desired expression for ψGOE .
In particular, we see that |ψGOE | is the density of a GOEppq{4 distribution,
or in other words that GOEppq˚ “ GOEppq{4. In particular the Gaussian orthogonal ensemble is its own G-conjugate, up to a constant factor.
Let us now compute the G-transform of the normalized Wishart distribution.
Unlike the GOEppq case, the elements of the matrix are not independent, but
the elements of its Cholesky decomposition are. By being careful about complex
changes of variables, we can reduce the computation of the G-transform to the
computation of characteristic functions of the Cholesky elements.
Proposition?3. Let n ě p´2. Then the G-transform of the normalized Wishart
distribution nrWp pn, Ip {nq ´ Ip s density on Sp pRq is given by
ˇ´ n`p`1
"
*ˇ
2
ˇ
ˇ
?
4T
ψNW pT q “ Cn,p exp 2i n tr T ˇˇIp ` i ? ˇˇ
n
with
˘
`
ppn`2pq
Γ2p n`p`1
2 2
4˘
`
Cn,p “ ppp`1q ppp`1q
.
Γp n2
π 2 n 4
Proof. Recall the notation m “ n ´ p ´ 1 used throughout this article. The
density of a Y „ Wp pn, Ip {nq distribution is
np
! n
) m
n 2 1rY ą 0s
` n ˘ exp ´ tr Y |Y | 2 .
fW pY q “
np
2
2 2 Γp 2
?
?
If we do a change of variables X “ npY ´ Ip q, so that Y “ Ip ` X{ n and
p
p
ź
ź
1
dYij “ ppp`1q
dXij ,
n 4 iďj
iďj
?
we see that the normalized Wishart distribution nrWp pn, Ip {nq ´ Ip s has den-

Chételat and Wells/Mid-scale Wishart asymptotics

13

sity
„

ppn`mq
n 4
X
`
˘
?
1
I
`
ą
0
np
p
n
2 2 Γp n2

fNW pXq “

! n ”
X ˇˇ m2
X ı)ˇˇ
¨ exp ´ tr Ip `?
ˇIp `? ˇ . (3.10)
2
n
n
` 4 ´1 ˘
p`1
n`p`1
1{2
Notice that fW is proportional to expt´ tr r n Ip s Y u|Y | 4 ´ 2 , so it
`
˘
must be proportional to the density of a matrix gamma distribution Gp n`p`1
, n4 Ip
4
ą p´1
when n`p`1
4
2 , i.e. n ě p ´ 2. In particular, it must be integrable. As fNW
1{2
was obtained by a linear change of variables from fW , fNW must be integrable
too, that is

ż
1{2

fNW pXqdX ă 8.

(3.11)

Sp pRq

Therefore, we can apply Equation (3.3) to obtain

ż
1

1{2

ψN W pT q “

“

p

22 π

1
p

22 π

ppp`1q
4

ppn´2q
4

1

ppp`1q
4

!
)
1{2
exp ´ i trpT Xq fNW pXq dX
Sp pRq

”
!
)
ı
´1{2
E exp ´ i trpT Xq fNW pXq
`n˘

„
"
*

n ”
X ı ˇˇ
X ˇˇ´ m4
?
?
E
exp
´
i
trpT
Xq
`
tr
I
`
`
.
ˇI
ˇ
p
p
ppp`1q
ppn`mq
4
n
n
π 4 n 8
?
If we rewrite the expectation in terms of Y “ Ip ` X{ n, this last expression
equals
1 ` ˘
ppn´2q
*

„ "´
`?
˘
?
2 4 Γp2 n2
n ¯
´m
4
.
(3.12)
“ ppp`1q ppn`mq exp i n tr T E etr ´i nT ` Ip Y |Y |
4
π 4 n 8
Since T is real symmetric, there must be a spectral decomposition T “ ODOt
with O real orthogonal and D real diagonal. As Ot Y O has the same distribution
as Y , namely Wp pn, Ip {nq, we can rewrite Equation (3.12) as
1 ` ˘
ppn´2q
„ "´
*

˘
`?
?
2 4 Γp2 n2
n ¯
´m
4
“ ppp`1q ppn`mq exp i n tr T E etr ´i nD ` Ip Y |Y |
.
(3.13)
4
π 4 n 8
Now, since Y is positive-definite it has a Cholesky decomposition Y “ U t U
with U upper-triangular. According to Bartlett’s theorem [see Muirhead, 1982,
Theorem 3.2.14], all the elements of U are independent, the diagonal elements
2
have the distribution Ukk
„ χ2n´k`1 {n and the upper diagonal elements have
Ukl „ Np0, 1{nq for k ă l. Since
„´

p ´
ÿ
?
?
Ip ¯
n ¯
1
tr
´ i nD `
Y “
´ i nD ` Ip
Ukl
Ulj
4
4
jk
j,k,l
“

2

Γp2

2

Chételat and Wells/Mid-scale Wishart asymptotics
p ´
ÿ

“
lďk

14

?
n¯ 2
´ i nDkk `
U
4 lk

śp

2
k“1 Ukk , we have by independence and Equation (3.13) that
#
+ p „
1 ` ˘
ppn´2q
*
"´
ź
?
?
2 4 Γp2 n2
n¯ 2
“ ppp`1q ppn`mq exp i n tr T
U
E exp
´ i nDkk `
4 lk
π 4 n 8
lăk
*

„
"´
p
ź
?
n¯ 2
2 ´m
E exp
´ i nDkk `
¨
Ukk pUkk
q 4 . (3.14)
4
k“1

and |Y | “

We will now compute these expected values in several steps. For a given 1 ď
k ď p, let
„
"´

*
?
n¯ 2
´m
A “ E exp
´ i nDkk `
(3.15)
Tkk Tkk 2 .
4
2
Since Tkk
„ χ2n´k`1 {n and m “ n ´ p ´ 1, we have
ż8 "
*
n´k`1
´
2
?
m
n
n¯
¯
exp
´ i nDkk `
x x´ 4
A “ n´k`1 ´
n´k`1
4
2 2 Γ
0
2
! n )
n´k`1
¨ x 2 ´1 exp ´ x dx
2
ż8 "
*
n´k`1
´
¯
n´2k`p`3
n 2
n ?
´1
4
¯
“ n´k`1 ´
dx
(3.16)
exp ´
` nDkk i x x
n´k`1
4
2
Γ
2
0
2
Consider the truncated integrands
" ´
¯ * n´2k`p`3
n ?
´1
4
hM pxq “ exp ´
` nDkk i x x
1r0 ă x ă M s .
4
Clearly this sequence is dominated by the integrable positive function h,
"
*
n´2k`p`3
n
´1
4
,
|hM pxq| ď hpxq “ exp ´ x x
4
ˆ ˙ n´2k`p`3
˙
ˆ
ż8
4
4
n ´ 2k ` p ` 3
hpxqdx “
ă 8.
Γ
n
4
0
Therefore, by the Dominated Convergence Theorem and Equation (3.16),
żM "
n´k`1
´n ?
¯ * n´2k`p`3
2
n
´1
4
¯ lim
A “ n´k`1 ´
exp ´
` nDkk i x x
dx.
M
Ñ8
n´k`1
4
2 2 Γ
0
2
`n ?
˘
By the change of variables z “ 4 ` nDkk i x, this can be rewritten
“

n

n´k`1
2

?

˘´ n´2k`p`3
4
nDkk i
4 `
´
¯
lim
n´k`1
M Ñ8
2 2 Γ n´k`1
2

`n

ż

nM
4

?
` nDkk M i

e´z z

n´2k`p`3
´1
4

dz.

(3.17)

0

To compute this integral, we use a contour argument. Consider the closed path
nM
C “ C1 ` C2 ` C3 given by C1 a path from 0 to nM
4 , C2 a path from 4 to

Chételat and Wells/Mid-scale Wishart asymptotics
y

nM
4

15

?
` nDkk M i

C3
C2

C1
0

nM
4

x

Fig 2. Contour C “ C1 ` C2 ` C3 when Dkk ě 0. The diagram is mirrored around the x axis
when Dkk ă 0.

?
?
` nDkk M i and finally C3 a path from nM
nDkk M i to 0. A diagram
4 `
is provided as Figure 2.
n´2k`p`3
´1
4
As k ď p, z ÞÑ e´z z
is entire and its integral over C must be zero.
Therefore
?
ˇ
ˇ
nM
ˇ
ˇ
4 ` nDkk M i
ˇ
˙ˇ
ˆ
ˇ
ˇ
n´2k`p`3
n
´
2k
`
p
`
3
´1
ˇ lim
ˇ
4
e´z z
dz ´ Γ
ˇM Ñ8
ˇ
4
ˇ
ˇ
ˇ
ˇ
0
?
ˇ
ˇ
nM
nM
ˇ
ˇ
4 ` nDkk M i
4
ˇ
ˇ
ˇ
ˇ
n´2k`p`3
n´2k`p`3
´1
´1
4
4
“ lim ˇˇ
e´z z
dz ´
e´x x
dxˇˇ
M Ñ8 ˇ
ˇ
ˇ 0
ˇ
0
?
ˇ
ˇ
ˇ
ˇ
nM
ˇ
ˇ
ˇ
ˇ
4 ` nDkk M i
ˇ
ˇ
ˇ
ˇ
ˇ
ˇ
ˇ
ˇ
n´2k`p`3
n´2k`p`3
´z
´z
´1
´1
ˇ
ˇ
ˇ
4
4
e z
“ lim ˇ
e z
dz ˇ “ lim ˇ
dz ˇˇ .
M Ñ8 ˇ
M
Ñ8
ˇ
ˇ
ˇ
ˇ C2
ˇ
ˇ nM
ˇ
4
nM
?
z´ 4
?
Do a change of variables z “ nM
4 ` y nDkk M i, so that y “ nDkk M i is real on
the path. It yields
ˇ
ˇ
1
ˇ
ˇ
ˇ
„
 n´2k`p`3
´1 ˇ
4
?
ˇ
?
?
M ˇˇ
nM
“ n|Dkk | lim nM ˇ e´y nDkk M i
`y nDkk M i
dy ˇˇ
M Ñ8 e 4 ˇ
4
ˇ
ˇ 0
ˇ
nM
4

ż

ż

ż

ż

ż

ż

ż
ď

?

1

ˇ
ˇ n´2k`p`3
n´2k`p`3
´1
4
ˇn
ˇ
4
?
M
ˇ ` y nDkk iˇ
dy
¨
lim
.
nM
ˇ4
ˇ
M Ñ8
e 4

n|Dkk |
0

This last integral is finite, since it is continuous on a bounded interval. Therefore

Chételat and Wells/Mid-scale Wishart asymptotics

the limit is zero and by Equation (3.17) and the previous expression,
´
¯
n´k`1
´n ?
¯´ n´2k`p`3
n 2 Γ n´2k`p`3
4
4
´
¯
A “
.
` nDkk i
n´k`1
n´k`1
4
2 2 Γ
2

16

(3.18)

Going back to (3.14), let us now consider the expectations in the second
products. For fixed 1 ď l ă k ď p, let
„
"´
*
?
n¯ 2
´ i nDkk `
B “ E exp
T
.
(3.19)
4 lk
2
Since Tlk
„ χ21 {n,

ż8
c
B “

"´

n
2π

exp

* n
?
n ¯ e´ 2 x
´ i nDkk `
x ? dx
4
x

0

ż8
c
“

" ´
¯ * 1
?
n
exp ´
` i nDkk x ? dx.
4
x

n
2π

(3.20)

0

Consider the truncated integrands

" ´
¯ * 1 „ 1
?
n
` i nDkk x ? 1
ăxăM .
(3.21)
hM pxq “ exp ´
4
x M
We see that they are dominated by a positive, integrable function hpxq,
c
ż8
n
e´ 4 x
π
|hM pxq| ď hpxq “ ? ,
ă 8.
hpxqdx “ 2
x
n
0
Therefore, by the Dominated Convergence Theorem and Equation (3.20), we
conclude that
M
c
" ´
¯ * 1
?
n
n
B “
lim
exp ´
` i nDkk x ? dx.
2π M Ñ8
4
x
1{M
`
˘
?
A complex change of variables z “ n4 ` i nDkk x yields

ż

ż
c
“

¯´ 12
?
n ´n
` i nDkk
lim
M Ñ8
2π 4

nM
4

?
` nDkk M i

?
e´z { z dz.
?
nDkk
n
4M `
M

(3.22)

i

Let’s compute this integral again using a contour integration argument. Consider
n
the contour C “ C1 ` C2 ` C3 ` C4 given by C1 a line from 4M
to nM
2 a line
4 , C?
?
?
nM
nM
nM
n
kk
from 4 to 4 ` nDkk M i, C3 a line from 4 ` nDkk M i to 4M ` nD
i
M
?
nD
n
n
kk
and C4 a line from?4M ` M i to 4M . A diagram is provided as Figure 3.
Since z ÞÑ e´z { z is holomorphic away from zero,
?
?
ˇ
ˇ
ˇ
ˇ
nM
nM
nM
ˇ
ˇ
ˇ
ˇ
4 ` nDkk M i
4 ` nDkk M i
4
ˇ
ˇ
ˇ
ˇ
ˇ
ˇ
ˇ
ˇ
?
?
?
?
´z
´z
´x
ˇ lim
ˇ “ lim ˇ
ˇ
e
{
z
dz
´
π
e
{
z
dz
´
e
{
x
dx
ˇM Ñ8
ˇ M Ñ8 ˇ
ˇ
ˇ
ˇ
ˇ
ˇ
?
?
nDkk
nDkk
ˇ
ˇ
ˇ
ˇ
n
n
n
`
i
`
i

ż

4M

ż

ż

M

4M

M

4M

Chételat and Wells/Mid-scale Wishart asymptotics
y

nM
4

17

?
` nDkk M i

C3
n `
4M
?
nDkk
M

C2
i

C4
0

n
4M

C1

nM
4

x

Fig 3. Contour C “ C1 ` C2 ` C3 ` C4 when Dkk ě 0. The diagram is mirrored around the
x axis when Dkk ă 0.

ˇ
ˇ
ˇ
ˇ
“ lim ˇ
M Ñ8 ˇ
ˇ
ˇ
ˇ
ˇ
ˇ
ď lim ˇˇ
M Ñ8 ˇ
ˇ

ż

ż

ˇ
ˇ
ˇ
?
?
ˇ
e´z { z dz `
e´z { z dz ˇ
ˇ
ˇ
C2
C4
?
ˇ
ˇ
ˇ
nM
n
ˇ ˇ
ˇ
4 ` nDkk M i
4M
ˇ ˇ
ˇ
ˇ
ˇ
ˇ
?
?
e´z { z dz ˇˇ ` ˇˇ
e´z { z dz ˇˇ .
ˇ ˇ
ˇ
ˇ ˇ n ` ?nDkk i ˇ
nM
4
4M
M
?
?
nDkk
n
By changes of variables z “ nM
i on the
4 ` y nDkk M i and z “ 4M ` y M
two respective integrals, we get
ˇ 1
ˇ
ˇ
ˇ
?
ˇ
ˇ
nD
M
i
y
kk
?
M ˇ
e
ˇ
“ n|Dkk | lim nM ˇ
1 dy ˇ
ı
”
M Ñ8 e 4 ˇ
ˇ
?
2
ˇ 0 nM
ˇ
4 ` y nDkk M i
ˇ 1
ˇ
?
ˇ
ˇ
nDkk
ˇ
ˇ
i
y
?
1 ˇ
e M
ˇ
` n|Dkk | lim
dy
ˇ
ˇ
n
1
”
ı
?
M Ñ8 M e 4M ˇ
ˇ
2
n
kk
ˇ 0 4M
ˇ
` y nD
i
M

ż

ż

ż

ż

ż1
ď

?

n|Dkk |
0

?
1
M
dy
¨
lim
nM
ˇ
ˇ 12
M
Ñ8
?
ˇn
ˇ
e 4
ˇ 4 ` y nDkk iˇ

ż1
1
1
?
n .
ˇ
ˇ 12 dy ¨ Mlim
Ñ8
?
M e 4M
ˇn
ˇ
0 ˇ 4 ` y nDkk iˇ
1
ˇ
ˇ
´ 2
¯´ 14
?
ˇ
ˇ´ 2
2
“ n16 ` y 2 nDkk
is continuous on r0, 1s, a bounded
Since ˇ n4 ` y nDkk iˇ
interval, we conclude that the integrals are finite and that the limits are zero.
`

?

n|Dkk |

Chételat and Wells/Mid-scale Wishart asymptotics

Therefore, by Equation (3.22),
c ´
¯´ 12
?
n n
B “
` i nDkk
.
2 4

18

(3.23)

Recall the definitions of A and B at Equations (3.15) and (3.19). Combining
1{2

both Equations (3.18) and (3.23) into the expression for ψNW at Equation (3.14)
provides
1 ` ˘
ppn´2q
"
*ź
p c ´
¯´ 12
?
?
2 4 Γp2 n2
n n
1{2
ψNW pT q “ ppp`1q ppn`mq exp i n tr T
` i nDkk
2 4
π 4 n 8
lăk
¯
´
n´k`1
p n 2
¯´ n´2k`p`3
´n ?
Γ n´2k`p`3
ź
4
4
´
¯
¨
nD
i
`
kk
n´k`1
n´k`1
4
Γ
k“1 2 2
2
´
¯
1 ` ˘
ppn´2q
"
*
n´k`1
k´1
p Γ n´2k`p`3
p
n
2
`
ź
ź
4
2
?
2
Γp 2
4
n 2
`
˘
“ ppp`1q ppn`mq
exp
i
n
tr
T
n´k`1
` k´1
Γ n´k`1
2
2
π 4 n 8
2
k“1
k“1 2
p ´
n´2k`p`3
¯
ź
´
´ k´1
n ?
4
2
´ nDkk i
¨
4
k“1
´
¯
n´2k`p`3
ˇ n`p`1
"
*ˇ
pp2n`p`1q
p
Γ
ź
ˇn
8
?
? ˇ´ 4
4
n
` n´k`1 ˘ exp i n tr T ˇˇ Ip ` i nT ˇˇ
“ ppn`2q ppp`1q
.
4
2 4 π 4 k“1 Γ
2
But by Muirhead [1982, Theorem 2.1.12],
p
`
˘
´
¯
ppp´1q ś
˘
`
π 4
Γ n`p`1
` 1´k
p Γ n´2k`p`3
4
2
ź
Γp n`p`1
4
k“1
4
` ˘ ,
`
˘ “
p
˘ “ Γp n
`
ppp´1q ś
Γ n´k`1
2
2
k“1
π 4
Γ n2 ` 1´k
2
k“1

so by taking a n{4 factor out of the determinant, we find that
`
˘
ˇ´ n`p`1
"
*ˇ
ppn`2pq
4
ˇ
?
Γp n`p`1
2 4
4T ˇˇ
1{2
4
ˇ
?
ψNW pT q “ ppp`1q ppp`1q
exp
i
n
tr
T
I
`
i
.
`
˘
p
ˇ
1{2
nˇ
Γp n2
π 4 n 8
Squaring this result yields the desired expression for ψNW .
By Proposition 3, when n ě p ´ 2 the G-conjugate of a normalized Wishart
distribution must have a density on Sp pRq given by
`
˘ˇ
ˇ´ n`p`1
ppn`2pq
4
ˇ
Γ2p n`p`1
2 2
16T 2 ˇˇ
4
` n ˘ ˇˇIp `
|ψNW |pT q “ ppp`1q ppp`1q
.
(3.24)
ˇ
n
Γp 2
π 2 n 4
As mentioned in the paragraph following Equation (3.4), the G-conjugate of
a χ2n {n distribution is a scaled tn{2 . Thus, by analogy, Equation (3.24) should
be represent some kind of generalization of the t distribution to the real symmetric matrices. Matrix-variate generalizations of the t distribution have been
investigated in the past, but not for symmetric matrices. Hence it appears the
concept is new.

Chételat and Wells/Mid-scale Wishart asymptotics

19

This motivates us to propose in Section 4 a candidate for a symmetric matrix
variate t distribution. Using that definition, the G-conjugate to the normalized
Wishart could then be regarded as the t distribution with n{2 degrees of freedom
and scale matrix Ip {8, which we denote Tn{2 pIp {8q. But regardless of its name,
this distribution will play a key role in our results about the middle-scale regime
asymptotics of Wishart matrices, and will be investigated in depth in Section
4.
4. The symmetric matrix variate t distribution
In Section 3, Equation (3.24), we proved ?
that when n ě p ´ 2, the G-conjugate
of the normalized Wishart distribution nrWp pn, Ip {nq ´ Ip s has density on
Sp pRq given by
`
˘ˇ
ˇ n`p`1
ppn`2pq
2 ˇ´ 4
ˇ
Γ2p n`p`1
16T
2 2
4
ˇ
` ˘ ˇˇIp `
|ψNW |pT q “ ppp`1q ppp`1q
.
(4.1)
n ˇ
Γp n2
π 2 n 4
Two remarks are in order. First, we are unaware of any matrix calculus tools
that could let us integrate this expression directly. Thus, the mere fact that
this expression integrates to unity, a consequence of being the G-conjugate of
another distribution, seems remarkable.?
Second, when p “ 1, this is the tn{2 { 8 distribution. Thus, as we mentioned
while discussing Equation (3.24), it is natural to interpret this distribution as the
parametrization of some generalization of the t distribution to Sp pRq, the space
of real-valued symmetric matrices. The purpose of this section is to propose
a candidate definition for such generalization, as well as prove several results
concerning the normalized Wishart G-conjugate.
To the best of our knowledge, no extension of the t distribution to symmetric
matrices has ever been proposed. However, a non-symmetric matrix variate t
distribution has been thoroughly investigated in the literature – see Gupta and
Nagar [1999, Chapter 4] for a thorough summary. Several definitions exist. For
our purposes, we say that a p ˆ q real-valued random matrix T has the matrix
variate t distribution with ν degrees of freedom and q ˆ q positive-definite scale
matrix Ω if it has density
¯
´
ˇ
ˇ´ ν`p`q´1
Γp ν`p`q´1
2
pˇ
2
T Ω´1 T t ˇˇ
´
ˇ
2
´
¯
|Ω|
I
`
.
p
ˇ
ˇ
pq
pq
ν
ν 2 π 2 Γ ν`p´1
p

2

It is not exactly clear what should be the proper analog of this distribution
for symmetric matrices. But it would be elegant if the degrees of freedom of
Equation (4.1) were to be exactly n{2, as in the univariate case. Thus, the
following definition seems natural.
Definition 2 (Symmetric matrix variate t distribution). We say a real symmetric p ˆ p matrix T has the symmetric matrix variate t distribution with
ν ě p{2 ´ 1 degrees of freedom and p ˆ p positive-definite scale matrix Ω,

Chételat and Wells/Mid-scale Wishart asymptotics

20

denoted Tν pΩq, if it has density
ˇ´ ν`pp`1q{2
ˇ
2
ˇ
T Ω´1 T ˇˇ
ˇ
.
fTn pΩq pT q 9 ˇIp `
ˇ
ν
With this definition, the G-conjugate to the normalized Wishart distribution,
whose density is given by Equation (4.1), is the Tn{2 pIp {8q distribution on Sp pRq.
In fact, since Equation (4.1) integrates to one, we can deduce the normalization constant of Definition 2. For an arbitrary degrees of freedom parameter ν,
imagine the density |ψNW | of the G-conjugate of a normalized Wishart distribution with n “ 2ν ě p ´ 2. By virtue of being a G-conjugate,
it must integrate to
1 ?
1
unity. Then from the change of variables T “ Ω´ 4 SΩ´ 4 { 8 which has Jacobian
ppp`1q
p`1
dT “ 8´ 4 |Ω|´ 4 dS, we see that

ż
1“

ż
|ψNW |pT q dT “
Sp pRq

2ppν´1q Γ2p
“
π

ppp`1q
2

ν

´

1
2

3ppp`1q
4

ν`pp`1q{2
2

ppp`1q
4

|Ω|

p`1
4

ˆ ´1
1 ˙
Ω 4 SΩ´ 4
?
dS
|ψNW |
8
Sp pRq

ż

¯
|Ω|

´ p`1
4

Γp pνq

ˇ
ˇ ν`pp`1q{2
´1 ˇ´
2
ˇ
ˇIp ` SΩ S ˇ
dS.
ˇ
ˇ
ν
Sp pRq

Thus we must have
fTn pΩq pT q “

2ppν´1q Γ2p

´

ν`pp`1q{2
2

¯

ˇ
ˇ´ ν`pp`1q{2
2
ˇ
T Ω´1 T ˇˇ
I
`
.
p
ˇ
ˇ
ν

´ p`1
4 ˇ

|Ω|
(4.2)
ppp`1q
ppp`1q
π 2 ν 4 Γp pνq
It would be interesting to see if this distribution satisfies the properties we would
expect of a t distribution, to ensure our guess is the “correct” one. However, this
would take us too far away from the topic of this article. Instead, we will focus
in the rest of this section on proving results about Tn{2 pIp {8q, the G-conjugate
to the normalized Wishart distribution.
Our first result will concern the asymptotic expansion of its normalization
constant. We mention that this constant is the same as the Cn,p term appearing
in the expression of the G-transform of the normalized Wishart in Proposition
3.
Lemma 1. The normalization constant of the Tn{2 pIp {8q distribution
`
˘
ppn`2pq
Γ2p n`p`1
2 2
4˘
`
Cn,p “ ppp`1q ppp`1q
Γp n2
π 2 n 4
has, for every K P N, the asymptotic expansion
" K`1
pp3p`1q
2 4
1 ÿ 1rk evens pk`2
Cn,p “ ppp`1q exp ´
2 k“1 kpk`1qpk`2q nk
π 4

(4.3)

K`1
´ pK`3 ¯*
1 ÿ 1`21rk evens pk`1
´
` o K`1 .
4 k“1
kpk`1q
nk
n

as n Ñ 8 with p{n Ñ 0.

Chételat and Wells/Mid-scale Wishart asymptotics

21

Proof. By Stirling’s approximation applied to log Γ, as well as Muirhead [1982,
Theorem 2.1.12], we find that
˙
ˆ
p
ÿ
ppp ´ 1q
i´1
log Γp pxq “
log π `
log Γ x ´
4
2
i“1
„ˆ
˙
ˆ
˙
p
ÿ
ppp ´ 1q
i´1 1
i´1
“
log π `
x´
´
log x ´
4
2
2
2
i“1
˙
ˆ ˙
ˆ
1
1
i´1
` log 2π ` O
´ x´
2
2
x
ppp ` 1q
p
ppp ´ 1q
“
log π ` log 2 ´ px `
4
2
4
˙
ˆ
˙
p ˆ
´p¯
ÿ
i
i´1
`
log x ´
`O
x´
2
2
x
i“1
as x Ñ 8. Thus
˙
ˆ
´ n ¯ ppp ` 1q
p
n`p`1
´ log Γp
“
log π ` log 2
2 log Γp
4
2
4
2
ˆ
˙
ˆ
˙
p
ppp ` 3q ÿ n`p`1
n`p`1 i´1
´
`
´ i log
´
4
2
4
2
i“1
˙
ˆ
˙
ˆ
p
ÿ n
i
n i´1
´
log
´
` op1q
´
2
2
2
2
i“1
“

ppp ` 1q
pp2n ` p ´ 1q
ppp ` 1q
ppp ` 3q
log π ´
log 2 `
log n ´
4
4
4
4
ˆ
˙
p
¯
1 ÿ´
2i´p´3
`
n ´ r2i´p´1s log 1 ´
2 i“1
n
ˆ
˙
p
i´1
1ÿ
pn ´ iq log 1 ´
` op1q,
´
2 i“1
n

as n Ñ 8 with p{n Ñ 0, and so by Equation (4.3),
pp3p ` 1q
ppp ` 1q
ppp ` 3q
log Cn,p “
log 2 ´
log π ´
4
4
4
˙
p
¯ ˆ
1 ÿ´
2i´p´3
`
n ´ r2i´p´1s log 1 ´
2 i“1
n
ˆ
˙
p
1ÿ
i´1
´
pn ´ iq log 1 ´
` op1q.
(4.4)
2 i“1
n
Let us now focus on the two sums in this expression. Recall that for any k ě 1,
x2
x3
xk
´ logp1 ´ xq “ x `
`
` ¨¨¨ `
` Opxk`1 q as x Ñ 0,
(4.5)
2
3
k
even for negative x. Therefore,
ˆ
˙
p
1ÿ
i´1
´
pn ´ iq log 1 ´
2 i“1
n

Chételat and Wells/Mid-scale Wishart asymptotics

22

ˆ K`4 ˙
p
p
K`2
K`1
1 ÿ ÿ pi´1qk
1 ÿ ÿ pi´1qk
p
´
`
O
n
i
2 i“1 k“1 k nk
2 i“1 k“1 k nk
nK`2
ˆ K`3 ˙
˙
p
p K`1 ˆ
p
p
1ÿ
1 ÿ ÿ i´1
i pi´1qk
`O
“
pi´1q `
´
k
K`1
2 i“1
2 i“1 k“1 k´1 k
n
n
n
ˆ
˙
˙
ˆ
p
K`1
pK`3
ppp ´ 1q 1 ÿ ÿ
1 pi´1qk
i´1
`o
“
´
`
k
4
2 i“1 k“1 kpk`1q k
n
nK`1
ˆ K`3 ˙
p´1
p´1
K`1
K`1
ÿ
ppp ´ 1q 1 ÿ
1
1 ÿ 1 ÿ k
p
k`1
“
´
i
i
.
´
`
o
k
k
4
2 k“1 kpk`1qn i“1
2 k“1 k n i“1
nK`1
“

Now let Bk denote the Bernoulli numbers, with the convention B1 “ 12 . Faulhaber’s formula provides
K`1
k`2
ppp ´ 1q 1 ÿ
1
1 ÿ
“
´
¨
Bk`2´l pp´1ql
4
2 k“1 kpk`1qnk k`2 l“1
ˆ K`3 ˙
K`1
k`1
1 ÿ
p
1 ÿ 1
l
¨
Bk`1´l pp´1q ` o
.
´
2 k“1 k nk k`1 l“1
nK`1
But by the binomial theorem,
pk`1
n

` op1q and

pp´1ql
“ op1q
nk
K`1
ÿ „ˆ

pp´1qk`2
nk

“

pk`2
nk

k`1

´ pk ` 2q pnk ` op1q,

pp´1qk`1
nk

“

for any 1 ď l ď k. Thus
˙
k`2
B0
pk`2 ´ pk`2qpk`1
ppp ´ 1q 1
´
“
4
2 k“1
k`2 kpk`1qpk`2q
nk
ˆ
˙

k`2
B1
pk`1
`
` op1q
k`1 kpk`1qpk`2q nk
˙

ˆ K`3 ˙
K`1 „ˆ
1 ÿ
k`1
B0 pk`1
p
` op1q ` o
.
´
2 k“1
k`1 kpk`1q nk
nK`1

Using that B0 “ 1 and B1 “ 12 , we obtain
“

K`1
ppp ´ 1q 1 ÿ
1
pk`2
´
4
2 k“1 kpk`1qpk`2q nk

ˆ K`3 ˙
K`1
1 ÿ
1
pk`1
p
´
`o
.
4 k“1 kpk`1q nk
nK`1

(4.6)

The analysis of the other sum of Equation (4.4) is similar but more involved,
as we must distinguish the cases where p is even and where p is odd. We find,
from Equation (4.5) again, that
˙
p
¯ ˆ
1 ÿ´
2i´p´3
n ´ r2i´p´1s log 1 ´
2 i“1
n
ˆ
˙
ˆ
˙
p
p
1ÿ
1ÿ
2i´p´3
2i´p´3
“´
`
p2i´p´1q log 1 ´
n log 1 ´
2 i“1
n
2 i“1
n

Chételat and Wells/Mid-scale Wishart asymptotics

23

ˆ K`4 ˙
p
p K`2
K`1
ÿ p2i´p´3qk
1ÿ
1 ÿ ÿ p2i´p´3qk
p
´
`
O
p2i´p´1q
n
k
k
K`2
2 i“1
k
n
2
k
n
n
i“1 k“1
k“1
˙
p
p K`1ˆ
1ÿ
1 ÿ ÿ 2i´p´1 2i´p´3 p2i´p´3qk
“
p2i´p´3q `
´
2 i“1
2 k“1 k“1
k
k`1
nk
ˆ K`3 ˙
p
p
`O
K`1
n
n
ˆ K`3 ˙
ˆ
˙
p K`1
k
ÿ
ÿ
p
2i´p´3 2 p2i´p´3q
1
`o
`
“p`
2 i“1 k“1 kpk`1q
k
nk
nK`1
“

“p`

p
K`1
ÿ
1 ÿ
1
p2i´p´3qk`1
2 k“1 kpk`1qnk i“1

ˆ K`3 ˙
p
1 ÿ
p
k
`
p2i´p´3q ` o
k
K`1
kn
n
i“1
k“1
„

p
ÿ
k`1
k`1
k`1
p´p´1q
`
p´p`1q
`
p2i´p´3q
k
K`1
ÿ

K`1
1 ÿ
1
2 k“1 kpk`1qn
i“3

ˆ K`3 ˙
p
K`1
ÿ 1 „
ÿ
p
k
k
k
`
p´p´1q ` p´p`1q ` p2i´p´3q ` o
.
k
kn
nK`1
i“3
k“1

“p`

(4.7)

At this point, it is simpler to analyze the cases where p is even and odd separately. If p is odd, define q “ pp´3q{2 and observe that by Faulhaber’s formula,
p
q
”
ıÿ
ÿ
p2i´p´3ql “ 1`p´1ql
p2iql
i“3

i“1

“

1rl evens
l`1

By the binomial theorem,
op1q and

p2qql
“
nk
k`1

p´p´1q
nk

l`1
ÿ

ˆ

s“1

p2qqk`2
nk

“

˙
l`1
Bl`1´s 2l`1´s p2qqs .
s
pk`2
nk

k`1

´ 3pk`2q pnk ` op1q,

p2qqk`1
nk

“

pk`1
nk

`

op1q for 1 ď l ď k. Moreover,

“

p´p`1qk`1
p´pqk`1
“
` op1q
k
n
nk

(4.8)

and
p´p´1qk
p´p`1qk
“
“ op1q.
(4.9)
nk
nk
Thus, for odd p, Equation 4.7 equals
„
ˆ
˙
K`1
1 ÿ
1
k`2 B0 pk`2 ´ 3pk`2qpk`1
“p`
` 1rk odds
2 k“1 kpk`1q
k`2 k`2
nk
ˆ
˙

k`2 2B1 pk`1
p´pqk`1
` 1rk odds
`2
` op1q
nk
k`1 k`2 nk

Chételat and Wells/Mid-scale Wishart asymptotics

`

24


ˆ K`3 ˙
ˆ
˙
K`1 „
1 ÿ 1
p
k`1 B0 pk`1
`
op1q
`
o
1rk evens
.
2 k“1 k
k`1 k`1 nk
nK`1

Moreover,
2p´pqk`1 ´ 31rk odds pk`1 ` 1rk odds pk`1 ` 1rk evens pk`1
“ ´1rk evens pk`1 .
Thus,
ˆ K`3 ˙
K`1
K`1
p
1 ÿ 1rk odds pk`2 1 ÿ 1rk evens pk`1
´
`o K`1 .
“p`
2 k“1 kpk`1qpk`2q nk
2 k“1 kpk`1q nk
n

(4.10)

When p is even, let q “ pp ´ 2q{2 and observe that by Faulhaber’s formula,
ˆÿ
˙
p
q
2q
q
”
ıÿ
ÿ
ÿ
l
l
l
l
l
p2i´p´3q “ 1`p´1q
p2i ´ 1q “ 21rl evens
i ´
p2iq
i“3

i“1

“

1rl evens
l`1

But by the binomial theorem,
p2qql
nk

pk`1
nk

l`1
ÿ
s“1

ˆ

i“1

i“1

˙
l`1
Bl`1´s p2 ´ 2l`1´s q p2qqs .
s

p2qqk`2
nk

“

pk`2
nk

k`1

´ 2pk `2q pnk ` op1q,

p2qqk`1
nk

“

` op1q and
“ op1q for 1 ď l ď k. If we apply Equations (4.8)–(4.9),
then Equation (4.7) becomes
„
ˆ
˙
K`1
1
k`2 p2 ´ 1qB0 pk`2 ´ 2pk ` 2qpk`1
1 ÿ
1rk odds
“p`
2 k“1 kpk`1q
k`2
k`2
nk
ˆ
˙

k`2 p2 ´ 2qB1 pk`1
` 1rk odds
`
op1q
k`1
k`2
nk
„
ˆ
˙

ˆ K`3 ˙
K
k`1 p2 ´ 1qB0 pk`1
p
1 ÿ 1
1rk evens
` op1q `o K`1 .
`
k
2 k“1 k
k`1
k`1
n
n
Moreover,
2p´pqk`1 ´ 21rk odds pk`1 ` 1rk evens pk`1 “ ´1rk evens pk`1 .
Thus again,
“p`

ˆ K`3 ˙
K`1
K`1
1 ÿ 1rk odds pk`2 1 ÿ 1rk evens pk`1
p
´
`o
.
k
k
2 k“1 kpk`1qpk`2q n
2 k“1 kpk`1q n
nK`1

(4.11)

which is the exact same result as in the odd p case (see Equation 4.10). Plugging
Equations (4.6) and (4.10)–(4.11) in Equation (4.4), we obtain
K`1
ppp ` 1q
1 ÿ
pp3p ` 1q
1rk evens pk`2
log 2 ´
log π ´
log Cn,p “
4
2
2 k“1 kpk`1qpk`2q nk
ˆ K`3 ˙
K`1
1 ÿ 1 ` 21rk evens pk`1
p
´
`o K`1 ,
4 k“1
kpk`1q
nk
n
as desired.

Chételat and Wells/Mid-scale Wishart asymptotics

25

Thus the constant Cn,p is closely related to the normalization constant of the
GOEppq distribution 2pp3p`1q{4 {π ppp`1q{2 .
We now turn our attention to the study of the asymptotic moments of a
Tn{2 pIp {8q distribution. We first remind the reader of some classic results. For a
Gaussian Orthogonal Ensemble matrix Z „ GOEppq, a moment-based approach
to Wigner’s theorem states that for any k P N, its k th moment satisfy

„
1 ´ Z ¯k
“ Ck{2 1rk evens ,
lim E tr ?
pÑ8
p
p
`2k˘
1
where Ck “ k`1
is the k th Catalan number. In fact, Anderson et al.
k
[2010, section
2.1.4 on‰ p.17] show that the variance of the k th moment satisfies
“1
?
lim Var p trpZ{ pqk “ 0, so we really have
pÑ8

1 ´ Z ¯k L2
tr ?
ÝÑ Ck{2 1rk evens
p
p
as p Ñ 8.
“ Now,‰ what do we know about the moments of Tn{2 pIp {8q? By symmetry,
E tr T k “ 0 for odd k, but it is much less clear what happens for even k. It
turns out that in many ways, if T „ Tn{2 pIp {8q then 4T „ Tn{2 p2Ip q mimics the
Gaussian Orthogonal Ensemble results outlined above, especially when p{n Ñ 0
as n Ñ 8. We have the following result.
Theorem 2. Let k P N and T „ T“n{2 pIp {8q.
‰ If p{n Ñ c P r0,“1q, the‰moments of
T satisfy the asymptotic bounds E tr T 2k “ Oppk`1 q and E tr2 T k “ Oppk`2 q
as n Ñ 8. In fact, for any k P N,
1 ´ 4T ¯k L2
tr ?
ÝÑ Ck{2 1rk evens
p
p
`2k˘
1
th
Catalan number.
as n, p Ñ 8 with p{n Ñ 0, where Ck “ k`1
k is the k
Although our proof will rely on the close relationship between the Wishart
and the t distribution, it is worthwhile to step back and think why a Tn{2 p2Ip q
should behave like a GOEppq when p{n Ñ 0. One good reason might be the
classic result that as n Ñ 8, the density of a t distribution converges pointwise
to a standard normal density. Thus, we might think that as long as p does not
grow too fast, in some aspects the symmetric t distribution should behave like
a GOEppq.
In the context of the proof, it will prove useful to use the notion of power sum
symmetric polynomials. For any integer partition κ “ pκ1 , . . . , κq q in decreasing
order κ1 ě ¨ ¨ ¨ ě κq ą 0, define its associated power sum polynomial to be
q
ź
rκ pZq “
tr Z κi .
(4.12)
i“1

The norm of the partition κ is |κ| “ κ1 ` ¨ ¨ ¨ ` κq ą 0, which should not be
confused with its length qpκq “ q (number of elements).
By convention, we will assume there also exists an empty partition ∅ “ pq
with length qp∅q “ 0, norm |∅| “ 0 and power sum polynomial r∅ pZq “ 1.

Chételat and Wells/Mid-scale Wishart asymptotics

26

Let’s now turn to the proof of the theorem. The odd moments of the Tn{2 pIp {8q
moments
sense
to focus on the even moments
“
‰ are zero by symmetry, so it makes
“
‰
E tr T 2k and the square moments E tr2 T k . Our first step in the proof is to
express these in terms of expectations of power sum polynomials of an inverse
Wishart Y ´1 „ W´1
p pn, Ip {nq, where by power sum polynomials we mean expressions like at Equation (4.12). Recall the useful shorthand m “ n ´ p ´ 1.
Lemma 2. Let T „ Tn{2 pIp {8q. Then for any k P N, whenever n is large enough
so that n ě p ` 16k ` 6, we can compute the 2k th moment of T by

ż

np

“
‰ p´1qk
E tr T 2k “
nk
p
ÿ

¨
i1 ,...,i2k

! n
)ˇ ˇ m
n2
` n ˘ exp ´ tr Y ˇY ˇ 4
np
4
2 2 Γp 2
Y ą0

ˇ
)ˇ
! n
ˇ m4 ˇ
Bs
Bs
Bs
ˇ
ˇ
...
dY
exp ´ trpY `Xq Y `X ˇˇ
Bs Xi1 i2k
Bs Xi3 i2 Bs Xi2 i1
4
X“0
“

“
‰
p´1qk ÿ p1q
bκ pn, m, pq E rκ pY ´1 q
k
n

(4.13)

|κ|ď2k

and its squared k th moment by

ż
np

“
‰ p´1qk
E tr2 T k “
nk

! n
)ˇ ˇ m
n2
` n ˘ exp ´ tr Y ˇY ˇ 4
4
2 Γp 2
np
2

p
ÿ

Y ą0

Bs
Bs
Bs
Bs
Bs
Bs
...
...
B
X
B
X
B
X
B
X
B
X
B
X
s j3 j2 s j2 j1 s i1 ik
s i3 i2 s i2 i1
i1 ,...,ik s j1 jk
j1 ,...,jk
ˇ
)ˇ
! n
ˇ m4 ˇ
ˇ
ˇ
dY
exp ´ trpY `Xq Y `X ˇˇ
4
X“0
k
ÿ
“
‰
p´1q
“
bκp2q pn, m, pq E rκ pY ´1 q ,
(4.14)
k
n
¨

|κ|ď2k`1

for Y

´1

„

W´1
p pn, Ip {nq

p1q

p2q

p1q

p2q

and some bκ , bκ . These bκ , bκ are polynomials in
p1q

n, m, p, indexed by integer partitions κ, whose degrees satisfy deg bκ ď 2k ` 1 ´
p2q
qpκq and deg bκ ď 2k ` 2 ´ qpκq. The sums are taken over all partitions of the
integers κ satisfying |κ| ď 2k and |κ| ď 2k ` 1 respectively, including the empty
partition.
Proof. Let fNW and ψNW
? stand for the density and the G-transform of a normalized Wishart matrix nrWp pn, Ip {nq ´ Ip s. In the proof of Proposition 3, we
1{2
concluded at Equation (3.11) that fNW had to be integrable when n ě p ´ 2, as
its integral was proportional to a multivariate gamma function. Let RpXq “ ´X
1{2
1{2
be the flip operator. Since fNW is integrable, fNW ˝R must be integrable as well,

Chételat and Wells/Mid-scale Wishart asymptotics

27

` 1{2
˘
1{2
and so their convolution fNW ‹ fNW ˝ R is well-defined and integrable.
At Equation (3.1) at the start of Section 3, we defined our notion of Fourier
transform for integrable functions on Sp pRq. Define the map ι : Sp pRq Ñ
Rppp`1q{2 that maps a symmetric matrix to its vectorized upper triangle, and
let τ : Sp pRq Ñ Sp pRq be the map
#
2Xij if i ‰ j
τ pXqij “
Xjj
if i “ 1.
Then in terms of the usual Fourier transform on Rppp`1q{2 ,
(`
˘
ppp´1q
Ftf upT q “ 2 4 F f ˝ ι´1 ι ˝ τ pT q .
This close relationship transfer properties to our Fourier transform on Sp pRq.
We will need three.
(
1. For any integrable function f , we have F f ˝ R “ Ftf u.
2. (Convolution) For any two integrable functions f1 and f2 , we have Ftf1 ‹
ppp`1q
p
f2 u “ 2 2 π 4 Ftf1 uFtf2 u.
3. (Fourier inversion) For any continuous integrable f with integrable Fourier
transform φ, we have

ż
f pXq “

1
p
2

2 π

ppp`1q
4

ei trpT Xq φpT q dT “ Ftφup´Xq,
Sp pRq

for all X P Sp pRq.
These properties are important for the following. Since fNW is real-valued, properties 1 and 2 provide
!
)
` 1{2
˘)
(
ppp`1q
p
1{2
1{2
1{2
2´ 2 π ´ 4 F fNW ‹ fNW ˝ R “ F fNW F fNW ˝ R
ˇ 1{2 ˇ
1{2 (
1{2 1{2
1{2 (
“ F fNW F fNW “ ψNW ψNW “ ˇψNW ˇ.
But then, since |ψNW | is integrable (in fact, to unity), the Fourier inversion
formula yields that

ż

` 1{2
˘
1{2
fNW ‹ fNW ˝ R pXq “

ˇ
ˇ
ei trpT Xq ˇψNW pT qˇ dT.

(4.15)

Sp pRq

Thus we might say the characteristic function of the Tn{2 pIp {8q distribution is
˘
1{2 ` 1{2
given by fNW ‹ fNW ˝ R . It is well known that the derivatives of the characteristic function of a distribution evaluated at zero provide its moments, up to a
˘
1{2 ` 1{2
constant. This suggests
should
try“ to repeatedly
differentiate fNW ‹ fNW ˝R
“ we2k
‰
‰
at zero to compute E tr T
and E tr2 T k , our ultimate goal.
Unfortunately, the convolution is given by an integral whose domain makes
it difficult to directly interchange the differentiation and integration symbols.
Because the integrand is orthogonally invariant, we found it easier to compute
the derivatives at zero by taking a limit over a sequence of decreasing positive-

Chételat and Wells/Mid-scale Wishart asymptotics

28

definite matrices at both sides instead. In this spirit, define on the open set
t0 ă X ă Ip u Ă Sp pRq the real-valued functions
p´1qk
H1 pXq “
nk

p
ÿ

˘?
Bs
Bs
Bs
1{2 ` 1{2
...
fNW ‹ fNW ˝R p nXq
Bs Xi1 i2k
Bs Xi3 i2 Bs Xi2 i1

i1 ,...,i2k

and
p´1qk
H2 pXq “
nk

p
ÿ

Bs
Bs
Bs
Bs
Bs
Bs
...
...
Bs Xj1 jk
Bs Xj3 j2 Bs Xj2 j1 Bs Xi1 ik
Bs Xi3 i2 Bs Xi2 i1

i1 ,...,ik
j1 ,...,jk

for fixed k, p and n. Here

˘?
1{2 ` 1{2
fNW ‹ fNW ˝R p nXq
Bs
Bs Xij

stands for the symmetric differentiation operator
?
“
as defined in Section 2. The n scaling in the argument
helps link the convolution to an expectation with respect to an inverse Wishart
distribution.
Let’s first relate these functions to the moments of the Tn{2 pIp {8q distribution. The symmetric differentation operator has the pleasant property that
Bs
Bs Xij trpXT q “ Tij for any two symmetric matrices X, T . Thus, for any 1 ď
l ď 2k and indices 1 ď i1 , . . . , i2l ď p, we find that
ˇ
ˇ
?
ˇ
ˇ
ˇ
ˇ
Bs
Bs
Bs
i n trpT Xq ˇ
ˇ
ˇ
...
ψNW pT qˇˇ
e
ˇ Bs Xi i
Bs Xi4 i3 Bs Xi2 i1
2l 2l´1
ˇ
ˇˇ
ˇ
ˇ
ˇ
“ nl ˇTi2l i2l´1 ¨ ¨ ¨ Ti4 i3 Ti2 i1 ˇˇψNW ˇpT q (4.16)
1`δij B
2 BXij ,

Bs
Bs Xij

for all X P Sp pRq.
We now show that the right hand side (4.15) is integrable. This is not a mere
formality: when p “ 1, asking if this expression is integrable is the same as
asking if the t distribution with n{2 degrees of freedom has an lth moment, and
it is well-known that the t distribution only possesses moments of order smaller
than its degrees of freedom. So the answer is most likely to be positive, but only
for n large enough.
Let us see why. For any symmetric matrix T ,
g
f p ´
¯ bˇ
a
fź
ˇ
2
|Tij | ď λ1 pT q ď e
1 ` λi pT 2 q “ ˇIp ` T 2 ˇ,
i“1

where λ1 pT 2 q ě ¨ ¨ ¨ ě λp pT 2 q ě 0 are the ordered eigenvalues of the positivedefinite matrix T 2 . Thus

ż
ˇ
ˇˇ
ˇ
ˇ
ˇ
nl ˇTi2l i2l´1 ¨ ¨ ¨ Ti2 i1 ˇˇψNW ˇpT q dT
Sp pRq

ż
3l
2

n
“ l Cn,p
4

ˇ´ n`p`1
ˇ
ˇ 4T
4
4Ti2 i1 ˇˇˇˇ
16T 2 ˇˇ
ˇ i2l i2l´1
¨ ¨ ¨ ? ˇˇIp `
dT
ˇ ?
ˇ
n
n
n
Sp pRq

Chételat and Wells/Mid-scale Wishart asymptotics

ż

ˇ
ˇ pn´2lq`p`1
2 ˇ´
4
ˇ
ˇIp ` 16T ˇ
dT
ˇ
n ˇ

3l
2

n
ď l Cn,p
4

Sp pRq
3l
2

ˆ

n Cn,p
ď l
4 Cn´2l,p

29

ż

n
n ´ 2l

˙ppp`1q
4

ˇ
ˇ´ pn´2lq`p`1
4
ˇ
16T 2 ˇˇ
ˇ
Cn´2l,p ˇIp `
dT.
n ´ 2l ˇ

(4.17)

Sp pRq

When n ´ 2l ě p ´ 2, the last integrand is the density of a Tn{2´l pIp {8q
distribution, so integrates to unity. Thus, when n ě p ` 4k ´ 2, the right
hand side of Equation (4.16) is an integrable function for all 1 ď l ď 2k and
1 ď i1 , . . . , i2l ď p. By Equation (4.15), and repeated differentiation under the
integral sign justified by the integrability bounds given by Equations (4.16) and
(4.17), we find that
k

H1 pXq “

p´1q
nk

ż

p
ÿ

?

Bs
Bs
Bs
...
Bs Xi1 i2k
Bs Xi3 i2 Bs Xi2 i1

ei

i1 ,...,i2k

ˇ

n trpT Xq ˇ

ˇ
ψNW pT qˇdT

Sp pRq

ż
tr T 2k ei

“

?

ˇ

n trpT Xq ˇ

ˇ
ψNW pT qˇdT

(4.18)

Sp pRq

and
p´1qk
H2 pXq “
nk

p
ÿ

Bs
Bs
...
Bs Xj1 jk
Bs Xj2 j1

i1 ,...,ik
j1 ,...,jk

ż
Bs
Bs
Bs
...
Bs Xi1 ik
Bs Xi3 i2 Bs Xi2 i1

ei

ˇ
ψNW pT qˇdT

Sp pRq

ż
“

?
ˇ
n trpT Xq ˇ

tr2 T k ei

?

ˇ

n trpT Xq ˇ

ˇ
ψNW pT qˇdT.

(4.19)

Sp pRq

for any X P Sp pRq and any n ě p ` 4k ´ 2.
˘
1{2 ` 1{2
Now let’s relate H1 and H2 to the definition of fNW‹ fNW˝R as a convolution.
This is where restricting H1 and H2 to small positive-definite matrices becomes
useful. By Equation (3.10), the expression is

ż
˘?
1{2 ` 1{2
fNW ‹ fNW ˝R p nXq “

ż
np

“

n2
` ˘
np
2 2 Γp n2

1{2

1{2

fNW pZqfNW pZ ´

?

nXqdZ

Sp pRq

! n
)ˇ
! n
)ˇ ˇ m
ˇm
exp ´ trpY `Xq ˇY `X ˇ 4 exp ´ tr Y ˇY ˇ 4 dY
4
4
Y `Xą0, Y ą0

?
ppp`1q
using the change of variables Y “ Ip ` Z{ n ´ X with dZ “ n 4 dY . For

Chételat and Wells/Mid-scale Wishart asymptotics

30

X ą 0, we have 1rY ` X ą 0, Y ą 0s “ 1rY ą 0s, and thus H1 , H2 satisfy
H1 pXq “

p
ÿ
p´1qk

nk

ż

Bs
Bs
Bs
...
Bs Xi1 i2k
Bs Xi3 i2 Bs Xi2 i1

i1 ,...,i2k

and
p´1qk
H2 pXq “
nk

p
ÿ

Y ą0

np

ˇm
ˇ
¨ˇY `X ˇ 4

! n
)
exp ´ trpY `Xq
4

! n
)ˇ ˇ m
n2
` n ˘ exp ´ tr Y ˇY ˇ 4 dY
np
4
2 2 Γp 2

(4.20)

Bs
Bs
Bs
...
Bs Xj1 jk
Bs Xj3 j2 Bs Xj2 j1

i1 ,...,ik
j1 ,...,jk

ż

Bs
Bs
Bs
...
Bs Xi1 ik
Bs Xi3 i2 Bs Xi2 i1

! n
)ˇ
ˇm
exp ´ trpY `Xq ˇY `X ˇ 4
4
Y ą0

np
2

! n
)ˇ ˇ m
n
` n ˘ exp ´ tr Y ˇY ˇ 4 dY.
(4.21)
4
2 Γp 2
We would now like to interchange the integral and differentiation signs. )To do
m
so, we must understand what the repeated derivatives of expt´ n4 tr Y |Y | 4
look like. Differentiating once, we see that:
! n
)ˇ
ˇm
Bs
exp ´ trpY `Xq ˇY `X ˇ 4
Bs Xi2 i1
4

„
! n
)ˇ
ˇm
n
m
´1
pY `Xqi2 i1 ´ pIp qi2 i1 exp ´ trpY `Xq ˇY `X ˇ 4 .
“
4
4
4
Differentiating twice, we see that:
! n
)ˇ
ˇm
Bs
Bs
exp ´ trpY `Xq ˇY `X ˇ 4
Bs Xi4 i3 Bs Xi2 i1
4
„
m
m
´1
´1
“ ´ pY `Xq´1
pY `Xq´1
i2 i4 pY `Xqi3 i1 ´
i2 i3 pY `Xqi4 i1
8
8
mn
m2
´1
´1
`
pY `Xq´1
pY `Xq´1
i4 i3 pY `Xqi2 i1 ´
i4 i3 pIp qi2 i1
16
16

! n
)ˇ
ˇm
n2
mn
´1
´1
´1
´
pIp q´1
pY
`Xq
`
pI
q
pI
q
trpY `Xq ˇY `X ˇ 4.
p i4 i3 p i2 i1 exp ´
i4 i3
i2 i1
16
16
4
So in general, it is clear that the repeated derivatives are) given by some polym
nomial in entries of pY ` Xq´1 , times expt´ n4 trpY ` Xq |Y ` X| 4 . We won’t
investigate further the nature of these polynomials beyond remarking that for
any indices 1 ď l ď 2k and 1 ď i1 , . . . , i2l ď p, and any symmetric matrices
X, Y P Sp pRq, we must have some crude bound
ˇ
ˇ
*
"
ˇ
ˇ
ˇ m4 ˇ
Bs
Bs
Bs
n
ˇ
ˇ
ˇ
...
exp ´ trpY `Xq Y `X ˇˇ
ˇ Bs Xi i
Bs Xi i Bs Xi i
4
¨

2l 2l´1

4 3

2 1

np
2

Chételat and Wells/Mid-scale Wishart asymptotics
l
ÿ

ÿ

ď

31

s
l
ˇ ź
ˇ
ˇ
ˇ
ˇ
ˇź
ˇ
ˇ
ˇpY `Xq´1
ˇ
ˇ
ˇ
pI
q
pn,
mq
ˇ
ˇaJ,s
p j2t j2t´1
j2t j2t´1

"t“1
*
ˇ
ˇm
n
¨ exp ´ trpY `Xq ˇY `X ˇ 4
4
for some polynomials aJ,s that do not depend on X or Y . We relegate a proof
of this result as Lemma 3 in Section 8. This can be uniformly bounded for all
0 ď X ď Ip by
l
! n
)“
ÿ
‰ mp
ď C1 pn, m, pq
trs pY ´1 q exp ´ tr Y 1 ` tr Y 4
(4.22)
4
s“0
s“0

t“s`1

JP
t1,...,pu2l

for some constant C1 pn, m, pq that does not depend on X or Y . But for any
n ě p ´ 2 and l ě 0,

ż
l
ÿ

)“
! n
‰ mp
trs pY ´1 q exp ´ tr Y 1 ` tr Y 4
4
s“0
np
! n
)ˇ ˇ m
n2
¨ np ` n ˘ exp ´ tr Y ˇY ˇ 4 dY
4
2 2 Γp 2

C1 pn, m, pq
Y ą0

ż

“

n

np
2

l
C1 pn, m, pq ÿ
`n˘
np
2 2 Γp 2
s“0

“

1`tr Y

‰ mp
4

trs pY ´1 q

Y ą0

! n
)ˇ ˇ n`p`1 p`1
´ 2
exp ´ tr Y ˇY ˇ 4
dY
2
`
˘ ”
ı
´ n ¯ mp
Γp n`p`1
mp
4
4˘
`
E p1 ` tr Y q 4 trs pY ´1 q
“
n
2
Γp 2
`
˘
for a Y with a matrix gamma distribution Gp n`p`1
, n2 Ip . The Cauchy-Schwarz
4
inequality then entails the bound
˘ ”
`
ı 12
´ n ¯ mp
ı 21 ”
Γp n`p`1
mp
4
4
` n ˘ E p1 ` tr Y q 2
E tr2s pY ´1 q .
(4.23)
ď
2
Γp 2
The first expectation is always finite when n ě p ´ 2. Since tr2s pY ´1 q can be
written as a sum of zonal polynomials indexed by partitions of the integer 2s, the
results of Muirhead [1982, Theorem 7.2.13] imply that the second expectation
is finite whenever n`p`1
ą 2s ` p´1
4
2 ô n ě p ` 8s ´ 2. Thus, in Equation (4.20)
with l ď k and (4.21) with l ď 2k, whenever n ě p ` 16k ´ 2 we are justified
in repeatedly differentiating under the integral sign by the integrability bounds
given by Equations (4.22) and (4.23), and obtain in that case

ż

p´1qk
H1 pXq “
nk

p
ÿ

Y ą0 i1 ,...,i2k

Bs
Bs
Bs
...
Bs Xi1 i2k
Bs Xi3 i2 Bs Xi2 i1

(4.24)

Chételat and Wells/Mid-scale Wishart asymptotics

32

np
! n
)ˇ
! n
)ˇ ˇ m
ˇ m4
n2
ˇ
ˇ
` n ˘ exp ´ tr Y ˇY ˇ 4 dY
exp ´ trpY `Xq Y `X
np
4
4
2 2 Γp 2

(4.25)
and

ż

p
ÿ

p´1qk
H2 pXq “
nk

i ,...,ik
Y ą0 1

Bs
Bs
Bs
...
Bs Xj1 jk
Bs Xj3 j2 Bs Xj2 j1

j1 ,...,jk

! n
)ˇ
ˇm
Bs
Bs
Bs
...
exp ´ trpY `Xq ˇY `X ˇ 4
Bs Xi1 ik
Bs Xi3 i2 Bs Xi2 i1
4
np
! n
)ˇ ˇ m
n2
¨ np ` n ˘ exp ´ tr Y ˇY ˇ 4 dY.
4
2 2 Γp 2
(4.26)
Let us now look at how H1 pXq and H2 pXq behave
as X Ñ 0. On one hand,
a
?
for any symmetric matrix T we have | tr T k | ď p tr T 2k ď p|Ip ` T 2 |k{2 , so
we must have the bounds
ˇ
ˇ
ˇ
ˇ pn´4kq`p`1
k
2 ˇ´
4
?
ˇ
ˇ
ˇˇ
ˇ
n
C
16T
n,p
2k
i
n
trpT
Xq
ˇ tr T e
ˇψNW pT qˇˇ ď
ˇ
ˇ
ˇ 16k Cn´4k,p Cn´4k,p ˇIp ` n ˇ
ˇ
and
ˇ
ˇ
ˇ pn´4kq`p`1
ˇ
k
2 ˇ´
4
ˇˇ
ˇ 2 k i?n trpT Xq ˇ
ˇ
pn
C
16T
n,p
ˇψNW pT qˇˇ ď
ˇ tr T e
ˇ
ˇ
.
ˇ 16k Cn´4k,p Cn´4k,p ˇIp ` n ˇ
ˇ
holding uniformly in X. When n ´ 4k ě p ´ 2 ô n ě p ` 4k ´ 2, the right hand
sides are proportional to the density of the G-conjugates of the normalized
Wishart distributions with n ´ 4k degrees of freedom, so are integrable. Thus,
by the dominated convergence theorem and Equations (4.18) and (4.19),

ż
lim

XÑ0 H1 pXq
0ăXăIp

“

ˇ
ˇ
“
‰
lim ?
tr T 2k XÑ0 ei n trpT Xq ˇψNW pT qˇdT “ E tr T 2k

(4.27)

0ăXăIp

Sp pRq

and
lim

XÑ0 H2 pXq
0ăXăIp

ż
“

lim ei?n trpT Xq ˇˇψ pT qˇˇdT “ E“tr2 T k ‰
tr2 T k XÑ0
NW

(4.28)

0ăXăIp

Sp pRq

for T „ Tn{2 pI{ 8q.
On the other hand, the integrands at Equations (4.25) and (4.26) take a
particularly simple form. Lemma 4 establishes by induction that there must be
p1q
p2q
p1q
polynomials bκ and bκ in n, m and p with degrees deg bκ ď 2k ` 1 ´ qpκq

Chételat and Wells/Mid-scale Wishart asymptotics

33

p2q

and deg bκ ď 2k ` 2 ´ qpκq such that

ż
p´1qk
H1 pXq “
nk

ÿ

bp1q
κ pn, m, pqrκ prY

` Xs

|κ|ď2k
Y ą0

"
*
n
q exp ´ trpY ` Xq
4

np

ˇ
ˇm
¨ˇY ` X ˇ 4
and

´1

! n
)ˇ ˇ m
n2
` n ˘ exp ´ tr Y ˇY ˇ 4 dY
4
2 Γp 2
np
2

(4.29)

ż

p´1qk
H2 pXq “
nk

ÿ

bp2q
κ pn, m, pqrκ prY
|κ|ď2k`1

´1

` Xs

"
*
n
q exp ´ trpY ` Xq
4

Y ą0

np
! n
)ˇ ˇ m
n2
` n ˘ exp ´ tr Y ˇY ˇ 4 dY
(4.30)
4
2 Γp 2
for any 0 ă X ă Ip and n ě p ` 16k ´ 2. The sums are taken over all partitions
of the integers κ satisfying |κ| ď 2k and |κ| ď 2k ` 1 respectively, including the
empty partition. But for any integer partition κ, the bound
ˇm
ˇ
“
‰ mp
n
n
rκ prY ` Xs´1 qe´ 4 trpY `Xq ˇY ` X ˇ 4 ď tr|κ| pY ´1 qe´ 4 tr Y 1 ` tr Y 4

ˇ
ˇm
¨ˇY ` X ˇ 4

np
2

holds uniformly in 0 ď X ď Ip . Thus for |κ| ď 2k ` 1, the right hand side is
integrable for n ě p ` 16k ` 6, by the same argument as for Equation (4.23).
Thus for such n, by the dominated convergence theorem and Equations (4.29)
and (4.30), we obtain that

ż

lim

XÑ0 H1 pXq
0ăXăIp

“

p´1qk
nk

lim

XÑ0
0ăXăIp
Y ą0

ÿ

´1
bp1q
q
κ pn, m, pqrκ prY ` Xs

|κ|ď2k

"
*
np
! n
)ˇ ˇ m
ˇ
ˇm
n
n2
¨ exp ´ trpY ` Xq ˇY ` X ˇ 4 np ` n ˘ exp ´ tr Y ˇY ˇ 4 dY
4
4
2 2 Γp 2
k ÿ
“
‰
p´1q
´1
“
bp1q
q
(4.31)
κ pn, m, pq E rκ pY
nk
|κ|ď2k

and
lim

ż

XÑ0 H2 pXq
0ăXăIp

“

p´1qk
nk

lim

XÑ0
0ăXăIp
Y ą0

ÿ

´1
bp2q
q
κ pn, m, pqrκ prY ` Xs

|κ|ď2k`1

"
*
np
)ˇ ˇ m
! n
ˇ
ˇm
n2
n
¨ exp ´ trpY ` Xq ˇY ` X ˇ 4 np ` n ˘ exp ´ tr Y ˇY ˇ 4 dY
4
4
2 2 Γp 2
k
ÿ
“
‰
p´1q
´1
“
bp2q
q ,
(4.32)
κ pn, m, pq E rκ pY
k
n
|κ|ď2k`1

Chételat and Wells/Mid-scale Wishart asymptotics

34

where Y follows a Wp pn, Ip {nq distribution. Combining Equations (4.27)–(4.28)
with Equations (4.31)–(4.32) and Lemma 4 concludes the proof.
Something remarkable about Lemma 2 is that it provides us with an algorithm to compute the moments of a symmetric t distribution in terms of the
moments of an inverse Wishart matrix. For example, when k “ 1, repeated
differentiation yields that
p
! n
)ˇ
ÿ
ˇm
Bs
Bs
exp ´ trpY `Xq ˇY `X ˇ 4
BX
BX
4
i1 ,i2 s i1 i2 s i2 i1
„
mpm ´ 2q
mn
n2 p
“
trpY ` Xq´2 ´
trpY ` Xq´1 `
16
8
16

! n
)ˇ
ˇm
m 2
´ tr pY ` Xq´1 exp ´ trpY `Xq ˇY `X ˇ 4 . (4.33)
8
4
We can recognize trpY ` Xq´2 and tr2 pY ` Xq as the power sum polynomials
rp2q prY ` Xs´1 q and rp1,1q prY ` Xs´1 q in the sense of Equation (4.12), so we
p1q

p1q

p1q

p1q

must have bp2q “ mpm ´ 2q{16, bp1,1q “ ´m{8, bp1q “ mn{8 and b∅ “ n2 p{16 in
the
“ result
‰ of Lemma 4. Hence Lemma 2 really tells us that whenever n ě p ` 22,
E tr T 2 for T „ Tn{2 pIp {8q can be expressed as
“
‰ m “ 2 ´1 ‰
‰
mpm ´ 2q “
E tr Y ´2 `
E tr Y
E tr T 2 “ ´
16n
8n
‰ np
m “
E tr Y ´1 ´
(4.34)
`
8
16
where Y „ Wp pn, Ip {nq.
Of course, this also works with square moments and higher k. For example,
the same strategy
for,
“
‰ say, square moments with k “ 2 yields that whenever
n ě p ` 38, E tr2 T 2 for T „ Tn{2 pIp {8q can be expressed as
`
˘
“ 2 2‰
‰ m pm ´ 2q “
‰
m m2 ´ 5m ` 10 “
E tr T “ ´
E tr Y ´4 `
E tr Y ´3 tr Y ´1
2
2
16n
4n
`
˘
m m3 ´ 4m2 ` 20m ´ 32 “ 2 ´2 ‰
m2 “ 4 ´1 ‰
E tr Y
E tr Y
`
`
2
256n˘
64n2
` 2
‰ m pm ´ 2q “
‰
m m ´ 2m ` 16 “
E tr Y ´2 tr2 Y ´1 `
´
E tr Y ´3
2
64n
8n
`
˘
‰
m m2 ´ 2m ` 16 “
m2 “ 3 ´1 ‰
´
E tr Y ´2 tr Y ´1 `
E tr Y
64n
32n
‰ m p´m ` pq “ 2 ´1 ‰
m pmp ´ 2p ´ 8q “
`
E tr Y ´2 ´
E tr Y
128
64
‰ n2 p 2
mnp “
´
E tr Y ´1 `
(4.35)
64
256
again where Y „ Wp pn, Ip {nq.
Unfortunately, as we consider larger orders, the repeated differentiation of
expt´ n4 tr Zu|Z|m{4 quickly becomes too cumbersome to perform by hand. But

Chételat and Wells/Mid-scale Wishart asymptotics

35

at“ least ‰in theory,
“ we can
‰ compute expressions like (4.34) and (4.35) for any
E tr T 2k and E tr2 T k , and Lemma 2 summarize that fact. That is, using the
Fourier inversion theorem we have reduced the problem of computing moments
of the t-distribution Tn{2 pIp {8q to that of computing expected power sum polynomials of the inverse Wishart distribution W´1
p pn, Ip {nq, for large enough n.
How can we compute expected power sum polynomials of an inverse Wishart?
There are two approaches in the literature. Letac and Massam [2004] found an
expression in terms of a different basis, the zonal polynomials, which behave
particularly nicely with respect to the inverse Wishart distribution, and whose
expectations have a simple closed form. From this, they provided an algorithm
for computing expected power sum polynomials to arbitrary order. Matsumoto
[2012] found expressions of coordinate-wise moments in terms of modified Wiengarten orthogonal functions, from which expectations of power sum polynomials
can be computed. We follow the idea of Letac and Massam [2004] in our asymptotic analysis.
For any integer partition κ, there exist coefficients cκ,λ (which depend solely
on κ and λ) such that
ÿ
rκ pY ´1 q “
cλ Cλ pY ´1 q
|λ|“|κ|

for Cλ the so-called zonal polynomials. For an overview of the topic with a focus
on random matrix theory, see Muirhead [1982, Chapter 7]. The coefficients cκ,λ
are explicitly computable. If we follow the normalization of zonal polynomials
of Muirhead [1982], for example, we find that
“ ‰ “ ‰

„
 „
„
r∅ “ C∅ ,
Cp2q
rp2q
1 ´ 21
and
.
“
(4.36)
“
‰ “
‰
Cp1,1q
rp1,1q
1 1
rp1q “ Cp1q ,
As mentioned, expectations of zonal polynomials with respect to a Wishart or
inverse Wishart distribution take a particularly simple form. From Muirhead
[1982, Theorem 7.2.13 and Equation (18) on p.237], the expected zonal polynomials for Y ´1 „ W´1
p pn, Ip {nq are
“
‰
n|λ|
E Cλ pY ´1 q “
Cλ pIp q
ś
qpλq
2|λ| i“1 m´i`1
2
śqpλq
i ´1
ź λź
2|λ| |λ|! iăj p2λi ´ 2λj ´ i ` jq |λ| qpλq
p ` p1 ´ i ` 2lq
n
“
(4.37)
śqpλq
m
´ p1 ´ i ` 2lq
i“1 l“0
i“1 p2λi ` qpλq ´ iq!
“
‰
for λ ‰ ∅, and E C∅ pY ´1 q “ 1. For example, the first few expected zonal
polynomials are
“
‰
“
‰ np
,
E C∅ pY ´1 q “ 1,
E Cp1q pY ´1 q “
m
“
‰ 2n2 p pp ´ 1q
“
‰
n2 p pp ` 2q
E Cp1,1q pY ´1 q “
,
E Cp2q pY ´1 q “
.
3m pm ` 1q
3m pm ´ 2q
‰
“
‰
“
‰
“
From this, we can exactly compute E rκ pY ´1 q and thus E tr T 2k and E tr2 T k ,

Chételat and Wells/Mid-scale Wishart asymptotics

36

as a function of p and n (or m). For example, by Equation (4.36) we find that
“
‰
“
‰ np
E r∅ pY ´1 q “ 1,
E rp1q pY ´1 q “
,
m
“
‰
“
‰
n2 p pm ` pq
n2 p pmp ´ p ` 2q
, E rp2q pY ´1 q “
.
E rp1,1q pY ´1 q “
m pm ´ 2q pm ` 1q
m pm ´ 2q pm ` 1q
and thus, by Equation (4.34), whenever n ě p ` 38
‰
“
np pmp ` m ` 2q
E tr T 2 “
.
16 pm ´ 2q pm ` 1q

(4.38)

In a similar way, we can compute the expected zonal polynomials and hence,
the expected power sum polynomials of Y ´1 for |κ| “ 3, 4. So from Equation
(4.35), we obtain for n ě p ` 38 that
´
“
‰
n2 p
m3 p3 ` 2m3 p2
E tr2 T 2 “
256 pm ´ 6q pm ´ 2q pm ´ 1q pm ` 1q pm ` 3q
` 5m3 p ` 4m3 ´ 3m2 p3 ` 6m2 p2 ` 9m2 p ` 24m2 ´ 12mp3
¯
´ 36mp2 ` 36mp ´ 36p .
(4.39)
Of course, this reasoning also works for other k’s. In particular, we essentially derived a (potentially inefficient) algorithm to compute the moments of a
Tn{2 pIp {8q distribution to arbitrary order on our path to proving this theorem.
At this point, it is worthwhile to realize that Equations (4.38) and (4.39) are
already enough to prove the theorem for small moments. For example, when
n Ñ 8 such that p{n Ñ c P r0, 1q, then m „ p1 ´ cqn and
ˆ 2
˙
“
‰
m2
np
np
n
1
`
`
„
p2 ,
E tr T 2 “
pm ´ 2q pm ` 1q 16m 16m 8m2
16p1 ´ cq
“
‰
which proves that E tr T 2 “ Opp2 q. In fact,

ˆ
˙
„
m2
1
p
2
3
2
2
1 ´ 4T ¯2
´1“
`
`
`
` 2` 2
E tr ?
pm´2q pm`1q p m m mp m
p
p
m p
1
p
Ñ 0.
(4.40)
„ `
p m
Moreover, when n, p Ñ 8 such that p{n Ñ 0, then m „ n and
„ˆ
˙ 
2
‰ 32 “
‰2
1 ´ 4T ¯2
256 “
E
tr ?
´1
“ 4 E tr2 T 2 ´ 2 E tr T 2 ` 1
p
p
p
p
ˆ
5
m
5
4
2
22
“
` 3`
`
pm´6q pm´2q pm´1q pm`1q pm`3q p2
p
m mp
27
32
p2
6p
33
60
115
52
`
`
`
` 2` 2` 2 ` 2 2` 2 3
2
3
2
mp
mp
m
m
m
m p m p
m p
24
12p2
3p2
6p
3
118
93
96p
´ 3 ` 3` 3` 3 ` 3 2` 3 3´
´ 4
m
m
m
m p m p
m p
m4
m

Chételat and Wells/Mid-scale Wishart asymptotics

48
84
144
144
36
36
´ 4´ 4 ´ 4 2´ 5 ´ 5 ´ 5 2
m
m p m p
m
m p m p
2
p2
5
` 2 Ñ 0.
„ 2`
p
m m

37

˙

(4.41)

Thus
1 ´ 4T ¯2 L2
tr ?
ÝÑ 1 “ C1
p
p
and the theorem is proven for the second moment.
In theory, we could proceed in the same way for any moment of interest, but
naturally we could never conclude that the theorem holds for all moments that
way. Nonetheless, the calculations give us some hints about how to argue in the
general case.
The idea is to express the moments of the symmetric t distribution as polynomials of p and p{m. There are two regimes where random matrix theory is
well understood: the classical regime where p is held fixed as n Ñ 8, and the
linear, high-dimensional regime where p grows linearly with n. From this, we
can therefore conclude a few facts regarding the behavior of symmetric t moments in these regimes. But these moments are polynomials, and a polynomial
is a very rigid object: results from the two extreme cases where p is fixed and
p grows linearly will be enough to prove results for every regime in between,
yielding the first part of the theorem. Proving the second part will then be the
simple matter of applying the GOE approximation of Jiang and Li [2015] and
Bubeck et al. [2016] to the specific shape found for the symmetric t moments
while proving the first part, namely Equations (4.54) and (4.55).
Proof of Theorem 2. Recall the expected zonal polynomial of an inverse Wishart
W´1
p pn, Ip {nq is given by Equation (4.37). Based on the previous calculations,
it is tempting to define
śqpλq
qpλq
i ´1
ź λź
2|λ| |λ|! iăj p2λi ´ 2λj ´ i ` jq
m
1
cλ “
, Rλ pmq “
śqpλq
m
´
p1
´ i ` 2lq
i“1 l“0
i“1 p2λi ` qpλq ´ iq!
ˆ
˙
qpλq
i ´1
ź λź
p
1 ´ i ` 2l
and
Pλ pm, pq “
`
(4.42)
m
m
i“1 l“0
so that
“
‰
E Cλ pY ´1 q “ c1λ n|λ| Rλ pmqPλ pm, pq.
With these expressions the expected power sum polynomials can be written as
˜
¸
p
ÿ
ź
ź
“
‰
´1
1 |κ|
´1
E rκ pY q “
cκ,λ cλ n
Rµ pmq
Rµ pmq Pλ pm, pq
|λ|“|κ|

“

|µ|“|κ|

|µ|“|κ|
µ­“λ

ź
ÿ
ź
n|κ|
¨
Rµ pmq ¨ m|κ|
cκ,λ c1λ
Rµ´1 pmqPλ pm, pq.
|κ|
m
|µ|“|κ|
|λ|“|κ|
|µ|“|κ|
µ­“λ

Chételat and Wells/Mid-scale Wishart asymptotics

38

In other words, if we define
ź
1
R|µ|
“
Rµ pmq,
|µ|“|κ|
1
P|λ|
pm, pq

“ m|κ|

ÿ

cκ,λ c1λ

|λ|“|κ|

ź

Rµ´1 pmqPλ pm, pq

(4.43)

|µ|“|κ|
µ­“λ

then
“
‰
n|κ| 1
1
E rκ pY ´1 q “ |κ| R|κ|
pmqP|κ|
pm, pq.
(4.44)
m
`
˘
ś
ś
λi ´1
qpλq
1´ 1´i`2l
But Rµ´1 pmq “ i“1 l“0
is a polynomial in 1{m, while Pλ pm, pq “
m
˘
śqpλq śλi ´1 ` p
1´i`2l
is a polynomial in p{m and 1{m, both of degree at
i“1
i“1
m `
m
most |µ| “ |λ| “ |κ|. Thus
ÿ
ź
1
cκ,λ c1λ
Rµ´1 pmqPλ pm, pq
P|κ|
pm, pq ” m|κ|
|λ|“|κ|

“ m|κ|

|κ| ÿ
|κ|
ÿ
i“0 j“0

|µ|“|κ|
µ­“λ

bij

´ p ¯i 1
m mj

(4.45)

for some coefficients bij that don’t depend on m, p (or n). Define the polynomials
ř|κ|
fj pαq “ i“0 bij αi , so that
|κ| ´
ÿ
p ¯ ´j
1
P|κ|
pm, pq “ m|κ|
f
m .
(4.46)
m
j“0
Let us show that for all 0 ď`j ă |κ|´qpκq, the polynomial
fj must be identically
˘
zero over the interval α P 0, 1{ maxp|κ| ´ 2, 0q . Indeed, say this was not the
case, and let 0 ď j0 ă |κ|
˘ smallest j with the property that
` ´ qpκq1 be the
. As fj0 is a polynomial, by continuity
fj0 pα0 q ‰ 0 for some α0 P 0, maxp|κ|´2,0q
it must be non-zero in a neighborhood of α0 , so we may as well assume
α‰0 is
“
rational without loss of generality. Now look at what happens to E rκ pY ´1 q as
α0
p grows to infinity at the very specific linear rate p “ t 1`α
pn ´ 1qu. Since α0 is
0
rational, there must be a subsequence nl such that pl is exactly an integer (for
example, if α0 “ a{b with a, b integers, we can take nl “ pa ` bql ` 1). Then for
α0
pl “ 1`α
pnl ´ 1q, we have exactly pl “ α0 ml .
0
` α ˘´1
1
0
Since α0 ă maxp|κ|´2,0q
, then |κ| ă 1 ` 1`α
. Thus by Hölder’s inequality
0
and Lemma 5,
ı
“
‰
1
1
1 ”
lim |κ|´qpκq´j ¨ qpκq E rκ pY ´1 q ď 0 ¨ lim
E tr Y ´|κ| “ 0.
0
lÑ8 m
lÑ8 pl
p
l

l

On the other hand, by Equations (4.44) and (4.46), the definition of j0 and the
fact that R|κ| pmq Ñ 1 as m Ñ 8,
“
‰
1
1
lim
¨ qpκq E rκ pY ´1 q
lÑ8 m|κ|´qpκq´j0
pl
l

Chételat and Wells/Mid-scale Wishart asymptotics

ˆ
“ lim

lÑ8

“ p1 `

nl
ml

˙|κ|

ˆ
1
R|κ|
pml q

ml
pl

˙qpκq ÿ
|κ|

39

fj pα0 qmjl 0 ´j

j“j0

´qpκq
α0 q|κ| α0
fj0 pα0 q.

As α0 ą 0, fj0 pα0 q must therefore equal zero, a contradiction. Hence, as claimed,
the
˘ fj pαq for 0 ď j ă |κ| ´ qpκq all vanish over the interval
` polynomials
1
0, maxp|κ|´2,0q
.
But a polynomial can have an infinite number of zeros only if all its coefficients
are zero, so we conclude that
bij “ 0

for

0 ď j ă |κ| ´ qpκq.

Thus, from Equations (4.44) and (4.45) we have
“
‰ ´ n ¯|κ| qpκq 1
R|κ| pmqPκ2 pm, pq
m
E rκ pY ´1 q “
m
where
|κ|
|κ|
´ p ¯i
ÿ
ÿ
1
Pκ2 pm, pq “
bij
.
j´|κ|`qpκq
m m
i“0 j“|κ|´qpκq
Going back to equations (4.13) and (4.14) and plugging in the above yields
that as long as n ě p ` 16k ` 6,
“
‰ m2k`1 1
R2k pmqQp1q
(4.47)
E tr T 2k “
κ pm, pq,
nk
“
‰ m2k`2 1
E tr2 T k “
R2k`1 pmqQp2q
(4.48)
κ pm, pq
nk
where
1
p1q
ÿ ´
1 ¯|κ| R|κ| pmq bκ pn, m, pq 2
p
k
`
Qp1q
pm,
pq
“
p´1q
1`
P pm, pq,
κ
1 pmq
m m
R2k
m2k`1´qpκq κ
|κ|ď2k
k
Qp2q
κ pm, pq “ p´1q

1
p2q
ÿ ´
1 ¯|κ| R|κ| pmq bκ pn, m, pq 2
p
P pm, pq.
1` `
1
m m
R2k`1
pmq m2k`2´qpκq κ

|κ|ď2k`1

Now, for any a ď b, we can associate a partition µ of norm |µ| “ a with the
partition µ˚ “ pµ1 ` b ´ a, µ` 2, . . . , µqpµq q of norm |µ˚ | “ b, which satisfies
˚
˚
ˆ
qpµ
i ´1
ź q µź

i“1

j“0

1´i`2j
1´
m

ˆ
qpµq
i ´1
ź µź

˙
“

i“1 j“0

1´i`2j
1´
m

µ1 `b
˙ ´a´1
ź ˆ
j“µ1

˙
2j
1´
.
m

By definition for the Rµ pmq’s at Equation (4.42), this means that every factor
that appears in Rµ´1 pmq appears in Rµ´1
˚ pmq, so by definition of the R|µ| pmq’s
p1q

1
at Equation (4.43), Ra pmqR´1 pmq is a polynomial in m
. Moreover, as bκ and
p2q
bκ are polynomials of degrees d1 pκq ” 2k ` 1 ´ qpκq and d2 pκq ” 2k ` 2 ´ qpκq

Chételat and Wells/Mid-scale Wishart asymptotics
p1q

40

p2q

respectively, there exists coefficients cijl and cijl such that
p1q
bκ pn, m, pq
m2k`1´qpκq

“

1
md1 pκq

d1 pκq d1 pκq
dÿ
´i ´i´j
1 pκq ÿ
ÿ
i“0

j“0

d1 pκq d1 pκq
dÿ
´i ´i´j
1 pκq ÿ
ÿ

“
i“0

j“0

p1q

cijl mi nj pl

l“0

1

p1q

cijl

l“0

md1 pκq´i´j´l

´
p
1 ¯j ´ p ¯l
1` `
m m
m

and
p2q
bκ pn, m, pq
m2k`2´qpκq

“

1
md2 pκq

d2 pκq d2 pκq
dÿ
´i ´i´j
2 pκq ÿ
ÿ
i“0

j“0

d2 pκq d2 pκq
dÿ
´i ´i´j
2 pκq ÿ
ÿ

“
i“0

j“0

l“0

p2q

cijl

p2q

cijl mi nj pl

l“0

´
1 ¯j ´ p ¯l
p
`
.
1`
m m
m
md2 pκq´i´j´l
1

As d1 p´i ´ j ´ l ě 0q, j, l ě 0 in the first case and d2 pκq ´ i ´ j ´ l, j, l ě 0
in the second case, we conclude that these two expressions are polynomials
p
1
in m
and m
. Therefore, looking back at (4.47) and (4.48), we conclude that
p2q
p1q
p
1
the Qκ pm, pq’s and Qκ pm, pq’s are polynomials in m
and m
. Therefore, if
p1q
p2q
n ě p ` 16k ` 6 there must be coefficients aij , aij and large enough integers
D1 , D2 such that
D1 ÿ
i
j
ÿ
“
‰ m2k`1
p1q p
E tr T 2k “
R
pmq
a
2k
ij
nk
mi
i“0 j“0
D1
´ m ¯k
ÿ
p1q
“
R2k pmq
gi ppqmk`1´i
(4.49)
n
i“0
D2 ÿ
i
j
ÿ
“
‰ m2k`2
p2q p
E tr2 T k “
R
pmq
a
2k`1
ij
nk
mi
i“0 j“0
D2
´ m ¯k
ÿ
p2q
“
R2k`1 pmq
gi ppqmk`2´i
(4.50)
n
i“0
p1q

for polynomials gi ppq “

i
ř
j“0

p1q

p2q

p1q

i
ř

p2q
aij pj .
j“0
p2q
gi must vanish

aij pj and gi ppq “

We will now proceed to show that gi and
on N for 0 ď
i0 ă k ` 1 and 0 ď i0 ă k ` 2 respectively. Our argument relies on the analysis
of the asymptotic behavior of the moments of T in the classical regime where p
is held fixed while n grows
to‰infinity.“
“
‰
Observe first that E tr T 2k and E tr2 T k must have a finite limit as n Ñ 8
with p held fixed. Indeed, since 16T 2 {n is positive definite, |Ip ` 16T 2 {n| is

Chételat and Wells/Mid-scale Wishart asymptotics

41

greater than one and so we have the bound

ż
‰
“
E tr T 2k “ Cn,p

tr T

2k

ˇ
ˇ n`p`1
2 ˇ´ 4
ˇ
ˇIp ` 16T ˇ
dT
ˇ
n ˇ

Sp pRq

ż
ď Cn,p

ˇ
ˇ´ n
ˇ
16T 2 ˇˇ 4
tr T 2k ˇˇIp `
dT.
n ˇ
Sp pRq

pp3p`1q
4

ppp`1q
4

When p is held fixed, lim Cn,p “ 2
{π
by Lemma 1. Moreover,
nÑ8
n
n
ˇ
ˇ
ˆ
˙
p
´
´
2ˇ 4
4
ź
ˇ
λi p4T 2 q
ˇIp ` 16T ˇ
“
1`
ˇ
ˇ
n
n{4
i“1
for λ1 p4T 2 q ě ¨ ¨ ¨ ě λp p4T 2 q ě 0 the eigenvalues of 4T 2 , and p1 ` x{nq´n is
monotone decreasing towards exppxq. Therefore, for a fixed dimension p we can
apply the monotone convergence theorem to obtain that

ż
pp3p`1q

“
‰
2 4
lim E tr T 2k ď
ppp`1q
nÑ8
π 4

“
‰
2
tr T 2k e´4 tr T dT “ E tr Z 2k ă 8

(4.51)

Sp pRq

for Z „ GOEppq{4. Repeating the argument with tr2 T k yields similarly

ż
pp3p`1q
4

“
‰
2
lim E tr2 T k ď
ppp`1q
π 4

nÑ8

“
‰
2
tr2 T k e´4 tr T dT “ E tr2 Z k ă 8.

(4.52)

Sp pRq

“
‰
“
‰
Thus indeed E tr T 2k and E tr2 T k have finite limits when p is held fixed.
p1q

p2q

We can use this to show that gi and gi must vanish on N for 0 ď i0 ă k `1
and 0 ď i0 ă k ` 2 as follows. Say the first statement wasn’t true, and let
p1q
0 ď i0 ă k ` 1 be the smallest i such that for some p0 P N, gi0 pp0 q ‰ 0. Then
“
‰
by Equation (4.49) and the definition of i0 , the limit of E tr T 2k as n Ñ 8
with p fixed at p0 satisfies
“
‰
D1
ÿ
E tr T 2k
p1q
p1q
k
“
1
¨
1
¨
lim
gi pp0 qmi0 ´i “ gi0 pp0 q.
lim
nÑ8
nÑ8 mk`1´i0
i“i
0

But m “ n´p´1 tends to infinity as n tends to infinity, and since k `1´i0 ą 0,
“
‰
p1q
Equation (4.51) means that E tr T 2k {mk`1´i0 must tend to zero. Thus gi0 pp0 q
has to equal zero, which contradicts our assumption. Thus for every 0 ď i ă k`1,
p1q
the polynomial gi must vanish on N.
p2q
Similarly, for 0 ď i ă k ` 2, the polynomial gi must vanish on N, because
if it wasn’t the case, we could take 0 ď i0 ă k ` 2 as the smallest i with the
p2q
property that for some p0 P N, gi0 pp0 q ‰ 0, and then by Equation (4.50) with

Chételat and Wells/Mid-scale Wishart asymptotics

42

p fixed at p0 as n Ñ 8 we would get
“
‰
D2
ÿ
E tr2 T k
p2q
p2q
k
“
1
¨
1
¨
lim
gi pp0 qmi0 ´i “ gi0 pp0 q.
lim
nÑ8
nÑ8 mk`2´i0
i“i
0

But then by Equation (4.52), as m tends to infinity and k ` 2 ´ i0 ě 1 the
“
‰
p2q
ratio E tr2 T k {mk`2´i0 must tend to zero. Thus we must have gi0 pp0 q “ 0, a
p2q

contradiction. Hence indeed for 0 ď i ă k ` 2, the polynomial gi must vanish
on N.
But of course, a polynomial can only have an infinite number of zeroes if its
coefficients are all zero, so we must have
p1q

for

0 ď i ă k ` 1,

p2q

for

0 ď i ă k ` 2.

aij “ 0
aij “ 0

(4.53)

Now say that p varies with n in such a way that limnÑ8 p{n “ α ă 1. Then
for large enough n, n ě p ` 16k ` 6 so by Equations (4.49) and (4.53),
D1 ÿ
i
´ m ¯k
j´pk`1q
ÿ
“
‰
1
p1q p
R2k pmq
aij i´pk`1q
lim k`1 E tr T 2k “ lim
nÑ8 n
nÑ8 p
m
i“k`1 j“0
«
p1q
D1 ÿ
k
ÿ
aij
“ p1 ´ αqk ¨ 1 ¨ lim
i´pk`1q ppk`1q´j
nÑ8
m
i“k`1 j“0
ff
D1
i
´ p ¯j´pk`1q 1
ÿ
ÿ
p1q
`
aij
m
mi´j
i“k`1
j“k`1
ff
«
˙
ˆ
p1q
D1
i´pk`1q
k
ÿ
ÿ
ak`1
α
p1q
(4.54)
“ p1 ´ αqk
`
a
`
˘pk`1q´j
ii
1´α
j“0 lim p
i“k`1
pÑ8

and by Equations (4.50) and (4.53),
D2 ÿ
i
´ m ¯k
j´pk`2q
ÿ
“
‰
1
p2q p
lim k`2 E tr2 T k “ lim
R2k`1 pmq
aij i´pk`2q
nÑ8 p
nÑ8 n
m
i“k`2 j“0
«
p2q
D
k`1
ÿ2 ÿ
aij
“ p1 ´ αqk ¨ 1 ¨ lim
i´pk`2q
nÑ8
m
ppk`2q´j
i“k`2 j“0
ff
D1
i
´ p ¯j´pk`2q 1
ÿ
ÿ
p2q
`
aij
m
mi´j
i“k`2 j“k`2
«
ff
ˆ
˙
p2q
D2
i´pk`2q
k`1
ÿ
ÿ
ak`2
α
p2q
“ p1 ´ αqk
aii
.
(4.55)
`
˘pk`2q´j `
1´α
j“0 lim p
i“k`2
pÑ8

Although we might not know what aij coefficients are, this shows at least that
the limits are finite.
In‰ particular, from “Equations
(4.54) and (4.55) we can
“
‰
conclude that E tr T 2k “ Oppk`1 q and E tr2 T k “ Oppk`2 q, which shows the
first claim of the theorem.
For the second claim, let n, p Ñ 8 with p{n Ñ α “ 0. Then Equations (4.54)

Chételat and Wells/Mid-scale Wishart asymptotics

43

and (4.55) specialize to
lim

nÑ8

1
pk`1
1

lim

nÑ8 pk`2

“
‰
p1q
E tr T 2k “ apk`1qpk`1q ,
“
‰
p2q
E tr2 T k “ apk`2qpk`2q .

(4.56)

What is interesting about this result is that these limits must be the same
regardless of the way p grows! As long as p Ñ 8 with p{n Ñ 0, the limits are
?
p1q
p2q
apk`1qpk`1q and apk`2qpk`2q , regardless of whether p „ log n or p „ n or some
other growth.
Now, Bubeck et al. [2016, Theorem 7] and Jiang and Li [2015, Theorem 1]
have shown that when p Ñ 8?with p3 {n Ñ 0, the total variation distance
between a normalized Wishart npWp pn, Ip {nq ´ Ip q matrix and a Gaussian
Orthogonal Ensemble GOEppq matrix tends to zero as n Ñ 8. Therefore, the
Hellinger distance satisfies also H 2 pψNW , ψGOE q “ H 2 pfNW , fGOE q Ñ 0 as n Ñ
8.
But convergence in Hellinger distance has strong implications for real-valued
statistics. Indeed, for T1 „ Tn{2 pIp {8q “ |ψNW |, T2 „ GOEppq{4 “ |ψGOE | and
any function g : Sp pRq Ñ R such that gpT1 q, gpT2 q are square-integrable,
ˇ
ˇ ˇˇ
”
ı ˇˇ
ˇ
ˇ
gpT q |ψNW |pT q ´ |ψGOE |pT q dT ˇˇ
ˇ ErgpT1 qs ´ ErgpT2 qs ˇ “ ˇˇ

ż

Sp pRq

ż
ˇ
ˇ
ď ˇˇ

”
ı ˇˇ
1{2
1{2
1{2
gpT qψNW pT q ψNW pT q ´ ψGOE pT q dT ˇˇ
Sp pRq

ˇ
ˇ
` ˇˇ

ż
”
ı ˇˇ
1{2
1{2
1{2
gpT qψGOE pT q ψNW pT q ´ ψGOE pT q dT ˇˇ
Sp pRq

ż

ż

„


1

gpT q2 |ψNW |pT qdT 2 `

ď

Sp pRq

ż
¨

1

gpT q2 |ψGOE |pT qdT 2
Sp pRq

ˇ
ˇ2 1
ˇ 1{2
ˇ
1{2
ˇψNW pT q ´ ψGOE pT qˇ dT 2
S pRq

p
1
1
“
‰
“
‰
“ E gpT1 q2 2 ` E gpT2 q2 2 HpψNW , ψGOE q
„

(4.57)

by the Cauchy-Schwarz inequality.
Let’s consider applying this result to gpT q “ tr T 2k {pk`1 and gpT q “ tr2 T k {pk`2 .
What do we know about these statistics? In the case where T2 „ GOEppq{4,
results of Anderson et al. [2010, Lemma 2.1.6 and the equation above Equation
(2.1.21)] provide that
„

tr T22k
Ck
lim E k`1
“ 2k
nÑ8
p
4

Chételat and Wells/Mid-scale Wishart asymptotics

44

„ 2 k
˙2
´ „ tr T k  ¯2 ˆ C
tr T
k{2
2
lim E k`22
“
1
rk
evens
“ lim E k{2`1
nÑ8
nÑ8
p
4k
p
„´

„

ˆ
˙2
´
tr T22k ¯2
tr T22k ¯2
Ck
lim E
“ lim E k`1
ă8
“
nÑ8
nÑ8
pk`1
p
42k
„
„´ 2 k ¯ 

tr T24k
tr T2 2
C2k
ď
lim
E
lim E
“ 4k ă 8
nÑ8
nÑ8
pk`2
p2k`1
4
because these expressions only depend on p, and since p Ñ 8 as n Ñ 8, taking
a limit as n Ñ 8 is the same as taking a limit as p Ñ 8. Moreover, in the
T1 „ Tn{2 pIp {8q case, using Jensen’s inequality and Equations (4.54)–(4.55) we
can at least see that
„´

„

tr T12k ¯2
tr T14k
lim E
ď lim E 2k`1 ă 8
nÑ8
nÑ8
pk`1
p
„ 2 2k 
„´ 2 k ¯ 
2
tr T1
tr T1
ď lim E
ă 8.
lim E
nÑ8
nÑ8
pk`2
pk`2
Therefore, using Equation (4.57) with gpT q “ tr T 2k {pk`1 and gpT q “ tr2 T k {pk`2
we find that when n, p Ñ 8 with p3 {n Ñ 0,
ˇ
ˇ
„

ˆ
„´
1
˙
2k
ˇ
Ck ˇˇ
tr T12k ¯2 2
C2k
ˇ lim E tr T1
´ 2k ˇ ď
lim E
` 2k ¨ 0 “ 0,
ˇ nÑ8
nÑ8
pk`1
4
pk`1
4
ˇ
ˇ
„ 2 k ´
¯
2ˇ
ˇ
Ck{2
ˇ lim E tr T1 ´
1rk evens ˇˇ
ˇ nÑ8
k`2
k
p
4
ˆ
„´ 2 k ¯  12
1{2 ˙
C2k
tr T1 2
ď
lim E
` 2k ¨ 0 “ 0.
nÑ8
pk`2
4
Since p3 {n Ñ 0 implies p{n Ñ 0, we conclude from Equation (4.56) that
2
Ck{2
Ck
p1q
p2q
2
apk`1qpk`1q “ 2k ,
apk`2qpk`2q “ 2k 1rk evens .
4
4
But then, by that same equation, we conclude that when n, p Ñ 8, not only
when p3 {n Ñ 0 but for all p such that p{n Ñ 0, we have
“
‰
Ck
1
lim
E tr T 2k “ 2k ,
nÑ8 pk`1
4
(4.58)
2
“ 2 k ‰ Ck{2
1
2
lim k`2 E tr T “ 2k 1rk evens
nÑ8 p
4
for
“ T „‰Tn{2 pIp {8q. To finish the proof, use Equation (4.58) with the fact that
E tr T k “ 0 for odd k to find that
„

„ 2 k
„
2
tr T k
tr T
tr T k
lim Var k{2`1 “ lim E k`2 ´ E k{2`1
nÑ8
nÑ8
p
p
p
2
´
¯2
Ck{2
Ck{2
2
“ 2k 1rk evens ´
1
rk
evens
“ 0.
4
4k
L2

Thus tr T k {pk`1 ÝÑ Ck {42k , as desired. This proves the second claim and
concludes the proof.

Chételat and Wells/Mid-scale Wishart asymptotics

45

A pleasant consequence of this result is that when n, p Ñ 8 with p{n Ñ 0, we
can conclude a version of the semicircle law holds for the Tn{2 p2Ip q distribution.
This is interesting because the Tn{2 p2Ip q distribution has dependent entries with
heavy tails, whose distribution varies with n, p.
?
?
?
?
Let 4T { p „ 4Tn{2 pIp {8q{ p, with eigenvalues λ1 p4T { pq ě ¨ ¨ ¨ ě λp p4T { pq.
Then define its empirical spectral measure to be
p
?
1ÿ
L4T {?p pAq “
1rλi p4T { pq P As .
p i“1
Since L4T {?p depends on the random matrix T , it is a random measure on R.
Corollary 1 (Semicircle law for the t distribution). The empirical spectral
?
measure L4T {?p of a 4Tn{2 pIp {8q{ p random matrix converges weakly, in square
mean, to the semicircle distribution
ż ?
4 ´ x2
1r|x| ď 2s dx.
LpAq “
2π
A
Proof. Let f be any continuous function R Ñ R that vanishes at infinity. By
the Stone-Weierstrass theorem, there exists a sequence f1 , f2 , . . . of polynomials
such that for any  ą 0, supxPR |f pxq ´ fl pxq| ă . To fix some notation, write
fl pxq “

degf
ÿl

alk xk .

k“1

Then since L4T {?p and L are both probability measures,
«ˆ ż
˙2 ff 12
ż
f pxqdL4T {?p pxq ´ f pxqdLpxq
E
R

R

«ˆ ż
˙2 ff 21
“
‰
f pxq ´ fl pxq dL4T {?p pxq
ď E
R

«ˆ ż
˙2 ff 21
ż
fl pxqdL4T {?p pxq ´ fl pxqdLpxq
`E
R

R

«ˆ ż
˙2 ff 21
“
‰
`E
fl pxq ´ f pxq dLpxq
R

ď `

degf
ÿl
k“1

«ˆ
˙2 ff 21
1 ´ 4T ¯k
|alk | E
tr ?
´ Ck{2 1rk evens
` .
p
p

By Theorem 2, the expectation tends to zero as n, p Ñ 8 with p{n Ñ 0. Thus
«ˆ ż
˙2 ff 21
ż
lim E
f pxqdL4T {?p pxq ´ f pxqdLpxq
ď 2.
nÑ8

R

R

But this is true for every  ą 0, so the limit must
ş be zero. Hence for every continuous f that vanishes at infinity, the integral f dL4T {?p converges in square

Chételat and Wells/Mid-scale Wishart asymptotics

46

ş
mean to f dL. By Chung [2001, Theoremş 4.4.1 and 4.4.2], this implies that for
every bounded continuous f , the integral f dL4T {?p converges in square mean
ş
to f dL. Thus the empirical spectral distribution L4T {?p converges weakly, in
square mean, to the semicircle distribution L, as desired.
5. Wishart asymptotics: the G-transform point-of-view
We now turn our attention to the main objective of this paper, namely studying
the behavior of Wishart matrices in the various middle-scale regimes. To do
this, we exploit the close connection between the Wishart and the symmetric
t distributions and make use of the results found Section 4. The main result
of this section, Theorem 3, states that we can approximate for every middlescale regime the G-transform ψNW of a normalized Wishart by a degree-specific
function ψK . This can be seen as an analogue of Theorem 1 in the G-transform
domain.
The reasoning behind the approximations is as follows. We could imagine
writing ψNW from Proposition 3 in exponential form, and expanding the terms
as a Taylor series would yield
ˇ
ˇ*
"
ˇ
?
4T ˇˇ
n`p`1
ˇ
log ˇIp ` i ? ˇ
ψNW pT q “ Cn,p exp 2i n tr T ´
2
n
"
ˆ
˙k *
8
´
?
n ` p ` 1 ÿ p´iqk p ¯ k2
4T
tr ?
.
“ Cn,p exp 2i n tr T `
2
k
n
p
k“1
Now imagine that the T ’s appearing in the expression follow a Tn{2 pIp {8q dis` 4T ˘k
tribution. By Theorem 2, we know that tr ?
“ Θppq when k is even, in
p
` 4T ˘k
2
an L sense. When k is odd, the theorem merely proves that tr ?
“ oppq,
p
k

but for a GOEppq matrix Z, we know that tr ?Zp is asymptotically normal
for odd k by Anderson et al. [2010, Theorem 2.1.31]. This would suggest that
` 4T ˘k
tr ?
“ Θp1q when k is odd. Thus we would have, in some sense,
p
"ÿ
ˆ k{2`1 ˙ ÿ
ˆ k{2 ˙ *
8
8
p
p
ψNW pT q “ Cn,p exp
Θ
`
Θ
k{2´1
k{2´1
n
n
k“2
k“2
even

“ Cn,p exp

"ÿ
8
k“0

odd

ˆ
Θ

k`2

p
nk

8
ÿ

˙
`

k“1

˜c
Θ

pp2k´1q`2
n2k´1

¸*
.

In other words, terms in the power series would be associated with some degree
K, such that they would be non-negligible in any middle-scale regime of degree
up to K, and negligible in higher degrees. In fact, a similar phenomenon occurs
with Cn,p , by Lemma 1. This suggests we should try truncating these power
series to derive degree-specific approximations.
Definition 3 (G-transform approximations). For any K P N, define the K th

Chételat and Wells/Mid-scale Wishart asymptotics

47

degree approximation ψK : Sp pRq Ñ C as
$ 2K`3`
,
2K`2´
’
/
odds
1
rK odds
& 1rK
.
ÿ
ÿ
¯
¯
´
´
k
k
k
k
pKq
p `1
4i tr T
4i tr T
ψK pT q “ Cn,p exp n
?
?
`
’
/
%2
n
k
2
n
k k“2

with
pKq
Cn,p
“

2

pp3p`1q
4

π

ppp`1q
4

k“1

" K`1
1 ÿ 1rk evens pk`2
exp ´
2 k“1 kpk`1qpk`2q nk
*
K`1
1 ÿ 1`21rk evens pk`1
´
.
4 k“1
kpk`1q
nk

Just like the G-transform of a normalized Wishart matrix, these functions
implicitly depend on n. The first three are
"
*
pp3p`1q
32
p`1
p`1
2 4
tr T 2 ,
ψ0 pT q “ ppp`1q exp ´4 tr T 2 ´i ? tr T 3 `2i ? tr T ´4
3 n
n
n
π 4
" 2
pp3p`1q
2 4
32
32
512
p
ψ1 pT q “ ppp`1q exp ´
´4 tr T 2 ´i ? tr T 3 ` tr T 4 `i 3{2 tr T 5
8n
3 n
n
5n
π 4
*
p
`
1
p
`
1
32pp
` 1q
1024
6
2
3
tr T `2i ? tr T ´4
tr T ´i
tr T
´
3n2
n
n
3n3{2
and
"
pp3p`1q
2 4
p4
32
32
p2
ψ2 pT q “ ppp`1q exp ´
´4 tr T 2 ´i ? tr T 3 ` tr T 4
´
2
48n
8n
3
n
n
4
π
512
1024
8192
p`1
`i 3{2 tr T 5 ´
tr T 6 ´i 5{2 tr T 7 `2i ? tr T
2
3n
n
5n
7n
p`1
32pp
`
1q
p
`
1
´4
tr T 2 ´i
tr T 3 `32 2 tr T 4
3{2
n
n
3n
*
512pp ` 1q
1024pp ` 1q
5
6
`i
tr
T
´
tr
T
.
3n3
5n5{2
These functions have the pleasant property that their modulus is bounded, up
to a constant, by the G-conjugate density |ψK |. Indeed, on one hand we can
rewrite Definition 3 into
K`1`
$
1rK
odds
’
? 2k
ÿ
&
n
tr p4T { nq
pKq
ψK pT q “ exp log Cn,p `
p´1qk
’
2
2k
%
k“1

K`1´
1rK odds

p`1
`
2

ÿ

p´1q
k“1

?

k tr p4T {

2k

2k

nq

K`1

n
´i
2

ÿ
k“1

p´1qk

? 2k`1
tr p4T { nq
2k ` 1

Chételat and Wells/Mid-scale Wishart asymptotics

48

,
? 2k`1 /
.
tr p4T { nq
´i
p´1qk
. (5.1)
/
2
2k ` 1
k“1
?
`
˘
On the other hand, for any x P R, we can write 1 ` ix “ 1 ` x2 exp i atanpxq
with the arctangent function taking values in p´π{2, π{2q. Thus by Proposition
3 we can rewrite ψNW as
ˇ
ˇ
"
ˇ
16T 2 ˇˇ
n`p`1
ˇ
log ˇIp `
ψNW pT q “ exp log Cn,p ´
4
n ˇ
ˆ
˙
*
?
n`p`1
4T
´i
tr atan ?
` 2i n tr T ,
(5.2)
2
n
with the understanding that the matrix-variate arctangent function operates on
eigenvalues by functional calculus. Now, for any x P R and odd integer L, there
is an elementary inequality
L
ÿ
1
x2l
ď ´ logp1 ` x2 q.
p´1ql
2l
2
l“1
K
ÿ
p`1

Notice that K ` 1 ˘ 1rK odds is always an odd integer. Thus, from the above
inequality and Equations (5.1) and (5.2), we can derive the bound
ˇ
! ´n p ` 1¯
ˇ ˇ
16T 2 ˇˇ)
ˇ
pKq
ˇψK ˇpT q ď Cn,p
exp ´
`
log ˇIp `
ˇ
4
4
n
pKq
ˇ
Cn,p ˇˇ
“
ψNW ˇpT q.
(5.3)
Cn,p
In particular, since ψNW is integrable whenever n ě p ´ 2 by Proposition 3,
Equation (5.3) implies that every ψK must also be integrable whenever n ě p´2.
In particular, for large enough n it makes sense to talk about the asymptotic
total variation or Hellinger distance between ψNW and ψK .
We now state the main result, which is that each function ψK approximates
the G-transform of a normalized Wishart for all middle-scale regimes of degree
K or lower, but no other.
log p
Theorem 3.Let lim log
n ă 1 as n Ñ 8. For any K P N, the total varianÑ8
tion
distance
between
the
G-transform of the normalized Wishart distribution
?
nrWp pn, Ip {nq ´ Ip s and the K th approximating function ψK satisfies

ż
`

˘

dTV ψNW , ψK “

ˇ
ˇ
ˇψNW pT q ´ ψK pT qˇ dT Ñ 0

as n Ñ 8

Sp pRq

if and only if pK`3 {nK`1 Ñ 0.
Proof. If statement. For the first part of the theorem, remark that by Equation
(3.7) it is equivalent to show that the Hellinger distance tends to zero, i.e. that

ż
`
˘
H2 ψNW , ψK “

ˇ
ˇ2
ˇ 1{2
ˇ
1{2
ˇψNW pT q ´ ψK pT qˇ dT Ñ 0
Sp pRq

as n Ñ 8

Chételat and Wells/Mid-scale Wishart asymptotics

49

when pK`3 {nK`1 Ñ 0. To control this quantity, we use the Kullback-Leibler
inequality for G-transforms. Notice that for any x P R and L P N,
ˇ
ˇ
L´1
ˇ 1
2l ˇ
ÿ
x2L
ˇ
2
lx ˇ
p´1q
,
(5.4)
ˇ ď
ˇ´ logp1 ` x q ´
ˇ 2
2l ˇ
2L
l“1
ˇ
ˇ
L´1
ˇ
2l`1 ˇ
ÿ
x2L
x
ˇ
ˇ
p´1ql
.
(5.5)
ˇ ď
ˇatanpxq ´
ˇ
2l ` 1 ˇ
2L
l“1
Let Log stand for the principal branch of the complex logarithm, and let us
study Log ψNW {ψK . Its real part can be bounded by
ˇ
ˇ
ˇ ˇˇ
ˇ
ˇ
2ˇ
ˇ
ˇ
ˇ
pKq
ˇ< Log ψNW pT q ˇ “ ˇˇ log Cn,p ´ n ` p ` 1 log ˇIp ` 16T ˇ ´ log Cn,p
ˇ
ˇ
ψK pT q ˇ ˇ
4
n ˇ
ˇ
K`1`

K`1´

k“1

k“1

1rK
odds
1rK
odds
? 2k
? 2k
ÿ
ÿ
n
p`1
tr p4T { nq
k tr p4T { nq
´
´
p´1q
p´1qk
2
2k
2
2k

ˇ
ˇ
ˇ
ˇ
ˇ
ˇ
ˇ

K`1`
ˇ
ˇ
1rK
odds
ˇ
ˇ
ˇ
ˇ
ˇ
ˇ
?
ÿ
2k
ˇ
ˇ
2ˇ
ˇ
ˇ nˇ 1
ˇ
16T
tr
p4T
{
nq
ˇ
pKq
k
ˇ´
ď ˇˇ log Cn,p ´log Cn,p ˇˇ ` ˇ´ log ˇˇIp `
p´1q
ˇ
ˇ
2ˇ 2
n ˇ
2k
ˇ
ˇ
k“1
K`1´
ˇ
ˇ
odds
ˇ
ˇ
ˇ 1rK
? 2k ˇˇ
ÿ
ˇ
2
ˇ
16T ˇˇ
p`1ˇ 1
tr p4T { nq ˇ
´
`
p´1qk
ˇ´ log ˇˇIp `
ˇ.
ˇ
ˇ
2 ˇ 2
n
2k
ˇ
ˇ
k“1
By Equation (5.4), this can be bounded by
ˇ
ˇ
?
ˇ
ˇ n trp4T { nq2K`4`21rK odds
pKq ˇ
ˇ
ď ˇ log Cn,p ´log Cn,p ˇ `
2 2K ` 4 ` 21rK odds
?
p ` 1 trp4T { nq2K`4´21rK odds
`
2
2K ` 4 ´ 21rK odds
ˇ
ˇ
2K`3`21rK odds
ˇ
ˇ
4
tr T 2K`4`21rK odds
pKq ˇ
“ ˇˇ log Cn,p ´log Cn,p
ˇ ` K ` 2 ` 1rK odds nK`1`1rK odds

42K`3´21rK odds pp ` 1q tr T 2K`4´21rK odds
.
(5.6)
K ` 2 ´ 1rK odds
nK`2´1rK odds
We can bound the imaginary part of Log ψNW {ψK in a similar way. Define
x
Pp´π,πs : R Ñ p´π, πs to be the projection mapping Pp´π,πs x “ x ´ 2πr 2π
´ 12 s.
A plot is given as Figure 4.
It satisfies = Log z “ Pp´π,πs = log z for all branches of log z, as well as the
inequality |Pp´π,πs x| ď |x|. Using this mapping, we can see that the imaginary
`

Chételat and Wells/Mid-scale Wishart asymptotics

50

y
π
-5π

-3π

-π

π

3π

5π

x

-π
Fig 4. Plot of Pp´π,πq on p´5π, 5πs.

part of Log ψNW {ψK can be bounded as
ˇ
»
ˇ ˇˇ
ˇ
ˆ
˙
ˇ
ˇ
?
n`p`1
4T
ˇ= Log ψNW pT q ˇ “ ˇˇ Pp´π,πs—
tr atan ?
` 2 n tr T
–´
ˇ
ˇ
ˇ
ψK pT q
2
n
ˇ
fiˇ
K
? 2k`1
? 2k`1 ˇˇ
ÿ
p`1
tr p4T { nq
tr p4T { nq
ffiˇ
`
`
p´1qk
p´1qk
flˇ
ˇ
2
2k ` 1
2
2k ` 1
ˇ
k“1
k“1
ˇ
ˇ
ˇ
K`1
ˆ
˙
?
? 2k`1 ˇˇ
ÿ
4T
2 n tr T
n ˇˇ
ˇ
k tr p4T { nq
´
´
ď ˇ tr atan ?
p´1q
ˇ
ˇ
2ˇ
n
n{2
2k ` 1
ˇ
ˇ
k“1
ˇ
ˇ
ˇ
K
˙ ÿ
ˆ
? 2k`1 ˇˇ
ˇ
p`1ˇ
4T
tr p4T { nq
ˇ
`
´
p´1qk
ˇ tr atan ?
ˇ.
ˇ
2 ˇ
n
2k ` 1
ˇ
ˇ
k“1
K`1

ÿ
n

By Equation (5.5), this can be bounded by
?
?
p ` 1 trp4T { nq2K`2
n trp4T { nq2K`4
`
ď
2
2K ` 4
2
2K ` 2
2K`3
2K`4
2K`1
4
tr T
4
pp ` 1q tr T 2K`4
ď
`
.
(5.7)
K`1
K `2 n
K `1
nK`1
Recall that the G-conjugate of the normalized Wishart distribution is the t distribution with n{2 degrees of freedom and scale matrix Ip {8, denoted Tn{2 pIp {8q
– see Equation (3.24) and Section 4 for details. Let us bound the expectations of
these absolute real and imaginary parts under this distribution. By Equations
(5.10), (5.6), (5.7) and Theorem 2, we find that for T „ |ψNW | “ Tn{2 pIp {8q,
‰
“
ˇ
„ˇ
2K`3`21rK odds E tr T 2pK`2`1rK oddsq
ˇ
ˇ
ψ
pT
q
4
NW
ˇ ď
E ˇˇ< Log
ψK pT q ˇ
K ` 2 ` 1rK odds
nK`1`1rK odds
“
‰
´ pK`3 ¯
42K`3´21rK odds pp ` 1q E tr T 2pK`2´1rK oddsq
`
`
o
K ` 2 ´ 1rK odds
nK`1
nK`2´1rK odds
ˆ K`2`1rK odds`1 ˙
ˆ
˙
´
K`2´1rK odds`1
K`3 ¯
p
p¨p
p
ď O
`O
` o K`1
K`1`
1
odds
K`2´
1
odds
rK
rK
n
n
n

Chételat and Wells/Mid-scale Wishart asymptotics

“ O

51

´ pK`3 ¯

(5.8)

nK`1

and
“
‰
“
‰
ˇ
„ˇ
2K`3 E tr T 2pK`2q
2K`1 pp`1qE tr T 2pK`2q
ˇ
ˇ
ψ
pT
q
4
4
NW
ˇ ď
E ˇˇ= Log
`
ψK pT q ˇ
K `2
nK`1
K `1
nK`1
ˆ K`2`1 ˙
ˆ
˙
´ pK`3 ¯
p
p ¨ pK`2`1
ď O
(5.9)
`O
“ O K`1
K`1
K`2
n
n
n
as n Ñ 8 with p{n Ñ 0.
Moreover, from Lemma 1, we see that
" ˆ K`3 ˙*
p
p
pKq
Cn,p “ Cn,p exp o
as n Ñ 8 with Ñ 0.
(5.10)
K`1
n
n
Thus, from Equations (5.3) and (5.10), we see that when pK`3 {nK`1 Ñ 0, the
asymptotic L1 norm of ψK is bounded by

ż

ż
!

lim

nÑ8

ˇ
ˇ
ˇψK pT qˇ dT ď lim exp ´o
nÑ8

´ pK`3 ¯)
nK`1

Sp pRq

ˇ
ˇ
ˇψNW pT qˇ dT “ 1.

(5.11)

Sp pRq

In fact, at the end of this proof we will see that this bound is sharp and the
limit must be exactly 1.
Using Equations (5.8), (5.9) and (5.11) with Proposition (1) implies that
when pK`3 {nK`1 Ñ 0,
´
¯
0 ď lim H2 ψNW , ψK
nÑ8

ż

ż
„
ď

lim

nÑ8


|ψK |pT q dT ´ 1 ` 0 ` 2 lim

|ψK |pT q dT 2 ¨ 0 2

Sp pRq

Sp pRq

nÑ8

1
2

1
2

1

1

ď r1 ´ 1s ` 0 ` 2 ¨ 1 ¨ 0 “ 0.
`
˘
Thus` H2 ψNW˘, ψK Ñ 0, hence by Equation (3.7) we must have the limit
dTV ψNW , ψK Ñ 0, as desired.
Only if statement. For the second part of the theorem, assume
that˘ the total
`
variation distance satisfies dTV pψNW , ψK q Ñ 0, hence H ψNW , ψK Ñ 0 by
Equation (3.7), as n Ñ 8. We will show by contradiction this implies that
pK`3 {nK`1 Ñ 0.
log p
Assume this wasn’t the case. Since limnÑ8 log
n ă 1, there must be an L P N
such that pL`3 {nL`1 Ñ 0, and since pK`3 {nK`1 Û 0, we must have K ă L.
By Equation (5.8), we must have for T „ |ψNW | “ Tn{2 pIp {8q that
ˇ
„ˇ
´ pL`3 ¯
ˇ
ψNW pT q ˇˇ
ˇ
lim E ˇ< Log
ď lim O L`1 “ 0,
ˇ
nÑ8
nÑ8
ψL pT q
n
so
1
ψNW pT q L1
< Log
ÝÑ 0.
(5.12)
2
ψL pT q

Chételat and Wells/Mid-scale Wishart asymptotics

52

Now write, by Equation (5.1) and Definition 3,
L`1`

RpT q ”

1
ψL pT q
1
< Log
“ log
2
ψK pT q
2

1rLÿ
odds
ˆ
˙2k
n
p´1qk
4T
`
tr ?
4
2k
n

pLq
Cn,p
pKq
Cn,p

K`2`
k“1rK odds
L`1´

1rLÿ
odds
ˆ
˙2k
p`1
p´1qk
4T
tr ?
`
4
2k
n
K`2´
k“1rK odds

“ ´

L`1
ÿ 1`21rk evens pk`1
1 ÿ
1rk evens pk`2 1 L`1
´
4 k“K`2 kpk`1qpk`2q nk
8 k“K`2 kpk`1q
nk

L`1`

`

1rLÿ
odds

n
4

L`1´

p´1qk ´ p ¯k
tr
2k
n

ˆ

4T
?
p

˙2k
`

p`1
4

K`2`
k“1rK odds

1rLÿ
odds

ˆ
˙2k
4T
p´1qk ´ p ¯k
tr ?
.
2k
n
p

K`2´
k“1rK odds

But as pL`3 {nL`1 , we must have p{n Ñ 0, so by Theorem 2, we have

1
p

L2

4T 2k
trp ?
Ñ
pq

Ck . Moreover, as we assumed that pK`3 {nK`1 Û 0, we must have p Ñ 8. Thus
L`1
1 ÿ
1rk evens pk´K´1
nK`1
RpT
q
“
´
pK`3
4 k“K`2 kpk`1qpk`2q nk´K´1
L`1`

1
1`21rk evens 1 pk´K´1
1
`
8 k“K`2 kpk`1q
p nk´K´1
4
L`1
ÿ

´

1rLÿ
odds

p´1qk pk´K´2 1
tr
2k nk´K´2 p

ˆ

4T
?
p

˙2k

ˆ

4T
?
p

˙2k

K`2`
k“1rK odds
L`1´

1
`
4
L2

ÝÑ 0 ` 0 `
“

1rK evens
8pK ` 2q

ˆ

1
1`
p

odds
˙ 1rL
ÿ

p´1qk pk´K´1 1
tr
2k nk´K´1 p

K`2´
k“1rK odds

CK`2 `

1rK odds
8pK ` 1q

CK`1

CK`1`1rK evens
ą 0.
8pK ` 1 ` 1rK evensq

Then by the reverse triangle inequality,

ż
`
˘
0 “ lim H2 ψNW , ψK “ lim
nÑ8

nÑ8

ˇ
ˇ2
ˇ 1{2
ˇ
1{2
ˇψNW pT q ´ ψK pT qˇ dT
Sp pRq

(5.13)

Chételat and Wells/Mid-scale Wishart asymptotics

53

ż
ď

lim

nÑ8

ˇ2
ˇ
ˇ
ˇ
ˇ|ψNW |1{2 pT q ´ |ψK |1{2 pT qˇ dT
Sp pRq

«ˇ
ˇ2 ff
"
*
ˇ
ˇ
1
ψ
pT
q
K
“ lim E ˇˇ exp
< Log
´ 1ˇˇ
nÑ8
2
ψNW pT q
for a T „ |ψNW | “ Tn{2 pIp {8q, that is
"
*
1
ψK pT q
L2
exp
< Log
ÝÑ 1.
2
ψNW pT q
Since Lp convergence implies convergence in probability, by the continuous mapping theorem we must have
ψK pT q P
1
< Log
ÝÑ 0
2
ψNW pT q
as n Ñ 8, so by Equation (5.12)
ψK pT q
1
ψNW pT q P
1
´ < Log
ÝÑ ´0 ´ 0 “ 0.
RpT q “ ´ < Log
2
ψNW pT q 2
ψL pT q
But then, from Equation (5.13) and Slutsky’s lemma [van der Vaart, 2000,
Lemma 2.8 (iii)],
ˆ K`1
˙´1
pK`3
n
P 8pK ` 1 ` 1rK evensq
¨ 0 “ 0.
“
RpT
q
RpT q ÝÑ
K`1
K`3
n
p
CK`1`1rK evens
K`3
as n Ñ 8. As pK`3 {nK`1 is deterministic, this
{nK`1 Ñ 0
` implies ˘that p
2
as n Ñ 8, a contradiction. Thus whenever H ψNW , ψK Ñ 0 as n Ñ 8 with
lim log p ă 1, we must have pK`3 {nK`1 Ñ 0, as desired. This concludes the
nÑ8 log n
proof.

Although Theorem 3 states that the functions ψK approximate ψNW , there
is no guarantee that they are G-transforms of a probability density. In other
words, nothing guarantees that their inverse G-transforms f˜K “ G ´1 tψK u are
real-valued, non-negative and integrate to unity. However, the reverse triangle
inequality applied to the L2 -norm provides that

ż
ˇ
ˇ
ˇ

ż
|ψNW pT q| dT 2 ´

ˇ
`
˘
1ˇ
|ψK pT q| dT 2 ˇ ď H ψNW , ψK ,

Sp pRq

Sp pRq

1

so Theorem 3 and the Plancherel theorem implies that

ż
lim

nÑ8

ż
|f˜K pT q| dT “ lim

|ψK pXq| dX “ 1

Sp pRq

Sp pRq

nÑ8

(5.14)

when pK`3 {nK`1 Ñ 0. That is, the theorem at least guarantees that |f˜K | is
asymptotically a density in its associated regime. We discuss this further in
Section 6.

Chételat and Wells/Mid-scale Wishart asymptotics

54

We independently know, by the results of Jiang and Li [2015] and Bubeck
et al. [2016], that a Gaussian orthogonal ensemble approximation holds in the
classical regime. Although ψ0 is not the G-transform of a GOEppq, a simple
Kullback-Leibler argument is sufficient to prove that it approximates ψGOE for
0th degree regimes.
Proposition 4. The total variation distance between the 0th degree G-transform
approximation ψ0 and the Gaussian orthogonal ensemble G-transform ψGOE
satisfies dTV pψ0 , ψGOE q Ñ 0 as n Ñ 8 with p3 {n Ñ 0.
Proof. We use a similar strategy to Theorem 3: namely, by Equation (3.7) is it
equivalent to prove that Hpψ0 , ψGOE q Ñ 0 as n Ñ 8 with p3 {n Ñ 0, and to control that quantity we can use the Kullback-Leibler inequality for G-transforms.
Let T „ |ψGOE | “ GOEppq{4. Since the Gaussian orthogonal ensemble has
been extensively studied, we understand well its empirical moments.
“
‰For example, according to Anderson et al. [2010, Lemma 2.2.2], we have E tr T 2 “ Opp2 q,
while
(2.1.45) of the same book we have Er| tr T |s “ Opp1{2 q and
“ from
‰ Equation
3
3{2
E | tr T | “ Opp q. Then from Definition 3 and Proposition 2, and using the
projection map Pp´π,πq as in the proof of Theorem 3, we find that

ˆ 3˙
„
‰
p`1 “
p
ψGOE pT q
2
“ ´4
E tr T “ O
E < Log
ψ0 pT q
n
n
and
ˇ
„ˇ
„
„ˇ
ˇ
ˇ
ˇ
ˇ
ψGOE pT q ˇˇ
p`1
3
ˇPp´π,πq ´ 32
ˇ
?
?
E ˇˇ= Log
“
E
tr
T
`
tr
T
ˇ
ˇ
ˇ
ψ0 pT q
3 n
n
ˆc 3 ˙
‰
p`1
p
32 “
.
ď ? E | tr T 3 | `2i ? Er| tr T |s “ O
3 n
n
n
ş
Since Sp pRq |ψ0 |pT qdT Ñ 0 when n Ñ 8 with p3 {n Ñ 0 by Equation (5.14), if
we apply Proposition 1 we find thatt
1

1

lim H2 pψ0 , ψGOE q ď 0 ` 0 ` 2 ¨ 1 2 ¨ 0 2 “ 0

nÑ8

for p3 {n Ñ 0. By Equation (3.7), this concludes the proof.
As a consequence, HpfNW , fGOE q “ HpψNW , ψGOE q ď HpψNW , ψ0 q`Hpψ0 , ψGOE q Ñ
0 when n Ñ 8 with p3 {n Ñ 0 by Theorem 3 and Proposition 4. Hence
dTV pfNW , fGOE q Ñ 0 by Equation (2.1) in the classical setting. This provides an
alternative proof of the results of Jiang and Li [2015] and Bubeck et al. [2016].
6. Wishart asymptotics: the density point-of-view
In Section 5,?we studied the asymptotic behavior of the normalized Wishart
distribution nrWp pn, Ip {nq ´ Ip s using its G-transform ψNW . In particular,
we derived an approximation to ψNW for every middle-scale regime of a given
degree. But although it is equivalent to study a probability distribution from a
density or a G-transform point of view, it is still natural to wonder if we can find

Chételat and Wells/Mid-scale Wishart asymptotics

55

approximations to the density of a normalized Wishart for every middle-scale
regime of a given degree.
Recall from Theorem 3 that dTV pψNW , ψK q Ñ 0 when pK`3 {nK`1 Ñ 0.
Define f˜K “ G ´1 tψK u. In general, there is no guarantee that these should be
real-valued. On the other hand, we know from Equation (5.3) that whenever
n ě p ´ 2, ψK must be integrable, and since the G-transform maps integrable
functions to integrable functions, f˜K must also be integrable. In fact, according to Equation (5.14), we know |f˜K | must be asymptotically a density when
pK`1 {nK`1 Ñ 0. This suggests we define the following densities.
Definition 4 (Density approximations). For any K P N and n ě p ´ 2, define
the K th degree density approximation as
|f˜K |pXq
,
fK pXq “ ş
|f˜K pY q|dY
Sp pRq

where f˜K “ G ´1 tψK u and ψK is as in Definition 3. The distribution on the real
symmetric matrices with density fK will be denoted FK .
The main interest is that we can asymptotically approximate the density fNW
of a normalized Wishart by the bona fide densities fK . This was the content
of Theorem 1 from Section 1, which we now prove as a simple corollary of its
G-transform analogue Theorem 3 from Section 5.
Proof of Theorem 1. As in the rest
? of this paper, we write the density of the
normalized Wishart distribution nrWp pn, Ip {nq ´ Ip s by fNW , and by Definition 4 the density of FK is fK . Notice that by Equation (2.1), to prove
dTV pfNW , fK q Ñ 0 it is equivalent to prove that HpfNW , fK q Ñ 0. From the
triangle inequality, the reverse triangle inequality, Theorem 3 and Equation
(5.14),
´
¯
´
¯
´
¯
lim H fNW , fK ď lim H fNW , |f˜K | ` lim H |f˜K |, fK
nÑ8

nÑ8

nÑ8

ż
ˇ
ˇ2
1
ˇ 1{2
ˇ
1{2
ˇfNW pXq ´ |f˜K |pXqˇ dX 2

“ lim

nÑ8

Sp pRq

ż
ˇ
ˇ 1{2
ˇ|f˜ |pXq ´ ş
ˇ K

` lim

nÑ8

Sp pRq

ˇ2
1{2
ˇ
|f˜K |pXq
ˇ dX 12
ˇ
1{2
˜
|fK |pY qdY

Sp pRq

ż
ď lim

nÑ8

ˇ
ˇ2
1
ˇ 1{2
ˇ
1{2
ˇfNW pXq ´ f˜K pXqˇ dX 2
Sp pRq

ż

ˇ
ˇ
` lim ˇˇ1 ´ ş
nÑ8

ˇ
ˇ
1
ˇ
ˇ
1{2
˜
|
f
|pY
qdY
K
Sp pRq

1

|f˜K |pXqdX 2
Sp pRq

Chételat and Wells/Mid-scale Wishart asymptotics

56

ˇ
1 ˇˇ
ˇ
“ 01{2 ` ˇ1 ´ 1{2 ˇ ¨ 11{2 “ 0.
1
when pK`3 {nK`1 Ñ 0. Thus HpfNW , fK q, hence dTV pfNW , fK q, tends to zero
when n Ñ 8 with pK`3 {nK`1 Ñ 0, as desired.
We defined fK in terms of the inverse G-transform of the ψK functions given
by Definition 3. How can we express this explicitely? By Equation (5.3), we see
that |ψK |pT q is asymptotically bounded by the Tn{2 pIp {8q density |ψNW |pT q,
which is integrable whenever n ´ p ` 2 ě 0. But |ψNW |1{2 pT q is proportional
n
Ip q density in the sense of Definition 2, which is integrable for
to a Tm{4 p 4m
m{4 ě p{2 ´ 1, that is whenever n ´ 3p ` 3 ě 0. Thus |ψNW |1{2 pT q and therefore
|ψK |1{2 pT q is integrable whenever n ´ 3p ` 3 ě 0. Hence we can use the Fourer
inversion theorem to conclude that fK is proportional to the integral
ˇ
ˇ
ˇ
ˇ
fK pXq 9 ˇG ´1 tψK uˇpXq
ˇ
ˇ
ˇ
9 ˇ
ˇ

ż

2K`3`

1rK
odds
ÿ
´ 4i ¯k tr T k
n
?
exp i trpXT q `
4
n
k
#

k“2

Sp pRq

2K`2´

1rK
odds
+ ˇ2
ÿ
ˇ
´ 4i ¯k tr T k
p `1
ˇ
?
`
dT ˇ
ˇ
4
n
k

(6.1)

k“1

?
whenever n ´ 3p ` 3 ě 0. In particular, if we do a change of variables Z “ 8T ,
we obtain Equation (1.3) from Section 1 whenever n ě 3p ´ 3, from which we
can derive Equations (1.1) and (1.2).
It would be quite pleasant if there was a way to solve the integral in Equation
(6.1) or (1.3) and obtain a (potentially quite complicated) closed form expression
for fK up to its normalization constant. So far, our efforts have been unfruitful.
We close our discussion with a final remark. At the end of Section 5, we
showed that ψ0 approximates ψGOE in 0th degree middle-scale regimes, from
which the classical asymptotic normality follows. It is natural to wonder if f0
approximates fGOE in the same context. An argument similar to that of Theorem 1 shows this is the case.
Proposition 5. The total variation distance between the 0th degree density approximation f0 and the Gaussian orthogonal ensemble G-transform fGOE satisfies dTV pf0 , fGOE q Ñ 0 as n Ñ 8 with p3 {n Ñ 0.
Proof. The Hellinger distance between f0 and fGOE satisfies
´
¯
´
¯
´
¯
lim H f0 , fGOE ď lim H f0 , |f˜0 | ` lim H |f˜0 |, fGOE
nÑ8

nÑ8

nÑ8

ż
ˇ
ˇ
ď lim ˇˇ1 ´ ş
nÑ8
Sp pRq

ˇ
ˇ
1
ˇ
ˇ
1{2
˜
|f |pY qdY

´
¯
1
|f˜0 |pXqdX 2 ` lim H f˜0 , fGOE
nÑ8

0

Sp pRq

Chételat and Wells/Mid-scale Wishart asymptotics

57

ˇ
´
¯
1 ˇˇ
ˇ
“ ˇ1 ´ 1{2 ˇ ¨ 11{2 ` lim H ψ0 , ψGOE “ 0.
nÑ8
1
By Equation (2.1), the result follows.
Of course, we could conclude from this that HpfNW , fGOE q ď HpfNW , f0 q `
Hpf0 , fGOE q Ñ 0 when n Ñ 8 with p3 {n Ñ 0, offering yet again another
proof that a Gaussian orthogonal ensemble approximation holds in the classical
setting.
7. The effect of phase transitions
Although we have established the existence of phase transitions, it does not shed
much light on how the behavior of a normalized Wishart distribution might differ
across phase transitions. To do this, it can be very illuminating to study the
asymptotics of some of its statistics. For example, we could study its empirical
moments.
?
For a normalized Wishart matrix X „ nrWp pn, Ip {8q ´ Ip s, a direct computation yields
„ˆ
˙ 
ˆ ˙
2
1 ´ X ¯2
5
4
8
20
20
5
1
E
tr ?
´1
“ 2` 3`
` 2 ` 3 “ 2 `o
p
p
p
p
np np
np
p
p2
so in every middle-scale regime, that is whenever n, p Ñ 8 with p{n Ñ 0,
?
›
›
´
¯
›
› 1 ´ X ¯2
5
all middle-scale
›
› tr ?
„
´
1
.
› 2
›p
regimes
p
p
L
Thus we have L2 convergence of the second empirical moment to 1, but otherwise nothing very interesting. There doesn’t seem to be any change of behavior
across the different middle-scale regimes. In contrast, the situation with the
symmetric t distribution is striking, and illustrates yet again that middle-scale
regime behavior becomes clearer under a G-transform. Indeed, we know from
4T 2
Theorem 2 that for a T „ Tn{2 pIp {8q, the quantity p1 trp ?
p q also converges to
1, but we know more. At Equation (4.41), we computed the exact L2 distance
4T 2
between p1 trp ?
p q and 1, and found that
˙ 
„ˆ
2
m5
1 ´ 4T ¯2
E
tr ?
´1
“
p
p
pm´6q pm´2q pm´1q pm`1q pm`3q
˙
„
ˆ
2
p2
2
p2
5
1
¨ 2`
` 2 `o
`
`
.
p
m m
p2
m m2
Thus the L2 distance must have middle-scale asymptotics
›
›
› 1 ´ 4T ¯2
›
› tr ?
´ 1››
›p
p
L2

Chételat and Wells/Mid-scale Wishart asymptotics

?
$
5
’
’
’
’
’
p
’
’
’
’
’
’
?
1
’
2
’
’
& 5 ` 2α ` α p
?
„
1
’
“ 5α´1 ` 2 ` α ?
’
’
n
’
’
?
’
’
´2
´1
’ “ 5α ` 2α ` 1 p
’
’
n
’
’
’
p
’
%
n

58

for

p2
Ñ0
n

for

p2
Ñα
n

for

¯
´
p2
second or higher
Ñ8
degree
n

´

`

classical or first
degree

second degree

¯

˘

(7.1)
4T 2
as n, p Ñ 8 with p{n Ñ 0. Thus there is a sharp change in behavior of p1 trp ?
pq
?
when p grows like n, and despite the symmetric t distribution satisfying a
semicircle law according to Corollary 1, it must ultimately behave differently
than a Gaussian orthogonal ensemble matrix. The first-order asymptotics look
the same: it is rather in the rate of this convergence that they differ.
This matters for both the symmetric t and the Wishart distribution because
rates of convergence can be distinguished in the strong topology. As a simple
example, consider the sequence of one-dimensional distributions
`
˘
`
˘
Fp “ N 0, 1{p ,
and Gp “ N 0, 1{p2 .

In the weak topology, these are asymptotically the same, since they converge to
the same distribution – namely Fp , Gp ñ δ0 as p Ñ 8, for δ0 the Dirac measure
at 0. In other words, in a metric that induces the weak topology such as the
Lévy metric,
dLévy pFp , Gp q Ñ 0.
Yet, by a direct computation of the Hellinger distance, which induces the strong
topology,
d
´
¯1{4
?
4p
dHellinger pFp , Gp q “ HpFp , Gp q “ 2 1 ´ 2
p ` 2p ` 1
?
Ñ 2ą0
as p Ñ 8. Thus it is clear that the strong topology captures rates of convergence in a way that the weak topology
can’t. But then, we should expect
?
a phase transition when p grows like n for the Tn{2 pIp {8q distribution. And
since the symmetric t is the G-conjugate
of the Wishart, this should imply a
?
phase transition when p grows like n for the Wishart distribution as well. This
is consistent with Theorem 3, and provides an alternative explanation for the
existence of the second phase transition.
A natural question then is to ask whether we can find symmetric t statistics
that exemplify all the middle-scale regime phase transitions. It is tempting to
look at the L2 error of the other empirical moments of the symmetric t distribution, because we can use the methodology developed in Section 4 to compute
their asymptotics to arbitrary order. As a reference, we compiled a table of the

Chételat and Wells/Mid-scale Wishart asymptotics

Normalized empirical moment
1 ´ 4T ¯
tr ?
p
p
1 ´ 4T ¯2
tr ?
p
p
1 ´ 4T ¯3
tr ?
p
p
1 ´ 4T ¯4
tr ?
p
p

Limit

Asymptotics of its squared L2 error

0

2
p2

C1 “ 1

5
2
p2
`
` 2
2
p
m
m

0

24
p2

C2 “ 2

50
25p2
97
`
`
2
p
m
m2

59

Table 1
Asymptotics of small normalized empirical moments of T „ Tn{2 pIp {8q.

few first few moments as Table 1.
As can be seen from the table, the odd moments seem to have uniform behavior across all middle-scale regimes. In contrast, the even ?
moments seem to
all change asymptotics at the second phase transition p “ Θp nq, but nowhere
else. Hence finding statistics that “flag” the other phase transitions remain an
open question.
8. Auxiliary results
This section compiles several lemmas used elsewhere in the article.
Lemma 3 (First derivatives lemma).For any indices 1 ď i1 , . . . , i2l ď p and real
symmetric matrix Z, there exist polynomials aJ,s pn, mq in n and m “ n ´ p ´ 1,
indexed by 0 ď s ď l and J “ pj1 , . . . , j2l q, such that
"
*
ˇ ˇm
Bs
Bs
Bs
n
...
exp ´ tr Z ˇZ ˇ 4
Bs Zi2l i2l´1
Bs Zi4 i3 Bs Zi2 i1
4
"
*
s
l
l
ź
ÿ
ÿ
ź
ˇ ˇm
n
´1
Zj2t j2t´1 exp ´ tr Z ˇZ ˇ 4 .
“
aJ,s pn, mq
pIp qj2t j2t´1
4
t“1
s“0
t“s`1
JP
t1,...,pu2l

Proof. To simplify notation, let
*
"
l
s
ź
ź
ˇ ˇm
n
´1
MJ,s pZq “
pIp qj2t j2t´1
Zj2t j2t´1 exp ´ tr Z ˇZ ˇ 4 ,
4
t“s`1
t“1
and let Ml “ tMJ,s | J P t1, . . . , pu2l , s ď lu be the set of all such terms “on 2l
indices”. Let xMl y denote the linear span of Ml , that is, the space of all linear
combinations of elements of Ml , with as coefficients real polynomials in n and
m. Then we are really claiming that
"
*
ˇ ˇm
Bs
n
Bs
Bs
...
exp ´ tr Z ˇZ ˇ 4
P xMl y.
(8.1)
Bs Zi2l i2l´1
Bs Zi4 i3 Bs Zi2 i1
4

Chételat and Wells/Mid-scale Wishart asymptotics

60

To see this, let J “ pj1 , . . . , j2l´2 q P t1, . . . , pu2l´2 and define the extension
q
Ja,b
“ pj1 , . . . , jq´1 , a, b, jq`1 , . . . , j2l´2 q P t1, . . . , pu2l to be J with indices a, b
inserted (in this order) at the q th position. Then using that
ı
1 ” ´1 ´1
Bs
´1
´1
Zab
“ ´ Zai
Zi2l´1 b ` Zai
Zi´1
2l
2l´1
2l b
Bs Zi2l i2l´1
2
and
"
*
ˇ ˇm
n
exp ´ tr Z ˇZ ˇ 4
Bs Zi2l i2l´1
4
"
*
”m
ı
ˇ ˇm
n
n
“
Zi2l i2l´1 ´ pIp qi2l i2l´1 exp ´ tr Z ˇZ ˇ 4 ,
4
4
4
we conclude that
s
s
Bs
1 ÿ
1 ÿ
MJi2r i
´
M 2r
MJ,s pZq “ ´
,s`1
2l 2l´1
Bs Zi2l i2l´1
2 r“1
2 r“1 Ji2l´1 i2l ,s`1
m
n
` MJ 2s`1 ,s`1 ´ MJ 2s`1 ,s
P xMl y.
i
i
i2l i2l´1
4
4
2l 2l´1
Thus, by linearity, Bs {Bs Zi2l i2l´1 maps xMl´1 y to xMl y. But naturally we have
expt´ n4 tr Zu|Z|m{4 P xM0 y, so by induction Equation (8.1) must then hold, as
desired.
Bs

Lemma 4 (Second derivatives lemma).For any k P N and any Z P Sp pRq,
p
ÿ

i1 ,...,i2k

ˇ ˇm
n
Bs
Bs
Bs
...
e´ 4 tr Z ˇZ ˇ 4
Bs Zi1 i2k
B s Z i3 i2 B s Z i2 i1
n

“ e´ 4

ˇ ˇ m4 ÿ p1q
Zˇ
bκ pn, m, pqrκ pZ ´1 q

tr Z ˇ

|κ|ď2k

and
p
ÿ
i1 ,...,ik
j1 ,...,jk

ˇ ˇm
n
Bs
Bs
Bs
Bs
Bs
Bs
e´ 4 tr Z ˇZ ˇ 4
...
...
Bs Zj1 jk
Bs Zj3 j2 Bs Zj2 j1 Bs Zi1 i2k
Bs Zi3 i2 Bs Zi2 i1
n

“ e´ 4

ˇ ˇ m4 ÿ p2q
Zˇ
bκ pn, m, pqrκ pZ ´1 q

tr Z ˇ

|κ|ď2k`1
p1q
bκ pn, m, pq

p2q
bκ pn, m, pq

p1q

with degrees deg bκ ď 2k `
and
for some polynomials
p2q
1 ´ qpκq and deg bκ ď 2k ` 2 ´ qpκq. The sums at the right hand sides are taken
over all integer partitions κ of norm at most 2k and 2k ` 1, including the empty
partition.
Proof. We give a spectral proof. Let OLOt be the spectral decomposition of Z,
with eigenvalues λ1 ě ¨ ¨ ¨ ě λp , and notice that
p
t
t
Bs Ohl
1 ÿ Oha Oai
1 ÿ Oha Oaj
“
Ojl `
Oil ,
Bs Zij
2 a‰l λl ´ λa
2 a‰l λl ´ λa

Bs λ h
“ Oih Ojh
Bs Zij

for any 1 ď i, j, h, l ď p. As a consequence, for any differentiable real-valued

Chételat and Wells/Mid-scale Wishart asymptotics

61

functions F1 pLq, . . . , Fp pLq, we have
ˆÿ
˙
p
p
p
p
ÿ
ÿ
Fb ´ Fa t
Bs
1 ÿ
BFa t
t
Oha
Oja Fa Oai
“
Oai `
Oha
O .
B
Z
2
λ
´
λ
Bλa ai
s
hj
b
a
a“1
a“1
j“1
a,b
b‰a

This suggests we define a new operator DL that would map the space of diagonal
matrices F pLq “ diagpF1 pLq, . . . , Fp pLqq that differentially depends on L, to
itself, by
p
p
ÿ
BFa
1 ÿ Fb ´ Fa
Bs
t
t
`
so that
OF Oji
“ OF Oki
.
DL tF ua “
2 b‰a λb ´ λa
Bλa
B
Z
j“1 s hj
In particular,
p
ÿ
i1 ,...,i2k

ˇ ˇm
n
Bs
Bs
Bs
...
e´ 4 tr Z ˇZ ˇ 4
Bs Zi1 i2k
Bs Zi3 i2 Bs Zi2 i1
p
ÿ

“
i1 ,...,i2k

Bs ” ´ n tr Z ˇˇ ˇˇ m4 ı
Bs
Bs
Z Ip
e 4
...
Bs Zi0 i2k
Bs Zi3 i2 Bs Zi2 i1
i1 i0

! n
ˇ ˇm )
2k
“ tr DL
e´ 4 tr Z ˇZ ˇ 4 Ip ,

(8.2)

and similarly
p
ÿ
i1 ,...,ik
j1 ,...,jk

ˇ ˇm
n
Bs
Bs
Bs
Bs
Bs
Bs
e´ 4 tr Z ˇZ ˇ 4
...
...
Bs Zj1 jk
Bs Zj3 j2 Bs Zj2 j1 Bs Zi1 ik
B s Z i3 i2 B s Z i2 i1

!
! n
ˇ ˇ m ))
k
k
“ tr DL
tr DL
e´ 4 tr Z ˇZ ˇ 4 Ip .

(8.3)

Let us look more closely at this operator DL . It satisfies the following.
(i) DL is linear, in the sense that for diagonals F pLq, GpLq and constants a,
b with respect to L,
DL taF ` bGu “ aDL tF u ` bDL tGu.
(ii) DL satisfies a restricted product rule, in the sense that for a diagonal
F pLq of the form F pLq “ f pLqIp for some function f pLq, and any diagonal
GpLq,
DL tF Gu “ DL tF uG ` F DL tGu.
Moreover, from the definition of DL ,
(
n
n n
DL e´ 4 tr L Ip “ ´ e´ 4 tr L Ip ,
4
(
DL trpL´s qIp “ ´sL´ps`1q ,
and

( m m
m
DL |L| 4 Ip “ |L| 4 Ip ,
4

s
(
s
1ÿ
trpL´rs`1´ts qL´t .
DL L´s “ ´ L´ps`1q ´
2
2 t“1

Chételat and Wells/Mid-scale Wishart asymptotics

62

Now define the spaces
#
Ml “

´n
4

bpn, m, pqe

ˇ
+
ˇ bpn, m, pq is a polynomial with
m
ˇ
ˇ
ˇ
tr L ˇ ˇ 4
L rκ pL´1 qL´s ˇ degree at most l ´ qpκq, and κ
ˇ and s satisfy |κ| ď l ´ s.

for l “ 1, . . . , 2k, and let xMl y denote the linear span of Ml , i.e. the space of
all real linear combinations of elements of Ml . Moreover, for a partition κ, let
κ ˘ i denote κ with the integer i added or removed, respectively. For example,
p3, 1, 1, 1q`2 “ p3, 2, 1, 1, 1q and p3, 2, 1, 1, 1q´1 “ p3, 2, 1, 1q. Note that |κ˘i| “
|κ| ˘ i. Then, for any F P Ml ,
)
!
ˇ ˇm
n
DL tF u “ DL bpn, m, pqe´ 4 tr L ˇLˇ 4 rκ pL´1 qL´s
! n
)ˇ ˇ m
“ bpn, m, pqDL e´ 4 tr L Ip ˇLˇ 4 rκ pL´1 qL´s
!ˇ ˇ m )
n
` bpn, m, pqe´ 4 tr L DL ˇLˇ 4 Ip rκ pL´1 qL´s
!
)
ˇ ˇm
n
` bpn, m, pqe´ 4 tr L ˇLˇ 4 DL rκ pL´1 qIp L´s
!
)
ˇ ˇm
n
` bpn, m, pqe´ 4 tr L ˇLˇ 4 rκ pL´1 qDL L´s
ı n
” n
ˇ ˇm
“ ´ bpn, m, pq e´ 4 tr L ˇLˇ 4 rκ pL´1 qL´s
” m4
ı n
ˇ ˇm
`
bpn, m, pq e´ 4 tr L ˇLˇ 4 rκ pL´1 qL´s
4
qpκq
ı n
ÿ”
ˇ ˇm
`
´ κi bpn, m, pq e´ 4 tr L ˇLˇ 4 rκ´κi pL´1 qL´s
i“1

ı n
ˇ ˇm
s
` ´ bpn, m, pq e´ 4 tr L ˇLˇ 4 rκ pL´1 qL´ps`1q
2
s ”
ı n
ÿ
ˇ ˇm
1
`
´ bpn, m, pq e´ 4 tr L ˇLˇ 4 rκ`ps`1´tq pL´1 qL´t
2
t“1
”

Thus DL tF u P xMl`1 y. It follows by linearity that DL maps xMl y to xMl`1 y.
ˇ ˇm
ˇ ˇm
n
2k ´ n
te 4 tr L ˇLˇ 4 Ip u P xM2k y.
Now, e´ 4 tr L ˇLˇ 4 Ip P M0 , so by induction DL
p1q

Hence, for some polynomials bκ,s pn, m, pq of degree at most 2k ´ qpκq,
! n
ÿ
ˇ ˇm )
ˇ ˇm
2k
´n
4 tr L ˇLˇ 4 r pL´1 q trpL´s q
tr DL
e´ 4 tr L ˇLˇ 4 Ip “
bp1q
κ
κ,s pn, m, pqe
|κ|`sď2k

ÿ
“

p1q

n

bκ1 pn, m, pqe´ 4

ˇ ˇ m4
Lˇ rκ1 pL´1 q

tr L ˇ

(8.4)

|κ1 |ď2k
p1q

p1q

p1q

p1q

for κ1 “ κ ` s, bκ1 “ bκ,s when s ‰ 0, while κ1 “ κ, bκ1 “ pbκ,s when
s “ 0. Notice that when s ‰ 0, the degree of the bκ1 ’s is at most 2k ´ qpκq “
2k ´ pqpκ1 q ´ 1q, while when s “ 0 it is at most 2k ´ qpκq ` 1 “ 2k ´ qpκ1 q ` 1.
p1q
Thus in both cases, deg bκ1 ď 2k ´ qpκ1 q ` 1, which by Equation (8.2) shows
the first statement of the lemma.
For the second statement of the lemma, by an argument analoguous to Equa-

Chételat and Wells/Mid-scale Wishart asymptotics

63

ˇ ˇ m4
Lˇ Ip u P xMk`1 y. Thus by induction
ˇ ˇm
n
k
k
again, we must have DL
ttr DL
te´ 4 tr L ˇLˇ 4 Ip uu P xM2k`1 y. Hence for some
n

k
tion (8.4) we find that tr DL
te´ 4

tr L ˇ

p2q

polynomials bκ,s pn, m, pq of degree at most 2k ` 1 ´ qpκq,
!
! n
ˇ ˇ m ))
k
k
DL
tr DL
e´ 4 tr L ˇLˇ 4 Ip
ÿ
ˇ ˇm
´n
4 tr L ˇLˇ 4 r pL´1 q trpL´s q
“
bp2q
κ
κ,s pn, m, pqe
|κ|`sď2k`1

ÿ
“

n

p2q

bκ1 pn, m, pqe´ 4

ˇ ˇ m4
Lˇ rκ1 pL´1 q

tr L ˇ

|κ1 |ď2k`1
p2q

p2q

p2q

p2q

for again κ1 “ κ ` s, bκ1 “ bκ,s when s ‰ 0, while κ1 “ κ, bκ1 “ pbκ,s when
p2q
s “ 0. By the same argument as before, deg bκ1 ď 2k ´ qpκ1 q ` 2, which by
Equation (8.3) shows the second statement of the lemma. This concludes the
proof.
We will also need in our proof a result about the asymptotics of inverse
moments of the Wishart distribution. Because we couldn’t find anything like it
in the literature, we think it is worthwhile to provide some context.
Let f : p0, 4q Ñ R be the restriction to the positive reals of a complex
function analytic in a neighborhood of p0, 4q. We are often interested in the
linear spectral statistic p1 tr f pY q for Y „ Wp pn, Ip {nq. Much is known about
its distributional properties in the high-dimensional regime where p Ñ 8 such
that lim np “ α ă 1. For example, if 0 ă α ă 1, there must be an  ą 0 such
nÑ8

that p{n P r, 1 ´ s for all n large enough, so Bai and Silverstein [2010, Theorem
9.10] and the dominated convergence theorem yield that

ż

rp{ns`

a
prp{ns` ´ tqpt ´ rp{ns´ q
f ptq
dt
2πrp{nst

1
P
tr f pY q ÝÑ lim
nÑ8
p

ż α rp{ns

´

`

“

a
pα` ´ tqpt ´ α´ q
dt
f ptq
2παt

(8.5)

α´

?
P
as n Ñ 8. Here, Ñ stands for convergence in probability and x˘ for p1 ˘ xq2 .
In fact, the theorem states more, namely a central limit theorem, but what we
want to draw to attention is the class of functions for which this result was
proven.
This is sometimes enough, but often we would like to understand the expectation of this linear spectral statistic. If f is bounded, then Equation (8.5)
implies that
α`
a
ı
pα` ´ tqpt ´ α´ q
1 ”
f ptq
dt.
(8.6)
lim E tr f pY q “
nÑ8 p
2παt

ż

α´

This is nice for a function f pzq like ez or sin z that happens to be bounded

Chételat and Wells/Mid-scale Wishart asymptotics

64

on a neighborhood of p0, 4q, but it unfortunately excludes many interesting
unbounded functions, such as log z or 1{z. In fact, for unbounded f , it is in
general not even clear if lim p1 Ertr f pY qs will be finite!
nÑ8

The following result shows that, at least in the case f pzq “ 1{z s for s P N,
we can use Stein’s lemma to obtain Equation (8.6) and its α “ 0 analogue.
Lemma 5. Let for Y „ Wp pn, Ip {nq and s be any integer s ě 1. Then as long
as n ě p ` 4s ` 2, the sth inverse moment satisfies the recursive bound
ˆ
˙ ”
ı
”
ı
pp ` 1qs
1´
E tr Y ´s ď E tr Y ´ps´1q .
n
In particular, as p Ñ 8 such that lim np “ α ă 1, if s ă α´1 ´ 1 then
nÑ8
$ α`
a
’
’
pα` ´ tqpt ´ α´ q
&
ı ’
1 ”
dt
if 0 ă α ă 1,
´s
lim E tr Y
“
2παts`1
nÑ8 p
’
α
´
’
’
%
1
if α “ 0.
? 2
for α˘ “ p1 ˘ αq .

ż

Proof. The classical
lemma
that for any differentiable function f :
ˇstates
“ˇ`Stein’s
˘
‰
2
B
R Ñ R such that E ˇ BZ
´ Z f pZqˇ ă 8 for Z „ Np0, 1q and lim f pzqe´z {2 “
zÑ˘8

0, we must have

„´
¯
B
´ Z f pZq “ 0.
E
BZ
Let Z „ Nnˆp p0, In b Ip q be an n ˆ p matrix of i.i.d. standard normal random
variables, and let Y “ n1 Z t Z „ Wp pn, Ip {nq. For any 1 ď α ď n and 1 ď
β, i, j ď p,
´s
p
s
ı
Bs Yβj
B
2 ÿ
1 ÿ ” ´l ´ps´l`1q
Bs
´l ´ps´l`1q
Yβi Yβj
` Yββ
“
Zαi
and
“´
Yij
,
BZαβ n i“1
Bs Yiβ
Bs Yiβ
2 l“1
so for δ the Kronecker delta,
˙
ˆ
`
˘
B
´ Zαβ ZY ´s αβ
BZαβ

p „
p
ÿ
2 ÿ
Bs
´s
´s
´s
“
δβj Yβj `
Zαj Zαi
Y ´ Zαβ Zαj Yβj
n i“1
Bs Yiβ βj
j“1
´s
“ Yββ
´

s
s
˘ `
˘
˘
1 ÿ`
1 ÿ ´l `
ZY l αβ ZY ´ps´l`1q αβ ´
Yββ ZY ´ps´l`1q Z t αα
n l“1
n l“1
`
˘
´Zαβ ZY ´s αβ
(8.7)

Let us first show that this expression is integrable. For any matrix X, |Xij | ď
1{2
}X}2 “ }X t X}2 . Thus by Equation (8.7),
ˇ
„ˇˆ
˙
ˇ
`
˘ ˇ
B
´s
ˇ
ˇ
E ˇ
´ Zαβ ZY
αβ ˇ
BZαβ

Chételat and Wells/Mid-scale Wishart asymptotics

«
ď E }Y ´s }2 `

s
ÿ

1

65

1

}Y 2l`1 }22 }Y 2s`2l´1 }22

l“1

ff

s
ÿ

}Y

`

´2l

1
2

}2 }Y

´2s`2l

1
2

1
2

}2 ` n}Y }2 }Y

´2s`1

1
2

}2

l“1

As Y is positive definite, }Y ˘a }2 ď tr Y ˘a for any a P N, so by the CauchySchwarz inequality,
„
 ÿ
s
”
ı 12 ”
ı 12
E tr Y ´2l`1 E tr Y ´2s`2l´1
ď E tr Y ´s `
l“1
s
”
ı 12 ”
ı 12
”
ı 12 ”
ı 12
ÿ
E tr Y ´2l E tr Y ´2s`2l ` n E tr Y E tr Y ´2s`1 ,
`
l“1

which is finite for n ě p ` 4s ` 2.
Moreover, pZY ´s qαβ can be expressed using minors and determinants as a
rational function of the entries of Z, so
lim

Zαβ Ñ˘8

2

pZY ´s qαβ e´Zαβ {2 “ 0.

So all conditions are fulfilled to apply Stein’s lemma to Equation 8.7 and obtain
ff
«
˙
p ˆ
n
`
˘
B
1 ÿ ÿ
´s
´ Zαβ ZY
0“E
αβ
n α“1 β“1 BZαβ
«
ff
s
s
1 ÿ
´s
´s
´l
´ps´lq
´ps´1q
“ E tr Y ´ tr Y ´
trpY q trpY
q ´ tr Y
n
n l“1
As trpY ´l q trpY ´ps´lq q ď p tr Y ´s for any 1 ď l ď s, and every term is integrable
as n ě p ` 4s ` 2, this means that
ˆ
˙ ”
ı
”
ı
pp ` 1qs
1´
E tr Y ´s ď E tr Y ´ps´1q .
(8.8)
n
This shows the first part of the proof.
For the second part, let S P N. If we let n Ñ 8 such that lim np “ α ă 1,
nÑ8

then any S ă α´1 we will have n ě p ` 4S ` 2 and n ě pp ` 1qS for n large
enough. So by repeatedly applying Equation (8.8) for s “ S, . . . , 1 and dividing
by p, we obtain
˙
S ˆ
ı
ı
ź
pp ` 1ql
1 ”
1 ”
1´
¨ E tr Y ´S ď
E tr Y ´0 “ 1.
n
p
p
l“1
Taking a limit in the above yields
S
ı
ź
1 ”
p1 ´ αlq lim E tr Y ´S ď 1.
nÑ8 p
l“1
Thus for any S ă α´1 , we have
S
ı
ź
1
1 ”
lim E tr Y ´S ď
ă 8.
nÑ8 p
1
´
αl
l“1

(8.9)

Chételat and Wells/Mid-scale Wishart asymptotics

66

In the case 0 ă α ă 1, if s`1 ă α´1 then by Jensen’s inequality and Equation
(8.9) applied to S “ s ` 1, we have
ı
˘1` 1s ı
1 ”
1 ”´
ď lim E tr Y ´ps`1q ă 8.
lim E tr Y ´s
nÑ8 p
nÑ8 p
1
´s
Thus p tr Y
is uniformly integrable, and by Equation (8.5),

żα

`
a
”
ı
pα` ´ tqpt ´ α´ q
1
dt
lim E tr T ´s “
nÑ8 p
2παts`1
α´
?
for α˘ “ p1 ˘ αq2 .
In contrast, by applying Jensen’s inequality twice,
„´
„
´s
„

¯´s 
1
1
1
ďE
1 “ E tr Y
tr Y
ď E tr Y ´s
p
p
p

so when α “ 0, by applying Equation (8.9) with S “ s, we obtain that
lim p1 Ertr Y ´s s “ 1, as desired.
nÑ8

9. Conclusion
The results of this paper raise more questions than they answer. We enumerate
some that we found particularly interesting.
(1) The univariate t distribution
with ν degrees of freedom is often defined as
?
the distribution of Z{ s, for Z „ Np0, 1q and s „ χ2ν {ν independent. In
the real symmetric matrix case, we could imagine studying the distribution
of S 1{4 ZS 1{4 , for Z „ GOEppq and S „ Wp pν, Ip {νq independent. Is this
the Tν p2Ip q distribution in the sense of Section 4?
(2) By Theorem 2 and Corollary 1, it is clear the empirical moments of a
symmetric t distribution are quite similar to those of a Gaussian orthogonal ensemble matrix, except perhaps in their rates of convergence. From
Anderson et al. [2010, Theorem 2.1.31], we know the empirical moments
of a Gaussian orthogonal ensemble are asymptotically normal. Are the
empirical moments of the symmetric t distribution also asymptotically
normal?
(3) In Section 4, we showed that the rate of convergence of the even normalized
empirical
moments of a symmetric t distribution change when p grows like
?
n. Can we find analogue symmetric t statistics that change their rates of
convergence when p grows like npK`1q{pK`3q for every K P N? This would
establish phase transitions for the symmetric t distribution. If so, can we
find approximating densities between every two transitions, just like in the
Wishart case?
(4) As a counterpart of Theorem 1, could we prove that dTV pfNW , fK q Û 0
whenever pK`3 {nK`1 Û 0 as n Ñ 8? This is delicate because we have no
guarantee that the L1 norm of ψK is asymptotically bounded for regimes
of degree K ` 1 or higher.

Chételat and Wells/Mid-scale Wishart asymptotics

67

(5) Can we find the normalization constant or, better, solve the expectation
of Equation (1.3) in closed form?
(6) What asymptotics hold for the symmetric t or the Wishart distribution in
a middle-scale regime of infinite degree? How do these asymptotics differ
from the other middle-scale regimes, or the high-dimensional regime?
(7) The symmetric t distribution was discovered as the G-conjugate of the
Wishart distribution. What other distributions can be realized as the Gconjugate of some well-known distribution?
(8) In Lemma 2, we expressed the characteristic function of the G-conjugate
F ˚ of a distribution F as f 1{2 ‹ pf 1{2 ˝ Rq, for f the density of F and R the
flip operator. To obtain the moments, we then repeatedly differentiated
under the convolution integral at zero, and obtained an expression of the
moments as an expectation with respect to f . The argument worked when
F ˚ was the symmetric t distribution. Can this argument be generalized to
other F ˚ ? If F ˚ is a well-known distribution, does this give rise to novel
and nontrivial expressions for its moments?
(9) The G-transform of a distribution encodes all the information relative to
that distribution. However, taking a modulus removes some information,
and so in some sense the G-conjugate distribution is “less informative”
than the original distribution. What happens when we repeatedly apply
the G-conjugation operator, destroying information every time? For example, is there an attractor distribution G that is the limit of this process
regardless of the initial distribution?
(10) Can we find distinct random operators which can be regarded, in some
sense, as the total variation limit of a normalized Wishart matrix between
every two phase transitions?
It appears to us that some of these questions might be very difficult to answer.
We would be pleased if future work were able to shed light on any of them.
References
Greg W. Anderson, Alice Guionnet, and Ofer Zeitouni. An Introduction to
Random Matrices. Cambridge University Press, 2010.
Zhidong Bai and Jack W. Silverstein. Spectral Analysis of Large Dimensional
Random Matrices, volume 20. Springer, 2010.
Maurice S. Bartlett. On the theory of statistical regression. Proceedings of the
Royal Society of Edinburgh, 53:260–283, 1933.
Sébastien Bubeck and Shirshendu Ganguly. Entropic CLT and phase transition
in high-dimensional Wishart matrices. International Mathematics Research
Notices, pages 243–258, 2016.
Sébastien Bubeck, Jian Ding, Ronen Eldan, and Miklós Z. Rácz. Testing for
high-dimensional geometry in random graphs. Random Structures & Algorithms, 49:503–532, 2016.
Kai-Lai Chung. A Course in Probability Theory. Academic Press, 2001.

Chételat and Wells/Mid-scale Wishart asymptotics

68

Michel Gaudin. Sur la loi limite de l’espacement des valeurs propres d’une
matrice alé atoire. Nuclear Physics, 25:447–458, 1961.
Arjun K. Gupta and Daya K. Nagar. Matrix Variate Distributions. CRC Press,
1999.
Tiefeng Jiang and Danning Li. Approximation of rectangular beta-Laguerre
ensembles and large deviations. Journal of Theoretical Probability, 28:804–
847, 2015.
Gérard Letac and Hélène Massam. All invariant moments of the wishart distribution. Scandinavian Journal of Statistics, 31:295–318, 2004.
Vladimir A. Marchenko and Leonid A. Pastur. Distribution of eigenvalues for
some sets of random matrices. Matematicheskii Sbornik, 114:507–536, 1967.
Sho Matsumoto. General moments of the inverse real wishart distribution and
orthogonal weingarten functions. Journal of Theoretical Probability, 25:798–
822, 2012.
Madan L. Mehta. On the statistical properties of the level-spacings in nuclear
spectra. Nuclear Physics, 18:395–419, 1960a.
Madan L. Mehta. Lemploi des polynômes orthogonaux pour calculer certain
déterminants. Rapport S PH (Saclay), 658, 1960b.
Robb J. Muirhead. Aspects of Multivariate Statistical Theory. Wiley Series in
Probability and Mathematical Statistics. Wiley, 1982.
Charles E. Porter and Norbert Rosenzweig. Statistical properties of atomic and
nuclear spectra. Technical report, Univ. of Minnesota, Minneapolis, 1960.
Mohsen Pourahmadi. High-dimensional Covariance Estimation: with HighDimensional Data. John Wiley & Sons, 2013.
Miklós Z. Rácz and Jacob Richey. A smooth transition from Wishart to GOE.
arXiv:1611.05838, 2016.
Aad W. van der Vaart. Asymptotic statistics, volume 3. Cambridge University
Press, 2000.
John Von Neumann and Herman H. Goldstine. Numerical inverting of matrices
of high order. Bulletin of the American Mathematical Society, 53:1021–1099,
1947.
Eugene P. Wigner. On the statistical distribution of the widths and spacings
of nuclear resonance levels. In Mathematical Proceedings of the Cambridge
Philosophical Society, volume 47, pages 790–798, 1951.
Eugene P. Wigner. Characteristic vectors of bordered matrices with infinite
dimensions. Annals of Mathematics, 62:548–564, 1955.
Eugene P. Wigner. Statistical properties of real symmetric matrices with many
dimensions. In Proceedings of the Fourth Canadian Mathematical Congress.
University of Toronto Press, 1957.
John Wishart. The generalised product moment distribution in samples from a
normal multivariate population. Biometrika, 20A:32–52, 1928.

