QMUL-PH-15-25

Permutation Centralizer Algebras
and Multi-Matrix Invariants

arXiv:1601.06086v1 [hep-th] 22 Jan 2016

Paolo Mattioli

1

and Sanjaye Ramgoolam

2

Centre for Research in String Theory, School of Physics and Astronomy,
Queen Mary University of London,
Mile End Road, London E1 4NS, UK

ABSTRACT
We introduce a class of permutation centralizer algebras which underly the combinatorics of multi-matrix gauge invariant observables. One family of such non-commutative
algebras is parametrised by two integers. Its Wedderburn-Artin decomposition explains
the counting of restricted Schur operators, which were introduced in the physics literature to describe open strings attached to giant gravitons and were subsequently used
to diagonalize the Gaussian inner product for gauge invariants of 2-matrix models. The
structure of the algebra, notably its dimension, its centre and its maximally commuting
sub-algebra, is related to Littlewood-Richardson numbers for composing Young diagrams.
It gives a precise characterization of the minimal set of charges needed to distinguish arbitrary matrix gauge invariants, which are related to enhanced symmetries in gauge theory.
The algebra also gives a star product for matrix invariants. The centre of the algebra
allows efficient computation of a sector of multi-matrix correlators. These generate the
counting of a certain class of bi-coloured ribbon graphs with arbitrary genus.

1
2

p.mattioli@qmul.ac.uk
s.ramgoolam@qmul.ac.uk

Contents
1 Introduction

3

2 Definitions and Key examples

7

3 Structure of the A(m, n) algebra
8
3.1 Symmetric Group characters and the pairing on the centre Z(m, n) . . . . 11
3.2 Maximal commuting subalgebra . . . . . . . . . . . . . . . . . . . . . . . . 14
4 Star product for composite operators
15
4.1 One matrix problem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
4.2 Two matrix problem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
4.3 Cartan subalgebra and the minimal set of charges . . . . . . . . . . . . . . 21
5 Computation of the finite N correlator
22
5.1 Coloured Ribbon graphs . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
6 Conclusions and future directions
6.1 Structure of the centre . . . . . . . . . . . . . . . . . . . . . . . . . . . .
6.2 Construction of quarter-BPS operators beyond zero coupling and the structure constants of A(m, n). . . . . . . . . . . . . . . . . . . . . . . . . . .
6.2.1 Central quarter BPS sector . . . . . . . . . . . . . . . . . . . . .
6.3 Non-commutative geometry and topological field theory . . . . . . . . .
6.4 Other examples of permutation centralizer algebras and correlators . . .

24
. 25
.
.
.
.

26
27
27
28

A Analytic formula for the dimension of M(m, n)
30
A.1 Multi-variable polynomials . . . . . . . . . . . . . . . . . . . . . . . . . . 31
B LR
B.1
B.2
B.3
B.4
B.5

rule for hook
(Y,Y) case . .
(Y,N) case . .
(N,N) case . .
(N,Y) case . .
A summary .

representations
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .

.
.
.
.
.

.
.
.
.
.

C Deriving the two point correlator

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

.
.
.
.
.

33
34
35
36
36
37
39

2

1

Introduction

A number of questions on gauge invariant functions and correlators of multiple-matrices
have been studied in the context of N = 4 Super Yang-Mills (SYM). The impetus for these
developments in physics has come from the AdS/CFT correspondence [1, 2, 3], notably
the duality between the N = 4 SYM theory with U (N ) gauge group and AdS5 ×S 5 . Local
composite operators are U (N ) gauge invariants. CFT gives extra motivation because of
the operator-state correspondence. Quantum states correspond to local operators, which
are composite fields. These can be matrix-valued fields which are space-time scalars,
fermions, field strengths or covariant derivatives of these. A generic problem is to understand U (N ) invariants constructed from a number n of such fields
j1
jn
F1,i
· · · Fn,i
n
1

(1.1)

This is subsequently used to understand their correlation functions. The n upper indices
each transform in the fundamental of U (N ) while the lower indices transform in the
anti-fundamental. Hence, an important ingredient is the nature of the invariants in
V ⊗n ⊗ V̄ ⊗n

(1.2)

The number of linearly independent invariants is n!. They are obtained by multiplying
(1.1) with a product of n Kronecker delta functions, contracted with a permutation σ ∈ Sn .
As n varies, we are interested in all possible values of n, so the properties of
C[S∞ ] =

∞
M

C[Sn ]

n=0

become important. If all the n operators are the same e.g. a complex matrix X = φ1 +iφ2
where φ1 , φ2 are two of the six hermitian matrices transforming in the vector of SO(6),
then the invariants are multi-traces, of which there are p(n), the number of partitions of
n. In terms of the permutations, the composite operators are
1
n
Oσ (X) = Xiiσ(1)
· · · Xiiσ(n)

(1.3)

Distinct σ related by conjugation, i.e. σ and γσγ −1 for some γ ∈ Sn give the same
operator
Oσ (X) = Oγσγ −1 (X)

(1.4)

When we consider invariants built from two types of matrices, say m copies of X and n
copies of Y , then we encounter equivalence classes
σ ∼ γσγ −1

(1.5)

where σ ∈ Sm+n and γ ∈ Sm × Sn .
The fact that the enumeration of gauge invariant operators can be effectively done by
using a formulation in terms of equivalence classes of permutations has driven significant
progress in the construction of operators and computation of correlators for the half-BPS
sector, the perturbations of the half-BPS operators as well as quarter BPS operators. Two
key facts have been used. One is that, by using the Fourier transformation which relates
functions on a group to matrix elements of irreducible representations, nice orthogonal
3

bases of functions on these equivalence classes can be found. In mathematics, in the
context of compact groups this is known as the Peter-Weyl theorem. In the context
of finite groups, this follows from the Schur orthogonality relations. This leads to the
construction of operators in the half-BPS sector parametrised by Young diagrams [4,
5]. For the two-matrix sector, one application of this thinking leads to restricted Schur
operators. These are labelled by three young diagrams and a pair of multiplicity labels: a
Young diagram R1 with m boxes, a Young diagram R2 with n boxes and a third diagram
R with m + n boxes. The two multiplicity labels each run over a space of dimension
equal to g(R1 , R2 ; R), which is equal to the Littlewood-Richardson (LR) coefficient for
the number of times R appears in the tensor product of R1 ⊗ R2 [6, 7, 8, 9, 10]. LR
coefficients will be reviewed as needed in this paper (see Appendix B).
One reason for the efficacy of permutation groups in enumeration of gauge invariant
operators is Schur-Weyl duality. This states that the tensor product of n copies of the
fundamental of U (N ) decomposes into a direct sum of irreps of Sn × U (N )
M
U (N )
VN⊗n =
VRSn ⊗ VR
(1.6)
R`N
c1 (R)≤N

Each summand is labelled by a Young diagram, and the Young diagrams are constrained
to have no more than N rows, equivalently the first column c1 (R) is no greater than
N . This uses the fact that Young diagrams are used to classify representations of Sn as
well as representations of U (N ). This is useful in the permutation approach to gauge
invariant operators, because it says that once we have organised operators according to
representation data for Sn , it is easy to implement finite N constraints. In the one-matrix
problem, the single Young diagram label R is cut-off at N , c1 (R) ≤ N . This leads
directly to the connection between the stringy exclusion principle for giant gravitons and
Young diagrams [11, 12, 13, 4]. In the two-matrix problem, the Young diagram R is
cut-off at c1 (R) ≤ N , which implies cut-offs for R1 , R2 . The 2-matrix problem can also
be approached using the walled Brauer algebra BN (m, n) and its representation theory
[14]. A third way to enumerate two-matrix invariants, also based on permutations but
involving Clebsch-Gordan multiplicities of Sn , keeps the U (2) global symmetry manifest
[15, 16].
Aside from enumerating gauge invariant operators, the permutation structures have
been used to compute correlators. Correlators in free field theory are obtained by sums
over Wick contractions. These sums are themselves parametrised by permutations. Correlators of gauge invariant operators are thus given in terms of these Wick permutations
and the permutations which enumerate the operators. Hence there are elegant formulae
for the correlation functions in terms of permutations. It can be shown that the 2-point
functions of gauge invariant operators in the 2-matrix sector are diagonalized by operators
constructed using representation bases. This was done with the Brauer basis in [14], with
the U (2) covariant basis in [15, 16] and with the restricted Schur basis in [17, 18]. The
restricted Schur and covariant basis results have been extended beyond N = 4 SYM to
the sector of holomorphic operators in general quiver gauge theories [19, 20, 21, 22, 23, 24]
which have been shown to include sectors related to generalized oscillators [25]. Aspects
involving Frobenius algebras have been studied in [26]. Within N = 4 SYM, perturbations of half-BPS giant graviton operators have been studied and integrability at one-loop
[27, 28, 29, 30, 27, 24] and beyond has been established.
As a way to understand the existence of the different bases in the multi-matrix problems, the paper [31] conducted a detailed study of enhanced symmetries in the free limit
4

of Yang Mills theories. The authors showed that Casimir-like elements constructed from
Noether charges of these enhanced symmetries can be used to understand these different
bases. Different sets of these Casimir-like charges each consist of mutually commuting simultaneously diagonalizable operators, which associate the labels of the basis with
eigenvalues of Casimir-like charges. Thus there is a set of Casimir-like elements for the
restricted Schur basis, another set for the covariant basis and yet another set for the
Brauer basis. The enhanced symmetries themselves take the form of products of unitary
groups, but the action of these Casimirs on gauge invariant operators can be related,
through applications of Schur-Weyl duality, to the algebraic structure of certain algebras
constructed from the equivalence classes of permutations or of Brauer algebra elements
discussed above. The discussion of charges which identify matrix invariants for general
classical groups has been given using a different approach in [32]. While a uniform treatment of the Young diagram labels has been achieved, a treatment of the multiplicity labels
running over Littlewood-Richardson coefficients in that approach remains an interesting
open problem.
This paper was motivated by the goal of obtaining a systematic understanding of the
algebraic structures involved in the construction of charges in [31]. To be more precise,
we will define the notion of permutation centralizer algebras. A particular class of these,
denoted as A(m, n), will be our main focus. Many of the important formulae we will use
have already appeared in the physics literature. Nevertheless the A(m, n), as associative
algebras with non-degenerate pairing, have not been made fully explicit. This paper
proposes that these algebras are interesting to study intrinsically, disentangled from the
contingencies of being embedded in a bigger symmetric group algebra, their simplicity
hidden among the application to matrix correlators for matrices of size N . Here we
define the algebras A(m, n), study their structure, and subsequently describe how they
are relevant to matrix theory invariants. We expect that a deeper study of this algebraic
structure has the potential to give a lot of information about correlators in free Yang-Mills
theory, in the loop corrected theory, at all orders in the 1/N expansion. This paper is a
step in this direction. Much as it is valuable to abstract Riemannian geometry from the
study of submanifolds of Euclidean spaces, abstracting a family of algebras intrinsic to
permutations hidden in the mathematics of matrix theory should be fruitful.
We describe the organization of the paper. In section 2 we introduce the definition
of permutation centralizer algebras. We consider four key examples of these algebras,
which are useful in the context of gauge-invariant operators. In section 3, we focus on the
algebras A(m, n) formed by equivalence classes of permutations in Sm+n , with equivalence
generated by conjugation with permutations in Sm × Sn . The dimension of this algebra is
X
|A(m, n)| =
g(R1 , R2 ; R)2
(1.7)
R1 `m,R2 `n
R`m+n

where g(R1 , R2 ; R) is the LR coefficient for the triplet of Young diagram (R1 , R2 , R) made
with (m, n, m + n) boxes respectively. We will show that this is an associative algebra
with a non-degenerate pairing. As a result, we know from the Wedderburn-Artin theorem
that it is isomorphic to a direct sum of matrix algebras Mat [33, 34]:
M
A(m, n) =
Mata
(1.8)
a

In eq. (3.5) we give a more precise version of this formula, where the index a is identified
with triplets (R1 , R2 , R) with non-vanishing LR coefficient g(R1 , R2 ; R). The construc5

tion of restricted Schur operators in gauge theory is used to give the Wedderburn-Artin
decomposition of A(m, n). Two sub-algebras will be of interest. The centre of the algebra
Z(m, n) is the subspace of the algebra which commutes with any element of A(m, n). The
dimension of this centre is equal to the number of triples (R1 , R2 , R) of Young diagrams,
with numbers of boxes equal to (m, n, m + n), for which the LR coefficient is non-zero.
It is useful to develop some formulae for the non-degenerate pairing on the centre, using characters of Sm+n , Sm , Sn . The Wedderburn-Artin decomposition also highlights
the importance of a maximally commuting sub-algebra M(m, n). The dimension of this
sub-algebra is the sum of Littlewood-Richardson coefficients g(R1 , R2 ; R). Appendix A
gives a multi-variable generating function for this sum of LR coefficients. We explain the
relevance of the this sub-algebra to the enhanced symmetry charges studied in [35]. In
particular we give a precise algebraic characterization (4.45) for the minimal number of
charges needed to identify all 2-matrix gauge-invariant operators. The evaluation of this
number is an open problem for the future.
In section 4, we explain some further physical implications of the permutation centralizer algebras. The simplest of these algebras is the algebra of class sums of permutations.
Given the one-to-one correspondence between matrix operators and conjugacy classes of
permutations given in (1.3), this means that there is a corresponding product on halfBPS operators. This is not the usual product obtained by multiplying the gauge invariant
operator built from X under which the dimension of the operator adds. The product on
the class sums rather gives a product for the BPS operators of fixed dimension, a product
which is associative and admits a non-degenerate pairing. We will refer to this as a star
product for half-BPS operators. We explain the relevance of this star product for the
computation of correlators. Similarly the product on the algebra A(m, n) gives a star
product for gauge invariant polynomials in two matrices, with degree m in the X’s and
degree n in the Y ’s. In the physics application, there is a closed associative star product
on the space of quarter-BPS operators at zero Yang-Mills coupling. Conversely the usual
product of gauge invariants gives a product on A(∞, ∞)
A(∞, ∞) =

∞
M

A(m, n)

(1.9)

m,n=0

which is the direct sum over all m, n. Thus A(∞, ∞) has two products one of which closes
at fixed m, n. This generalizes a structure seen in the study of symmetric polynomials.
In section 5, we show that the study of the structure of the algebra A(m, n) we
developed in section 3 is useful for the computation of correlators of 2-matrix gauge
invariants. In particular, we identify an efficiently computable sector of central gauge
invariant operators whose correlators can be computed using the knowledge of characters
of Sm+n , Sm , Sn . It does not require the knowledge of more detailed data such as matrix
R
elements Dij
(σ) or branching coefficients for Sm+n → Sm ×Sn . To illustrate the simplicity
of this central sector, we compute the two-point function
Tr(X m Y n )Tr((X † )m (Y † )n )

(1.10)

at finite N . The computation requires a calculation of Littlewood-Richardson coefficients
g(R1 , R2 ; R) where R1 , R2 are hook-shaped Young diagrams. This computation is given in
Appendix B. Further technical aspects of the computation are given in Appendix C. The
computation agrees with the one in [36] which was done with explicit Young-Yamanouchi
symbols which can be used to construct states in irreps R and describe their reduction to
R1 , R2 .
6

In section 6, we outline some future research directions related to the present results.

2

Definitions and Key examples

When studying the representation theory of a group G, it is useful to introduce the algebra
C[G] which consists of formal linear combinations of group elements, equipped with the
multiplication inherited from the group. In the group algebra C[G], for each conjugacy
class, we can form a sum over all the elements in the conjugacy class of G. Such class
sums commute with any element of G and form the central sub-algebra of C[G], i.e. the
sub-algebra which commutes with all C[G]. We will refer to Z[C[G]] as the centre of C[G].
Conjugacy classes are in 1-1 correspondence with irreducible representations and there is
a basis of the centre consisting of projectors of the form
PR =

dR X
χR (g) g −1
|G| g∈G

(2.1)

Of primary interest to us is the group algebra of C[Sn ] and its centre Z[C[Sn ]]. The
elements in Z[C[Sn ]] are sums over conjugacy classes t of Sn
X
T =
σ
(2.2)
σ∈t

Given any σ ∈ Sn , we can generate an element of this subalgebra by summing over γ ∈ H.
X
γσγ −1
(2.3)
γ∈Sn

Some properties of group algebras and their centre can be found in [37, 33]. In the context
of AdS/CFT , group algebras C[Sn ] and associated representation theory play a role in
the half-BPS sector of N = 4 SYM in 4D [4, 5] and also in the symmetric orbifolds
in AdS3/CFT2 [11, 38]. Motivated by developments in AdS/CFT we will introduce a
generalization of this construction.
Definition: Consider an associative algebra A containing a sub-algebra B = C[H], the
group algebra of a finite group H. Now define the sub-space of A of elements which are
invariant under conjugation by H. This subspace will contain group averages of the form
X
γσγ −1 ,
σ∈A
(2.4)
γ∈H

which commute with elements of B. It is easy to verify that these sub-spaces are subalgebras. We have
!
!
!
X
X
X
X
γ1 σγ1−1
γ2 σγ2−1 =
γ1
σ1 γ3 σ2 γ3−1 γ1−1
(2.5)
γ1 ∈H

γ2 ∈H

γ1 ∈H

γ3 ∈H

where we set γ3 = γ1−1 γ2 . This shows that the product of two group averages is still a group
average. This sub-algebra of A commuting with B, in cases where H is a permutation
group, will be called a permutation centralizer algebra.
Three cases of primary interest will be
7

• Example 1 The algebra A = C[Sn ]. The algebra B = C[Sn ]. The centralizer of B
is Z[C[Sn ]] .
• Example 2 A = C[Sm+n ] ; B = C[Sm × Sn ]. We will call this algebra A(m, n).
• Example 3 A = BN (m, n) - the walled Brauer algebra ; B = C[Sm × Sn ]. This
algebra is called BN (m, n).
• Example 4 A = C[Sn × Sn ] ; B = C[Sn ] where the latter is the Sn diagonally
embedded in the product group. This should be called K(n).
The case where A is itself a group algebra has been studied in mathematics, for example,
in [39].
Our primary interest in this paper will be in A(m, n) of example 2. Z[C[Sn ]] of Example 1 will be a useful guide and a source of analogies in our investigations. Fourier transformation on A(m, n) will be related to restricted Schur operators studied in AdS/CFT.
These are parametrised by representation theory data (R, R1 , R2 , i, j) consisting of Young
diagrams R1 , R2 , R with m, n, m + n boxes as well as multiplicity indices i, j. The latter
take values 1 ≤ i, j ≤ g(R1 , R2 ; R) where g(R1 , R2 ; R) is the LR multiplicity for the triple
of Young diagrams computed with the LR combinatoric rule (see for example [40]). Unlike Z[C[Sn ]], the algebra A(m, n) is not commutative. The central sub-algebra Z(m, n),
consisting of the subspace Z(m, n) ⊂ A(m, n) which commutes with all of A(m, n) will
play a predominant role. Likewise the algebras BN (m, n) and K(n) in Examples 3 and 4
are non-commutative.

3

Structure of the A(m, n) algebra

The algebra A(m, n) is constructed by taking all the elements in C[Sm+n ] which are
invariant under C[Sm ×Sn ]. Any element of σ ∈ C[Sm+n ] can be mapped to a σ̄ ∈ A(m, n)
by the group averaging
X
σ̄ =
γ −1 σγ
(3.1)
γ∈Sm ×Sn

The σ̄ are formal sums of permutations τ lying in the same orbit of σ under the Sm × Sn
action. Each τ has a stabiliser group, given by those γ ∈ Sm × Sn for which
γ −1 τ γ = τ

(3.2)

The stabilisers of two permutations τ1 , τ2 in the same orbit are generally different (they
are conjugate to each other), but they have the same dimension. By the Orbit-Stabiliser
theorem, σ̄ is then a sum of permutations weighted by the same coefficient:
X
σ̄ = |AutSm ×Sn (σ)|
τ
(3.3)
τ ∈Orbit(σ, Sm ×Sn )

A(m, n) is a finite-dimensional associative algebra (the associativity follows from the
associativity of C[Sm+n ]), which we can equip with the non-degenerate symmetric bilinear
form
hσ̄1 , σ̄2 i = δ(σ̄1 σ̄2 ) ,
8

σ̄1,2 ∈ A(m, n)

(3.4)

Here the delta function on the group algebra C[Sm+n ] is a linear function which obeys
δ(σ) = 1 for σ = 1 and δ(σ) = 0 otherwise.
The non-degeneracy of the bilinear form (3.4) implies that A(m, n) is semi-simple.
According to the Wedderburn-Artin theorem, it can then be decomposed into a direct
sum of matrix algebras:
M
A(m, n) =
Span{QR
(3.5)
R1 ,R2 ,i,j ; i, j}
R`m+n
R1 `m, R2 `n

In this equation R, R1 and R2 are representations of Sm+n , Sm and Sn respectively.
The integers i, j run over the multiplicity g(R1 , R2 ; R) of the branching R → R1 ⊗ R2 :
0 ≤ i, j ≤ g(R1 , R2 ; R). An explicit expression for QR
R1 ,R2 ,i,j is given in terms of the
restricted Schur characters [17, 31, 22], defined as
R→R1 ,R2 ;i
R→R1 ,R2 ;j
R
Bm→l
χR
R1 ,R2 ,i,j (σ) = Dm,m0 (σ) Bm0 →l1 ,l2
1 ,l2

(3.6)

R→R1 ,R2 ;j
R
Here Dm,m
0 (σ) are the matrix elements of σ in the irreducible representation R. Bm→l ,l
1 2
is the branching coefficient for the representation branching R → R1 ⊗ R2 , in the j-th
copy of R1 ⊗ R2 ⊂ R. l1,2 are states in R1,2 . The restricted Schur characters χR
R1 ,R2 ,i,j (σ)
are invariant under conjugation by C[Sm × Sn ] elements. With these definitions we can
write
X
χR
(3.7)
QR
=
R1 ,R2 ,i,j (σ) σ
R1 ,R2 ,i,j
σ

which is manifestly invariant under the action of C[Sm × Sn ]. It follows that
S
R,S
QR
δR1 ,S1 δR2 ,S2 (δjk QR
R1 ,R2 ,i,j QS1 ,S2 ,k,l = δ
R1 ,R2 ,i,l )

(3.8)

This is in accordance with the decomposition (3.5). Consequently it is useful to write
QR
R1 ,R2 ,i,j as
X
QR
=
|R → R1 , R2 , m1 , m2 , iihR → R1 , R2 , m1 , m2 , j|
(3.9)
R1 ,R2 ,i,j
m1 ,m2



Moreover, the basis QR
R1 ,R2 ,i,j is complete as we now explain. The number of distinct
is equal to the number of restricted Schur characters, which is in turn equal
to R1 ,R2 ,R g(R1 , R2 ; R)2 . On the other hand the dimension of A(m, n) is by definition
equal to the number of elements of C[Sm+n ] invariant under the C[Sm × Sn ] action. Using
the Burnside lemma, it is possible to show that this dimension |A(m, n)| is given as
X
|A(m, n)| =
g(R1 , R2 , R)2
(3.10)
’s
QR
RP
1 ,R2 ,i,j

R1 `m,R2 `n
R`m+n

P
In each of the blocks in (3.5) there is a projector of the form PRR1 ,R2 = i QR
R1 ,R2 ,i,i .
Let now PR , PR1 and PR2 be the projectors onto the irreps R, R1 and R2 of Sm+n , Sm and
Sn respectively. Since
hR → R1 , R2 , m1 , m2 , i|PR PR1 PR2 |R → R1 , R2 , m01 , m02 , ji
= R → R1 , R2 , m1 , m2 , i|PRR1 ,R2 |R → R1 , R2 , m01 , m02 , j = δm1 ,m01 δm2 ,m02 δi,j
9

(3.11)

for all triplets R, R1 , R2 , we can write
PRR1 ,R2 = PR PR1 PR2

(3.12)

so that the projectors PRR1 ,R2 are just products of ordinary Sm+n , Sm and Sn projectors.
The set {PRR1 ,R2 } forms a basis for the centre of A(m, n), which we call Z(m, n). Its
dimension is then given by the number of non vanishing LR coefficients g(R1 , R2 ; R), or
X
|Z(m, n)| =
(1 − δ(g(R1 , R2 ; R)))
(3.13)
R1 `m,R2 `n
R`m+n

Here δ(g(R1 , R2 ; R)) = 1 if g(R1 , R2 ; R) = 0 and δ(g(R1 , R2 ; R)) = 0 otherwise. The
generating function for the dimension of the centre is [41]
Z(x, y) =

Y
i

1
(1 − xi − y i )

(3.14)

We will now argue that the collection of the generators of the centres of C[Sm+n ],
(m+n)
(m)
(n)
C[Sm ] and C[Sn ], that we denote as {Tp
}, {Tq1 } and {Tq2 } respectively, is a set
of generators for Z(m, n). Here p, q1 and q2 are integer partitions of m + n, m and n
(m+n)
respectively. For example, for the partition p = (p1 , p2 , ...) of m + n, the operator Tp
consists of a sum over permutations belonging to the conjugacy class p = (p1 , p2 , ...):
X
Tp(m+n) =
(i1 i2 · · · ip1 )(ip1 +1 ip1 +2 · · · ip1 +p2 ) · · ·
(3.15)
i1 ,··· ,ip1 +p2 +··· ∈[m+n]
(m+n)

(m)

(n)

Tp
are sums of conjugates by elements of Sm+n , whereas Tq1 and Tq2 are sums over
(m+n)
(m)
(n)
Sm ⊂ Sm+n and Sn ⊂ Sm+n respectively. To show that {Tp
, Tq1 , Tq2 } generate the
whole centre Z(m, n) we can use the following argument. Using the Wedderburn-Artin
decomposition (3.5), we see that the centre of A(m, n) is the direct sum of the centres of
the matrix algebras Span{QR
R1 ,R2 ,i,j ; i, j}. For each of these matrix blocks, that is for any
fixed representations R, R1 , R2 for which g(R1 , R2 ; R) 6= 0, the centre is one-dimensional,
and is spanned by
X
PRR1 ,R2 =
QR
(3.16)
R1 ,R2 ,i,i
i=1

Using the equation (3.8), it is immediate to check that
[PRR1 ,R2 , QR
R1 ,R2 ,i,j ] = 0 ,

∀ i, j

(3.17)

We know that PRR1 ,R2 = PR PR1 PR2 , with PR , PR1 and PR2 projectors on the representations
R, R1 and R2 . Therefore every central element of A(m, n) can be generated with the
collection of projectors {PR , PR1 , PR2 }. For an R irrep of Sn , the projector is
PR =

1
1 X
χR (σ) σ =
n! σ∈S
n!
n

X

χR (σp ) Tp(n)

(3.18)

p∈Partitions(n)

where σp is a representative permutation belonging to the conjugacy class p ` n. This
means that every projector PR can be written as a linear combination of the central
10

(n)

elements {Tp }. We can then write the set {PR , PR1 , PR2 } in terms of the central elements
(m+n)
(m)
(n)
{Tp
, Tq1 , Tq2 }. Since we know that the former generates the whole Z(m, n), we
can now conclude that the latter is a complete set of generators for the centre Z(m, n)
as well. The basis thus obtained will be useful in the following sections. However, it is
important to point out that such a basis is overcomplete. An easy way to see it is to
note that, given (3.12), PR PR1 PR2 = 0 if g(R1 , R2 , R) = 0. Therefore, taking a triplet
(R1 , R2 , R) for which g(R1 , R2 , R) = 0 we have, using (3.18):
X
1
χR (σp )χR1 (σq1 )χR2 (σq2 ) Tp(m+n) Tq(m)
Tq(n)
=0
(3.19)
1
2
(m + n)!m!n! p`(m+n)
q1 `m, q2 `n

(m+n)

(m)

(n)

This shows that {Tp
, Tq1 , Tq2 } is indeed an overcomplete basis.
(m+n)
(m)
(n)
We can also argue that {Tp
, Tq1 , Tq2 } generate Z(m, n) just by using the SchurWeyl duality as in [31]. The T (m) elements are Schur-Weyl dual to U (N ) Casimirs of acting
on the upper m indices of X-type matrices. This action is generated by
∂
(Ex )ij = (Dx )ij = Xli
(3.20)
∂Xlj
The T (n) elements are Schur-Weyl dual to U (N ) Casimirs acting on the upper n indices
of Y -type matrices. We have
∂
(Ey )ij = (Dy )ij = Yli j
(3.21)
∂Yl
Finally, the T (m+n) elements are Schur-Weyl dual to U (N ) Casimirs acting on the upper
n and m indices of both X- and Y -type matrices, and the generator is
Eji = (Ex )ij + (Ey )ij

(3.22)

We then have three distinct types of Casimirs:
(m+n)

Ck

(m)

= (Ex )ii12 (Ex )ii23 · · · (Ex )iik1

(n)

= (Ey )ii12 (Ey )ii23 · · · (Ey )iik1

Ck

Ck
(m+n)

(m)

= Eii21 Eii32 · · · Eii1k

(3.23)

(n)

But the Ck
, the Ck and the Ck operators measure respectively the R, R1 and
R2 labels of the restricted Schurs χR
R1 ,R2 ,i,j . Therefore they can be used to isolate every
subspace R1 ⊗ R2 ⊆ R, and to build all the correspondent projectors PRR1 ,R2 . Since we
know that each of these projectors is in a 1-1 correspondence with an element of Z(m, n),
the whole centre Z(m, n) is obtained.
On the other hand, non-central elements are needed to measure the multiplicity labels
i, j. This observation will be developed in section 4.

3.1

Symmetric Group characters and the pairing on the centre
Z(m, n)

A central element Za ∈ Z(m, n) can be expanded in terms of the projectors PRR1 ,R2 as
X
Za =
ZaR,R1 ,R2 PRR1 ,R2
(3.24)
R,R1 ,R2

11

We can then define
X
hR → R1 , R2 , m1 , m2 , i|Za |R → R1 , R2 , m1 , m2 , ji
χR
R1 ,R2 ;i,j (Za ) =
m1 ,m2

= δij

X X

ZaS,S1 ,S2 hR → R1 , R2 , m1 , m2 , i|PSS1 ,S2 |R → R1 , R2 , m1 , m2 , ji

S,S1 ,S2 m1 ,m2

= δij ZaR,R1 ,R2 dR1 dR2

(3.25)

and
χR
R1 ,R2 (Za ) =

X

R,R1 ,R2
g(R1 , R2 , R)dR1 dR2
χR
R1 ,R2 ;i,i (Za ) = Za

(3.26)

i

From these equations it also follows that for any central element Za
χR
R1 ,R2 ,i,j (Za ) =

δi,j
χR (Za )
g(R1 , R2 ; R) R1 ,R2
(m+n)

Another useful expansion is in terms of {Tp
elements generate the centre, we can write

(m)

(3.27)
(n)

}, {Tq1 } and {Tq2 }. Since these

Za = Zap,q1 ,q2 Tp(m,n) Tq(m)
Tq(n)
1
2

(3.28)
(m+n)

(m)

(n)

for some Zap,q1 ,q2 coefficients. However, since the basis generated by {Tp
, Tq1 , Tq2 }
is overcomplete, such coefficients are not unique. Using the expansion (3.28), we can write
(m+n)

p,q1 ,q2
χR
R1 ,R2 ,i,j (Za ) = δij Za

χR (Tp
dR

)

χR1 (Tq(m)
)χR2 (Tq(n)
)
1
2

(3.29)

and
χR
R1 ,R2 (Za )

=

X

χR
R1 ,R2 ,i,i (Za )

i

=

(m+n)
χR (Tp
)
p,q1 ,q2
χR1 (Tq(m)
)χR2 (Tq(n)
Za
)
g(R1 , R2 , R)
1
2
dR

(3.30)
From these equation we see that all the restricted characters of central elements are
determined by characters of Sm+n , Sm , Sn . Just as the centre of Sn is generated by class
sums, which are dual to irreducible characters of Sn , the centre Z(m, n) of A(m, n) is
dual to the characters χR
R1 ,R2 which are nothing by products of characters. Therefore, to
compute restricted characters of elements in Z(m, n) we only need the ordinary symmetric
group character theory.
We will now use some of the known equations for the character of symmetric group
and use them to compute restricted characters in Z(m, n). Our aim will be to compute
the dual pairing (3.4) for central elements. Equation (B.12) in [22] reads
X
dR
(m + n)! X
R
χR
δ(σγτ γ −1 ) =
R1 ,R2 ,i,j (σ)χR1 ,R2 ,i,j (τ )
m!n! γ∈S ×S
dR 1 dR 2
R,R ,R ,i,j
m

n

1

(3.31)

2

By setting τ = 1 this equation simplifies to
(m + n)!δ(σ) =

X
R

12

dR χ R
R1 ,R2 ,i,i (σ)

(3.32)

where we used
χR
R1 ,R2 ,i,j (1) = δij dR1 dR2

(3.33)

We can immediately use this result to show that δ(QR
R1 ,R2 ,i,j ) = δij dR1 dR2 . This is because,
using (3.7)
 X R
δ QR
χR1 ,R2 ,i,j δ (σ) = χR
(3.34)
R1 ,R2 ,i,j =
R1 ,R2 ,i,j (1) = δij dR1 dR2
σ

It is also worthwhile to notice that, for O ∈ A(m, n), Tr(O) = δ(O). Therefore we could
have obtained the same result by considering
X X X
)
=
Tr(QR
hS → S1 , S2 , m01 , m02 , k|R → R1 , R2 , m1 , m2 , ii
R1 ,R2 ,i,j
S,S1 ,S2

m1 ,m2
m01 ,m02

k

× hR → R1 , R2 , m1 , m2 , j|S → S1 , S2 , m01 , m02 , ki
= δij dR1 dR2

(3.35)

where we used the definition (3.9).
Let us now go back to eq. (3.32). If we replace σ by a central element Za , using the
expansion (3.28) and eq. (3.30), we find
X
(m + n)! δ(Za ) =
Zap,q1 ,q2 g(R1 , R2 , R)χR (Tp(m+n) )χR1 (Tq(m)
)χR2 (Tq(n)
)
(3.36)
1
2
R,R1 ,R2

By further replacing σ → Za , τ → Zb in (3.31) we get, in a similar fashion
X
dR
R
(m + n)! δ(Za Zb ) =
χR
R1 ,R2 ,i,j (Za )χR1 ,R2 ,i,j (Zb )
d
d
R1 R2
R,R ,R ,i,j
1

2

p0 ,q10 ,q20

= Zap,q1 ,q2 Zb

X g(R1 , R2 , R)
×
dR dR1 dR2
R,R ,R
1

(3.37)

2

(m+n)

× χR (Tp(m+n) ) χR1 (Tq(m)
) χR2 (Tq(n)
)χR (Tp0
1
2

(m)

(n)

) χR1 (Tq0 ) χR2 (Tq0 )
1

2

Comparing the LHS above with eq. (3.4) we find that for central elements Za , Zb
X g(R1 , R2 , R)
1
p0 ,q 0 ,q 0
×
(3.38)
hZa , Zb i = Zap,q1 ,q2 Zb 1 2
(m + n)! R,R ,R dR dR1 dR2
1

2

(m+n)

)χR (Tp0
× χR (Tp(m+n) ) χR1 (Tq(m)
) χR2 (Tq(n)
1
2

(m)

(n)

) χR1 (Tq0 ) χR2 (Tq0 )
1

2

Thus we have an explicit way of computing the dual paring on the centre Z(m, n) in
terms of ordinary Sn characters.
Similarly, there is a character expansion for δ(Za Zb Zc ). We begin by writing
X
dR
R
χR
(m + n)! δ(Za Zb Zc ) =
R1 ,R2 ,i,j (Za Zb )χR1 ,R2 ,i,j (Zc )
dR1 dR2
R,R ,R ,i,j
1

=

X
R,R1 ,R2

2

dR
χR (Za Zb )χR
R1 ,R2 (Zc )
dR1 dR2 g(R1 , R2 ; R) R1 ,R2
13

(3.39)

R
Since Za is central, Za = (ca )R
R1 ,R2 1, where (ca )R1 ,R2 is a constant. This constant can be
obtained by considering:
R
R
R
χR
R1 ,R2 ,i,j (Za ) = (ca )R1 ,R2 χR1 ,R2 (1) = (ca )R1 ,R2 dR1 dR2 g(R1 , R2 ; R)

(3.40)

We therefore have that
χR
R1 ,R2 (Za Zb ) =

R
χR
R1 ,R2 (Za ) χR1 ,R2 (Zb )
dR1 dR2 g(R1 , R2 ; R)

(3.41)

Using (3.41) in (3.39), and then exploiting (3.30), we obtain
(m + n)! δ(Za Zb Zc ) =
p0 ,q10 ,q20

= Zap,q1 ,q2 Zb

Zcp

X
d2 d2
R,R1 ,R2 R1 R2

00 ,q 00 ,q 00
1 2

X
R,R1 ,R2

(m+n)

× χR (Tp0

(m)

dR
R
χR (Za )χR
R1 ,R2 (Zb )χR1 ,R2 (Zc )
g(R1 , R2 ; R)2 R1 ,R2

g(R1 , R2 ; R)
χR (Tp(m+n) ) χR1 (Tq(m)
) χR2 (Tq(n)
)×
1
2
2 2
2
dR dR 1 dR 2
(m+n)

(n)

) χR1 (Tq0 ) χR2 (Tq0 ) χR (Tp00
1

2

(m)

(n)

) χR1 (Tq00 ) χR2 (Tq00 )
1

2

(3.42)

More generally, we can use (3.41) to compute the identity coefficient of an arbitrary large
products of central elements, δ(Za Zb · · · Zk ), just by using ordinary symmetric group
characters.

3.2

Maximal commuting subalgebra

In this section we describe the Maximal commuting subalgebra M(m, n) of A(m, n):
Z(m, n) ⊆ M(m, n) ⊆ A(m, n)

(3.43)

We often refer to M(m, n) as the Cartan subalgebra of A(m, n). M(m, n) is spanned by
elements of the form QR
R1 ,R2 ,i,i (no sum over i). For fixed R1 , R2 and R, the total number
of basis elements is g(R1 , R2 ; R), so that its dimension is
X
|M(m, n)| =
g(R1 , R2 ; R)
(3.44)
R1 `m,R2 `n
R`m+n

In Appendix A we derived the dimension formula
XX
|M(m, n)| =
Fp Fq Fp+q Sym(p + q)

(3.45)

p`m q`n

where p, q are partitions of m and n, Fp , Fq , Fp+q are combinatorial quantities
Q pi +qi dependent
only on the partitions p, q and p + q respectively, and Sym(p + q) = i i
(pi + qi )! is
a symmetry factor.
We now turn to the problem of constructing a basis for M(m, n). According to
the definition (3.9), to write the basis elements QR
R1 ,R2 ,i,i we first need to compute the
branching coefficients for the branching R → R1 ⊗ R2 . These quantities are in general
computationally hard to obtain 1 , and require a choice of a basis in Sm+n representations
1

see for example a discussion of the difficulty and the simplifications in a “distant corners approximation” in [30]

14

adapted to Sm × Sn . However, using the correspondence with matrix algebras given
by the Wedderburn-Artin decomposition, we can construct the Cartan by solving, in
each block, the following equations for (g(R1 , R2 ; R) − 1) linearly independent elements
QR
R1 ,R2 ,a ∈ A(m, n)
R
PRR1 ,R2 QR
R1 ,R2 ,a = QR1 ,R2 ,a

(3.46a)

PRR1 ,R2 , QR
R1 ,R2 ,a = 0

 R
QR1 ,R2 ,a , QR
R1 ,R2 ,b = 0

(3.46b)
(3.46c)

In the second equation, we are using the pairing defined in (3.4).

4

Star product for composite operators

In the previous sections we discussed the algebra A(m, n) and its centre Z(m, n). We
noted that central elements are special, as all their properties only depend on ordinary
symmetric group character theory. An example of this is eq. (3.42). In this section we will
take advantage of this fact to compute physically relevant quantities, in particular two
and three point functions of BPS operators in N = 4 SYM. To do so, we will first start by
discussing the one matrix sector in N = 4 SYM, reviewing the permutation description of
U (N ) matrix invariants which are Gauge Invariants Operators (GIOs) in the conformal
field theory. We will stress that for this case there is an underlying Z[C[Sn ]] algebra. The
one matrix problem will be used as a guide to extend to the two matrix problem, that we
treat in subsection 4.2. Here the underlying algebra will be A(m, n).

4.1

One matrix problem

Let us consider a matrix invariant constructed with n copies of the same matrix Z. Any
such invariant can be written in terms of a contraction

Oσ (Z) = tr Z ⊗n σ ,
σ ∈ Sn
(4.1)
subject to the equivalence relation
Oσ (Z) = Oγ −1 σγ (Z) ,

γ ∈ Sn

(4.2)

Polynomials in Z like the one in (4.1) can be multiplied together. Set σ1 ∈ Sn1 ,
σ2 ∈ Sn2 . By multiplying together Oσ1 (Z) and Oσ2 (Z) we get
Oσ1 (Z) Oσ2 (Z) = Oσ1 ◦σ2 (Z)

(4.3)

where σ1 ◦ σ2 ∈ Sn1 × Sn2 ⊂ Sn1 +n2 . Therefore for the usual product of matrix invariants,
σ1 ◦ σ2 lives in the symmetric group of degree n1 + n2 . We can define
M
C[S∞ ] =
C[Sn ]
(4.4)
n

which is closed under the circle product
◦ : C[S∞ ] ⊗ C[S∞ ] → C[S∞ ]
15

(4.5)

However, we can define another associative product, that we call star product, which closes
on the operators of fixed degree:
Oσ1 (Z) ∗ Oσ2 (Z) = Oσ1 σ2 (Z) ,

σ1,2 ∈ Sn

(4.6)

It is immediate to see how this product is different from the ordinary GIO multiplication
product (4.3): σ1 , σ2 and σ1 σ2 are all permutations of n elements, and the star product
is generally non-commutative. Let [σ] be the conjugacy class of σ. We now define a map
from the multi-trace GIOs to the class-algebra
Oσ (Z) →

X
1
Tσ
τ≡
size of [σ]
|Tσ |

(4.7)

τ ∈[σ]

This map is 1-1 at large N . Let us focus on this case. We can expand the product of
Ti , Tj ∈ Z[C[Sn ]] as
Ti Tj = Cijk Tk

(4.8)

Here the Cijk are the class algebra structure constants. By multiplying both sides above
by Tl and taking the coefficient of the identity we get
δ (Ti Tj Tl ) = Cijk δ (Tk Tl ) = δk,l Cijk |Tl | = Cijk |Tk |

(4.9)

Now we expand the star product Oσ1 (Z) ∗ Oσ2 (Z) as
Oσ1 (Z) ∗ Oσ2 (Z) =

X
p

X δ(Tσ1 Tσ2 Tσp )
|Tσp |
p
(Z)
=
C[σ
O
Oσp (Z)
σ
|Tσ1 | |Tσ2 | 1 ][σ2 ] p
|Tσ1 | |Tσ2 |
p

(4.10)

where the sum is over the conjugacy classes p of Sn . σp is a representative element of the
conjugacy class p. This equation will lead to a new expression for the two point functions
of GIOs built from Z, Z † in N = 4 SYM. First observe that setting Z to the identity
N × N matrix
Oσ (Z = 1N ) = N Cσ

(4.11)

where Cσ is the number of cycles in the permutation σ. Now consider taking the star
product of Oσ1 (Z), Oσ2 (Z) and then setting Z = 1N . We have, according to (4.10)
Oσ1 (Z) ∗ Oσ2 (Z)

=

=
Z=1N

X
1
δ(Tσ1 Tσ2 Tσp ) Oσp (1N )
|Tσ1 | |Tσ2 | p

X
X
1
1
δ(Tσ1 Tσ2 Tσp ) N Cσp =
δ(γTσ1 γ −1 Tσ2 Ω)
|Tσ1 | |Tσ2 | p
n!|Tσ1 | |Tσ2 | γ∈S

(4.12)

n

where we set Ω =
be [4]

P

p

Tσp N Cσp . On the other hand the free field correlator is known to

Oσ1 (Z)Oσ† 2 (Z) =

X
1
δ(γTσ1 γ −1 Tσ2 Ω)
|Tσ1 | |Tσ2 | γ∈S
n

16

(4.13)

so that
Oσ1 (Z)Oσ† 2 (Z) = n! Oσ1 (Z) ∗ Oσ2 (Z)

(4.14)
Z=1N

The two point function Oσ1 (Z)Oσ† 2 (Z) is therefore proportional to the star product
Oσ1 (Z) ∗ Oσ2 (Z) followed by the evaluation Z → 1N .
Similar considerations lead to the following expression for the extremal three point
function. In this case, we find that Oσ1 (Z)Oσ2 (Z)Oσ† 3 (Z) is proportional to the usual
product Oσ1 (Z) Oσ2 (Z), followed by the star product with Oσ3 (Z), followed by the evaluation Z → 1N . To see this, take σ1 ∈ Sn1 , σ2 ∈ Sn2 and consider
1
δ (Tσ1 ◦σ2 Tσ3 Ω)
(4.15)
|Tσ1 ◦σ2 | |Tσ3 |
Z=1N
P
where Tσ1 ◦σ2 ∈ Z[C[Sn1 +n2 ]], Tσ3 ∈ Z[C[Sn1 +n2 ]] and Ω = σ∈Sn +n σN Cσ . On the other
1
2
hand the correlator in N = 4 SYM [4] is
(Oσ1 (Z) Oσ2 (Z)) ∗ Oσ3 (Z)

Oσ1 (Z)Oσ2 (Z)Oσ† 3 (Z) =

X
γ∈Sn1 +n2

=


(n1 + n2 )!
δ γ(σ1 ◦ σ2 )γ −1 σ3−1 Ω =
δ (Tσ1 ◦σ2 Tσ3 Ω)
|Tσ1 ◦σ2 | |Tσ3 |
(4.16)

so that
Oσ1 (Z)Oσ2 (Z)Oσ† 3 (Z) = (n1 + n2 )! (Oσ1 ◦σ2 (Z)) ∗ Oσ3 (Z)

(4.17)
Z=1N

Given that these correlators are neatly expressed in terms of the star product, it would
be interesting to give an interpretation of the latter in the dual AdS5 × S5 side.
We will now write similar equations for the two matrix problem.

4.2

Two matrix problem

For the two matrix problem, the GIOs are polynomials in the X, Y matrices. Formally,
we can write them in terms of a permutation σ ∈ Sm+n as

Oσ (X, Y ) = Tr X ⊗m ⊗ Y ⊗n σ
(4.18)
As in the one matrix problem, there is an equivalence relation
γ ∈ Sm × Sn

Oσ (X, Y ) = Oγσγ −1 (X, Y ) ,

(4.19)

To each of these GIO Oσ we can associate a specific element Nσ of A(m, n) that we
call a necklace. We define a necklace Nσ as
X
1
Nσ =
γσγ −1
(4.20)
|AutSm ×Sn (σ)| γ∈S ×S
m

n

or equivalently as
Nσ =

X
τ ∈Orbit(σ, Sm ×Sn )

17

τ

(4.21)

where the sum is restricted to the permutations τ in the group orbit of σ under Sm × Sn .
We can think of the necklaces as the normalised version of the σ̄ elements defined in (3.3).
The set of necklaces form a basis for A(m, n). We associate a GIO to a necklace simply
by mapping
Oσ (X, Y ) →

1
Nσ
|Nσ |

(4.22)

For example, for the GIO corresponding to the permutation σ̃ = (1, 2, 4, 5)(3, 6) ∈ S6 :
Oσ̃ (X, Y ) = Tr(X 2 Y 2 )Tr(XY )
we associate, through the map (4.22), the A(3, 3) element
X
X
Nσ̃ =
γ σ̃γ −1 =
(a1 , a2 , b̄1 , b̄2 )(a3 , b̄3 )
S3 ×S3

(4.23)

(4.24)

a1 6=a2 6=a3 ∈{1,2,3}
b̄1 6=b̄2 6=b̄3 ∈{4,5,6}

Similarly, for the GIO specified by σ̃ = (1, 2, 3) ∈ S6
Oσ̃ (X, Y ) = Tr(X 2 Y )Tr(Y )3

(4.25)

we associate the A(2, 4) necklace
Nσ̃ =

X

(a1 , a2 , b̄1 )

(4.26)

a1 6=a2 ∈{1,2}
b̄1 ∈{3,4,5,6}

Notice that in the necklaces we do not explicitly write the single cycle permutations, but
rather we leave them implicit. In the last example, these single cycle permutations would
account for the multi-trace Tr(Y )3 component of Oσ̃ = Tr(X 2 Y )Tr(Y )3 .
From these examples it is clear how these necklaces are built by taking products of
cyclic objects, which in turn are constructed using two different types of beads. Such
cyclic objects are well studied in Polya theory. They can be related to the single cycle
permutations in Sm+n with equivalences generated by Sm × Sn . These equivalence classes
form the algebra A(m, n). We can imagine having blue beads corresponding to integers
[1, 2, ..m] and red beads corresponding to integers [m + 1, m + 2, ..., m + n]. Therefore,
we can pictorially depict the necklaces of examples (4.24) and (4.26) as in figure 1. The
same structure is present in the GIO Oσ corresponding to the necklace Nσ . In this case
the single-traces are the cyclic objects, and the role of the blue and red beads is played
by the X and Y type fields respectively.

Figure 1: Pictorial interpretations of the necklaces in the examples (4.24) and (4.26).

18

The map (4.22) is 1-1 at large N : as in the 1-matrix problem, we now focus on this
case. There is a natural product on the space of two matrix GIOs coming from multiplying
the multi-traces. For such a product, the degrees of the permutations add:
Oσ1 (X, Y )Oσ2 (X, Y ) = Oσ1 ◦σ2 (X, Y )

(4.27)

Here σ1 ∈ Sm1 +n1 is a representative of a class in A(m1 , n1 ) and σ2 ∈ Sm2 +n2 represents a
class in A(m2 , n2 ), while σ1 ◦ σ2 ∈ Sm1 +n1 × Sm2 +n2 ⊂ Sm1 +m2 +n1 +n2 represents a class in
A(m1 + m2 , n1 + n2 ). Continuing the analogy with (4.4), we can define
M
A(∞, ∞) =
A(m, n)
(4.28)
m,n

and for σ̄1 ∈ A(m1 , n1 ) and σ̄2 ∈ A(m2 , n2 ) we have
◦ : A(∞, ∞) ⊗ A(∞, ∞) → A(∞, ∞)

(4.29)

As in the one matrix case, there is however a second type of product of GIOs that we can
construct. The product on A(m, n) can in fact be used to define a closed and associative
star product on the space of the multi-trace operators with fixed numbers (m, n) of (X, Y ),
in the same fashion as (4.6):
Oσ̄1 (X, Y ) ∗ Oσ̄2 (X, Y ) = Oσ̄1 σ̄2 (X, Y ) ,

σ̄1,2 ∈ A(m, n)

(4.30)

Notice that here σ̄1 , σ̄2 and σ̄1 σ̄2 are all of the same degree, and that the star product
is non-commutative. We will use this star product to express the two point function of
GIOs built from X, Y .
Since the set of necklaces {Na } forms a basis for A(m, n), we can expand the product
Na Nb as
c
Na Nb = Ca,b
Nc

(4.31)

c
for some structure constants Ca,b
. Moreover, the necklaces are orthogonal in the metric
(3.4):

hNa , Nb i = δ(Na Nb ) = δa,b |Nb |

(4.32)

Here |Na | is the number of permutations in the necklace Na . We can write
c
δ(Na Nb Nc ) = |Nc |Ca,b

(4.33)

Now use the map (4.22) to map the two matrix invariants Oa (X, Y ) and Ob (X, Y ) to the
necklaces Na and Nb respectively. Then
Oa (X, Y ) ∗ Ob (X, Y ) =

X
c

c
Ca,b

X
1
|Nc |
Oc (X, Y ) =
δ(Na Nb Nc ) Oc (X, Y )
|Na | |Nb |
|N
|
|N
|
a
b
c
(4.34)

As for the one matrix problem case, by setting X = Y = 1N we get
Oa (X, Y ) ∗ Ob (X, Y )

=
X=Y =1N

19

1
δ(Na Nb Ω)
|Na | |Nb |

(4.35)

where Ω =
D

P

Oa (X, Y

σ∈Sm+n

σN Cσ . On the other hand the free field correlator [5, 17] is

)Ob† (X, Y

E

) =

X

δ(γaγ −1 b−1 Ω) =

γ∈Sm ×Sn

=

X
1
δ(γNa γ −1 Nb Ω)
|Na | |Nb | γ∈S ×S
m

n

m! n!
δ(Na Nb Ω)
|Na | |Nb |

(4.36)

Therefore, in analogy with (4.14) and (4.17), we can write the two point function as
D

Oa (X, Y

)Ob† (X, Y

E
) = m! n! Oa (X, Y ) ∗ Ob (X, Y )

(4.37)
X=Y =1N

and the extremal three point function as
Oa (X, Y )Ob (X, Y )Oc† (X, Y ) = (m1 + m2 )!(n1 + n2 )! Oa◦b (X, Y ) ∗ Oc (X, Y )

X=Y =1N

(4.38)
where a ∈ Sm1 +n1 , b ∈ Sm2 +n2 and c ∈ Sm1 +n1 +m2 +n2 . Finally, notice that the pairing (3.4)
is proportional to the planar correlator [42, 43, 44] of BPS operators: given Oa (X, Y ) and
Ob (X, Y ), we have
D
E
Oa (X, Y )Ob† (X, Y )
= m!n! ha, bi
(4.39)
planar

where the pairing on the RHS is the one in eq. (3.4).
Let us now focus on the centre of A(m, n). In section 3 we argued that the centre
(m+n)
(m)
(n)
(m+n)
(m)
is generated by {Tp
, Tq1 , Tq2 }. We remind the reader that {Tp
}, {Tq1 } and
(n)
{Tq2 } are the generators of the centres of C[Sm+n ], C[Sm ] and C[Sn ] respectively, and
that p, q1 and q2 are integer partitions of m + n, m and n. A GIO OTp(m+n) (X, Y ) can be
understood as a descendant of a single matrix 1/2 BPS state OTp(m+n) (X) under the U (2)
∂
internal symmetry that mixes the X and Y fields. In fact, given (D− )ij = Yki ∂X
j : we can
k
write

OTp(m+n) (X, Y ) ∼ (D− )n OTp(m+n) (X)

(4.40)

This means that central elements (and their corresponding matrix gauge invariants), de(m,n) (m) (n)
scribed in terms of the over-complete basis {Tp
Tq1 Tq2 }, are formed from composites
which employ both the usual product and the star product :
[Descendant Operators] ∗ {(X-Operators ) (Y -Operators )}
(m+n)

(4.41)
(m)

The descendant GIOs are associated to Tp
elements, X- and Y - GIOs to Tq1 and
(n)
Tq2 elements respectively. In terms of the permutations we are taking the product in
A(m, n) along with the circle product ◦ : A(m, 0) ⊗ A(0, n) → A(m, n).
Single-trace symmetrised traces are U (2) descendants of single-trace operators built
from a single matrix. In terms of the permutation language, they correspond to singlecycle permutations which are invariant under any reshuffling2 . On the other hand, U (2)
2

Further details of symmetrised traces in terms of an operation on the permutations in the Oσ (X, Y )
can be found in [44].

20

descendants of multi-trace operators built from one matrix form a subspace of the space
spanned by products of symmetrised single-trace states. In other words, not all products
of single-trace descendants are themselves descendants. One way to see this explicitly is
the following. Let STm,n be the space of symmetrised traces with m copies of X and n
copies of Y matrices. The generating function for the dimension Dim(STm,n ) is
X
1
=
Dim(STm,n ) xm y n
i
j
1−xy
m,n
i,j∈Ω
Y

(4.42)

where Ω = {0 ≤ i ≤ ∞} ∪ {0 ≤ j ≤ ∞} \ {i = j = 0}. Let STm+n be the space of
symmetrised traces with a total of m + n matrices, with any number of X or Y . We have
Dim(STm+n ) =

m+n
X

Dim(STi,m+n−i )

(4.43)

i=0

On the other hand, the total number of U (2) descendants obtained from a multi-trace
operator with m + n copies of X is
(m + n + 1) p(m + n)

(4.44)

p(m + n) is the number of partitions of m + n (the number of highest weight states), while
m + n + 1 is the number of descendants for a fixed highest weight. It can now be checked
that Dim(STm+n ) > p(m + n)(m + n + 1). This indeed proves our original claim.

4.3

Cartan subalgebra and the minimal set of charges

In [31], it was observed that, in the free limit, multi-matrix gauge theories have enhanced
symmetries including products of unitary groups. There are Noether charges for these
enhanced symmetries. Casimirs constructed from these charges have eigenvalues which
can distinguish all the labels R, R1 , R2 , i, j of restricted Schur operators. Because of SchurWeyl duality, these charges are also expressible in terms of permutations. Given the
definitions in this paper, this action of permutations amounts to the action of A(m, n) on
itself by the left or right regular representation. We can now characterize more precisely
what is a minimal set of charges which can measure all the labels. In section 3.2 we
introduced the Cartan subalgebra M(m, n), and gave a prescription to build a basis for
it. We need to find a subspace Cm,n of M(m, n) such that polynomials in some basis
elements ca ∈ Cm,n with coefficients taking values in the centre Z(m, n) span M(m, n).
In other words Cm,n contains a minimal set of generators for M(m, n) as a polynomial
algebra over Z(m, n). A minimal set of generators for Z(m, n), along with the basis
elements of the subspace Cm,n , provide a complete set of charges, which can measure all
min
the labels of the QR
(Z(m, n)) be the
R1 ,R2 ,i,j by left and right multiplication. Let N
minimal number of elements of Z(m, n) which generate Z(m, n) as a polynomial algebra.
min
Also, let NZ(m,n)
(M(m, n)) be the minimal number of elements of M(m, n) which generate
M(m, n) as a polynomial algebra over Z(m, n). Left multiplication by these generators
correspond to enhanced symmetry charges which measure the multiplicity index i of
restricted Schur operators. Right multiplication by the same generators correspond to
other enhanced symmetry charges which measure the multiplicity index j of restricted
Schur operators. Hence the minimal number of charges is
min
N min (Z(m, n)) + 2NZ(m,n)
(M(m, n))

21

(4.45)

An important open problem is to determine this function of (m, n) in general. This will
tell us how many bits of information completely specify all the operators in a multi-matrix
set-up.
The above discussion is complete for the case where m + n < N , which is adequate
for a treatment of the physics at all orders in the 1/N expansion. For finite N effects,
where we consider m + n > N , the charges given by the above still determine all the
multi-matrix invariants, but they are not a minimal set any more. The discussion can be
easily adapted to this case. Define
M
M
(m,
n)
=
(4.46)
Anull
Span{QR
N
R1 ,R2 ,i,j ; i, j}
R`m+n:c1 (R)>N R1 `m, R2 `n

The quotient
AN (m, n) = A(m, n)/Anull
N (m, n)

(4.47)

is a closed sub-algebra of blocks surviving the finite N cut. It has a centre ZN (m, n) and
a Cartan MN (m, n) which are simply related to Z(m, n) and M(m, n) by quotienting out
min
the parts belonging to Anull
(ZN (m, n)) be the number of generators in a
N (m, n). Let N
(MN (m, n))
minimal generating set for ZN (m, n) as a polynomial algebra. Let NZmin
N (m,n)
be the number of generators in a minimal generating set for MN (m, n) as a polynomial
algebra over ZN (m, n). The minimal number of charges needed is
N min (ZN (m, n)) + 2NZmin
(MN (m, n))
N (m,n)

(4.48)

We expect (4.45),(4.48) will have implications for information theoretic discussions of
AdS/CFT such as [45, 46].

5

Computation of the finite N correlator

In this section we will derive a finite N generating function for the two point function of
operators of the form
O = Tr(X m Y n )

(5.1)

in the free field metric. Operators like the one in (5.1) correspond to A(m, n) elements
1
(X) (Y )
T1̄,1 T[m] T[n]
m!n!
(X,Y )

(X)

(Y )

(X,Y )

(X)

(5.2)
(Y )

where T1̄,1 = T2
−T2 −T2 . Here T2
, T2 and T2 are the sum of transpositions
in Sm+n , Sm and Sn respectively. T1̄,1 can be understood as a joining operator, merging
the (1 · · · m) type cycles with the (m + 1 · · · m + n) type cycles.
The two point function (4.36) therefore reads, with O = Tr(X m Y n )
hOO† i =

=

1
2
m! n!2

X

X



(X) (Y )
(X) (Y )
δ γ T1̄,1 T[m] T[n] γ −1 T1̄,1 T[m] T[n] σ N Cσ

γ∈Sm ×Sn σ∈Sm+n


1 
(X) (Y )
(X) (Y )
δ T1̄,1 T[m] T[n] T1̄,1 T[m] T[n] Ω
m!n!
22

(5.3)

P
where we set Ω = σ∈Sm+n σN Cσ . This quantity can be computed using only ordinary
character theory. Using eq. (3.42) and using the shorthand notation g = g(R1 , R2 ; R) we
write


2
X X
1
dR
(X) (Y )
R
hOO† i =
χ
T
T
T
χR
1̄,1 [m] [n]
R1 ,R2
R1 ,R2 (Ω)
(m + n)!m!n! R `m R`m+n d2R1 d2R2 g 2
1
R2 `n

(5.4)
(X,Y )

We now expand T1̄,1 = T2
(X)

(Y )

(X)

− T2
(X,Y )

T1̄,1 T[m] T[n] = T2

(Y )

− T2

(X)

so that

(Y )

(X)

T[m] T[n] − T2

(X)

(Y )

(X)

(Y )

(Y )

T[m] T[n] − T[m] T2 T[n]

(5.5)

g dR1 dR2
(n + m)! DimN (R)
dR

(5.6)

We also have (see e.g. [22])



R

χR
R1 ,R2 (Ω) = χR1 ,R2

X

σN Cσ  =

σ∈Sm+n

Eq. (5.4) simplifies then to
hOO† i =

1 X
m!n! R `m
1
R2 `n

X
R`m+n

2


1
(X) (Y )
T
DimN (R) χR
T
T
1̄,1 [m] [n]
R1 ,R2
dR1 dR2 g

On the other hand, as shown in Appendix C


(X) (Y )
=
T
T
T
χR
1̄,1
R1 ,R2
[m] [n]


(X,Y )
)

c
+c
R
R
 (−1) 1 2 g (m − 1)!(n − 1)! χR (Td2
−
R
=


0 otherwise

(X)

χR1 (T2
dR1

)

(Y )

−

χR2 (T2
dR2

)

(5.7)


;

R1 , R2 hooks

Here cRi is the number of boxes in the first column of the Young diagram associated
with the representation Ri . This expression restricts the sums over representations R1 `
m, R2 ` n in (5.7) to a sum over hook representations h1 ` m, h2 ` n.
We now need an equation for g(h1 , h2 ; R), with h1 and h2 hook representations of Sm
and Sn respectively. We specify any representation R by the sequence of pairs of integers
R = ((a1 , b1 ), (a2 , b2 ), ...(ad , bd )). In a Young diagram interpretation, aj (1 ≤ j ≤ d) is
the number of boxes to the right of the j-th diagonal box, and bj is the number of boxes
below the j-th diagonal box. We refer to d as the ‘depth’ of the representation R. Let us
write h1 = (k1 , l1 ), h2 = (k2 , l2 ) and R = ((a1 , b1 ), (a2 , b2 )). In Appendix B we show that
g(h1 , h2 ; R) = δk1 +k2 ,a1 δl1 +l2 +1,b1 δ−1,a2 δ0,b2 + δk1 +k2 +1,a1 δl1 +l2 ,b1 δ0,a2 δ−1,b2
+

1
X

min(k1 −¯
1 ¯2 ,k2 −1 2 )

min(l1 −¯
1 ¯2 ,l2 −1 2 )

X

X

1 ,2 =0

i=1 ¯2

j=¯
1 2

δk1 +k2 −i+¯1 2 ,a1 δl1 +l2 −j+1 ¯2 ,b1 δi−1 ¯2 ,a2 δj−¯1 2 ,b2
(5.8)

23

where ¯1,2 = 1 − 1,2 . Using this identity, in Appendix C we derive the formula
hTr(X m Y n )Tr(X m Y n )† i
=

m
X

n
X

k1 ,l1 =0 k2 ,l2 =0

n+m
X

g δ(k1 + l1 − m) δ(k2 + l2 − n) F (a1 , b1 , a2 , b2 , k1 , l1 , k2 , l2 ) (5.9)

a1 ,b1 =0
a2 ,b2 =0

where we defined the function
k1 !k2 !l1 !l2 ! (a1 − a2 )(b1 − b2 )
4(a1 + b2 + 1)(a2 + b1 + 1)(k1 + l1 + 1)(k2 + l2 + 1)





a1 + b 1
a2 + b 2
N + a1
N + a2
×
×
b1
b2
a1 + b 1 + 1 a2 + b 2 + 1

F (a1 , b1 , a2 , b2 , k1 , l1 , k2 , l2 ) =

× ((a1 + b1 + 1)(a1 − b1 ) + (a2 + b2 + 1)(a2 − b2 )+
− (k1 + l1 + 1)(k1 − l1 ) − (k2 + l2 + 1)(k2 − l2 ))2

(5.10)

In [36] a closed form for the two point function has been given by using a different
approach based on Young-Yamanouchi symbols. We have checked agreement of (5.9)
with that closed form for up to n = m = 10. It is an interesting exercise to simplify
(5.9) into the closed form obtained in [36]. It will also be interesting to apply the present
franework to obtain formulae analogous to (5.9) for more general GIOs corresponding to
central elements of A(m, n).
In this section we have shown how to calculate a particular two point function of
a central operator, without explicitly constructing projectors. The result rather follows
from knowing how central operators of interest are generated via the star product of pure
X gauge invariants, pure Y gauge invariants and descendants of half-BPS operators.

5.1

Coloured Ribbon graphs

The correlator computations above can be expressed in terms of ribbon graphs, equivalently the usual double-line graphs of large N expansions, but with edges coming in two
colors, as explained for example in [47]. The graphs can be organised by the minimum
genus of the surface they can be embedded in and these graphs of a given genus contribute to a fixed power of N . For small m, n, we have checked with GAP that directly
computing the permutation sums for a given genus agree with the analytic result (5.9) we
have derived.

6

Conclusions and future directions

In this paper, we initiated a systematic study of permutation centralizer algebras, in
connection with gauge invariant operators. We focused our attention on the algebras
A(m, n) which are related to restricted Schur operators studied in the context of giant
gravitons in AdS/CFT. Other closely related algebras are related to the Brauer basis for
multi-matrix invariants, the covariant basis and to tensor models.
While many of the key formulae we have used were already understood in the literature
on giant gravitons, we have emphasized the intrinsic structure of A(m, n) as an associative algebra with a non-degenerate pairing. This means that it has a Wedderburn-Artin
24

decomposition, which gives a basis for the algebra in terms of matrix-like linear combinations. The construction of these matrix-units in terms of representation theory data from
Sm+n , Sm , Sn has already been extensively used in the context of giant gravitons, although
the link to the Wedderburn-Artin decomposition has not been made explicit before. In
addition to explicating this link, the new emphasis in this paper has been on the structure
of the centre Z(m, n) and the maximally commuting sub-algebra M(m, n).
We have used the structure of M(m, n) as a polynomial algebra over Z(m, n) to characterize the minimal number of charges needed to identify any 2-matrix gauge invariant
(section 4.3). It will be interesting to generalize this discussion to gauge invariants for
more general gauge groups.
Two key structural facts about A(m, n) have played a role in the computation of
correlators in Section 5. The first is that (xm ) ∗ (y n ) = (xm y n ) and the second is that
(xm y n ) is part of Z(m, n). The non-degenerate pairing on A(m, n), when restricted to
elements in the centre, can be expressed in terms of characters of Sn , Sm , Sn+m without
requiring more detailed representation theory data such as matrix elements and branching
coefficients. These are in general computationally hard to calculate, although there has
been progress in the context of “perturbations of half-BPS giants”. This makes it very
interesting to understand the structure of the centre A(m, n). A special case is Z[C[Sn ]],
which is the algebra of class sums in Sn .

6.1

Structure of the centre

A number of questions about A(n), A(m, n) and the centre A(m, n) can be explored
experimentally, with the help of group theory software, notably GAP. In particular, since
Z(m, n) is generated by the centre of Sm , the centre of Sn and that of Sn+m it is a useful
first step to know about these centres.
Since Sn is generated by transpositions, one might naively expect that the sum of
permutations T2 will generate A(n). This is actually not true. We know that T2 obeys a
relation of degree p(n)

Y
χR (T2 )
=0
(6.1)
T2 −
d
R
R`n
If this is the only relation, then we know that T2 alone generates Z[C[Sn ]]. However
simpler relations occur when there are coincidences in the normalized characters, e.g.
two different irreps have the same normalized character. In fact the the failure of T2
to generate centre is always correctly predicted by the degeneracies of the normalized
characters. If we take
Y0
χR (T2 )
)=0
(6.2)
(T2 −
dR
R
where the product is taken over a maximal set of irreps with distinct normalized characters, we are getting an element in C[Sn ] which vanishes in all irreps. It is a central element,
so the matrix elements in any irrep are proportional to the identity. We conclude that the
above element vanishes. Given that the Peter-Weyl theorem gives an isomorphism between C[Sn ] and matrix elements of irreps, it follows that something which has vanishing
matrix elements in all irreps should be identically zero.
Even for large n, it is possible to check that the centre of C[Sn ] is generated by a small
number of Tp ’s. Using GAP we tested that T[2,1n−2 ] and T[3,1n−3 ] are enough to generate
25

the centre for C[Sn ] up to n = 14. The procedure we used to perform these checks is
the following. We know that the set of projectors {PR }, with R integer partition of n,
generate the centre of Sn . We can compute the overlap of PR with the k-th power of Tp ,
that we simply write as Tpk :
Tpk , PR

=

δ(Tpa


= dR

1 X
1 X
PR ) =
χS (Tpk )χS (PR ) =
dS
n! S`n
n! S`n
χR (Tp )
dR



χS (Tp )
dS

k
χS (PR )

k
(6.3)

Similarly, we can derive
Tpk Tql , PR


= dR

χR (Tp )
dR

k 

χR (Tq )
dR

l
(6.4)

Now we construct the AB ×p(n) matrix M (A, B), whose matrix elements are the overlaps
(6.4):
k 
l

χR (Tq )
χR (Tp )
(6.5)
M (A, B)|(k,l),R = dR
dR
dR
with 0 ≤ k < A and 0 ≤ l < B. By computing the rank of this matrix we obtain the
number of independent central elements in C[Sn ] that are obtained by taking at most
A − 1 powers of Tp and B − 1 powers of Tq . This method can be easily generalised to
obtain the number of central elements generated by the string of operators Tpk11 Tpk22 · · · TpkNN .
These studies on the centre of C[Sn ] inspire a similar analysis for centre of A(m, n).
The task is to find a minimal set of generators for Z(m, n) as a polynomial algebra.
The importance of this problem is discussed in section 4.3. Concretely, we would like to
determine N min (Z(m, n)). There are many approaches one can take in this case, which
would be interesting to investigate in the future. For example, using GAP we checked
(m+n)
(m+n)
that low powers of the sum of two- and three-cycles permutations, T2
and T3
,
together with the generators of the centres of C[Sm ] and C[Sn ], generate the whole centre
Z(m, n). We leave a more systematic discussion of this problem for future work.

6.2

Construction of quarter-BPS operators beyond zero coupling and the structure constants of A(m, n).

The centre of C[Sm+n ] is denoted by Z[C[Sm+n ]]. Z[C[Sm+n ]] is a commutative subalgebra of A(m, n). The A(m, n) algebra is a module over Z[C[Sm+n ]]. We can write
k
Tp Ni = C̃pj
Nk

(6.6)

for some coefficients C̃. The Tp are themselves linear combinations of necklaces:
Tp = Tpi Ni

(6.7)

l
Nl
Tp Ni = Tpj Nj Ni = Tpj Cjk

(6.8)

Hence

26

Another subspace in A(m, n) is the subspace of symmetrised traces. A symmetrised
trace Sv can be parametrised by a vector partition v of (m, n). We can expand Sv on the
basis of necklaces {Nk } as
Sv = Svk Nk

(6.9)

Symmetrised traces and their products are quarter-BPS at weak coupling in the large N
limit. One can get the complete set of 1/N corrected BPS states at large N by acting
on Sv with Ω−1 which belongs to Z[C[Sm+n ]] ⊗ C(1/N ) [15, 16, 43, 22]. The coefficients
of Tp are easily computable. The expansion of Tp in terms of necklaces is also easily
computable. The non-trivial part of the calculation is the Cijk of the necklace algebra
A(m, n). For any symmetrised trace Sv , the corrected operator is
−1
j
−1 j k
−1 j l k
Ω−1
k Sv = Ωp Tp Sv Nj = Ωp Sv C̃pj Nk = Ωp Sv Tp Clj Nk

6.2.1

(6.10)

Central quarter BPS sector

A subspace of symmetrised trace elements is central. The symmetrised trace elements give
a subspace of A(m, n) and the central elements form another subspace. The intersection is
the space of central symmetrised traces. The dimension of this subspace can be computed
for small m, n using GAP. Suppose S C is an element in this subspace. Then elements
Ω−1 S C in A(m, n) are very interesting. They are quarter-BPS beyond zero coupling and
they are central, so computations of their correlators have the simplicity of the centre.
The computations can be done using knowledge of the characters of Sm , Sn , Sm+n , without
knowing branching coefficients. From AdS/CFT this central quarter BPS sector should
have a dual in the space-time theory, e.g. some sub-class of states in the tensor product
of super-graviton states. An interesting question is to compute their correlators in spacetime and verify the matching with the gauge theory computations.

6.3

Non-commutative geometry and topological field theory

Studies in non-commutative geometry in string theory suggest that open strings can be
associated to non-commutative algebras and the centre is related to closed strings [48].
If we apply this thinking to A(m, n) and Z(m, n), how do we interpret these emergent open and closed strings? The traditional view is that Yang-Mills theory is the
open-string picture in AdS/CFT with the closed string picture given by the AdS description, so this is an intriguing question. Non-commutative algebras and their centre
have also been discussed in non-commutative geometry in [49]. The study of the pair
{A(m, n), Z(m, n)} should form an interesting example of this discussion. Additionally
we have the Cartan M(m, n) here, with physical relevance in distinguishing the multiplicity labels. So a more complete picture of strings and non-commutative geometry for the
triple {A(m, n), M(m, n), Z(m, n)} looks desirable. Given that the infinite direct sum
A(∞, ∞) comes up in connection with matrix invariants, it would also be interesting to
study the triple {A(∞, ∞), M(∞, ∞), Z(∞, ∞)} from this point of view. Some relevant
work in this direction is in [26].

27

6.4

Other examples of permutation centralizer algebras and
correlators

Based on our study of A(m, n), we outline some properties of the other examples of permutation centralizer algebras given in section 2 and sketch the connection to correlators.
We leave a more detailed development for the future.
Consider BN (m, n), which is the subspace of the Brauer algebra BN (m, n) invariant
under C[Sm × Sn ]. This is Example 3 in Section 2. Brauer algebras were used to construct gauge invariant operators in [14] from tensor products of a complex matrix and its
conjugate. For an element b in the walled Brauer algebra BN (m, n), we use

trm,n Z ⊗m ⊗ Z̄ ⊗n b
(6.11)
where the trace is taken in V ⊗m ⊗ V̄ ⊗n , a tensor product of fundamentals and antifundamentals of U (N ). We focus here on the case m + n ≤ N . The number of gauge
invariant operators is
X
γ
)2
(6.12)
(Mα,β
γ,α,β
γ
where γ labels an irrep of BN (m, n), while α, β are irreps of Sm and Sn respectively. Mα,β
is a multiplicity with which (α, β) appears in the reduction of γ from BN (m, n) to its
C[Sm × Sn ] sub-algebra. The sum of squared dimensions in (6.12) is the dimension of
the algebra BN (m, n). This is a non-commutative algebra. The dimension of its centre
γ
is the number of triples (γ, α, β) for which Mα,β
is non-vanishing. There is a maximally
commuting sub-algebra of dimension equal to the sum
X γ
Mα,β
(6.13)
γ,α,β

This follows since the (γ, α, β) give a Wedderburn-Artin decomposition of BN (m, n). A
tractable sector of correlators should be given by the centre of BN (m, n) and more detailed
study of the structure of this centre will be useful.
The next algebra of interest is the sub-algebra K(n) of C[Sn ]×C[Sn ] which is invariant
under conjugation by Diag(C[Sn ]). Let us denote this as Adiag (n, n). We can generate
elements in this algebra by summing over the elements of the sub-group
X
σ1 ⊗ σ2 →
γσ1 γ −1 ⊗ γσ2 γ −1
(6.14)
γ∈Sn

The dimension of this algebra is
X

C(R, S, T )2

(6.15)

R,S,T

where C(R, S, T ) is the Kronecker coefficient, i.e. the number of times the irrep T of Sn
appears in the tensor product R ⊗ S. The dimension of the centre is the number of triples
(R, S, T ) for which the C(R, S, T ) is non-zero. A maximal commuting sub-algebra has
dimension
X
C(R, S, T )
(6.16)
R,S,T

28

These properties follow from the fact the Wedderburn-Artin decomposition of the algebra
K(n) has blocks labelled by triples (R, S, T ) with non-vanishing C(R, S, T ). An explicit
formula for this decomposition is
X X
1 R,S,T,τ2
QR,S,T
SiR,S,T,τ
Sj1 ,j2 ,i3 DiR1 j1 (σ1 )DiS2 j2 (σ2 ) σ1 ⊗ σ2
(6.17)
τ1 ,τ2 =
1 ,i2 ,i3
σ1 ,σ2 i1 ,i2 ,i3 ,j1 ,j2

The D’s are representation matrices for Sn irreps. The S’s are Clebsch-Gordan coefficients.
One verifies, using equivariance properties of the Clebsch’s that these are invariant under
conjugation by the diagonal Sn .
There is another definition of K(n) which is more symmetric in (R, S, T ). C(R, S, T )
is also the multiplicity of invariants of the diagonal Sn acting on R ⊗ S ⊗ T . K(n) can be
defined as the subalgebra of C[Sn ] ⊗ C[Sn ] ⊗ C[Sn ] which is invariant under left action by
the diagonal C[Sn ] and right action by the diagonal C[Sn ]. These invariant elements can
again be constructed by averaging
X
(γ1 σ1 γ2 , γ1 σ2 γ2 , γ1 σ3 γ2 )
(6.18)
γ1 ,γ2

A representation basis is given by
1 R,S,T,τ2
σ1 ⊗ σ2 ⊗ σ3 DiR1 ,j1 (σ1 )DiS2 ,j2 (σ2 )DiT3 ,j3 (σ3 )SiR,S,T,τ
Sj1 ,j2 ,j3
1 ,i2 ,i3

(6.19)

labelled by R, S, T, τ1 , τ2 .
These triples of permutations (σ1 , σ2 , σ3 ), with equivalences given by left and right
diagonal action have appeared in the enumeration invariants for tensor models built from
3-index tensors [50]. The simplification from a description in terms of permutation triples
to one in terms of permutation pairs was also described there, which lead to a connection
between 3-index tensor invariants and Belyi maps. By analogy with the discussion in
this paper, we expect that the centre of K(n) will lead to a class of simpler correlators in
tensor models. The discussion of A(∞, ∞) will analogously lead to
K(∞) =

∞
M

K(n)

(6.20)

n=0

This space will have two products: one related to the algebra structure of K(n) and one
related to the multiplication of tensor invariants. Somewhat related algebraic structures
appear in [51] and it would be useful to better understand these relations. As a last
remark, consider the Kronecker multiplicities C(R, R, T ), i.e. in the special case where
R = S. These have also appeared in the construction of gauge-invariant multi-matrix
operators in a basis which is covariant under the global symmetries [15, 16]. The structure
of K(n) can thus also be expected to have implications for multi-matrix correlators in the
covariant basis.

Acknowledgements
We thank David Berenstein, Robert de Mello Koch and Edward Hughes for useful
discussions, and Robert de Mello Koch for comments on an earlier draft. SR is supported
29

by STFC consolidated grant ST/L000415/1 “String Theory, Gauge Theory & Duality.”
SR thanks the Simons Summer workshop 2015 for hospitality while part of this work was
done. PM is supported by a Queen Mary University of London studentship.

Analytic formula for the dimension of M(m, n)

A

In this section we derive a formula for the dimension of M(m, n). This dimension is equal
to the sum of Littlewood-Richardson coefficients
X
X
Dim M(m, n) =
g(R1 , R2 , R)
(A.1)
R1 `m,R2 `n R`m+n

The sum of squares of the Littlewood-Richardson coefficients is the dimension of A(m, n)
and has a simple 2-variable generating function. It is natural to ask if we can write a nice
generating function for the dimension of M(m, n). While we have not been able to derive
something of comparable simplicity, we will derive two interesting expressions (A.11) and
(A.28) in terms of multi-variable polynomials.
Let Tp denote a conjugacy class of permutations with cycle structure determined by a
vector (p1 , p2 , · · · ), i.e. permutations with pi cycles of length i. Let now σp be an element
in Tp . For σp ∈ Tp , it is known that [52]


X
Y
tpi i
(A.2)
χR (σp ) =
Coeff fi (ti ),
pi !
i
R
where
fi (ti ) = e

it2
(1−(−1)i )
ti + 2i
2

(A.3)

We can define
F (t1 , t2 , · · · ) =

Y

fi (ti )

(A.4)

i

and write
X

χR (σp ) = Coeff

F (t1 , t2 , · · · ),

Y tpi
i

i

R

!

pi !

(A.5)

It is also useful to define
 
ti
f˜i (ti ) = fi
i



t2 t3
ti
e
F (t1 , t2 , · · · ) = F t1 , , · · · = F { }
2 3
i
∞
Y ti Y
it2
i
=
ei
e 2i
i=1

i:odd

30

(A.6)

We can write the LR coefficients in terms of Tp ’s as
1 X X
χR1 (σ1 )χR2 (σ2 )χR (σ1 ◦ σ2 )
m!n! σ ∈S σ ∈S
n
m 2
1
XX
Y
1
=
χR1 (Tp )χR2 (Tq )χR (Tp ◦ Tq )
ipi +qi pi !qi !
i
p`m q`n

g(R1 , R2 , R) =

(A.7)
This uses the fact that thePnumber of permutations in the class Tp is n!/
use the above formula for R χR (Tp ), to obtain
X
g(R1 , R2 , R)

Q

i

ipi pi !. Now

R1 ,R2 ,R

=
=
=

XXY
p`m q`n
X
X
p`m q`n
X
X

Coeff(f˜i (si ), tpi i ) Coeff (f˜i (ti ), tqi i ) Coeff (fi (ui ), upi i +qi )(pi + qi )!

i

Coeff (Fe(~s)Fe(~t)F (~u),
Coeff (Fe(~s)Fe(~t)Fe(~u),

Y
i
Y

spi i tqi i upi i +qi (pi + qi )!−1 )
spi i tqi i upi i +qi )ipi +qi (pi + qi )!

(A.8)

i

p`m q`n

It is useful to make the substitutions si → si zi , ti → ti zi , ui → z̄i and to introduce a
pairing 3
hzjk , z̄il i = δij δkl k! ik

(A.9)

F(zi , s) = Fe(ti → si zi )

(A.10)

With these substitutions define

Then we can write
Dim (M(m, n)) = h Coeff (F(zi , s)F(zi , t)F(zi , u = 1), sm tn )i

(A.11)

This has been checked for very simple cases, e.g. up to (m, n) = (3, 3)

A.1

Multi-variable polynomials

It is useful to isolate the multi-variable polynomials in the zi variables at each order in
the s, t variables. Let us introduce the quantities
 2i 2 
Y
s zi
A(~z, s) =
exp
2i
i


Y
si zi
B(~z, s) =
exp
(A.12)
i
i=1,3,..
It follows from previous formulae (A.6) and (A.10) that
F(~z, s) = A(~z, s)B(~z, s)

(A.13)

Alternatively we can think about expectation values in a Fock space with zi → ai , z̄i → a†i . This
would allow us to write the subsequent formulae in terms of quantities in a 2D field theory. This
perspective could be fruitful, but we will leave its exploration for the future
3

31

Introducing polynomials Fm (~z) for each order in s we can rewrite the latter quantity as
X
F(~z, s) =
Fm (~z)sm
(A.14)
m=0

We will now write formulae for the coefficients of sm in A and B. For A(~z, s) we derive
A(~z, s) =

∞
X

∞
X

A2m (~z)s2m =

m=0

∞
Y
s2ipi z 2ipi

p1 ,p2 ,···=0 i=1

i
p
i
(2i) p

i!

(A.15)

so that
A2m (~z) =

zi2ipi
(2i)pi pi !
p`m

X

(A.16)

We can also define Am (~z) to be zero for odd m and equal to the above for the even values.
It is useful to define the coefficients of z12p1 z24p2 . . . zi2ipi in the A(~z, s = 1) as
A[p] = A[p1 ,p2 ··· ] =

Y
i

1
pi !(2i)pi

(A.17)

zi2ipi

(A.18)

so that we may write
A2m =

X

A[p]

∞
Y
i=1

p`m

Similarly, for B(~z, s) we obtain
∞
Y

s(2i+1) z2i+1
B(~z, s) =
exp
(2i + 1)
i=0



(A.19)

and
Bm (~z) =

X

ziipi
(i)pi pi !
odd

Y

{p1 ,p3 ··· }`m i

(A.20)

Therefore it is natural to define
Y 1
B[p1 ,p3 ,··· ] =
ipi pi !
Xi
Y ip
Bm (~z) =
B[p1 ,p3 ,··· ]
zi i
p`m

(A.21)

i odd

Going back to (A.14) we get, using the formulae just derived
Fm (~z) =

m
X

m

Ak (~z)Bm−k (~z) =

k=0

b2c
X

A2k (~z)Bm−2k (~z)

k=0

m

=

b2c
X
X X
k=0 r`k

A[r] B[q]

Y
i

q`m−2k
q odd

32

i(2ri +qi )

zi

(A.22)

Grouping terms with the same power of zi we obtain
X
Y ip
F(~z, s = 1) =
F[p1 ,p2 ...]
zi i

(A.23)

i

[p1 ,p2 ...]

with
F[p] =

X

X

A[r1 ,r2 ··· ] B[q1 ,q2 ,··· ]

Y

δ(pi , 2ri )

i even

[r1 ,r2 ,...] [q1 ,q2 ··· ]

Y

δ(pi , 2ri + qi ) (A.24)

i odd

Note that the function F(~z, s) is closely related to the generating function for the cycle
indices of Sn which is
#
"∞
X ti zi
Z(~z, t) = exp
i
i=1
1/2
Ã(~z, s) = Z(zi → zi2 , s → s2 )
B̃(~z, s) = (Z(z2i+1 → z2i+1 , z2i → 0))1/2
(A.25)
We can work with the same function if we change the pairing. With the pairing
E
D
kj
ki
= δi,j δki ,kj ki ! iki
(A.26)
zi , zj
we can write the above formulae as
Dim(M(m, n)) = hFm (~z)Fn (~z), Fm+n (~z)i

(A.27)

or, equivalently,
Dim(M(m, n)) =

XX

Fp1 ,p2 ··· Fq1 ,q2 ,··· Fp1 +p2 ,q1 +q2 ,···

=

ipi +qi (pi + qi )!

i

p`m q`n

XX

Y

Fp Fq Fp+q Sym(p + q)

(A.28)

p`m q`n

This is eq. (3.45).

B

LR rule for hook representations

Here we derive the LR decomposition rule for the tensor product of two hook representations. Let us consider three representations R, R1 and R2 of Sm+n , Sm and Sn
respectively. The LR coefficient g(R1 , R2 ; R) gives the multiplicity with which the representation R1 ⊗ R2 appear in the representation R upon its restriction to Sm × Sn . There
is a systematic procedure to obtain such coefficients [40], that we now briefly review. We
take the Young diagrams corresponding to R1 and R2 , and we start by decorating the
latter as follows. We write ‘1’ in all the boxes of the first row, ‘2’ in all the boxes of the
second row and so on in a similar fashion until the last row. Then we proceed to move
all the ‘1’ boxes from R2 to R1 , ensuring that that we produce legal Young diagrams and
no two copies of ‘1’ appear in the same column. We then move the ‘2’ boxes following
the same rules, and so on. In doing so, we also require a reading condition. At any step,
reading from right lo left along the first row and then subsequent rows, the number of ‘1’
33

boxes must be greater or equal to the number of ‘2’ boxes. Similarly, the number of ‘2’
boxes must be greater or equal to the number of ‘3’ boxes, and so on.
At the end of this procedure we are left with a collection of Young diagrams, made
with m + n boxes. If two or more of the resulting diagrams are identical (that is, the
not only match in shape but also in the numbering of their boxes), we only retain one of
them. Otherwise, if k diagrams R appear with the same shape but different numbering,
we can say that g(R1 , R2 ; R) = k. These will be the prescriptions that we will follow to
derive our LR formula.
We specify any representation R by the sequence of pairs of integers R = ((a1 , b1 ), (a2 , b2 ),
...(ad , bd )). In a Young diagram interpretation, aj (1 ≤ j ≤ d) is the number of boxes
to the right of the j-th diagonal box, and bj is the number of boxes below the j-th diagonal box. We refer to d as the ‘depth’ of the representation R. Hooks therefore are
representations of depth 1. Schematically, in this appendix we will obtain the RHS of
M
(k1 , l1 ) ⊗ (k2 , l2 ) =
((a1 , b1 ), (a2 , b2 ))
(B.1)
In our derivation we imagine to keep the first hook fixed, and to add to it boxes coming
from the second diagram. In doing so we are careful to follow the LR prescription. The
boxes of the second diagram are decorated by a ‘1’ or a ‘v’, depending whether they
come from the first row of the diagram or not. The tensor product (k1 , l1 ) ⊗ (k2 , l2 ) will
decompose into a direct sum of a varying number of depth 2 representation and precisely
two hooks (regardless of the actual value of k1,2 , l1,2 ). These hooks are
Hook 1: (k1 + k2 + 1, l1 + l2 )
Hook 2: (k1 + k2 , l1 + l2 + 1)

(B.2)

Notice that we can rewrite them using the notation we use for the depth two diagram as
Hook 1: ((k1 + k2 + 1, l1 + l2 ), (0, −1))

(B.3)

Hook 2: ((k1 + k2 , l1 + l2 + 1), (−1, 0))

(B.4)

This notation will be helpful at a later stage.
We now turn to the depth two representations. We proceed systematically, grouping
them into four categories according to the two yes/no questions:
1) Is there a 1 in the first column of the resulting diagram?
2) Is there a v in the first row of the inner hook of the resulting diagram?
We now analyse these four possibilities.

B.1

(Y,Y) case

The diagrams in this class are of the form

34

Figure 2: (Y,Y) case

They can be described by the expression
(Y, Y ) :

((k1 + k2 − i, l1 + l2 − j), (i, j))

(B.5)

where i and j are constrained by the boundaries
0 ≤ i ≤ min(k1 , k2 − 1)
0 ≤ j ≤ min(l1 , l2 − 1)

(B.6)

The upper bound on i is min(k1 , k2 − 1) because, if k1 ≥ k2 , we cannot remove all the
k2 1 type boxes from the first row. This has to be avoided since by construction the
rightmost box in the second row has to be a v type box. A diagram with no 1 type
boxes on the first row and a v type box at the end of the second row would violate the
LR reading condition.

B.2

(Y,N) case

The diagrams in this class are of the form

Figure 3: (Y,N) case

35

They can be described by the expression
(Y, N ) :

((k1 + k2 − i, l1 + l2 − j + 1), (i − 1, j))

(B.7)

with the boundaries
1 ≤ i ≤ min(k1 , k2 )
0 ≤ j ≤ min(l1 , l2 )

B.3

(B.8)

(N,N) case

The depth two diagrams in this class are of the form

Figure 4: (N,N) case

They can be described by the expression
(N, N ) :

((k1 + k2 − i, l1 + l2 − j), (i, j))

(B.9)

with the boundaries
0 ≤ i ≤ min(k1 − 1, k2 )
0 ≤ j ≤ min(l1 − 1, l2 )

B.4

(N,Y) case

The diagrams in this class are of the form

36

(B.10)

Figure 5: (N,Y) case

These can be described by the equation
(N, Y ) :

((k1 + k2 − i + 1, l1 + l2 − j), (i, j − 1))

(B.11)

The boundary for i is
0 ≤ i ≤ min(k1 , k2 )

(B.12)

The upper bound is k2 and not k2 + 1 because we cannot remove all the 1 from the first
row, as the rightmost box in the second row has to be a v type box. In this way, we are
enforcing the LR reading condition. On the other hand, the boundary for j is
1 ≤ j ≤ min(l1 , l2 )

(B.13)

The lower bound is a 1 as by construction there has to be a v box in the first row of the
inner hook.

B.5

A summary

These four cases comprise all possible valid depth two diagrams. Summarising our result,
we have
• (Y, Y ) case:

((k1 + k2 − i, l1 + l2 − j), (i, j))
0 ≤ i ≤ min(k1 , k2 − 1)
0 ≤ j ≤ min(l1 , l2 − 1)

• (Y, N ) case:

(B.14)

((k1 + k2 − i, l1 + l2 − j + 1), (i − 1, j))
1 ≤ i ≤ min(k1 , k2 )
0 ≤ j ≤ min(l1 , l2 )

37

(B.15)

• (N, N ) case:

((k1 + k2 − i, l1 + l2 − j), (i, j))
0 ≤ i ≤ min(k1 − 1, k2 )
0 ≤ j ≤ min(l1 − 1, l2 )

• (N, Y ) case:

(B.16)

((k1 + k2 − i + 1, l1 + l2 − j), (i, j − 1))
0 ≤ i ≤ min(k1 , k2 )
1 ≤ j ≤ min(l1 , l2 )

We now introduce the boolean parameters

 0 If the answer to the first question is no
1 =
 1 If the answer to the first question is yes
and


 0
2 =
 1

If the answer to the second question is no

(B.17)

(B.18)

(B.19)

If the answer to the second question is yes

With this notation we can compactly rewrite (B.14) - (B.17) as
((k1 + k2 − i + ¯1 2 , l1 + l2 − j + 1 ¯2 ), (i − 1 ¯2 , j − ¯1 2 ))

(B.20)

where the sign¯denotes the logical negation of a boolean variable, so that ¯1,2 = 1 − 1,2 .
In this notation, i and j have the boundaries
1 ¯2 ≤ i ≤ min(k1 − ¯1 ¯2 , k2 − 1 2 )
¯1 2 ≤ j ≤ min(l1 − ¯1 ¯2 , l2 − 1 2 )

(B.21)

By denoting h1 = (k1 , l1 ) and h2 = (k2 , l2 ), together with R = ((a1 , b1 ), (a2 , b2 )) we
can then write
g(h1 , h2 ; R) = δk1 +k2 ,a1 δl1 +l2 +1,b1 δ−1,a2 δ0,b2 + δk1 +k2 +1,a1 δl1 +l2 ,b1 δ0,a2 δ−1,b2
+

1
X

min(k1 −¯
1 ¯2 ,k2 −1 2 )

min(l1 −¯
1 ¯2 ,l2 −1 2 )

X

X

1 ,2 =0

i=1 ¯2

j=¯
1 2

δk1 +k2 −i+¯1 2 ,a1 δl1 +l2 −j+1 ¯2 ,b1 δi−1 ¯2 ,a2 δj−¯1 2 ,b2
(B.22)

38

where we also added the two hooks in the depth two notation, (B.3) and (B.4). Explicitly,
summing over the 1,2 parameters, we get the lengthier expression
g(h1 , h2 ; R) =
min(k1 ,k2 ) min(l1 ,l2 )

= δk1 +k2 ,a1 δl1 +l2 +1,b1 δ−1,a2 δ0,b2 +

X

X

i=1

j=0

δk1 +k2 −i,a1 δl1 +l2 −j+1,b1 δi−1,a2 δj,b2

min(k1 ,k2 ) min(l1 ,l2 )

+ δk1 +k2 +1,a1 δl1 +l2 ,b1 δ0,a2 δ−1,b2 +

+

X

i=0

j=0

X

i=0

j=1

δk1 +k2 −i+1,a1 δl1 +l2 −j,b1 δi,a2 δj−1,b2


min(k1 ,k2 −1) min(l1 ,l2 −1)

X

X

min(k1 −1,k2 ) min(l1 −1,l2 )

+

X

X

i=0

j=0

 δk1 +k2 −i,a1 δl1 +l2 −j,b1 δi,a2 δj,b2 (B.23)

From this equation it is clear that g(h1 , h2 ; R) can be either 0, 1 or 2. In particular,
g(h1 , h2 ; R) = 2 only if R = ((k1 + k2 − i, l1 + l2 − j), (i, j)) and 0 ≤ i < min(k1 , k2 ),
0 ≤ j < min(l1 , l2 ).

C

Deriving the two point correlator

In this Appendix we will derive eq. (5.9) from eq. (5.7). Let us start by considering the
quantity


(X) (Y )
χR
T
T
T
(C.1)
1̄,1 [m] [n]
R1 ,R2
where we remind the reader that R1 , R2 and R are irreps of Sm , Sn and Sm+n respectively.
(Y )
(X)
(X,Y )
, T2 and T2 as the sum of transpositions in Sm+n , Sm and Sn
Let us define T2
respectively. We can expand (C.1) as


(X) (Y )
χR
T
T
T
1̄,1 [m] [n]
R1 ,R2
=

χR
R1 ,R2



(X,Y ) (X) (Y )
T2
T[m] T[n]

(X,Y )

χR (T2
=g
dR

−

)

(X)



−

χR
R1 ,R2



(X) (X) (Y )
T2 T[m] T[n]

(Y )

χR1 (T[m] ) χR2 (T[n] ) −
1

g dR1 dR2

1

−

χR
R1 ,R2
(X)

χR
R1 ,R2 (T2

g d R 1 dR 2


(Y )
(X) (Y )
R
R
χR1 ,R2 (T2 ) χR1 ,R2 T[m] T[n]

39





(X) (Y ) (Y )
T[m] T2 T[n]





(X) (Y )
T
) χR
T
+
R1 ,R2
[m] [n]

(X,Y )

χR (T2
=g
dR

)

(X)

(Y )
(X)
χR1 (T[m] ) χR2 (T[n] )

χR (T
− 1 2
dR1

)

(Y )

(X)

χR
R1 ,R2 (T[m] T[n] )+

(Y )

χR (T ) R
(X) (Y )
χR1 ,R2 (T[m] T[n] )
− 2 2
dR 2
"
#
(X,Y )
(X)
(Y )
χR (T2
) χR1 (T2 ) χR2 (T2 )
(Y )
(X)
−
−
= g χR1 (T[m] ) χR2 (T[n] )
dR
dR1
dR2

(C.2)

But now

 (−1)cR1 +1 (m − 1)! if R is a hook representation
1
χR1 (T[m] ) =

0
otherwise

(C.3)

where cR1 is the number of boxes in the firs column of the Young diagram associated with
the representation R1 . A similar equation holds for χR2 (T[n] ). We then have
χR
R1 ,R2

=



(X) (Y )
T1̄,1 T[m] T[n]



=



(X,Y )
)

c
+c
R
R
 (−1) 1 2 g (m − 1)!(n − 1)! χR (Td2
−
R



(X)

χR1 (T2
dR1

)

(Y )

−

χR2 (T2
dR2

)


;

R1 , R2 hooks

0 otherwise

(C.4)
this is eq. (5.8). Let us now restrict to the case in which both R1 , R2 are hooks representations. We will denote there representations as h1 = R1 = (k1 , l1 ) and h2 = R2 = (k2 , l2 ).
This also forces the representation R to be at most of depth two, as we derived in Appendix
B. We now consider such a representation. With the notation given at the beginning of
this section, R = ((a1 , b1 ), (a2 , b2 )), it is immediate to write an equation for the normalised
2)
character χRd(T
R
1X
χR (T2 )
=
ri (ri − 2i + 1) = a1 (a1 + 1) + (a2 + 2)(a2 − 1)+
dR
2 i
+2

bX
2 +2

bX
1 +1

i=3

i=b2 +3

(3 − 2i) + 2

(C.5)

(1 − i)

1
1
= (a21 + a22 + a1 + a2 ) − (b21 + b22 + b1 + b2 )
2
2
1
1
= (a1 + b1 + 1)(a1 − b1 ) + (a2 + b2 + 1)(a2 − b2 )
2
2

(C.6)

We now need the equivalent of this formula for the depth one representations h1 and h2 ,
i.e. the hooks. Such an equation can be directly obtained by setting (a2 , b2 ) = (−1, 0) or

40

(a2 , b2 ) = (0, −1) in (C.5). We can then write (C.4) as

 (−1)ch1 +ch2
(X) (Y )
g (m − 1)!(n − 1)!×
χR
T
T
T
=
1̄,1 [m] [n]
h1 ,h2
2
h
× (a1 + b1 + 1)(a1 − b1 ) + (a2 + b2 + 1)(a2 − b2 )+
− (k1 + l1 + 1)(k1 − l1 ) − (k2 + l2 + 1)(k2 − l2 )

(C.7)
i

where R = ((a1 , b1 ), (a2 , b2 )) and h1 = (k1 , l1 ), h2 = (k2 , l2 ).
The last piece we need is an equation for the U (N ) dimension of a depth two representation R = ((a1 , b1 ), (a2 , b2 )). It is straightforward to write





(a1 − a2 )(b1 − b2 )
a1 + b 1 a2 + b 2
N + a1
N + a2
DimN (R) =
(a1 + b2 + 1)(a2 + b1 + 1)
b1
b2
a1 + b1 + 1 a2 + b2 + 1
(C.8)
This equation reduces to its depth 1 equivalent by imposing (a2 , b2 ) = (−1, 0) or (a2 , b2 ) =
(0, −1). It is also helpful to recall the dimension formula for a Sl+k+1 hook representation
(k, l):


k+l
dR =
(C.9)
k
Let us now consider eq. (5.7):
hOO† i =

1 X
m!n! R `m
1
R2 `n

X
R`m+n

2


1
(X) (Y )
T
DimN (R) χR
T
T
1̄,1
R1 ,R2
[m] [n]
dR1 dR2 g

(C.10)

Inserting eq. (C.7), (C.8) and (C.9) into the above equation gives
hTr(X m Y n )Tr(X m Y n )† i
=

m
X

n
X

k1 ,l1 =0 k2 ,l2 =0

n+m
X

g δ(k1 + l1 − m) δ(k2 + l2 − n) F (a1 , b1 , a2 , b2 , k1 , l1 , k2 , l2 )

a1 ,b1 =0
a2 ,b2 =0

(C.11)
where we defined the function
k1 !k2 !l1 !l2 ! (a1 − a2 )(b1 − b2 )
4(a1 + b2 + 1)(a2 + b1 + 1)(k1 + l1 + 1)(k2 + l2 + 1)





a1 + b 1
a2 + b 2
N + a1
N + a2
×
×
b1
b2
a1 + b 1 + 1 a2 + b 2 + 1

F (a1 , b1 , a2 , b2 , k1 , l1 , k2 , l2 ) =

× ((a1 + b1 + 1)(a1 − b1 ) + (a2 + b2 + 1)(a2 − b2 )+
− (k1 + l1 + 1)(k1 − l1 ) − (k2 + l2 + 1)(k2 − l2 ))2

41

(C.12)

References
[1] J. M. Maldacena, “The Large N limit of superconformal field theories and
supergravity,” Int. J. Theor. Phys. 38 (1999) 1113–1133, arXiv:hep-th/9711200
[hep-th]. [Adv. Theor. Math. Phys.2,231(1998)].
[2] E. Witten, “Anti-de Sitter space and holography,” Adv. Theor. Math. Phys. 2
(1998) 253–291, arXiv:hep-th/9802150 [hep-th].
[3] S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, “Gauge theory correlators from
noncritical string theory,” Phys. Lett. B428 (1998) 105–114,
arXiv:hep-th/9802109 [hep-th].
[4] S. Corley, A. Jevicki, and S. Ramgoolam, “Exact correlators of giant gravitons from
dual N=4 SYM theory,” Adv.Theor.Math.Phys. 5 (2002) 809–839,
arXiv:hep-th/0111222 [hep-th].
[5] S. Corley and S. Ramgoolam, “Finite factorization equations and sum rules for
BPS correlators in N=4 SYM theory,” Nucl.Phys. B641 (2002) 131–187,
arXiv:hep-th/0205221 [hep-th].
[6] V. Balasubramanian, M.-x. Huang, T. S. Levi, and A. Naqvi, “Open strings from
N=4 superYang-Mills,” JHEP 08 (2002) 037, arXiv:hep-th/0204196 [hep-th].
[7] V. Balasubramanian, D. Berenstein, B. Feng, and M.-x. Huang, “D-branes in
Yang-Mills theory and emergent gauge symmetry,” JHEP 03 (2005) 006,
arXiv:hep-th/0411205 [hep-th].
[8] R. de Mello Koch, J. Smolic, and M. Smolic, “Giant Gravitons - with Strings
Attached (I),” JHEP 06 (2007) 074, arXiv:hep-th/0701066 [hep-th].
[9] R. de Mello Koch, J. Smolic, and M. Smolic, “Giant Gravitons - with Strings
Attached (II),” JHEP 09 (2007) 049, arXiv:hep-th/0701067 [hep-th].
[10] D. Bekker, R. de Mello Koch, and M. Stephanou, “Giant Gravitons - with Strings
Attached. III.,” JHEP 02 (2008) 029, arXiv:0710.5372 [hep-th].
[11] J. M. Maldacena and A. Strominger, “AdS(3) black holes and a stringy exclusion
principle,” JHEP 12 (1998) 005, arXiv:hep-th/9804085 [hep-th].
[12] J. McGreevy, L. Susskind, and N. Toumbas, “Invasion of the giant gravitons from
Anti-de Sitter space,” JHEP 06 (2000) 008, arXiv:hep-th/0003075 [hep-th].
[13] V. Balasubramanian, M. Berkooz, A. Naqvi, and M. J. Strassler, “Giant gravitons
in conformal field theory,” JHEP 04 (2002) 034, arXiv:hep-th/0107119
[hep-th].
[14] Y. Kimura and S. Ramgoolam, “Branes, anti-branes and brauer algebras in
gauge-gravity duality,” JHEP 0711 (2007) 078, arXiv:0709.2158 [hep-th].
[15] T. W. Brown, P. Heslop, and S. Ramgoolam, “Diagonal multi-matrix correlators
and BPS operators in N=4 SYM,” JHEP 0802 (2008) 030, arXiv:0711.0176
[hep-th].
42

[16] T. W. Brown, P. Heslop, and S. Ramgoolam, “Diagonal free field matrix
correlators, global symmetries and giant gravitons,” JHEP 0904 (2009) 089,
arXiv:0806.1911 [hep-th].
[17] R. Bhattacharyya, S. Collins, and R. d. M. Koch, “Exact Multi-Matrix
Correlators,” JHEP 0803 (2008) 044, arXiv:0801.2061 [hep-th].
[18] S. Collins, “Restricted Schur Polynomials and Finite N Counting,” Phys.Rev. D79
(2009) 026002, arXiv:0810.4217 [hep-th].
[19] P. Caputa, R. d. M. Koch, and P. Diaz, “Operators, Correlators and Free Fermions
for SO(N) and Sp(N),” JHEP 06 (2013) 018, arXiv:1303.7252 [hep-th].
[20] R. de Mello Koch, B. A. E. Mohammed, J. Murugan, and A. Prinsloo, “Beyond the
Planar Limit in ABJM,” JHEP 05 (2012) 037, arXiv:1202.4925 [hep-th].
[21] J. Pasukonis and S. Ramgoolam, “From counting to construction of BPS states in
N=4 SYM,” JHEP 02 (2011) 078, arXiv:1010.1683 [hep-th].
[22] J. Pasukonis and S. Ramgoolam, “Quivers as Calculators: Counting, Correlators
and Riemann Surfaces,” JHEP 1304 (2013) 094, arXiv:1301.1980 [hep-th].
[23] P. Mattioli and S. Ramgoolam, “Quivers, Words and Fundamentals,” JHEP 03
(2015) 105, arXiv:1412.5991 [hep-th].
[24] R. de Mello Koch, R. Kreyfelt, and N. Nokwara, “Finite N Quiver Gauge Theory,”
Phys. Rev. D89 no. 12, (2014) 126004, arXiv:1403.7592 [hep-th].
[25] D. Berenstein, “Extremal chiral ring states in the AdS/CFT correspondence are
described by free fermions for a generalized oscillator algebra,” Phys. Rev. D92
no. 4, (2015) 046006, arXiv:1504.05389 [hep-th].
[26] Y. Kimura, “Multi-matrix models and Noncommutative Frobenius algebras
obtained from symmetric groups and Brauer algebras,” Commun. Math. Phys. 337
no. 1, (2015) 1–40, arXiv:1403.6572 [hep-th].
[27] R. de Mello Koch and S. Ramgoolam, “A double coset ansatz for integrability in
AdS/CFT,” JHEP 06 (2012) 083, arXiv:1204.2153 [hep-th].
[28] W. Carlson, R. d. M. Koch, and H. Lin, “Nonplanar Integrability,” JHEP 03
(2011) 105, arXiv:1101.5404 [hep-th].
[29] R. d. M. Koch, B. A. E. Mohammed, and S. Smith, “Nonplanar Integrability:
Beyond the SU(2) Sector,” Int. J. Mod. Phys. A26 (2011) 4553–4583,
arXiv:1106.2483 [hep-th].
[30] R. d. M. Koch, M. Dessein, D. Giataganas, and C. Mathwin, “Giant Graviton
Oscillators,” JHEP 10 (2011) 009, arXiv:1108.2761 [hep-th].
[31] Y. Kimura and S. Ramgoolam, “Enhanced symmetries of gauge theory and
resolving the spectrum of local operators,” Phys.Rev. D78 (2008) 126003,
arXiv:0807.3696 [hep-th].

43

[32] P. Diaz, “Novel charges in CFT‘s,” JHEP 09 (2014) 031, arXiv:1406.7671
[hep-th].
[33] R. Goodman and N. Wallach, Representations and Invariants of the Classical
Groups. Encyclopedia of Mathematics and its Applications. Cambridge University
Press, 2000.
[34] A. Ram, “Representation theory; dissertation, chapter 1.”
http://www.ms.unimelb.edu.au/ ram/Preprints/dissertationChapt1.pdf, 1990.
Accessed: 18 January 2016.
[35] Y. Kimura and S. Ramgoolam, “Enhanced symmetries of gauge theory and
resolving the spectrum of local operators,” Phys. Rev. D78 (2008) 126003,
arXiv:0807.3696 [hep-th].
[36] R. Bhattacharyya, R. de Mello Koch, and M. Stephanou, “Exact Multi-Restricted
Schur Polynomial Correlators,” JHEP 06 (2008) 101, arXiv:0805.3025 [hep-th].
[37] M. Hamermesh, Group Theory and its Application to Physical Problems. Dover
Publications, 1989.
[38] F. Larsen and E. J. Martinec, “U(1) charges and moduli in the D1 - D5 system,”
JHEP 06 (1999) 019, arXiv:hep-th/9905064 [hep-th].
[39] S. Danz, J. Ellers, and J. Murray, “The centralizer of a subgroup in a group
algebra,” Proceedings of the Edinburgh Mathematical Society (Series 2) (2013)
49–56.
[40] W. Fulton and J. Harris, Representation Theory: A First Course. Graduate Texts
in Mathematics / Readings in Mathematics. Springer New York, 1991.
[41] M. Bianchi, F. A. Dolan, P. J. Heslop, and H. Osborn, “N=4 superconformal
characters and partition functions,” Nucl. Phys. B767 (2007) 163–226,
arXiv:hep-th/0609179 [hep-th].
[42] D. Vaman and H. L. Verlinde, “Bit strings from N=4 gauge theory,” JHEP 11
(2003) 041, arXiv:hep-th/0209215 [hep-th].
[43] T. Brown, “Cut-and-join operators and N=4 super Yang-Mills,” JHEP 1005 (2010)
058, arXiv:1002.2099 [hep-th].
[44] J. Pasukonis and S. Ramgoolam, “From counting to construction of BPS states in
N=4 SYM,” JHEP 1102 (2011) 078, arXiv:1010.1683 [hep-th].
[45] V. Balasubramanian, J. de Boer, V. Jejjala, and J. Simon, “The Library of Babel:
On the origin of gravitational thermodynamics,” JHEP 12 (2005) 006,
arXiv:hep-th/0508023 [hep-th].
[46] V. Balasubramanian, B. Czech, K. Larjo, and J. Simon, “Integrability versus
information loss: A Simple example,” JHEP 11 (2006) 001,
arXiv:hep-th/0602263 [hep-th].
[47] R. d. M. Koch and S. Ramgoolam, “From Matrix Models and Quantum Fields to
Hurwitz Space and the absolute Galois Group,” arXiv:1002.1634 [hep-th].
44

[48] G. W. Moore and G. Segal, “D-branes and K-theory in 2D topological field
theory,” arXiv:hep-th/0609042 [hep-th].
[49] D. Berenstein and R. G. Leigh, “Resolution of stringy singularities by
noncommutative algebras,” JHEP 06 (2001) 030, arXiv:hep-th/0105229
[hep-th].
[50] J. Ben Geloun and S. Ramgoolam, “Counting Tensor Model Observables and
Branched Covers of the 2-Sphere,” arXiv:1307.6490 [hep-th].
[51] A. Mironov, A. Morozov, and S. Natanzon, “A Hurwitz theory avatar of open-closed
strings,” Eur. Phys. J. C73 no. 2, (2013) 2324, arXiv:1208.5057 [hep-th].
[52] I. MacDonald, The Theory of Groups. Clarendon Press, 1968.

45

