Cauchy Identities for the Characters of the Compact Classical
Groups
Amir SepehriI

arXiv:1610.04659v2 [math.RT] 17 Nov 2016

Department of Statistics, Sequoia Hall, Stanford, CA 94305, USA

Abstract
Motivated by statistical applications, this paper introduces Cauchy identities for characters of the
compact classical groups. These identities generalize the well-known Cauchy identity for characters of the unitary group, which are Schur functions of symmetric function theory. Application
to statistical hypothesis testing is briefly sketched.
Keywords: Cauchy identity, Character theory, Classical groups, Hypothesis testing
2000 MSC: 43A75, 43A45, 62G10

1. Introduction
Consider generating uniformly random three dimensional rotations, i.e. elements of SO(3).
The standard way of generating a random three dimensional rotation is to rotate uniformly along
the north pole, and then move the north pole to a uniform point on the sphere. This is a special
case of the famous subgroup algorithm of Diaconis and Shahshahani (1987). However, a very
common and easy to make mistake is to perform this algorithm in the reverse order; that is,
first move the north pole to a random point on the sphere, and then rotate uniformly along that
point. This can be rephrased in terms of axis-angle representation of three dimensional rotations.
A well known fact, due to Euler, asserts that any three dimensional rotation is a pure rotation
along a single fixed axis. It is common to think that the uniform distribution of the rotation
induces uniform distributions on the axis and the angle, independently. However, this leads to
the naive algorithm above. For instance, a paper with sampling as its main concern mentions this
as an ‘intuitively correct’ way of sampling uniformly from SO(3); even though they use another
correct algorithm in a different parameterization of rotations (Kuffner, 2004). Figure 1 below
shows samples generated using the two algorithms above. As it can be seen, it is not easy to tell
the difference between the two only by looking at this plot. However, a more careful inspection
of Figure 1 may suggest that the naive sampling has more rotations with a small angle, that is,
rotations with the red and black arrows ‘close’ to each other. For a more formal inspection,
consider only the rotation angles, as represented in Figure 2. Figure 2 illustrates the scatter-plot
of the angles for the two samplers as well as a histogram for each. It is clear that the underlying
I The author is supported by a Weiland Graduate Fellowship in the School of Humanities and Sciences, Stanford
University.
Email address: asepehri@stanford.edu (Amir SepehriI )
Preprint submitted to Elsevier
November 18, 2016

(a) Uniform sampling

(b) Naive sampling

Figure 1: Samples from the uniform distribution (a) and naive sampling (b). Each rotation is
represented by its axis and angle. The point on the sphere represents the axis, the angle between
the red and black arrows show the angle (black arrow corresponds to the rotation with angle zero,
i.e. the identity). Each plot is based on N = 1500 samples.
distributions are different. The scatter-plot shows that the angle for uniform sampler has a higher
density around θ = π and a lower density around θ = 0. This makes intuitive sense because
all the rotations with small angles are similar to the identity regardless of the axis. Therefore,
a uniform angle distribution results in an oversampling of rotations closer to the identity. This
is the case for the naive sampler. For a more rigorous assessment of the differences between
the two samplers, formal statistical hypothesis tests can be used. If only the angle of rotation is
of concern, the problem becomes a one dimensional goodness of fit testing problem. This is a
classical problem with an enormous literature around it, including several volumes of textbooks.
In many problems, it may not be true that the angle distribution is different between two samplers;
even if it is, it might be hard to guess that by looking at the data. A more natural problem is
to test if the full distributions on SO(3) are equal. The problem of goodness of fit for rotation
data has been studied extensively, and there are various tests available. Two commonly used tests
are Downs’ generalization of the Rayleigh test (Downs, 1972) and Prentice’s generalization of
Giné’s Gn test (Prentice, 1978; Giné, 1975). For more details and references to the literature,
see Mardia and Jupp (2000, Section 13.2.2). Given a sample x1 , . . . , xn ∈ SO(3), the generalized
Rayleigh statistic is given by
T R = 3N Tr(X̄ T X̄),

X̄ =

N
1 X
xi .
N i=1

The value of T R for the naive sampler based on the sample of size N = 1500 (used to generate
figures) is 1379.452 which corresponds to p-value of virtually 0 based on 1000 simulations of
the test value under the uniform distribution. Note that T R is asymptotically distributed as a χ29
which yields an asymptotic p-value consistent with the simulation-based p-value. This is not
surprising because the difference between the two distributions is very significant. In fact, only
N = 20 observations yield an asymptotic p-value of less than 0.002 and simulated p-value of
2

100

Frequency

0 40

3.0
2.0

0.5

1.0

1.5

2.0

2.5

3.0

2.0

2.5

3.0

1.0

1.5

2.0

2.5

3.0

Uniform sampling

30 60

0.5

0

Frequency

1.0
0.0

Naive sampling

0.0

0.0

0.5

1.0

1.5

Figure 2: Rotation angles. Left: the scatter-plot of rotation angles for the uniform versus the
naive sampler. Top right: the histogram of rotation angles for the uniform sampler. Bottom right:
the histogram of rotation angles for the naive sampler. Each plot is based on N = 1500 samples.
0.003.
More generally, data from manifolds, including SO(3), appear naturally in many applications. The general problem of testing for uniformity on manifolds has been studied by various
researchers. Giné (1975), in one of the most notable contributions in this area, introduced a
method to construct non-parametric statistical tests of uniformity on compact Riemannian manifolds. These tests are called Sobolev tests, corresponding to each Sobolev norm k.k s . The case
s = 0 is briefly sketched below. Let M be a compact Riemannian manifold, µ the uniform
measure on M, and ∆ the Laplace-Beltrami operator (Laplacian) of M acting on the space of
Schwartz distributions by duality. Denote by Ek the k-th invariant eigenspace of ∆ with eigenvalue σk . Let { fi }∞
i=0 be an orthonormal basis of eigenfunctions for the Laplacian. Consider the
following hypothesis testing problem. Given observations x1 , . . . , xn ∈ M, drawn independently
P
from a probability measure ν, test whether ν = µ. Let νn be the empirical measure νn = n1 δ xi .
For a sequence of positive weights α = (α1 , α2 , . . .), with sup |αk σks | < ∞ for some s > 21 dim M,
define the test statistic T nα as
#2
∞
X
X "Z
α
fi d(νn (x) − µ) ,
T n (x) = n
αk
k=1

fi ∈Ek

M

i.e. a weighted sum of the Fourier coefficients of νn (x) − µ with respect to the orthonormal system
{ fi } with weights depending only on the eigenspaces. See Giné (1975) for more details, including
statistical properties of the Sobolev tests.
Although the program introduced in Giné (1975) constructs statistical tests with desirable
statistical properties, there remains substantial work to be done to carry the program out for any
particular example. Several authors have studied and derived Sobolev tests for different examples
including circular and directional data, tests of symmetry, and unitary eigenvalues; see Prentice
(1978); Wellner (1979); Jupp and Spurr (1983, 1985); Hermans and Rasson (1985); Baringhaus
(1991); Sengupta and Pal (2001); Coram and Diaconis (2003).
Coram and Diaconis (2003) introduce statistical tests of correspondence between zeros of
the Riemann zeta function and eigenvalues of random unitary matrices. They incorporate Giné’s
3

program for the eigenvalue distribution, µ, induced by the Haar measure on the unitary group
U(n). The orthonormal basis for L2 (µ) is given by the characters χλ of the unitary group. The
weights are chosen as αλ = z|λ| for a parameter 0 < z < 1. Given observations x1 , . . . , xN , the
test statistic is given by
T Nz =

X

z|λ|

λ,0

=

2
N
1 X
χλ (xi )
N i=1

N
1 X X |λ|
z χλ (xi )χλ (x j ).
N 2 i, j=1 λ,0

To compute T Nz , a closed form for the inner sum on the right hand side is needed. The Cauchy
identity provides a closed form for
X
z|λ| χλ (g)χλ (h),
(1)
λ,0

as given below in (4).
Sepehri (2016) considers a similar goodness-of-fit testing problem on SO(2m + 1). Let µ be
the eigenvalue distribution induced by the Haar measure on SO(2m + 1). Given data x1 , . . . , xN
drawn independently from a measure ν, test if ν = µ (note that each xi is a vector of eigenvalues).
Following the Giné’s program, the author considers the test based on
S zN

=

X
λ,0

=

z

|λ|

2
N
1 X
soλ (gi ) ,
N i=1

N
1 X X |λ|
z soλ (gi )soλ (g j ),
N 2 i, j=1 λ,0

where soλ is the character of SO(2m + 1) corresponding to the partition λ, and gi ∈ SO(2m + 1)
is a rotation with eigenvalue xi .
To use the test based on S zN , a closed form for the inner sum is needed. The Cauchy identity
(7) for the orthogonal group is derived to address this problem. It suggests a tractable method to
compute S zN . We refer the reader to Sepehri (2016) for more details and statistical properties of
the test based on S zN , as well as numerical results on synthetic and real data.
There is a sizable literature on generalizations of the Cauchy identity to larger classes of
orthogonal polynomials. A major line of research focuses on generalizing the Cauchy identity
to more general families of orthogonal polynomials, such as Hall-Littlewood polynomials Macdonald (1995); Betea and Wheeler (2016); Warnaar (2008); Kirillov and Noumi (1999); Vuletić
(2009). Although the identities introduced in these works generalize the Cauchy identity to extended families of orthogonal polynomials, which specialize to Schur function, they do not offer
any identity that results is a closed form for (1).
There is also substantial work concerning generalization of the Cauchy identity to generalized Schur functions, which include the characters of the compact groups Sundaram (1986,
1990a,b); Fu and Lascoux (2009); Brent et al. (2016); Okada (1998); Krattenthaler (2001). The
most related works among these to our case of interest seem to be Sundaram (1986, 1990a,b).
4

They introduced Cauchy-type identities for the orthogonal and symplectic groups which present
closed form expressions for
X
X
sλ (x)spλ (y± ) and
sλ (x)so2n+1,λ (y± ).
λ

λ

All of the identities introduced in this part of the literature involves series where each term
is the product of a Schur function and a generalized Schur function; whereas, what is needed in
the application of interest in this paper is an identity that involves series in which each term is a
product of two generalized Schur functions of the same type, as in (1) .
None of the aforementioned work gives closed form expressions for (1). However, it was
brought to our attention that a closed form expression for a similar quantity for the symplectic
group has recently appeared in an independent work (Wheeler and Zinn-Justin, 2016, Remark
7). In particular, it proves our identity (13) concerning the symplectic group. The author was not
aware of that work while writing this paper, and proof techniques are different. To the best of our
knowledge, closed form expressions for (1) have not appeared in the literature for the orthogonal
groups of odd and even dimensions (Type B and D root systems). The main contribution of the
current paper is to derive and prove such a generalization for other root systems.
2. Background
The Cauchy identity is a well-known result about characters of the unitary group. For two
sets of variables, x1 , . . . , xn and y1 , . . . , ym , the Cauchy identity asserts that
X

sλ (x1 , . . . , xn )sλ (y1 , . . . , ym ) =

λ

Y
i, j

1
,
1 − xi y j

(2)

where the sum is over all partitions λ of all non-negative integers and sλ is the Schur polynomial.
The special case of m = n is of particular interest here. Schur polynomials are related to characters of unitary groups in the following way: Let g ∈ U(n) with eigenvalues α1 , . . . , αn and λ be a
partition with at most n parts. Then,
χλ (g) = sλ (α1 , . . . , αn ),
where χλ is the irreducible character corresponding to λ. Thus, the representation theoretical
version of the Cauchy identity can be written as
X

χλ (g)χλ (h) =

n Y
n
Y

λ

i=1 j=1

1
,
1 − αi β j

(3)

where g ∈ U(n) has eigenvalues α1 , . . . , αn and h ∈ U(n) has eigenvalues β1 , . . . , βn . Homogeneity of Schur polynomials yields

X

z|λ| χλ (g)χλ (h) =

λ

n Y
n
Y
i=1 j=1

for any 0 < z < 1.
5

1
1 − zαi β j

(4)

The Cauchy identity (3) has a conceptual interpretation and generalization from a representation theoretical point of view. This is sketched briefly following Bump (2004, ch. 38). Let
Matn (C) be the set of all n × n complex matrices and Σ◦ , the ring of polynomials on Matn (C).
Define the action Π◦ of U(n) × U(n) on Σ◦ as follows. For (g, h) ∈ U(n) × U(n) and f ∈ Σ◦ define
(Π◦ (g, h) f ) (x) = f (gt xh).
Then, as a U(n) × U(n) representation, Π◦ decomposes into irreducible representations as
X
Π◦ 
πλ ⊗ πλ ,

(5)

λ

where λ runs through all partitions of length ≤ n, and πλ is the corresponding irreducible representation of U(n). Taking traces on both sides of (5) (in a proper sense of trace for operators
on infinite dimensional spaces) yields the Cauchy identity (3). For details of this and an excellent textbook treatment of the Cauchy identity, see Bump (2004, ch. 38). A symmetric function
theoretical point of view has been detailed in Macdonald (1995, ch. 1).
3. Generalization of the Cauchy Identity
The identities for classical groups are presented based on the type of root systems.
3.1. Type B
This section gives a Cauchy identity for the compact classical group of type B, i.e. special
orthogonal group in odd dimensions. The result for the full orthogonal group, O(2m + 1), will
follow from this. Before stating the result, recall that irreducible representations of SO(2m+1) are
labeled by partitions λ, of non-negative integers, with at most m parts, see e.g. Proposition 3.1.20
in Goodman and Wallach (2009, ch. 3). Any element g ∈ SO(2m + 1) has an eigenvalue equal
to one, and the rest of eigenvalues come in conjugate pairs. The character soλ corresponding to
partition λ is given explicitly by the Weyl Character Formula (Goodman and Wallach, 2009, ch.
7) as follows:


−(λ +m− j+ 12 )
λ +m− j+ 12
− xi j
det xi j

 ,
(6)
soλ (g) =
m− j+ 21
−(m− j+ 21 )
det xi
− xi
−1
where x1 , x1−1 , . . . , xm , xm
, 1 are eigenvalues of g. We are ready to state our first theorem:

Theorem 3.1 (Cauchy identity for SO(2m + 1)). Let soλ be the character of SO(2m + 1) corre−1
sponding to the partition λ, and g, h ∈ SO(2m+1) with eigenvalues equal to x1 , x1−1 , . . . , xm , xm
,1
−1
−1
and y1 , y1 , . . . , ym , ym , 1, respectively. Then,
X

z|λ| soλ (g)soλ (h) =

λ

m Q
z( 2 ) i< j



(1 − z)m det (C)



−1 Q
−1
−1
yi + y−1
i< j xi + xi − (x j + x j )
i − (y j + y j )

(7)

where λ runs over all partitions of non-negative integers with at most m parts, 0 < z < 1 is
an arbitrary parameter and C is an m × m matrix defined as
Ci j =

(1 + z)2 + z(xi + xi−1 + y j + y−1
j )
−1 −1
(1 − zxi y j )(1 − zxi−1 y j )(1 − zxi y−1
j )(1 − zxi y j )

6

.

Proof. The proof builds on the classical Cauchy identity in Macdonald (1995, pg. 34). Let
x1 , . . . , xm and y1 , . . . , ym be two sets of variables. Consider the matrix X defined as
X
−(k+ 1 )
−(k+ 1 )
k+ 1
k+ 1
Xi j =
zk (xi 2 − xi 2 )(y j 2 − y j 2 ).
k≥0

The (i, j) entry, Xi j , can be simplified as follows
Xi j =

X

−(k+ 12 )

k+ 21

zk (xi

− xi

k+ 12

)(y j

−(k+ 12 )

− yj

)

k
1

1

= xi2 y j2

X

−1

1

zk xik ykj − xi 2 y j2

X

k
1
2

=

1

−1

zk xi−k ykj − xi2 y j2

X

k
−1
2

1
2

xi y j

1 − zxi y j

−

1
2

zxi−1 y j

1−

k

1
2

xi y j

−

−1

−1
2

−1
2

xi y j
1−

zxi y−1
j

+

−1

2
2
zk xik y−k
j + xi y j

X

zk xi−k y−k
j

k
−1
2

xi y j

1 − zxi−1 y−1
j

,

which simplifies to
−1

1

Xi j =

1

−1

(xi2 − xi 2 )(y j2 − y j2 )(1 − z)[(1 + z)2 + z(xi + xi−1 + y j + y−1
j )]
−1 −1
(1 − zxi y j )(1 − zxi−1 y j )(1 − zxi y−1
j )(1 − zxi y j )

.

(8)

This shows the relation between X and the determinant on the right hand side of identity (7).
Now expand det(X) as follows



X
−(k+ 12 )
k+ 12
−(k+ 12 ) 
k+ 21
k

)(y j − y j
)
det(X) = det  z (xi − xi
k≥0


Y X
X
1
1
1
1 

k+
−(k+
)
k+
)
−(k+
k
 z (xi 2 − xi 2 )(yπ(i)2 − yπ(i) 2 ) .
sign(π)
=
π∈S m

i≤m

k≥0

Expanding the product on the right hand side leads to

X
Y  a +1
X
−(a + 1 )
ai + 12
−(a + 1 )
zai (xi i 2 − xi i 2 )(yπ(i)
det(X) =
sign(π)
− yπ(i)i 2 ) .
π∈S m

a1 ,a2 ,...,am ≥0 i≤m

Changing the order of summation and using the definition of m × m determinants yields
det(X) =

X
a1 ,a2 ,...,am ≥0

=

X

P

z

i

ai

X
π∈S m

sign(π)


Y  a +1
−(a + 1 )
ai + 21
−(a + 1 )
(xi i 2 − xi i 2 )(yπ(i)
− yπ(i)i 2 )
i≤m



−(a + 1 )
a +1
−(a + 1 )
a +1
z det (xi i 2 − xi i 2 )(y j i 2 − y j i 2 )
|a|

a∈Zm
+

=

X
a∈Zm
+


Y
a +1
−(a + 1 )
a +1
−(a + 1 )
z|a| det y j i 2 − y j i 2
(xi i 2 − xi i 2 ).
i

The last equality can be rewritten as follows by ordering the m-tuple a ∈ Zm .
7

X

det(X) =

X

a1 ≥...≥am ≥0 σ∈S m

i

X

Y a +1
−(a + 1 )
−(a + 1 )
a +1
z|a| det y j i 2 − y j i 2
sign(σ)
(xi σ(i) 2 − xi σ(i) 2 )

X

=



Y a +1
−(a + 1 )
−(a + 1 )
a +1
z|a| det y j i 2 − y j i 2 sign(σ)
(xi σ(i) 2 − xi σ(i) 2 )

σ∈S m

a1 ≥...≥am ≥0

=

i





−(a + 1 )
−(a + 1 )
a +1
a +1
z det y j i 2 − y j i 2 det x j i 2 − x j i 2 .

X

|a|

a1 ≥...≥am ≥0

Note that if two of ai ’s are equal then both determinants on the right hand side are zero. Thus,
it may be assumed that a1 > a2 > . . . > am ≥ 0. This is equivalent to the condition that ai −(m−i)
is a non-increasing sequence. Define λi = ai − m + i, then λ1 ≥ . . . ≥ λm ≥ 0. The last equality
translates to the following in terms of λ:
X
 
 
m
det(X) =
z|λ|+( 2 ) det Aλ det Bλ ,
λ1 ≥...≥λm ≥0
λ +m−i+ 21

Aλi, j = y j i

λ +m−i+ 12
x ji

Bλi, j =

−(λi +m−i+ 12 )

− yj
−

,

−(λ +m−i+ 12 )
xj i
.

Assume x1 , . . . , xm and y1 , . . . , ym are such that there exist g, h ∈ SO(2m + 1) with eigenvalues
−1
−1
, 1 and y1 , y−1
equal to x1 , x1−1 , . . . , xm , xm
1 , . . . , ym , ym , 1, respectively. Dividing by a factor of the
following form




m
m−i+ 1
−(m−i+ 12 )
m−i+ 1
−(m−i+ 12 )
z( 2 ) det y j 2 − y j
det x j 2 − x j
,
and using the Weyl character formula yields
m

X

z soλ (g)soλ (h) =
|λ|

λ


m−i+ 1
det y j 2

z−( 2 ) det(X)


,
−(m−i+ 21 )
m−i+ 1
−(m−i+ 21 )
− yj
det x j 2 − x j

(9)

where λ runs over all partitions of an arbitrary integer with at most m parts. To finish the proof,
the right side of (9) is simplified; from (8),
 1

−1
1
−1
 (x 2 − x 2 )(y 2 − y 2 )(1 − z)[(1 + z)2 + z(xi + x−1 + y j + y−1 )] 
i
j
i
j
j
 i

det(X) = det 

−1
−1
−1 −1

(1 − zx y )(1 − zx y )(1 − zx y )(1 − zx y )
i j

i

j

i j

i

j

Y 1
−1
1
−1
= (1 − z)m
(xi2 − xi 2 )(yi2 − yi 2 ) det (C) ,
i

where C is defined in the statement of the theorem. To get to the final form, one more simplification is needed. Note that
m−i+ 12

xj

−(m−i+ 12 )

− xj

1

−1

= (x j2 − x j2 )(xm−i
+ xm−i−1
+ . . . + x−(m−i)
)
j
j
j

leads to



 Y 1
−1
m−i+ 1
−(m−i+ 12 )
det x j 2 − x j
= det xm−i
+ xm−i−1
+ . . . + x−(m−i)
(x j2 − x j2 ) .
j
j
j
1≤ j≤m

8

Subtracting (i + 1)th row from the ith row for i = 1, 2, . . . , m − 1 yields



 Y 1
−1
−(m−i+ 21 )
m−i+ 1
det x j 2 − x j
= det xm−i
+ x−(m−i)
(x j2 − x j2 )
j
j
1≤ j≤m



m−i
= det (x j + x−1
j )

 Y

−1

1

(x j2 − x j2 )

1≤ j≤m

=

Y

(xi +

i< j

xi−1

− (x j + x−1
j ))

Y

1

−1

(x j2 − x j2 ),

1≤ j≤m

where the second equality can be proved by simple row operations, and the last equality is the
Vandermonde identity. Simplification of the determinant involving y is identical to the calculations above. Substituting these simplifications in (9) proves the theorem.
Remark 3.2. It is noteworthy that unlike the unitary group, the characters oλ of the orthogonal
group are not homogeneous functions of the eigenvalues. This perhaps is the reason that the
Cauchy identity is much more complicated for the orthogonal group than for the unitary group.
With Theorem 3.1 in hand, a Cauchy identity for the orthogonal group O(2m+1) is a straightforward corollary .
Corollary 3.3 (Cauchy identity for O(2m+1)). Let oλ be the character of O(2m+1) correspond−1
ing to the partition λ, and g, h ∈ O(2m + 1) with eigenvalues equal to x1 , x1−1 , . . . , xm , xm
, det(g)
−1
−1
and y1 , y1 , . . . , ym , ym , det(h), respectively, then
X
λ

 (1+det(gh)z)2 +z(det(h)(x +x−1 )+det(g)(y +y−1 )) 
i
j
(1 − det(gh)z)m det (1−zx y )(1−zx−1 y )(1−zxi y−1 )(1−zx−1 y−1 )j
j
i j
i j
j
i
i
Q 
,
z|λ| oλ (g)oλ (h) = m Q 
−1
−1
−1
−1
(
)
z2
i< j yi + yi − (y j + y j )
i< j xi + xi − (x j + x j )

(10)

where λ runs over all partitions of arbitrary integers with at most m parts.
Proof. Let n = 2m+1. Recall that the irreducible representations of O(n) are labeled by partitions
λ of arbitrary non-negative integers for which λ01 + λ02 ≤ n, where λ0 is the transpose of λ defined
as λ0i = #{ j : λ j ≥ i} (see Goodman and Wallach (2009, ch. 10)). For a partition λ with λ01 ≤ m,
define λ̃ as the partition with λ̃0 = (n − λ01 , λ02 , ..., λ0l ). For each partition λ with λ01 ≤ m, there exist
an irreducible representation of O(n) labeled by λ and one labeled with λ̃. Moreover, the Weyl
character formula asserts that


λ +m− j+ 21
−(λ +m− j+ 12 )
det xi j
− det(g)xi j

 ,
oλ (g) =
m− j+ 21
−(m− j+ 12 )
det xi
− det(g)xi
oλ̃ (g) = det(g)oλ (g).
The last equality justifies summing only over λ with at most m parts. The rest of the proof is
identical to that in Theorem 3.1.
3.1.1. Limit as z → 1
The Identity (7) exhibits interesting limiting behavior as z tends to 1. Assume throughout this
section that g and h share no common eigenvalues. The right hand side is equal to zero at z = 1,
9

but the left hand side does not converge at z = 1, because the series defining Xi, j do not converge.
However, under an extended notion of series convergence, i.e. Cesaro Cα -summability, the left
hand side converges to a limit at z = 1 which coincides with the value of the right hand side.
Definition 3.4 (Ck -summability). Given a sequence {an }, define
A(−1)
= an ,
n

A(k)
n =

n
X

A(k−1)
k = 1, 2, . . . .
i

i=1

Let
be the corresponding sequence
starting from the initial sequence e1 = 1, ei = 0 for
P
i > 1 . The series n an is called Cα -summable to the value s if
En(α)

A(α)
n

lim

n→∞

A(α)
n
En(α)

= s.

It is known that the geometric series 1 + z + z2 + . . . is C1 -summable for |z| = 1, z , 1
1
with C1 -limit 1−z
. This ensures that the entries Xi, j of the matrix X are C1 -summable to their
corresponding closed forms given in (8). Then, the following lemma implies that all the terms
in the expansion of det(X) are C2m−1 -summable. Therefore, the left hand side of (7) is C2m−1 summable at z = 1.
P
P
Lemma 3.5 (Theorem 277 in Knopp (1948)). If an is Cα -summable to the value A and bn is
Cβ -summable to the value B, then their Cauchy product,
X
X
cn =
(a0 bn + a1 bn−1 + . . . + an b0 ),
is certainly Cγ -summable to the value C = AB, where γ = α + β + 1.
The following lemma justifies taking the limit z → 1 on both sides of (7). This, together with
the discussion above, imply that, for z → 1, the left hand side converges to the value of the right
hand side which is equal to zero.
P
Lemma 3.6 (Theorem 278 in Knopp (1948)). If a power series f (z) = an zn has convergence
radius 1 and is Ck -summable to the value s at the point z = 1, then
f (z) → s
for every mode of approach of z to 1, in which z remains within an angle of vertex +1, bounded
by two fixed chords of the unit circle. This certainly includes the case z → 1 for z ∈ [0, 1].
P
MacPhail (1941) proved that the series n P(n)einθ is Ck -summable for any polynomial P of
degree less than k. This allows differentiating both sides of identity (7) with respect to z. More
precisely, write identity (7) as
X m
(1 − z)m det (C)


.
z( 2 )+|λ| soλ (g)soλ (h) = Q 
−1
−1 Q
−1
−1
i< j yi + yi − (y j + y j )
i< j xi + xi − (x j + x j )
λ
Differentiating m times with respect to z yields
m

∂m
m
X
m
m [(1 − z) det (C)]
2 + |λ| !
∂z
+|λ|−m
(
)
m
 z2
Q 
.
soλ (g)soλ (h) = Q 
−1
−1
−1
−1
i< j yi + yi − (y j + y j )
i< j xi + xi − (x j + x j )
λ
2 + |λ| − m !
10

Note that the left hand side is Cesaro summable at z = 1 to the value of the right hand side at
z = 1. This follows from the Leibniz rule, Cesaro summablity of the derivatives of geometric
series, and lemma 3.5 using an argument similar to the one above. Then, applying lemma 3.6,
and simplification of the right hand side, yields


u +v
 
(−1)m det (uii−v jj)2
X m + |λ|! m
2
Q 
,
z( 2 )+|λ|−m soλ (g)soλ (h) = Q 
(11)
lim
z→1
m
i< j ui − u j
i< j vi − v j
λ
where ui = xi + xi−1 + 2 and v j = y j + y−1
j + 2, and the limit is taken over real number z ∈ (0, 1)
approaching 1. This suggests using the right hand side of (11) for computational purposes.
3.2. Type C
The type C root system corresponds to the compact (real) symplectic group. The statement
and proof of the type C Cauchy identity is very similar to that of type B. Irreducible representations of the symplectic group of order m, Sp(2m), are labeled by partitions λ, of non-negative
integers, with at most m parts (see Proposition 3.1.20 in Goodman and Wallach (2009, ch. 3)).
−1
Let g ∈ Sp(2m) have eigenvalues x1 , x1−1 , . . . , xm , xm
. The character spλ corresponding to partition λ is given explicitly by the Weyl Character Formula:
 λ +m− j+1

−(λ +m− j+1)
det xi j
− xi j

 .
spλ (g) =
(12)
det xim− j+1 − xi−(m− j+1)
With this notation, the main result of this section is:
Theorem 3.7 (Cauchy identity for Sp(2m)). Let spλ be the character of Sp(2m) correspond−1
ing to the partition λ, and g, h ∈ Sp(2m) with eigenvalues equal to x1 , x1−1 , . . . , xm , xm
and
−1
−1
y1 , y1 , . . . , ym , ym , respectively. Then,


1
2 m
(1
−
z
)
det
−1 −1
X
(1−zxi y j )(1−zxi−1 y j )(1−zxi y−1
j )(1−zxi y j )
Q 

z|λ| spλ (g) spλ (h) = m Q 
−1
−1
−1
z( 2 ) i< j yi + y−1
i< j xi + xi − (x j + x j )
λ
i − (y j + y j )

(13)

where λ runs over all partitions of non-negative integers with at most m parts and 0 < z < 1
is an arbitrary parameter.
Proof. Let x1 , . . . , xm and y1 , . . . , ym be two sets of variables. Consider the matrix X defined as
X
Xi j =
zk (xik+1 − xi−(k+1) )(yk+1
− y−(k+1)
).
j
j
k≥0

Write
Xi j =

(1 − z2 )(xi − xi−1 )(y j − y−1
j )
−1 −1
(1 − zxi y j )(1 − zxi−1 y j )(1 − zxi y−1
j )(1 − zxi y j )

11

.

(14)

Expand det(X) as follows


X

−(k+1)
−(k+1)
det(X) = det  zk (xik+1 − xi
)(yk+1
− yj
)
j
k≥0


X
Y X

−(k+1)
−(k+1) 
k
k+1
k+1

 z (xi − xi
=
)(yπ(i) − yπ(i) ) .
sign(π)
π∈S m

i≤m

k≥0

Expanding the product on the right hand side leads to
X
X
Y 

ai +1
i +1)
det(X) =
sign(π)
zai (xiai +1 − xi−(ai +1) )(yπ(i)
− y−(a
) .
π(i)
π∈S m

a1 ,a2 ,...,am ≥0 i≤m

Changing the order of summation and using the definition of m × m determinant yields
det(X) =

X

z

P

i

=

sign(π)

π∈S m

a1 ,a2 ,...,am ≥0

X

X

ai

Y

i +1)
i +1
(xiai +1 − xi−(ai +1) )(yaπ(i)
− y−(a
)
π(i)



i≤m


 Y a +1
i +1)
(xi i − xi−(ai +1) ).
z det yaj i +1 − y−(a
j
|a|

a∈Zm
+

i

Rewrite the last equality as follows by ordering the m-tuple a ∈ Zm .
det(X) =

X

X

a1 ≥...≥am ≥0 σ∈S m

=

X

|a|

Y a +1


−(a +1)
i +1)
(xi σ(i) − xi σ(i) )
z|a| det yaj i +1 − y−(a
sign(σ)
j
i

z det



yaj i +1

−

i +1)
y−(a
j



det



xaj i +1

−

i +1)
x−(a
j



.

a1 ≥...≥am ≥0

If two of ai ’s are equal then both determinants on RHS are zero. Thus, assume a1 > a2 >
. . . > am ≥ 0. Define λi = ai − m + i, then λ1 ≥ . . . ≥ λm ≥ 0. The last equality translates to the
following in terms of λ:
X
 
 
m
det(X) =
z|λ|+( 2 ) det Aλ det Bλ ,
λ1 ≥...≥λm ≥0

Aλi, j

i +m−i+1)
= yλj i +m−i+1 − y−(λ
,
j

i +m−i+1)
Bλi, j = xλj i +m−i+1 − x−(λ
.
j

Assume x1 , . . . , xm and y1 , . . . , ym are such that there exist g, h ∈ Sp(2m) with eigenvalues equal
−1
−1
to x1 , x1−1 , . . . , xm , xm
and y1 , y−1
1 , . . . , ym , ym , respectively. Dividing by a factor of the following
form




m
det xm−i+1
− x−(m−i+1)
,
z( 2 ) det ym−i+1
− y−(m−i+1)
j
j
j
j
and using the Weyl character formula yields
m

X
λ

z|λ| spλ (g)spλ (h) =

z−( 2 ) det(X)


,
−(m−i+1)
det ym−i+1
− yj
det xm−i+1
− x−(m−i+1)
j
j
j


12

(15)

where λ runs over all partitions of an arbitrary integer with at most m parts. Using equation (14)


2
(xi − xi−1 )(y j − y−1


j )(1 − z )
det(X) = det 

−1
−1
−1 −1 
(1 − zxi y j )(1 − zxi y j )(1 − zxi y j )(1 − zxi y j )


Y


1
2 m
−1
−1
 .

= (1 − z )
(xi − xi )(yi − yi ) det 
(1 − zx y )(1 − zx−1 y )(1 − zx y−1 )(1 − zx−1 y−1 ) 
i j

i

i

j

i j

i

j

To get to the final form, one more simplification is needed. Note that
m−i
xm−i+1
− x−(m−i+1)
= (x j − x−1
+ xm−i−2
+ . . . + x−(m−i)
)
j )(x j
j
j
j
j

leads to



 Y
det xm−i+1
− x−(m−i+1)
= det xm−i
+ xm−i−2
+ . . . + x−(m−i)
(x j − x−1
j ).
j
j
j
j
j
1≤ j≤m

Subtracting (i + 2)th row from the ith row for i = 1, 2, . . . , m − 2 yields



 Y
(x j − x−1
det xm−i+1
− x−(m−i+1)
= det xm−i
+ x−(m−i)
j )
j
j
j
j
1≤ j≤m


 Y
m−i
(x j − x−1
= det (x j + x−1
j )
j )
1≤ j≤m

Y
Y
=
(xi + xi−1 − (x j + x−1
(x j − x−1
j ))
j ).
i< j

1≤ j≤m

Substituting these simplifications in (15) proves the theorem.
3.3. Type D
The situation is more subtle for type D, which corresponds to the special orthogonal groups
in even dimensions. Focus on SO(2m). The irreducible representations of SO(2m) are indexed by
sequences of integers λ = (λ1 , λ2 , ..., λm ) and λ− = (λ1 , λ2 , ..., λm−1 , −λm ), where λ is a partition
−1
with at most m parts. The Weyl character formula is as follows. Let x1 , x1−1 , . . . , xm , xm
are
eigenvalues of g ∈ SO(2m). If λ̃1 < m , i.e. λm = 0, then
 λ +m− j

−(λ +m− j)
det xi j
+ xi j

 .
soλ (g) =
(16)
det xim− j + xi−(m− j)
If λ̃1 = m, i.e. λm > 0, then the character of the irreducible representation corresponding to λ is
given by
 λ +m− j

 λ +m− j

−(λ +m− j)
−(λ +m− j)
det xi j
+ xi j
− det xi j
− xi j


soλ (g) =
,
(17)
2 det xim− j + xi−(m− j)
and the one corresponding to λ− is given as
 λ +m− j

 λ +m− j

−(λ +m− j)
−(λ +m− j)
det xi j
+ xi j
+ det xi j
− xi j


soλ− (g) =
.
2 det xim− j + xi−(m− j)
13

(18)

Define χλ as soλ if λ̃1 < m, and as soλ + soλ− if λ̃1 = m. Note that χλ is the restriction of the
irreducible character of O(2m) to SO(2m); it is an irreducible character of SO(2m) if and only if
λ̃1 < m. We will state the Cauchy identity in terms of {χλ }.
Theorem 3.8 (Cauchy identity for SO(2m)). Let χλ be as above, and g, h ∈ SO(2m) with eigen−1
−1
values equal to x1 , x1−1 , . . . , xm , xm
and y1 , y−1
1 , . . . , ym , ym , respectively. Then,

X
λ



det 1−zx1 i y j + 1−zx1−1 y + 1−zx1 y−1 + 1−zx1−1 y−1
j
i j
i
i
j
Q 

z|λ| χλ (g)χλ (h) = m Q 
−1
−1
−1
−1
(
)
z2
i< j yi + yi − (y j + y j )
i< j xi + xi − (x j + x j )

(19)

where λ runs over all partitions of non-negative integers with at most m parts and 0 < z < 1
is an arbitrary parameter.
Proof. Let x1 , . . . , xm and y1 , . . . , ym be two sets of variables. Consider the matrix X defined as
X
Xi j =
zk (xik + xi−k )(ykj + y−k
j )
k≥0

=

1
1
1
1
+
+
+
.
−1
−1
1 − zxi y j 1 − zxi y j 1 − zxi y j
1 − zxi−1 y−1
j

Expand det(X) to get (similar to proof of Theorem 3.1)


X

k
k
−k
k
−k

det(X) = det  z (xi + xi )(y j + y j )
k≥0


X
Y X

k
k
−k
k
−k
 z (xi + xi )(yπ(i) + yπ(i) )
=
sign(π)
π∈S m
i≤m k≥0
X
Y 
X

i)
i
zai (xiai + xi−(ai ) )(yaπ(i)
+ y−(a
sign(π)
=
π(i) )
π∈S m

=

a1 ,a2 ,...,am ≥0 i≤m
X
P
i ai

X

z

sign(π)

π∈S m

a1 ,a2 ,...,am ≥0

=

X

|a|

z det



yaj i

+

i)
y−(a
j

a∈Zm
+

Y

i)
i
(xiai + xi−(ai ) )(yaπ(i)
+ y−(a
π(i) )
i≤m

Y a
(xi i + xi−(ai ) ).
i

Rewrite the last equality as follows by ordering the m-tuple a ∈ Zm .
X

det(X) =





i)
i)
z|a| det yaj i + y−(a
det xaj i + x−(a
.
j
j

a1 ≥...≥am ≥0

Clearly a1 > a2 > . . . > am ≥ 0. Define λi = ai − m + i, then λ1 ≥ . . . ≥ λm ≥ 0. The last
equality translates to the following in terms of λ:
X
 
 
m
det(X) =
z|λ|+( 2 ) det Aλ det Bλ ,
λ1 ≥...≥λm ≥0

Aλi, j

i +m−i)
= yλj i +m−i + y−(λ
,
j

i +m−i)
.
Bλi, j = xλj i +m−i + x−(λ
j

14

Assume x1 , . . . , xm and y1 , . . . , ym are such that there exist g, h ∈ SO(2m) with eigenvalues equal
−1
−1
to x1 , x1−1 , . . . , xm , xm
and y1 , y−1
1 , . . . , ym , ym , respectively. Dividing by a factor of the following
form




m
−(m−i)
m−i
z( 2 ) det ym−i
+ y−(m−i)
det
x
+
x
,
j
j
j
j
and using definition of χλ yields
m

X
λ

z−( 2 ) det(X)



,
z χλ (g)χλ (h) =
det ym−i
+ y−(m−i)
det xm−i
+ x−(m−i)
j
j
j
j
|λ|

(20)

where λ runs over all partitions of an arbitrary integer with at most m parts. To get to the final
form note that




m−i
det xm−i
+ x−(m−i)
= det (x j + x−1
j )
j
j
Y
=
(xi + xi−1 − (x j + x−1
j )).
i< j

Substituting these simplifications in (20) proves the theorem.
Remark 3.9. A similar z → 1 analysis of Type C and Type D Cauchy identities is possible and
gives rise to similar expressions as in Type B case. Carrying out the details is straightforward
and omitted here.
4. Discussion
Although the identities introduced in the current paper provide a closed formula for an otherwise intractable infinite sum, they are still computationally expensive as they require evaluation
of an m × m determinant for m dimensional problems. Proving any alternative representations of
these determinants which result in faster computations would be of significant practical importance. This is an interesting future direction to pursue. Another interesting direction to pursue
would be the following. Find a conceptual proof of the Cauchy identity for compact groups
similar to the one sketched for the unitary group. Note that an analogue of (5) holds for other
compact groups. The main challenge is finding an appropriate basis for the coordinate ring of
the group in order to compute the traces on both sides. Of course, any new applications in other
parts of mathematics or applied problems are very much appreciated.
5. Acknowledgments
The author owes a great deal of gratitude to his doctoral advisor, Professor Persi Diaconis,
for introducing him to the subject, suggesting the problem, and for his constant encouragement
and support. The author is grateful to Professor Daniel Bump for illuminating discussions, and
to Professor Arun Ram for sharing his notes on characters of classical groups and symmetric
functions. We would like to thank Dr. Michael Wheeler for bringing to our attention their
recently published work (Wheeler and Zinn-Justin, 2016).

15

References
Baringhaus, L., 1991. Testing for spherical symmetry of a multivariate distribution. The Annals of Statistics, 899–917.
Betea, D., Wheeler, M., 2016. Refined cauchy and littlewood identities, plane partitions and symmetry classes of alternating sign matrices. Journal of Combinatorial Theory, Series A 137, 126–165.
Brent, R. P., Krattenthaler, C., Warnaar, S. O., 2016. Discrete analogues of macdonald-mehta integrals. arXiv preprint
arXiv:1601.06536.
Bump, D., 2004. Lie groups. Springer.
Coram, M., Diaconis, P., 2003. New tests of the correspondence between unitary eigenvalues and the zeros of riemann’s
zeta function. Journal of Physics A: Mathematical and General 36 (12), 2883.
Diaconis, P., Shahshahani, M., 1987. The subgroup algorithm for generating uniform random variables. Probability in
the engineering and informational sciences 1 (01), 15–32.
Downs, T. D., 1972. Orientation statistics. Biometrika 59 (3), 665–676.
Fu, A. M., Lascoux, A., 2009. Non-symmetric cauchy kernels for the classical groups. Journal of Combinatorial Theory,
Series A 116 (4), 903–917.
Giné, E., 1975. Invariant tests for uniformity on compact riemannian manifolds based on sobolev norms. The Annals of
statistics, 1243–1266.
Goodman, R., Wallach, N. R., 2009. Symmetry, representations, and invariants. Vol. 66. Springer.
Hermans, M., Rasson, J., 1985. A new sobolev test for uniformity on the circle. Biometrika 72 (3), 698–702.
Jupp, P., Spurr, B., 1983. Sobolev tests for symmetry of directional data. The Annals of Statistics, 1225–1231.
Jupp, P., Spurr, B., 1985. Sobolev tests for independence of directions. The Annals of Statistics, 1140–1155.
Kirillov, A. N., Noumi, M., 1999. q-difference raising operators for macdonald polynomials and the integrality of transition coefficients. In: Algebraic Methods and q-Special Functions, CRM Proceedings and Lecture Notes. Vol. 22. pp.
227–243.
Knopp, K., 1948. Theory and application of infinite series. Hafner New York.
Krattenthaler, C., 2001. Advanced determinant calculus. In: The Andrews Festschrift. Springer, pp. 349–426.
Kuffner, J. J., 2004. Effective sampling and distance metrics for 3d rigid body path planning. In: Robotics and Automation, 2004. Proceedings. ICRA’04. 2004 IEEE International Conference on. Vol. 4. IEEE, pp. 3993–3998.
Macdonald, I. G., 1995. Symmetric functions and Hall polynomials. Vol. 354. Clarendon press Oxford.
MacPhail, M., 1941. Cesàro summability of a class of series. Bulletin of the American Mathematical Society 47 (6),
483–487.
Mardia, K., Jupp, P., 2000. Directional statistics.
Okada, S., 1998. Applications of minor summation formulas to rectangular-shaped representations of classical groups.
Journal of Algebra 205 (2), 337–367.
Prentice, M., 1978. On invariant tests of uniformity for directions and orientations. The Annals of Statistics, 169–176.
Sengupta, A., Pal, C., 2001. On optimal tests for isotropy against the symmetric wrapped stable-circular uniform mixture
family. Journal of Applied Statistics 28 (1), 129–143.
Sepehri, A., 2016. New tests of uniformity on the compact classical groups as diagnostics for weak-∗ mixing of markov
chains. manuscript under preparation.
Sundaram, S., 1986. On the combinatorics of representations of sp(2n, c). Ph.D. thesis, Massachusetts Institute of Technology.
Sundaram, S., 1990a. The cauchy identity for sp(2n). Journal of Combinatorial Theory, Series A 53 (2), 209–238.
Sundaram, S., 1990b. Orthogonal tableaux and an insertion algorithm for so(2n + 1). Journal of Combinatorial Theory,
Series A 53 (2), 239–256.
Vuletić, M., 2009. A generalization of macmahons formula. Transactions of the American Mathematical Society 361 (5),
2789–2804.
Warnaar, S. O., 2008. Bisymmetric functions, macdonald polynomials and sl3 basic hypergeometric series. Compositio
Mathematica 144 (02), 271–303.
Wellner, J. A., 1979. Permutation tests for directional data. The Annals of Statistics, 929–943.
Wheeler, M., Zinn-Justin, P., 2016. Refined cauchy/littlewood identities and six-vertex model partition functions: Iii.
deformed bosons. Advances in Mathematics 299, 543–600.

16

