Notes
Notes
Predrag Cvitanović
The fastest way to watch any week’s lecture videos is by letting YouTube run the
course playlist
and skipping ‘extras’ manually.
Contents
1 Linear algebra 7
Homework HW1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.1 Week’s videos, reading . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.2 Other sources (optional) . . . . . . . . . . . . . . . . . . . . . . . . 9
1.3 Special projects . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.4 Matrix-valued functions . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.5 A linear diversion . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.6 Eigenvalues and eigenvectors . . . . . . . . . . . . . . . . . . . . . . 14
1.6.1 Yes, but how do you really do it? . . . . . . . . . . . . . . . . 21
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
exercises 23
3 Group representations 43
Homework HW3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3.1 Week’s videos, reading . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.1.1 Matrix representations, Schur’s Lemma . . . . . . . . . . . . 44
3.1.2 Wonderful Orthogonality Theorem . . . . . . . . . . . . . . . 44
3
CONTENTS
5 It takes class 59
Homework HW5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
5.1 It’s all about class . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
5.1.1 Dirac characters, Burnside’s method (optional) . . . . . . . . 61
5.1.2 William G. Harter (optional) . . . . . . . . . . . . . . . . . . 62
5.2 Other sources (optional) . . . . . . . . . . . . . . . . . . . . . . . . 63
5.3 Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
exercises 65
6 For fundamentalists 67
Homework HW6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
6.1 Other sources (optional) . . . . . . . . . . . . . . . . . . . . . . . . 68
6.2 Thoughts (optional) . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
6.3 ChaosBook notes . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
6.4 Chaotic 3-spring, integrable 3-vortex systems (optional) . . . . . . . . 73
6.5 Eigenfunctions (optional) . . . . . . . . . . . . . . . . . . . . . . . . 76
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
exercises 80
8 Space groups 87
Homework HW8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
8.1 Other sources (optional) . . . . . . . . . . . . . . . . . . . . . . . . 88
8.2 Thoughts (optional) . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
8.3 Space groups . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
8.3.1 Wallpaper groups . . . . . . . . . . . . . . . . . . . . . . . . 91
8.3.2 One-dimensional line groups . . . . . . . . . . . . . . . . . . 93
8.3.3 Time reversal symmetry . . . . . . . . . . . . . . . . . . . . 93
birdtracks.eu/course3 4 2020-12-23
CONTENTS
2020-12-23 5 birdtracks.eu/course3
CONTENTS
birdtracks.eu/course3 6 2020-12-23
group theory - week 1
Linear algebra
Homework HW1
7
GROUP THEORY - WEEK 1. LINEAR ALGEBRA
◦ The fastest way to watch any week’s lecture videos is by letting YouTube run the
course playlist
Downside is that the playlist plays also all ‘extra’ videos - you can skip through
those, if you are short on time. Or patience.
Matrices : 2 kinds
Derivative of a matrix function
Exponential, logarithm of a matrix
Determinant is a volume
log det = tr log (updated Aug 18, 2020)
Multi-matrix functions (optional, for the QM inclined)
The subject of linear algebra is a vast and very alive research area, generates innumer-
able tomes of its own, and is way beyond what we can exhaustively cover here. Some
linear operators and matrices reading (optional reading for week 1, not required for
this course; whenever the text is colored, you can click on the live hyperlink in the pdf
version of these notes):
Stone and Goldbart [15], Mathematics for Physics: A Guided Tour for Graduate
Students, Appendix A. This is an advanced summary where you will find almost
everything one needs to know.
In sect. 1.4 I make matrix functions appear easier than they really are. For an in-
exercise 1.3
depth discussion, consult Golub and Van Loan [8] Matrix Computations, chap. 9
Functions of Matrices (click here).
Much more than you ever wanted to know about linear algebra: Axler [3] Down
with determinants! (click here).
Steve Trettel Linear Algebra and the Periodic Table is a gentle 53 min tour from
vectors to function spaces to quantum mechanics. True, what they teach you
as QM is 95% linear algebra, but Trettel does not mention that QM is 95% one
amazing experimental fact: ~ is a nature-given constant. Mathematicians...
1. The talk by David Weitz on melting of crystal lattices. Can you do a calculation
on a Wigner lattice or a graphene, or a silicon carbide polytype used as a substrate
in our graphene lab (ask Claire Berger about it), using our group theory methods
as applied to space groups (2- or 3-D lattices)?
2. If you are really wild about string theory, then you can read Giles and Thorn [7]
Lattice approach to string theory, and write up what you have learned as the
project report. The Giles-Thorn (GT) discretization of the worldsheet begins
with a representation of the free closed or open string propagator as a light-
cone worldsheet path integral defined on a lattice. The sequel Papathanasiou and
Thorn [14] Worldsheet propagator on the lightcone worldsheet lattice gives in
Appendix B 2D lattice Neumann open string, Dirichlet open string, and closed
string propagators. Discrete Green’s functions are explained, for example, by
Chung and Yau [4] who give explicitly, in their Theorem 6, a 2-dimensional lat-
tice Green’s function for a rectangular region R[`1 ×`2 ] . The paper is cited over
100 times, maybe there is a better, more up-to-date one to read in that list.
I recommend that you take a final, as these are hard and time-consuming projects,
and the faculty does not want to overburden you with course work. However, if a
project dovetails with your research interests, it might be worth it. Fly it by me.
Why should a working physicist care about linear algebra? Physicists were bliss-
fully ignorant of group theory until 1920’s, but with Heisenberg’s sojourn in Helgoland,
everything changed. Quantum Mechanics was formulated as
i
φ(t) = Û t φ(0) , Û t = e− ~ tĤ , (1.1)
where φ(t) is the quantum wave function t, Û t is the unitary quantum evolution opera-
tor, and Ĥ is the Hamiltonian operator. Fine, but what does this equation mean? In the
first lecture we deconstruct it.
The matrices that have to be evaluated are very high-dimensional, in principle in-
finite dimensional, and the numerical challenges can quickly get out of hand. What
made it possible to solve these equations analytically in 1920’s for a few iconic prob-
lems, such as the hydrogen atom, are the symmetries, or in other words group theory,
which start sketching out in the second lecture (and fill in the details in the next 27
lectures).
Whenever you are confused about an “operator”, think “matrix”. Here we recapit-
ulate a few matrix algebra concepts
Q that we found essential. The punch line is (1.35):
Hamilton-Cayley equation (M − λi 1) = 0 associates with each distinct root λi of a
matrix M a projection onto ith vector subspace
Y M − λj 1
Pi = .
λi − λj
j6=i
What follows - for this week - is a jumble of Predrag’s notes. If you understand the
examples, we are on the roll. If not, ask :)
How are we to think of the quantum operator
d 2 dA dA
A = A+A . (1.5)
dx dx dx
d −1
The derivative of the inverse of a matrix, if the inverse exists, follows from dx (AA ) =
0:
d −1 1 dA 1
A =− . (1.6)
dx A dx A
The first equation follows from the second one by the binomial theorem, so these in-
deed are equivalent definitions. That the terms of order O(N −2 ) or smaller do not
matter for a function of a single variable follows from the bound
N N N
x− x + δxN x+
1+ < 1+ < 1+ ,
N N N
where |δxN | < . If lim δxN → 0 as N → ∞, the extra terms do not contribute. A
proof for matrices would probably require defining the norm of a matrix (and, more
generally, a norm of an operator acting on a Banach space) first. If you know an easy
proof, let us know.
1
det (1 + A/N ) = 1 + tr A + O(N −2 ) .
N
hence
N N
1 tr A
det e A
= lim −2
1 + tr A + O(N ) = lim 1 + = etr A (1.10)
N →∞ N N →∞ N
ln det M = tr ln M , (1.11)
a crazy useful identity that you will run into over and over again.
t2 t3
etA Be−tA = B + t[A, B] + [A, [A, B]] + [A, [A, [A, B]]] + · · · (1.12)
2 3!
sometimes used to establish the equivalence of the Heisenberg and Schrödinger pic-
tures of quantum mechanics, follows by recursive evaluation of t derivatives
d
etA Be−tA = etA [A, B]e−tA .
dt
Expanding exp(A + B), exp A, exp B to first few orders using (1.7) yields
1
e(A+B)/N = eA/N eB/N − [A, B] + O(N −3 ) , (1.13)
2N 2
and the Trotter product formula: if B, C and A = B + C are matrices, then
N
eA = lim eB/N eC/N (1.14)
N →∞
In particular, we can now make sense of the quantum evolution operator (1.1) as a
succession of short free flights (kinetic term) interspersed by small acceleration kicks
(potential term),
N
e−itĤ = lim e−i∆t T̂ e−i∆t V̂ , ∆t = t/N , (1.15)
N →∞
Linear fields are the simplest vector fields, described by linear differential equations
which can be solved explicitly, with solutions that are good for all times. The state
space for linear differential equations is M = Rd , and the equations of motion are
written in terms of a vector x and a constant stability matrix A as
ẋ = v(x) = Ax . (1.16)
Solving this equation means finding the state space trajectory
x(t) = (x1 (t), x2 (t), . . . , xd (t))
passing through a given initial point x0 . If x(t) is a solution with x(0) = x0 and
y(t) another solution with y(0) = y0 , then the linear combination ax(t) + by(t) with
a, b ∈ R is also a solution, but now starting at the point ax0 + by0 . At any instant in
time, the space of solutions is a d-dimensional vector space, spanned by a basis of d
linearly independent solutions.
How do we solve the linear differential equation (1.16)? If instead of a matrix
equation we have a scalar one, ẋ = λx , the solution is x(t) = etλ x0 . In order to solve
the d-dimensional matrix case, it is helpful to rederive this solution by studying what
happens for a short time step ∆t. If time t = 0 coincides with position x(0), then
x(∆t) − x(0)
= λx(0) , (1.17)
∆t
which we iterate m times to obtain Euler’s formula for compounding interest
m
t
x(t) ≈ 1 + λ x(0) ≈ etλ x(0) . (1.18)
m
The term in parentheses acts on the initial condition x(0) and evolves it to x(t) by
taking m small time steps ∆t = t/m. As m → ∞, the term in parentheses converges
to etλ . Consider now the matrix version of equation (1.17):
x(∆t) − x(0)
= Ax(0) . (1.19)
∆t
A representative point x is now a vector in Rd acted on by the matrix A, as in (1.16).
Denoting by 1 the identity matrix, and repeating the steps (1.17) and (1.18) we obtain
Euler’s formula for the exponential of a matrix:
m
t
x(t) = J t x(0) , J t = etA = lim 1 + A , (1.20)
m→∞ m
where J t = J(t) is a short hand for exp(tA).
However, r, the number of distinct eigenvalues, is in general smaller than the dimension
of the matrix, r ≤ d (see example 1.3).
From (1.22) it follows that
matrix (M − λi 1) annihilates e(i) , thus the product of all such factors annihilates any
vector, and the matrix M satisfies its characteristic equation
d
Y
(M − λi 1) = 0 . (1.24)
i=1
This humble fact has a name: the Hamilton-Cayley theorem. If we delete one term from
this product, we find that the remainder projects x from (1.23) onto the corresponding
eigenspace: Y Y
(M − λj 1)x = (λi − λj ) xi e(i) .
j6=i j6=i
with the dimension of the ith subspace given by di = tr Pi . For each distinct eigen-
value λi of M,
(M − λj 1)Pj = Pj (M − λj 1) = 0 , (1.27)
the colums/rows of Pj are the right/left eigenvectors e(j) , e(j) of M which (provided
M is not of Jordan type, see example 1.3) span the corresponding linearized subspace.
The main take-home is that once the distinct non-zero eigenvalues {λi } are com-
puted, projection operators are polynomials in M which need no further diagonaliza-
tions or orthogonalizations.
It follows from the characteristic equation (1.27) that λi is the eigenvalue of M on
Pi subspace:
M Pi = λi Pi (no sum on i) . (1.28)
Using M = M 1 and completeness relation (1.26) we can rewrite M as
M = λ1 P1 + λ2 P2 + · · · + λd Pd . (1.29)
Any matrix function f (M) takes the scalar value f (λi ) on the Pi subspace, f (M) Pi =
f (λi ) Pi , and is thus easily evaluated through its spectral decomposition
X
f (M) = f (λi )Pi . (1.30)
i
This, of course, is the reason why anyone but a fool works with irreducible reps: they
reduce matrix (AKA “operator”) evaluations to manipulations with numbers.
By (1.27) every column of Pi is proportional to a right eigenvector e(i) , and its
every row to a left eigenvector e(i) . In general, neither set is orthogonal, but by the
idempotence condition (1.26), they are mutually orthogonal,
The non-zero constant cj is convention dependent and not worth fixing, unless you feel
nostalgic about Clebsch-Gordan coefficients. We shall set cj = 1. Then it is convenient
to collect all left and right eigenvectors into a single matrix.
Figure 1.2: Qualitatively distinct types of exponents {λ(1) , λ(2) } of a [2×2] Jacobian
matrix. Here the eigenvalues of the Jacobian matrix are multipliers Λ(j) , and the expo-
nents are defined as the deformation rates λ(j) = log(Λ(j) )/t.
These two possibilities are refined further into sub-cases depending on the signs of
the real part. In the case of real λ1 > 0, λ2 < 0, x1 grows exponentially with time, and x2
contracts exponentially. This behavior, called a saddle, is sketched in figure 1.1, as are
the remaining possibilities: in/out nodes, inward/outward spirals, and the center. The
magnitude of out-spiral |x(t)| diverges exponentially when µ > 0, and in-spiral contracts
into (0, 0) when µ < 0; whereas, the phase velocity ω controls its oscillations.
If eigenvalues λ1 = λ2 = λ are degenerate, the matrix might have two linearly
independent eigenvectors, or only one eigenvector, see example 1.3. We distinguish
two cases: (a) A can be brought to diagonal form and (b) A can be brought to Jordan
form, which (in dimension 2 or higher) has zeros everywhere except for the repeating
eigenvalues on the diagonal and some 1’s directly above it. For every such Jordan
[dα ×dα ] block there is only one eigenvector per block.
We sketch the full set of possibilities in figures 1.1 and 1.2.
Example 1.2. Complex eigenvalues: in-out spirals. As M has only real entries, it
will in general have either real eigenvalues, or complex conjugate pairs of eigenvalues.
Also the corresponding eigenvectors can be either real or complex. All coordinates used
in defining a dynamical flow are real numbers, so what is the meaning of a complex
eigenvector?
If λk , λk+1 eigenvalues that lie within a diagonal [2 × 2] sub-block M0 ⊂ M form
a complex conjugate pair, {λk , λk+1 } = {µ + iω, µ − iω}, the corresponding com-
plex eigenvectors can be replaced by their real and imaginary parts, {e(k) , e(k+1) } →
{Re e(k) , Im e(k) }. In this 2-dimensional real representation, M0 → A, the block A is a
sum of the rescaling×identity and the generator of SO(2) rotations in the {Re e(1) , Im e(1) }
plane.
µ −ω 1 0 0 −1
A= =µ +ω .
ω µ 0 1 1 0
Trajectories of ẋ = A x, given by x(t) = J t x(0), where (omitting e(3) , e(4) , · · · eigen-
directions)
cos ωt − sin ωt
J t = etA = etµ , (1.32)
sin ωt cos ωt
spiral in/out around (x, y) = (0, 0), see figure 1.1, with the rotation period T and the
radial expansion /contraction multiplier along the e(j) eigen-direction per a turn of the
spiral:
exercise 1.4
T = 2π/ω , Λradial = eTµ . (1.33)
We learn that the typical turnover time scale in the neighborhood of the equilibrium
(x, y) = (0, 0) is of order ≈ T (and not, let us say, 1000 T, or 10−2 T).
Example 1.3. Degenerate eigenvalues. While for a matrix with generic real
elements all eigenvalues are distinct with probability 1, that is not true in presence of
symmetries, or spacial parameter values (bifurcation points). What can one say about
situation where dα eigenvalues are degenerate, λα = λi = λi+1 = · · · = λi+dα −1 ?
Hamilton-Cayley (1.24) now takes form
r
Y X
(M − λα 1)dα = 0 , dα = d . (1.34)
α=1 α
M can be brought to diagonal form. The characteristic equation (1.34) can be re-
placed by the minimal polynomial,
r
Y
(M − λα 1) = 0 , (1.35)
α=1
where the product includes each distinct eigenvalue only once. Matrix M acts multi-
plicatively
M e(α,k) = λi e(α,k) , (1.36)
on a dα -dimensional subspace spanned by a linearly independent set of basis eigen-
vectors {e(α,1) , e(α,2) , · · · , e(α,dα ) }. This is the easy case. Luckily, if the degeneracy is
due to a finite or compact symmetry group, relevant M matrices can always be brought
to such diagonalizable form.
M can only be brought to upper-triangular, Jordan form. This is the messy case,
so we only illustrate the key idea in example 1.4.
takes us a long way toward developing intuition about arbitrary finite-dimensional matri-
ces. The eigenvalues
1 1p
λ1,2 = tr M ± (tr M)2 − 4 det M (1.37)
2 2
are the roots of the characteristic (secular) equation (1.21):
Distinct eigenvalues case has already been described in full generality. The left/right
eigenvectors are the rows/columns of projection operators (see example 1.5)
M − λ2 1 M − λ1 1
P1 = , P2 = , λ1 6= λ2 . (1.38)
λ1 − λ2 λ2 − λ1
v(2) helps span the 2-dimensional space, (M − λ)2 v(2) = 0, but is not an eigenvector,
as Mv(2) = λv(2) + e(1) . For every such Jordan [dα × dα ] block there is only one
eigenvector per block. Noting that
m
mλm−1
λ
Mm = ,
0 λm
picks up a power-low correction. That spells trouble (logarithmic term ln t if we bring the
extra term into the exponent).
That M satisfies its secular equation (Hamilton-Cayley theorem) can be verified by ex-
plicit calculation:
2
4 1 4 1 1 0 0 0
−6 +5 = .
3 2 3 2 0 1 0 0
Matrices Pi are orthonormal and complete. The dimension of the ith subspace is given
by di = tr Pi ; in case at hand both subspaces are 1-dimensional. From the charac-
teristic equation it follows that Pi satisfies the eigenvalue equation M Pi = λi Pi . Two
consequences are immediate. First, we can easily evaluate any function of M by spec-
tral decomposition, for example
58591 19531
M7 − 3 · 1 = (57 − 3)P1 + (1 − 3)P2 = .
58593 19529
Second, as Pi satisfies the eigenvalue equation, its every column is a right eigenvector,
and every row a left eigenvector. Picking first row/column we get the eigenvectors:
1 1
{e(1) , e(2) } = { , }
1 −3
3 1
{e(1) , e(2) } = { , },
1 −1
with overall scale arbitrary. The matrix is not symmetric, so {e(j) } do not form an orthog-
onal basis. The left-right eigenvector dot products e(j) · e(k) , however, are orthogonal
as in (1.31), by inspection. (Continued in example ??.)
Example 1.6. Computing matrix exponentials. If A is diagonal (the system is un-
coupled), then etA is given by
λ1 t
e
λ1 t
λ t
λ2 t e 2
exp = .
. .. . ..
λd t
λd t e
If A is diagonalizable, A = F DF −1 , where D is the diagonal matrix of the eigen-
values of A and F is the matrix of corresponding eigenvectors, the result is simple:
An = (F DF −1 )(F DF −1 ) . . . (F DF −1 ) = F Dn F −1 . Inserting this into the Taylor se-
ries for ex gives eAt = F eDt F −1 .
But A may not have d linearly independant eigenvectors, making F singular and
forcing us to take a different route. To illustrate this, consider [2×2] matrices. For any
linear system in R2 , there is a similarity transformation
B = U −1 AU ,
where the columns of U consist of the generalized eigenvectors of A such that B has
one of the following forms:
λ 0 λ 1 µ −ω
B= , B= , B= .
0 µ 0 λ ω µ
These three cases, called normal forms, correspond to A having (1) distinct real eigen-
values, (2) degenerate real eigenvalues, or (3) a complex pair of eigenvalues. It follows
that
λt
Bt e 0 Bt λt 1 t Bt at cos bt − sin bt
e = , e =e , e =e ,
0 eµt 0 1 sin bt cos bt
and eAt = U eBt U −1 . What we have done is classify all [2×2] matrices as belonging to
one of three classes of geometrical transformations. The first case is scaling, the second
is a shear, and the third is a combination of rotation and scaling. The generalization of
these normal forms to Rd is called the Jordan normal form. (J. Halcrow)
Commentary
Remark 1.1. Projection operators. The construction of projection operators given in
sect. 1.6.1 is taken from refs. [5, 6]. Who wrote this down first we do not know, lineage cer-
tainly goes all the way back to Lagrange polynomials [13], but projection operators tend to get
drowned in sea of algebraic details. Arfken and Weber [1] ascribe spectral decomposition (1.30)
to Sylvester. Halmos [9] is a good early reference - but we like Harter’s exposition [10–12] best,
for its multitude of specific examples and physical illustrations. In particular, by the time we
get to (1.27) we have tacitly assumed full diagonalizability of matrix M. That is the case for
the compact groups we will study here (they are all subgroups of U(n)) but not necessarily in
other applications. A bit of what happens then (nilpotent blocks) is touched upon in example 1.4.
Harter in his lecture Harter’s lecture 5 (starts about min. 31 into the lecture) explains this in great
detail - its well worth your time.
References
[1] G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists: A Compre-
hensive Guide, 6th ed. (Academic, New York, 2005).
[2] G. B. Arfken, H. J. Weber, and F. E. Harris, Mathematical Methods for Physi-
cists: A Comprehensive Guide, 7th ed. (Academic, New York, 2013).
[3] S. Axler, “Down with determinants!”, Amer. Math. Monthly 102, 139–154 (1995).
[4] F. Chung and S.-T. Yau, “Discrete Green’s functions”, J. Combin. Theory A 91,
191–214 (2000).
[5] P. Cvitanović, “Group theory for Feynman diagrams in non-Abelian gauge the-
ories”, Phys. Rev. D 14, 1536–1553 (1976).
[6] P. Cvitanović, Classical and exceptional Lie algebras as invariance algebras, Ox-
ford Univ. preprint 40/77, unpublished., 1977.
[7] R. Giles and C. B. Thorn, “Lattice approach to string theory”, Phys. Rev. D 16,
366–386 (1977).
[8] G. H. Golub and C. F. Van Loan, Matrix Computations, 4th ed. (J. Hopkins Univ.
Press, Baltimore, MD, 2013).
[9] P. R. Halmos, Finite-Dimensional Vector Spaces (Princeton Univ. Press, Prince-
ton NJ, 1948).
[10] W. G. Harter, “Algebraic theory of ray representations of finite groups”, J. Math.
Phys. 10, 739–752 (1969).
[11] W. G. Harter, Principles of Symmetry, Dynamics, and Spectroscopy (Wiley, New
York, 1993).
[12] W. G. Harter and N. dos Santos, “Double-group theory on the half-shell and the
two-level system. I. Rotation and half-integral spin states”, Amer. J. Phys. 46,
251–263 (1978).
[13] K. Hoffman and R. Kunze, Linear Algebra, 2nd ed. (Prentice-Hall, Englewood
Cliffs NJ, 1971).
Exercises
1.1. Trace-log of a matrix. Prove that
det M = etr ln M .
Generalize the square root function f (x) = x1/2 to a square root f (A) = A1/2 of a
matrix A.
a) Which one(s) of these is/are the square root of A
2 2 −2 10 −2 −2 2 −10
, , , ?
0 3 0 3 0 −3 0 −3
b) Assume that the eigenvalues of a [d × d] matrix are all distinct. How many square root
matrices does such matrix have?
c) Given a [2×2] matrix A with a distinct pair of eigenvalues {λ1 , λ2 }, write down a
formula that generates all square root matrices A1/2 . Hint: one can do this using the 2
projection operators associates with the matrix A. 2 points
1.4. Real representation of complex eigenvalues. (Verification of example 1.2.) λk , λk+1
eigenvalues form a complex conjugate pair, {λk , λk+1 } = {µ + iω, µ − iω}. Show that
(a) corresponding projection operators are complex conjugates of each other,
P = Pk , P∗ = Pk+1 ,
(c)
Pk 1 1 i R
= .
Pk+1 2 1 −i Q
(d) The · · · + λk Pk + λ∗k Pk+1 + · · · complex eigenvalue pair in the spectral decom-
position (1.29) is now replaced by a real [2×2] matrix
µ −ω R
··· + + ···
ω µ Q
or whatever you find the clearest way to write this real representation.
Homework HW2
Bonus points
Exercise 2.6 (a), (b) and (c) Permutation of three objects 2 points
Exercise 2.7 Three masses on a loop 6 points
Extra points accumulate, can help you later if you miss a few problems.
25
GROUP THEORY - WEEK 2. FINITE GROUPS - DEFINITIONS
Please do not get intimated by the length of this week’s notes - they are here more
for me than for you, as notes on these topics for future reference. If you understand the
main sequence of video clips, and recommended reading, that should suffice to do the
problems. The rest is optional, you can quickly skim over...
Section playlist
Dresselhaus et al. [4] Chapter 1 Basic Mathematical Background: Introduction.
The MIT course 6.734 online version contains much of the same material.
ChaosBook Chapter 10. Flips, slides and turns
(extra) Discussion: There might be many examples of it, but a ‘group’ itself
is an abstract notion. (3 min)
(extra) Discussion: permutations, symmetric group, simple groups, Italian
renaissance, French revolution, Galois (5:23 min)
For a deep dive into this material, here is your rabbit hole.
Nathan Carter Visual Group Theory (read it online through your university
library) seems very good. The next two online courses are based on it:
For a typical (but for this course advanced) application see, for example, Stone
and Goldbart [10], Mathematics for Physics: A Guided Tour for Graduate Stu-
dents, Section 14.3.2 Vibrational spectrum of H2 O (click here).
◦ Glance through sect. 2.6 Group presentations and sect. 2.8 Literature, but I do
not expect you to understand this material.
Table 2.1: The dihedral group D3 group multiplication table. Actually, we prefer
cyclic and dihedral groups notation ‘rotations’ r` and ‘flips’ σm , as in table 4.1.
where {xi } are small deviations from the equilibrium, resting points of the molecules
labelled i. Vij is a symmetric matrix, so it can be brought to a diagonal form by an
orthogonal transformation, to a set of N uncoupled harmonic oscillators or normal
modes of frequencies {ωi }.
N
X mi
ẏi2 + ωi2 yi2 .
x → y = U x, H= (2.3)
i=1
2
Group multiplication (or Cayley) tables, such as Table 2.1, define each distinct
discrete group, but they can be hard to digest. A Cayley graph, with links labeled
by generators, and the vertices corresponding to the group elements, has the same
information as the group multiplication table, but is often a more insightful presentation
of the group.
For example, the Cayley graph figure 2.1 is a clear presentation of the dihedral
group D4 of order 8,
Quaternion group is also of order 8, but with a distinct multiplication table / Cayley
graph, see figure 2.2. For more of such, see, for example, mathoverflow Cayley graph
discussion.
In 1937 R. Brauer [2] introduced diagrammatic notation for the Kronecker δij op-
eration, in order to represent “Brauer algebra” permutations, index contractions, and
matrix multiplication diagrammatically. His equation (39)
(send index 1 to 2, 2 to 4, contract ingoing (3·4), outgoing (1·3)) is the earliest published
diagrammatic notation I know about. While in kindergarten (disclosure: we were too
poor to afford kindergarten) I sat out to revolutionize modern group theory [3]. But I
To draw this, Brauer style, it is convenient to turn his drawing on a side. For 2-index
tensors, there are two permutations:
σ(23) = , σ(13) =
Here group element labels refer to the standard permutation cycles notation. There is
really no need to indicate the “time direction" by arrows, so we omit them from now
on.
The symmetric sum of all permutations,
1 n b1 b2 o
Sa1 a2 ...ap ,bp ...b2 b1 = δa1 δa2 . . . δabpp + δab12 δab21 . . . δabpp + . . .
p!
1
S = = + + + ... , (2.8)
p!
...
...
...
...
yields the symmetrization operator S. In birdtrack notation, a white bar drawn across
p lines [6] will always denote symmetrization of the lines crossed. A factor of 1/p! has
been introduced in order for S to satisfy the projection operator normalization
S2 = S
= . (2.9)
...
...
Fourier space duality, but you are not likely to recognize it. There the average was not
over all permutations, but the zero-th Fourier mode φ̃0 was the average over only cyclic
permutations. Every finite discrete group has such fully-symmetric representation, and
in statistical mechanics and quantum mechanics this is often the most important state
(the ‘ground’ state).
A subset of indices a1 , a2 , . . . aq , q < p can be symmetrized by symmetrization
matrix S12...q
...
...
S12...q = q. (2.10)
...
Overall symmetrization also symmetrizes any subset of indices:
SS12...q = S
= . (2.11)
... ...
... ...
... ...
...
... ...
σS = S
= . (2.12)
...
...
Diagrammatically this means that legs can be crossed and uncrossed at will.
One can construct a projection operator onto the fully antisymmetric space in a
similar manner [3]. Other representations are trickier - that’s precisely what the theory
of finite groups is about.
From Emory Math Department: A pariah is real! The simple finite groups fit into
18 families, except for the 26 sporadic groups. 20 sporadic groups AKA the Happy
Family are parts of the Monster group. The remaining six loners are known as the
pariahs.
Hang in there! And relax. None of this will be on the test. As a matter of fact, there
will be no test.
Question 2.1. Henriette Roux asks
Q What did you do this weekend?
A The same as every other weekend - prepared week’s lecture, with my helpers Avi the Little,
Edvard the Nordman, and Malbec el Argentino, under Master Roger’s watchful eye, see here.
2.9 Examples
Example 2.1. Discrete symmetries in physics:
• Point groups i.e., subgroups of O(3).
• Point groups + discrete translations e.g., symmetry groups of crystals.
• Permutation groups
SΨ(x1 , x2 , . . . xn ) = Ψ(x2 , x1 , . . . xn ).
• Boson wave functions are symmetric while fermion wave functions are anti-symmetric
under exchange of variables.
(B. Gutkin)
Example 2.2. The group multiplication table for D3 : See table 4.1.
L R
M m M
(B. Gutkin)
Example 2.4. Vibrations of a classical CO2 molecule: Consider one carbon and
two oxygens constrained to the x-axis [1] and joined by springs of stiffness k, as shown
in figure 2.4. Newton’s second law says
k
ẍ1 = − (x1 − x2 )
M
k k
ẍ2 = − (x2 − x3 ) − (x2 − x1 )
m m
k
ẍ3 = − (x3 − x2 ) . (2.14)
M
The normal modes, with time dependence xj (t) = xj exp(itω) , are the common fre-
quency ω vibrations that satisfy (2.14),
A −A 0 x1 x1
2
Hx = −a 2 a −a x2 = ω x2 , (2.15)
0 −A A x3 x3
First thing to always ask yourself is: does the system have a symmetry? Yes! Note
that the CO2 molecule (2.14) of figure 2.4 is invariant under x1 ↔ x3 interchange, i.e.,
coordinate relabeling by matrix σ that commutes with our law of motion H,
0 0 1 0 −A A
σ = 0 1 0 , σH = Hσ = −a 2 a −a . (2.16)
1 0 0 A −A 0
d+ = 2 , d− = 1 . (2.18)
In this subspace the outer oxygens are moving in opposite directions, with the carbon
stationary.
On the 2-dimensional symmetric subspace, the trace yields the sum of the remain-
ing two eigenvalues
1 k k
λ+ + λ0 = tr HP+ = (tr H + tr Hσ) = (a + A) + a = +2 .
2 M m
We could disentangle the two eigenfrequencies by evaluating tr H2 P+ , for example, but
thinking helps again.
There is still another, translational symmetry, so obvious that we forgot it; if we
change the origin of the x-axis, the three coordinates xj → xj − δx change, for any
continuous translation δx, but the equations of motion (2.14) do not change their form,
H x = H x + H δx = ω 2 x ⇒ H δx = 0 . (2.20)
(0)
So any translation e = δx = (δx, δx, δx) is a nul, ‘zero mode’ eigenvector of H
in (2.16), with eigenvalue λ0 = ω02 = 0, and thus the remaining eigenfrequency is
2
ω+ = k/M + 2 k/m. As we can add any nul eigenvector e(0) to the corresponding
e(+) eigenvector, there is some freedom in choosing e(+) . One visualization of the
corresponding eigenvector is the carbon moving opposite to the two oxygens, with total
momentum set to zero.
(Taken from AWH Example 6.2.3 Degenerate eigenproblem, but done here using
symmetries.)
where {xi } are small deviations from the resting the equilibrium, resting points of the
molecules labelled i. Vij is a symmetric matrix, so it can be brought to a diagonal form
by an orthogonal transformation, a set of N uncoupled harmonic oscillators or normal
modes of frequencies {ωi }.
N
X mi 2
ẏi + ωi2 yi2 .
x → y = U x, H= (2.21)
i=1
2
Consider now the ring of pair-wise interactions of two kinds of molecules sketched in
figure 2.5 (a), given by the potential
N
1X xi
k1 (xi − yi )2 + k2 (xi+1 − yi )2 ,
V (z) = zi = , (2.22)
2 i=1 yi
whose [2N ×2N ] matrix form is (aside to the cognoscenti: this is a Toeplitz matrix):
This potential matrix is a holy mess. How do we find an orthogonal transformation (2.21)
that diagonalizes it? Look at figure 2.5 (a). Molecules lie on a circle, so that suggests
we should use a Fourier representation. As the i = 1 labelling of the starting molecule
on a ring is arbitrary, we are free to relabel them, for example use the next molecule
pair as the starting one. This relabelling is accomplished by the [2N ×2N ] permutation
matrix (or ‘one-step shift’, ‘stepping’ or ‘translation’ matrix) M of form
0 0 ... 0 I z1 zn
I 0 ... 0 0 z2 z1
1 0 xi
0 I ... 0 0 z3 z2
= , I= , zi = (2.23)
0 1 yi
. .. .. .. .. .. ..
.. . . . . . .
0 0 ... I 0 zn zn−1
| {z }
M
Yn X1
ω y1
Xn
x1
Y1
optical m−1
X2 acoustic
x3
Y2
k
y3
(a) (b) −n/2 n/2 (c) x2
y2
Figure 2.5: (a) Chain with circular symmetry. (b) Dependance of frequency on the
representation wavenumber k. (c) Molecule with D3 symmetry. (B. Gutkin)
Projection operators corresponding to M are worked out in example 2.6. They are N
distinct [2N ×2N ] matrices,
1 0
0 1
λ 0
1 1
η1 = √ ,
0 η2 = √ λ .
n
n
.. ..
.
.
λn−1 0
0 λn−1
ω2
k1 + k2 −(k1 + k2 λ)
0 = det − I =⇒
−(k1 + k2 λ̄) k1 + k2 2
1 2 k
ω± (k) = k1 + k2 ± |k1 + k2 λ | , (2.25)
2
one acoustic (ω(0) = 0), one optical, see figure 2.5 (b) and the acoustic and optical
phonons wiki. (B. Gutkin)
where we have multiplied all denominators and numerators by λ−n . The numerator is
now a matrix polynomial of form (x − λ)(x − λ2 ) · · · (x − λN −1 ) , with the zeroth root
(x − λ0 ) = (x − 1) quotiented out from the defining matrix equation M N − 1 = 0. Using
1 − xN
= 1 + x + · · · + xN −1 = (x − λ)(x − λ2 ) · · · (x − λN −1 )
1−x
we obtain the projection operator in form of a discrete Fourier sum (rather than the
product (1.25)),
N −1
1 X i 2π
Pn = e N nm M m .
N m=0
This form of the projection operator is the simplest example of the key group theory tool,
projection operator expressed as a sum over characters,
1 X
Pn = χ̄(g)D(g) ,
|G| g∈G
upon which stands all that follows in this course. (B. Gutkin and P. Cvitanović)
In ChaosBook Chapter 10. Flips, slides and turns cosets are absolutely essential. The
significance of the coset is that if a solution has a symmetry, then the elements in a coset act on
the solution the same way, and generate all equivalent copies of this solution. Example 10.7.
Subgroups, cosets of D3 should help you understand that.
Henriette Roux writes: When talking about the cosets of a subgroup we demonstrated
multiplication between cosets with a specific example, but this wasn’t leading to
something along the lines of that the set of all left cosets of a subgroup (or the
set of all the right cosets of a subgroup) form a group, correct? It didn’t appear
so in the example since the “unit” {E, A} we looked appears to only have the
properties of an identity with multiplication from one direction (the direction
depending on if it is the set of left cosets or the set of right cosets). In the context
of the lecture I think this point was related to Lagrange’s theorem (although we
didn’t call it that) and I vaguely remember cosets being used in the proof of
Lagrange’s theorem but I wasn’t connecting it today. Are we going to cover that
in a future lecture?
Predrag You are right - Lagrange’s theorem (see the wiki) simply says the order of a
subgroup has to be a divisor of the order of the group. We used cosets to partition
elements of G to prove that. But what we really need cosets for is to define (see
Dresselhaus et al. [4] Sect. 1.7) Factor Groups whose elements are cosets of a
self-conjugate subgroup (click here). I will not cover that in a subsequent lecture,
so please read up on it yourself.
Henriette Roux You talked about the period of an element X, and said that that period
is the set
{E, X, · · · , X n−1 } , (2.27)
where n is the order of the element X. I had thought that set was the subgroup
generated by the element X and that the period of the element X was a synonym
for the order of the element X? Is that incorrect?
Predrag To keep things as simple as possible, in Thursday’s lecture I followed Sect. 1.3
Basic Definitions of Dresselhaus et al. textbook [4], to the letter. In Def. 3 the
order of an element X is the smallest n such that X n = E, and they call the set
(2.27) the period of X. I do not like that usage (and do not remember seeing it
anywhere else). As you would do, in ChaosBook.org Chap. Flips, slides and
turns I also define the smallest n to be the period of X and refer to the set (2.27)
as the orbit generated by X. When we get to compact continuous groups, the
orbit will be a (great) circle generated by a given Lie algebra element, and look
more like what we usually think of as an orbit.
I am not using my own ChaosBook.org here, not to confuse things further by
discussing both time evolution and its discrete symmetries. Here we focus on the
discrete group only (typically spatial reflections and finite angle rotations).
References
[1] G. B. Arfken, H. J. Weber, and F. E. Harris, Mathematical Methods for Physi-
cists: A Comprehensive Guide, 7th ed. (Academic, New York, 2013).
[2] R. Brauer, “On algebras which are connected with the semisimple continuous
groups”, Ann. Math. 38, 857 (1937).
[3] P. Cvitanović, Group Theory: Birdtracks, Lie’s and Exceptional Groups (Prince-
ton Univ. Press, Princeton NJ, 2008).
[4] M. S. Dresselhaus, G. Dresselhaus, and A. Jorio, Group Theory: Application to
the Physics of Condensed Matter (Springer, New York, 2007).
[5] R. Penrose, “Angular momentum: An approach to combinatorical space-time”,
in Quantum Theory and Beyond, edited by T. Bastin (Cambridge Univ. Press,
Cambridge, 1971).
[6] R. Penrose, “Applications of negative dimensional tensors”, in Combinatorial
mathematics and its applications, edited by D. J. A. Welsh (Academic, New
York, 1971), pp. 221–244.
[7] R. Penrose, The Road to Reality - A Complete Guide to the Laws of the Universe
(A. A. Knopf, New York, 2005).
[8] R. Penrose and M. A. H. MacCallum, “Twistor theory: An approach to the quan-
tisation of fields and space-time”, Phys. Rep. 6, 241–315 (1973).
[9] D. S. Silver, “The new language of mathematics”, Amer. Sci. 105, 364 (2017).
[10] M. Stone and P. Goldbart, Mathematics for Physics: A Guided Tour for Graduate
Students (Cambridge Univ. Press, Cambridge UK, 2009).
[11] M. Tinkham, Group Theory and Quantum Mechanics (Dover, New York, 2003).
[12] E. P. Wigner, Group Theory and Its Application to the Quantum Mechanics of
Atomic Spectra (Academic, New York, 1931).
Exercises
2.1. Gx ⊂ G. The maximal set of group actions which maps a state space point x into itself,
Gx = {g ∈ G : gx = x} , (2.28)
is called the isotropy group (or stability subgroup or little group) of x. Prove that the set
Gx as defined in (2.28) is a subgroup of G.
2.2. Transitivity of conjugation. Assume that g1 , g2 , g3 ∈ G and both g1 and g2 are
conjugate to g3 . Prove that g1 is conjugate to g2 .
2.3. Isotropy subgroup of gx. Prove that for g ∈ G, x and gx have conjugate isotropy
subgroups:
Ggx = g Gx g −1
2 3.
(a) List the group elements and the corresponding geometric operations
(b) Find the subgroups of the group D3 .
(c) Find the classes of D3 and the number of elements in them, guided by the geometric
interpretation of group elements. Verify your answer using the definition of a class.
(d) List the conjugacy classes of subgroups of D3 . (continued as exercise 4.1)
2.5. C4 invariant potential. Consider the Schrödinger equation for a particle moving in a
two-dimensional bounding potential V , such that the spectrum is discrete. Assume that
V is CN -invariant (in some literature, ZN -invariant), i.e., V remains invariant under the
rotation R by the angle 2π/N . For N = 3 case, figure 2.6 (a), the spectrum of the system
can be split into two sectors: {En0 } non-degenerate levels corresponding to symmetric
eigenfunctions φn (Rx) = φn (x) and doubly degenerate levels {En± } corresponding to
non-symmetric eigenfunctions φn (Rx) = e±2πi/3 φn (x).
(a) What is the spectral structure in the case of N = 4, figure 2.6 (b)?
How many sectors appear and what are their degeneracies?
(b) What is the spectral structure for general N ?
(c) A constant magnetic field normal to the 2D plane is added to V .
How will it affect the spectral structure?
(d) (bonus question) Figure out the spectral structure if the symmetry group of potential
is D3 (also includes 3 reflections), figure 2.6 (c).
(Boris Gutkin)
Figure 2.6: Hard wall potential with (a) symmetry C3 , (b) symmetry C4 , and (c) symmetry
D3 .
Figure 2.7: Three identical masses are constrained to move on a hoop, connected by
three identical springs such that the system wraps completely around the hoop. Find
the normal modes.
Group representations
Homework HW3
Bonus points
Exercise 3.4 Abelian groups 1 point
Exercise 3.5 Representations of CN 1 point
Total of 10 points = 100 % score. Extra points accumulate, can help you later if you
miss a few problems.
43
GROUP THEORY - WEEK 3. GROUP REPRESENTATIONS
Dresselhaus et al. [1] Sect. 2.7 ‘Wonderful’ Orthogonality Theorem, sect. 2.8
Representations and vector spaces.
Dresselhaus et al. [1] Sects. 3.1 Characters and Class to 3.5 The number of
irreducible representations.
Tinkham [3] covers the same material in Chapter 3 Theory of Group Representa-
tions, in a more compact way.
The Group Orthogonality Theorem: Let Dµ , Dµ0 be two irreducible matrix rep-
resentations of a compact group G of dimensions dµ , dµ0 , where the sum is over all
elements of the group, G = {g}, and |G| is their number, or the order of the group:
G
1 X (µ) 0 0 1 0
D (g)a b D(µ ) (g −1 )a b0 = δµµ0 δaa δbb0 .
|G| g dµ
This is a remarkable formula, one relation for each of the d2µ + d2µ0 matrix entries. Still,
the explicit matrix entries reflect largely arbitrary coordinate choices - there should a
more compact statement of irreducibility, and there is: the “character orthogonality
theorem” (3.1).
Consider a reducible representation D(g), i.e., a representation of group element g
that after a suitable similarity transformation takes form
(a)
D (g) 0 0 0
0 (b)
D (g) 0 0
D(g) = 0 ,
0 D(c) (g) 0
..
0 0 0 .
where ca , the multiplicity of the ath irreducible representation (colloquially called “ir-
rep”), is determined by the character orthonormality relations,
class
1 X
ca = χ(a)∗ χ = Nk χ(a) (Ck−1 ) χ(Ck ) . (3.1)
h
k
Knowing characters is all that is needed to figure out what any reducible representation
decomposes into! Work out exercise 4.2 as an example.
Section playlist
This week’s Dresselhous exposition (or the corresponding chapter in Tinkham [3])
comes from Wigner’s classic [4] Group Theory and Its Application to the Quantum
Mechanics of Atomic Spectra, which is a harder going, but the more group theory you
learn the more you’ll appreciate it. Eugene Wigner got the 1963 Nobel Prize in Physics,
so by mid 60’s gruppenpest was accepted in finer social circles.
In this course, we learn about full reducibility of finite and compact continuous
groups in two parallel ways. On one hand, I personally find the multiplicative projec-
tion operators (1.25), coupled with the notion of class algebras (Harter [2] (click here)
appendix C) most intuitive - a block-diagonalized subspace for each distinct eigenvalue
of a given all-commuting matrix.
On the other hand, the character weighted sums (here related to the multiplica-
tive projection operators as in ChaosBook Example A24.2 Projection operators for
discrete Fourier transform) offer a deceptively ‘simple’ and elegant formulation of
full-reducibility theorems, preferred by all standard textbook expositions.
References
[1] M. S. Dresselhaus, G. Dresselhaus, and A. Jorio, Group Theory: Application to
the Physics of Condensed Matter (Springer, New York, 2007).
[2] W. G. Harter and N. dos Santos, “Double-group theory on the half-shell and the
two-level system. I. Rotation and half-integral spin states”, Amer. J. Phys. 46,
251–263 (1978).
[3] M. Tinkham, Group Theory and Quantum Mechanics (Dover, New York, 2003).
[4] E. P. Wigner, Group Theory and Its Application to the Quantum Mechanics of
Atomic Spectra (Academic, New York, 1931).
Exercises
3.1. 1-dimensional representation of anything. Let D(g) be a representation of a group
G. Show that d(g) = det D(g) is one-dimensional representation of G as well.
(B. Gutkin)
3.2. 2-dimensional representation of S3 .
(i) Show that the group S3 can be generated by two permutations:
1 2 3 1 2 3
a= , d= .
1 3 2 3 1 2
with z = ei2π/3 , provide proper (faithful) representation for these elements and
find representation for the remaining elements of the group.
(iii) Is this representation irreducible?
(B. Gutkin)
3.3. 3-dimensional representations of D3 . The dihedral group D3 is the symmetry group
of the equilateral triangle. It has 6 elements
where C is rotation by 2π/3 and σ (i) is reflection along one of the 3 symmetry axes.
(i) Prove that this group is isomorphic to S3
(ii) Show that matrices
1 0 0 z 0 0 0 0 1
(1)
D(E) = 0 1 0 , D(C) = 0 1 0 , D(σ ) = 0 −1 0 ,
0 0 1 0 0 z2 1 0 0
(3.2)
generate a 3-dimensional representation D(g) of D3 . Hint: Calculate products for
representations of group elements and compare with the group table (see lecture).
(iii) Show that this is a reducible representation which can be split into one dimensional
A and two-dimensional representation Γ. In other words find a matrix R such that
A(g) 0
RD(g)R−1 =
0 Γ(g)
for all elements g of D3 . (Might help: D3 has only one (non-equivalent) 2-dim
irreducible representation).
(B. Gutkin)
3.4. Abelian groups. Let G be a group with only one-dimensional irreducible representa-
tions. Show that G is Abelian.
(B. Gutkin)
3.5. Representations of CN . Find all irreducible representations of CN .
(B. Gutkin)
Homework HW4
Bonus points
Exercise 4.4 Irreducible representations of dihedral group Dn 2 points
Exercise 4.5 Perturbation of Td symmetry 6 points
Exercise 4.7 Two particles in a potential 4 points
Total of 10 points = 100 % score. Bonus points accumulate, can help you later if you
miss a few problems.
49
GROUP THEORY - WEEK 4. HARD WORK BUILDS CHARACTER
Table 4.1: The D3 group multiplication table. The same as table 2.1, but written as a
class operator multiplication table.
D3 1 r r2 σ1 σ2 σ3
1 1 r r2 σ1 σ2 σ3 D3 C1 C2 C3
r r r2 1 σ3 σ1 σ2 C1 C1 C2 C3
r2 r2 1 r σ2 σ3 σ1 C2 C2 2C1 +C2 2C3
σ1 σ1 σ2 σ3 1 r r2 C3 C3 2C3 3C1 +3C2
σ2 σ2 σ3 σ1 r2 1 r
σ3 σ3 σ1 σ2 r r2 1
2021-06-01 Lecture 5
Character orthogonality theorem
2. it is all about class and character. “Character", in particular, I find very surprising
- one complex number suffices to characterize a matrix!
For a continuous group version of the character orthogonality theorem, see sect. 9.4. In
particular, the replacement of an irrep matrix representation D(µ) (g)a b by its character
χ(µ) (g) (a single scalar quantity) leads to no loss of any of the matrix indices structure.
I enjoyed reading Mathews and Walker [6] Chap. 16 Introduction to groups. You
can download it from here. Goldbart writes that the book is “based on lectures by
Richard Feynman at Cornell University.” Very clever. Try working through the exam-
ple of fig. 16.2: deadly cute, you get explicit eigenmodes from group theory alone. The
main message is that if you think things through first, you never have to construct the
representation matrices in explicit form - recasting the calculation in terms of invari-
ants, such characters, will get you there much faster.
You might find Gutkin notes useful:
Lect. 4 Representation Theory II, up to Sect. 4.5 Three types of representations:
Character tables. Dual character orthogonality. Regular representation. Indicators for
real, pseudo-real and complex representations. See example 4.3 “Irreps for quaternion
multiplication table.”
Oliver Pierson ChaosBook.org chapter Discrete factorization - Character tables
(10:05 min)
Oliver Pierson ChaosBook.org chapter Discrete factorization - Projection into
invariant subspaces (5:31 min)
Lect. 5 Applications I. Vibration modes go through Wigner’s theorem, Cn symme-
try and D3 symmetry. Study Example 5.1. Cn symmetry. More quantum mechanics
applications follow in
sect. 6.2 Applications II. Quantum Mechanics, Sect. 2. Perturbation theory.
Does the proof in the Lect. 4 Representation Theory II Appendix that the number
of irreps equals the number of classes make sense to you? For an easy argument, see
Vedensky Theorem 5.2 The number of irreducible representations of a group is equal
to the number of conjugacy classes of that group. For a proof, work though Murnaghan
Theorem 7. If you prefer a proof that your professor cannot understand, click here.
For the record (I retract the heady claim I made in class):
Mathworld.Wolfram.com: “A character table often contains enough information to
identify a given abstract group and distinguish it from others. However, there exist
nonisomorphic groups which nevertheless have the same character table, for example
D4 (the symmetry group of the square) and Q8 (the quaternion group).”
exercise 4.3
Fun read along these lines: Hart and Segerman [2] discuss the distinction between
abstract groups and symmetry groups of objects. They exhibit two very different ob-
jects with
D4 = hr, σ | σrσ = r−1 , r4 = σ 2 = ei (4.1)
symmetry (describing the group this way is called a presentation of D4 ), and explain
the Cayley graph for D4 (its edges with arrows correspond to rotations, the other edges
Simon Berman You would think that the analysis of three masses connected by har-
monic strings, see figure 4.1, is a simple exercise finding irreps of D3 symmetry,
but no, it merits a 2019 Phys. Rev. Lett., see Katz and Efrati [3] Self-driven frac-
tional rotational diffusion of the harmonic three-mass system. The article even
starts with our figure 4.1. We continue the discussion in sect. 6.4.
Example 4.1. D3 symmetry: Reflections and rotations of a triangle, figure 2.5 (c)
0 0 0 0 1 0 −1 0 0 0 0 0
0 0 0 0 0 1 0 1 0 0 0 0
1 0 0 0 0 0 0 0 0 0 −1 0
D(T ) =
0
, D(σ1 ) = (4.2)
1 0 0 0 0
0
0 0 0 0 1
0 0 1 0 0 0 0 0 −1 0 0 0
0 0 0 1 0 0 0 0 0 1 0 0
0 0 0 0 −1 0 0 0 −1 0 0 0
0 0 0 0 0 1 0 0 0 1 0 0
0 0 −1 0 0 0 −1 0 0 0 0 0
D(σ2 ) =
0
, D(σ3 ) =
0 0 1 0 0
0
1 0 0 0 0
−1 0 0 0 0 0 0 0 0 0 −1 0
0 1 0 0 0 0 0 0 0 0 0 1
(4.3)
G = {[e]; [g, g 2 ]; [σ1 , σ2 , σ3 ]}, χ(1) = {1, 1, 1}, χ(2) = {1, 1, −1}, χ(3) = {2, −1, 0}
1 X (i)
Pi = χ (g)D(g)
3 g∈G
0 0 0 0 0 0 1 0 1 0 1 0
0 1 0 1 0 1 0 0 0 0 0 0
1 0 0 0 0 0 0 1 1 0 1 0 1 0
P1 = , P2 = (4.4)
30 1 0 1 0 1 30 0 0 0 0 0
0 0 0 0 0 0 1 0 1 0 1 0
0 1 0 1 0 1 0 0 0 0 0 0
The 3 equal masses connected by harmonic springs system of figure 4.1 is a text-
book example of such system, see for example problems 6.37 and 9.16 in Kotkin and
Serbo [4] Collection of Problems in Classical Mechanics.
ν1 ν2
(ν1 3 −ν)2 2
The vibrational modes associated with the two 1-dimensional representations are
given by
0 1
1 0
0 1
P1V = α
1
and P 2V = β ,
0
0 1
1 0
respectively. Here P1V represents symmetric mode shown in figure 4.1 (red). The sec-
ond mode P2V corresponds to the rotations of the whole system. The projection opera-
tor for the two-dimensional representation is
2 0 −1 0 −1 0
0
2 0 −1 0 −1
2 2 1 −1 0 2 0 −1 0
P3 = (2D(I) − D(T ) − D(T )) = (4.5)
6 30 −1 0 2 0 −1
−1 0 −1 0 2 0
0 −1 0 −1 0 2
From this we have to separate two vectors corresponding to shift in x and y directions.
1 0
0 1
√
−1/2 − 3/2
ηx = √
3/2 , ηy =
−1/2
√
−1/2 3/2
√
− 3/2 −1/2
2 0 0 0
0 0 2 0
1 −1 1 1 1 0 1 0
P3V = α √ +β √ +γ √ +δ √ ,
6 0 2 0 6 −1 2 1
−1 −1 0 0
0 0 −1 −1
| {z } | {z } | {z } | {z }
ξ1 ξ2 ξ3 ξ4
p p
where ηx = 3/2(ξ4 + ξ1 ), ηy = 3/2(ξ3 − ξ2 ). Vectors ξi are columns of P3 and their
(B. Gutkin)
Example 4.2. (Pseudo)real and complex representations. There are three types
of representation: real, pseudo-real and complex (see Montaldi for details). For real
representations matrices D(g) can be brought into real form such that Dij (g) = D̄ij (g).
This implies in particular that all the characters are real. For pseudo-real representation
the characters are also real but matrices D(g) cannot be brought into real form. Finally,
for complex representations the characters are complex. In the last case D(g) and
the conjugate D̄(g) constitute two different representation (since they characters are
different), while in the real and pseudo-real case both representations are equivalent,
i.e., D̄(g) = U D(g)U † .
Indicator. To distinguish between three types of representations one looks at the indi-
cator:
1 X (α) 2
Ind(α) = χ (g ) ∈ {1, 0, −1} , (4.6)
|G| g∈G
where 1, −1, 0 are obtained for real, complex and pseudo-real representations, respec-
tively.
Proof: For a general irreducible representation we have
where α 6= β for a complex representation (since χ(α) (g) 6= χ̄(α) (g)) and α = β for real
and pseudo-real representations. From D(α) (g 2 ) = D(α) (g)D(α) (g) follows
mα X mα X
X 1 X (α) (β) †
Ind(α) = Uk,j Di,k (g)D̄j,n (g)Uni ,
i,j=1 g∈G
|G| g∈G
k,n=1
with mα being dimension of α. By the orthogonality theorem this expression is zero for
α 6= β which is the case of complex α. For real and pseudo-real representations we
have
1
Ind(α) = tr U Ū .
mα
Now note, that for α = β eq. (4.7) yields
By the first Schur’s lemma it follows then that U Ū = γI, or U = γU > which also implies
γ 2 = 1. This leaves only two possibilities γ = 1 for real and γ = −1 for pseudo-real
representations. In the first case we have U U > = I and Ind(α) = 1, while in the second
one U U > = −I and Ind(α) = −1. Note finally, that 1 = det U Ū = γ mα . So γ = −1
might appear only if mα is even. In other words, a pseudo-real irreducible representation
must be of even dimension.
the last representation is pseudo-real. Note that this representation can be realized
using Pauli matrices:
{±I, ±σx , ±σy , ±σz }.
(B. Gutkin)
References
[1] M. S. Dresselhaus, G. Dresselhaus, and A. Jorio, Group Theory: Application to
the Physics of Condensed Matter (Springer, New York, 2007).
[2] V. Hart and H. Segerman, The quaternion group as a symmetry group, in Proc.
Bridges 2014: Mathematics, Music, Art, Architecture, Culture, edited by G. H.
G. Greenfield and R. Sarhangi (2014), pp. 143–150.
[3] O. Katz-Saporta and E. Efrati, “Self-driven fractional rotational diffusion of the
harmonic three-mass system”, Phys. Rev. Lett. 122, 024102 (2019).
[4] G. L. Kotkin and V. G. Serbo, Collection of Problems in Classical Mechanics
(Elsevier, 2013).
[5] L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory
(Pergamon Press, Oxford, 1959).
[6] J. Mathews and R. L. Walker, Mathematical Methods of Physics (W. A. Ben-
jamin, Reading, MA, 1970).
[7] M. Tinkham, Group Theory and Quantum Mechanics (Dover, New York, 2003).
Exercises
4.1. Characters of D3 . (continued from exercise 2.4) D3 , the group of symmetries of an
equilateral triangle: has three irreducible representations, two one-dimensional and the
other one of multiplicity 2.
(a) All finite discrete groups are isomorphic to a permutation group or one of its sub-
groups, and elements of the permutation group can be expressed as cycles. Express
the elements of the group D3 as cycles. For example, one of the rotations is (123),
meaning that vertex 1 maps to 2, 2 → 3, and 3 → 1.
(b) Use your representation from exercise 2.4 to compute the D3 character table.
(c) Use a more elegant method from the group-theory literature to verify your D3 char-
acter table.
(d) Two D3 irreducible representations are one dimensional and the third one of multi-
plicity 2 is formed by [2×2] matrices. Find the matrices for all six group elements
in this representation.
since
0 1 0 a b
1 0 0 b = a
0 0 1 c c
(a) What are the selection rules? Between which energy levels of the system are tran-
sitions possible?
(b) Would the answer be different if a magnetic field B 0 cos(ωt) is added instead?
Explain how and why.
It takes class
Homework HW5
Exercise 5.1 Vibration modes of CH4 , parts (a) (b) (c) i 8 points
Exercise 5.2 Keep it classy (a) 2 points
Bonus points
Exercise 5.1 Vibration modes of CH4 , part (c) ii 2 points
Exercise 5.2 Keep it classy (b) 2 points
Exercise 5.2 Keep it classy (c) 4 points
Total of 10 points = 100 % score. Extra points accumulate, can help you later if you
miss a few problems.
59
GROUP THEORY - WEEK 5. IT TAKES CLASS
Show class, have pride, and display character. If you do, win-
ning takes care of itself.
— Paul Bryant
constructs ‘class operators’, shows that they form a basis for the algebra of ‘central’
or ‘all-commuting’ operators, and uses their characteristic equations to construct the
projection operators (1.27) from the ‘structure constants’ of the finite group, i.e., its
class multiplication tables. Expanded, these projection operators are indeed the same
as the ones obtained from character orthogonality.
I find Harter’s Sect. 3.3 Second stage of non-Abelian symmetry analysis particu-
larly illuminating. It shows how physically different (but mathematically isomorphic)
higher-dimensional irreps are constructed corresponding to different subgroup embed-
dings. One chooses the irrep that corresponds to a particular sequence of physical
symmetry breakings.
You might want to have a look at Harter [9] Double group theory on the half-shell
(click here). Read appendices B and C on spectral decomposition and class algebras.
Article works out some interesting examples.
See also remark 1.1 Projection operators and perhaps watch Harter’s online lecture
from Harter’s online course.
There is more detail than what we have time to cover here, but I find Harter’s
Sect. 3.3 Second stage of non-Abelian symmetry analysis particularly illuminating. It
shows how physically different (but mathematically isomorphic) higher-dimensional
irreps are constructed corresponding to different subgroup embeddings. One chooses
the irrep that corresponds to a particular sequence of physical symmetry breakings.
[...] the evaluation of Dirac and similar character operators is all that is
required for the solution of the standard molecular problems in the spirit
of Dirac’s original program which avoids appeal to formal group theory.
Dirac characters use not only the abstract group information, but also account for
the symmetry information contained in the basis set used. The diagonalization of Dirac
characters has three main advantages:
1. It can be realized by means of a quite simple and general algorithm.
2. The projective irreps obtained are just the ones that are needed to reduce the
starting basis set into irreducible sets.
3. No tabulated quantities are required to construct the projective irreps.
The scheme is completely general, in the sense that it applies to all space groups.
Ananda Dasgupta had 1.68K followers on YouTube, now he has one more:
Lecture 15 (start at about 35 min into the lecture) has a nice discussion of Dirac
characters, their relation to characters, and motivates the algorithmic Burnside’s
method for computing characters via class multiplication tables (Hi )jk .
More generally, the whole course is of interest, if covers most topics of our course
in greater depth:
Cini [2] Topics and Methods in Condensed Matter Theory (2007) (click here)
Jacobs [10] Group Theory with Applications in Chemical Physics, (click here)
(2005)
El-Batanouny and Wooten [1] Symmetry and Condensed Matter Physics: A Com-
putational Approach (2008) (click here). In sect. 4.3 they describe the Burnside’s
method. They give an example of Mathematica code that constructs the charac-
ter table. If needed, on might use Dixon’s method, which is more clever for
numerical computations.
The CRYSTAL package performs ab initio calculations of the ground state energy,
energy gradient, electronic wave function and properties of periodic systems.
Uses Dirac characters.
physics library had been checked out in 1960-62 by Gell-Mann, Zweig and Glashow.
That only half of the entering students were meant to complete their PhD’s there led to
lots of ugly competition. Harter transferred to UC Irvine, and, upon graduation, got a
job at USC in LA. After a few years he suggested in a faculty meeting that the way they
could improve their quality as a department was “to get rid of all the old farts." These
same “old farts" soon voted to deny him tenure. He ended up in Campinas, Brazil. Fox
rescued him from there by bringing him for an interview at Georgia Tech, where he
was hired in late 1970’s. He was brilliant, an asset for teaching, making all sorts of
demonstration devices. He built a giant rotating table upon which he placed billiard
balls, a wonderful demonstration of mechanical analogues for charged particle motion
in crossed E and B fields. Everyone (except for one nefarious character) liked him, his
work, and especially his devices. The faculty unanimously supported his promotion
to tenure. He did not, however, think much of the Director of School of Physics, and
made that clear. After an argument with the Director, he stormed out, offended. So, he
was denied tenure and moved in 1985 to University of Arkansas where he is a professor
today.
In 1987 Harter and Weeks used Harter’s theory of the rotational dynamics of mole-
cules to calculate the rotational-vibrational spectra of the soccer ball-shaped molecule
Buckminsterfullerene, C60, or “buckyball." C60 had been proposed in 1985 by chemists,
who had seen a mass-spectra peak of atomic mass 720. By 1989 the Harter theory cal-
culations led to a realization that chemists had been making C60 since the early 1970s.
In 1992 Science named C60 “Molecule of the Year," and in 1996 Curl, Kroto and
Smalley were awarded the Nobel Prize in Chemistry for their discovery of fullerenes.
You can find here many Soft Elegant Educational tools developed by Harter, and
follow his lectures on line. He is a great teacher. Georgia Tech’s loss.
5.3 Discussion
Henriette Roux I have a few questions about the exercise 5.1 part (d) Vibration modes
of CH4 : Find all modes of the methane molecule.
1. When we use the angle of improper rotation, is it true that reflection equals
to the π improper rotation?
2. I assume it is π and it gives me other characters are zero. In the case of all
symmetry, this will give the , which we usually get non-negative integer.
As a result, I’m not perfectly sure that the character formulas you give are
correct.
References
[1] M. El-Batanouny and F. Wooten, Symmetry and Condensed Matter Physics: A
Computational Approach (Cambridge Univ. Press, Cambridge UK, 2008).
[2] M. Cini, Topics and Methods in Condensed Matter Theory - From Basic Quan-
tum Mechanics to the Frontiers of Research (Springer, Berlin, 2007).
[3] E. M. Corson, “Note on the Dirac character operators”, Phys. Rev. 73, 57–60
(1948).
[4] P. Cvitanović, Group Theory: Birdtracks, Lie’s and Exceptional Groups (Prince-
ton Univ. Press, Princeton NJ, 2008).
[5] P. A. M. Dirac, The Principles of Quantum Mechanics (Oxford Univ. Press,
1930).
[6] M. S. Dresselhaus, G. Dresselhaus, and A. Jorio, Group Theory: Application to
the Physics of Condensed Matter (Springer, New York, 2007).
[7] M. Hamermesh, Group Theory and Its Application to Physical Problems (Dover,
New York, 1962).
[8] W. G. Harter, Principles of Symmetry, Dynamics, and Spectroscopy (Wiley, New
York, 1993).
[9] W. G. Harter and N. dos Santos, “Double-group theory on the half-shell and the
two-level system. I. Rotation and half-integral spin states”, Amer. J. Phys. 46,
251–263 (1978).
[10] P. Jacobs, Group Theory with Applications in Chemical Physics (Cambridge
Univ. Press, 2005).
[11] J. Mathews and R. L. Walker, Mathematical Methods of Physics (W. A. Ben-
jamin, Reading, MA, 1970).
a) b)
Figure 5.1: a) Two classes of rotational symmetries, and a class of reflection symme-
tries of a tetrahedron. (left) Hold the Tetra Pak by a tip, turn it by a third. (middle)
Hold the Tetra Pak by the midpoints of a pair of opposing edges, make a half-turn.
(right) Exchange the vertices outside the reflection plane. b) Methane molecule with
the symmetry Td .
Exercises
5.1. Vibration modes of CH4 .
Tetrahedral group T describes rotational symmetries of a tetrahedron. The order of the
group is |T| = 12, and its conjugacy classes are:
• The identity mapping.
• Four rotations by ϕ = 2π/3, with each of the four rotation axes going through a
vertex, and piercing the midpoint of the triangle opposite.
• Four inverse rotations by ϕ = −2π/3.
• Three rotations by ϕ = π, one for each of the three rotation axes going through
midpoints of opposing edges.
The full group of tetrahedron symmetries Td includes also reflections. This is the sym-
metry group of molecules such as methane CH4 , see figure 5.1).
(a) What is the order of the group Td ? Show that the group is isomorphic to i) the
group of permutations S4 ; ii) to the group O of rotational symmetries of the cube.
iii) Show that T is normal subgroup of Td .
(b) Find all conjugacy classes of the group. Which of these classes correspond to proper
(det R(ϕ) = +1), improper (det R(ϕ) = −1) rotations ?
Information on T might help. Note that ϕ might be also 0.
(c) i) Find all irreducible representations of the group & build the character table.
A shortcut: find all one-dimensional representations, assume that characters are
integers, then use the orthogonality relationship between characters.
ii) Really compute the character table, without assuming that characters are integers
(2 bonus points).
One-dimensional representations + orthogonality of characters is not enough to
build the whole character table for Td . One needs more black magic, such as rep-
resentation of permutation group by matrices.
(d) Find all modes of the methane molecule. Which of them correspond to vibrations,
translations and rotations? What are the degeneracies?
Figure 5.2: Three identical masses are constrained to move on a hoop, connected by
three identical springs such that the system wraps completely around the hoop. Find
all symmetries of the equations of motion.
Path: Find characters of the full (reducible) representation by using formulas from
the lecture:
(
ng (1 + 2 cos(ϕ)) rotation,
χ(g) =
ng (−1 + 2 cos(ϕ)) improper rotation .
Here ng is the number of atoms staying at the same place under the action of g, ϕ is
the rotation angle corresponding to g = R(ϕ). Then decompose this representation
into irreducible representations. Identify the rotational and translational parts.
(e) To what representation corresponds the most symmetric ”breathing“ mode and
why? Is it infrared active, i.e., can this mode can be excited by electromagnetic
field?
(B. Gutkin)
5.2. Keep it classy. Check out Harter’s PowerPoint presentation :)
(a) Go through the derivation of the three projection operators for D3 .
(3) (3) (3)
(b) Decompose P(3) = P1 + P2 . Construct Pij . Verify that they are idempotent.
(d) Compute the [2 × 2] irreducible matrix representation D(3) (g)ij for every group
element g, in the spirit of Harter’s slides 13-8 and 13-9.
5.3. Three masses on a loop. (Exercise 2.7 revisited.) Three identical masses, connected
by three identical springs, are constrained to move on a circle hoop as shown in figure 5.2.
For fundamentalists
Homework HW6
Total of 10 points = 100 % score. Extra points accumulate, can help you later if you
miss a few problems.
67
GROUP THEORY - WEEK 6. FOR FUNDAMENTALISTS
Henriette Roux What do the parameters σ, ρ and b stand for in the Lorenz equations
(6.4)?
Predrag The short answer is the truncation of the Navier-Stokes that leads to Lorenz
equations is so drastic that they have no longer any physical meaning; in his
1963 paper [17] Lorenz played with the parameters until he empirically found
an interesting example of deterministic chaos. Since then, applied mathemati-
cians have reverse-engineered various physical systems to find situations where
parameters σ, ρ and b mean something, see remark 6.1 (copied to here from
ChaosBook.org). The discrete symmetry of the original Navier-Stokes system
(‘left’ is as good as ‘right’) happened to survive the drastic truncation from 105
Fourier modes (for physically accurate simulations) to 3. I prefer to teach non-
linear dynamics using the Rössler system, precisely because it has no discrete
symmetry, just chaos.
(a) (b)
Figure 6.1: (a) The pair of full-space 9-cycles, the counter-clockwise 121232313 and
the clockwise 131323212 correspond to (b) one fundamental domain 3-cycle 001.
Table 6.1: D3 correspondence between the binary labeled fundamental domain prime cycles
p̃ and the full 3-disk ternary labeled cycles p, together with the D3 transformation that maps
the end point of the p̃ cycle into the irreducible segment of the p cycle. White spaces in the
above ternary sequences mark repeats of the irreducible segment; for example, the full space
12-cycle 1212 3131 2323 consists of 1212 and its symmetry related segments 3131, 2323. The
multiplicity of p cycle is mp = 6np̃ /np . The shortest pair of fundamental domain cycles related
by time reversal (but no spatial symmetry) are the 6-cycles 001011 and 001101.
p̃ p gp̃ p̃ p gp̃
0 12 σ12 000001 121212 131313 σ23
1 123 C 000011 121212 313131 232323 C2
01 12 13 σ23 000101 121213 e
001 121 232 313 C 000111 121213 212123 σ12
011 121 323 σ13 001011 121232 131323 σ23
0001 1212 1313 σ23 001101 121231 323213 σ13
0011 1212 3131 2323 C2 001111 121231 232312 313123 C
0111 1213 2123 σ12 010111 121312 313231 232123 C2
00001 12121 23232 31313 C 011111 121321 323123 σ13
00011 12121 32323 σ13 0000001 1212121 2323232 3131313 C
00101 12123 21213 σ12 0000011 1212121 3232323 σ13
00111 12123 e 0000101 1212123 2121213 σ12
01011 12131 23212 31323 C 0000111 1212123 e
01111 12132 13123 σ23 ··· ··· ···
(r)2 = 1 condition decomposes the state space into two linearly irreducible subspaces
M = M+ ⊕ M− , the z-axis M+ and the [x, y] plane M− , with projection operators
onto the two subspaces given by
0 0 0 1 0 0
+ 1 − 1
P = (1 + r) = 0 0 0 , P = (1 − r) = 0 1 0 . (6.7)
2 2
0 0 1 0 0 0
so
ẋ− −σ σ x− 0
= +
ẏ− ρ −1 y− −z x−
1
ż+ = −b z+ + (x+ + x− )(y+ + y− ) , (6.8)
4
where z+ = z. As (ẋ+ , ẏ+ ) = (0, 0), values of (x+ , y+ ) are conserved parts of the initial
condition. We define the fundamental domain by the (arbitrary) condition x̂− ≥ 0, and
whenever exits the domain,we replace the function dependence by the corresponding
fundamental domain coordinates,
(a) (b)
Figure 6.2: (a) Lorenz flow cut by y = x Poincaré section plane P through the z axis
and both E1,2 equilibria. Points where flow pierces into section are marked by dots.
To aid visualization of the flow near the E0 equilibrium, the flow is cut by the second
Poincaré section, P 0 , through y = −x and the z axis. (b) Poincaré sections P and P 0
laid side-by-side. (E. Siminos)
and record that we have applied r (that is the ‘reconstruction equation’ in the case of a
discrete symmetry). When we integrate (6.8), the trajectory coordinates (x̂− (t), ŷ− (t))
are discontinuous whenever the trajectory crosses the fundamental domain border.
That, however, we do not care about - the only thing we need are the Poincaré sec-
tion points and the Poincaré return map in the fundamental domain.
Poincaré section hypersurface can be specified implicitly by a single condition, through
a function U (x) that is zero whenever a point x is on the Poincaré section,
In order that there is only one copy of the section in the fundamental domain, this con-
dition has to be invariant, U (gx̂) = U (x̂) for g ∈ G, or, equivalently, the normal to it has
to be equivariant
∂j U (gx̂) = g∂j U (x̂) for g ∈ G . (6.10)
There are two kinds of compact (finite-time) orbits. Periodic orbits x(Tp ) = x(Tp )
are either self dual under rotation r, or appear in pairs related by r; in the fundamental
domain there is only one copy x̂(Tp ) = x̂(Tp ) of each. Relative periodic orbits (or ‘pre-
periodic orbits’) x̂(Tp ) = rx(Tp ) they are periodic orbits.
As the flow is C2 -invariant, so is its linearization ẋ = Ax. Evaluated at E0 , A com-
mutes with r, and the E0 stability matrix A decomposes into [x, y] and z blocks.
The 1-dimensional M+ subspace is the fixed-point subspace, with the z-axis points
left point-wise invariant under the group action
(here x = (x, y, z) is a 3-dimensional vector, not the coordinate x). A C2 -fixed point x(t)
in Fix (C2 ) moves with time, but according to (6.2) remains within x(t) ∈ Fix (C2 ) for all
times; the subspace M+ = Fix (C2 ) is flow invariant. In case at hand this jargon is a bit
of an overkill: clearly for (x, y, z) = (0, 0, z) the full state space Lorenz equation (6.4) is
reduced to the exponential contraction to the E0 equilibrium,
ż = −b z . (6.12)
(x̂, ŷ, z) = (r cos 2θ, r sin 2θ, z) = ((x2 − y 2 )/r, 2xy/r, z) , (6.13)
as in figure 6.4 (a). In contrast to the original G-equivariant coordinates [x, y, z], the
Lorenz flow expressed in the new coordinates [x̂, ŷ, z] is G-invariant. In this representa-
tion the M̃ = M/C2 fundamental domain flow is a smooth, continuous flow, with (any
choice of) the fundamental domain stretched out to seamlessly cover the entire [x̂, ŷ]
plane.
(E. Siminos and J. Halcrow)
Simon Berman According to the 2019 Phys. Rev. Lett., of Katz-Saporta and Efrati [15],
Self-driven fractional rotational diffusion of the harmonic three-mass system, a
system of three masses connected by harmonic springs might be the simplest
mechanical system (homonuclear triatomic molecule, such as ozone, except the
three couplings are not the same) that exhibits a geometric phase. Away from
its resting configuration the system is nonlinear, and once its rotational SO(2)
symmetry is reduced, and as its energy is increased, it exhibits all kinds of shape-
dependent chaotic geometric phases. Katz and Efrati [15] mostly do numerical
simulations and plot displacement vs. time diffusion plots in its 6D phase space,
like this is still early 1960’s. The earlier arXiv:1706.09868 version has more
information than the PRL (no doubt thanks to impatience of the referees, plus
space constraints of a PRL). One suspects that a bit of thinking along periodic
orbit theory lines could yield some insight into the diffusive properties of its
shape-changing dynamics.
In the symmetry-reduced or the ‘shape’ state space there is a D3 symmetry. One
sees it in their [15] Hamiltonian (2): the bij vectors can be viewed as the three
coordinates of an equilateral triangle in the w1 −w2 plane. Since the Hamiltonian
only depends on |w| and in a symmetric way on w · bij , it has a D3 symmetry for
(w1 , w2 ) components of the w vector, and a reflection symmetry for w3 . So the
total symmetry group is D3 × C 1/2 .
By going to relative w’s coordinates, one has quotiented only the 2D Euclidean
translations and SO(2) rotations, no discrete symmetries, so D3 still remains.
Now, anyone who has taken ChaosBook.org/course1 knows that the next step
is to quotient D3 , and do the calculation in the 1/6th of the phase space, i.e., the
fundamental domain.
I’m curious whether I’m right, because soon we’ll look at space groups (infinite
lattices with discrete symmetries) and there I have confused understanding of
how to quotient the space group, but that is related to diffusion in space, rather
than the angular diffusion, as in this 3-springs system.
We can make this a course project for a student in this course (a project instead
of taking the final). To be especially pedagogical, we’ll ask them to do it in Julia
(there is one potential candidate on Piazza).
Predrag: N vortex system Went to hear Tomoki Ohsawa ( Google Scholar), talk
about Symplectic reduction and the Lie–Poisson dynamics of point vortices on
the plane, arXiv:1808.01769 .
I had previously written to Tomoki’s friend Molei about how much I had already
suffered through Weinstein, Marsden, etc. moment maps, for decades. We all
have to do symmetry reductions, but with Marsden it is always the moment map,
and then the climax is the rigid 3D body example which is the end-all of every
article and book. Perhaps due to my pleas, Ohsawa gave us a gentle, sensible
seminar, Weinstein-Marsden for humans, where he explained why moment map
is called ‘moment,’ etc.. As nice a birthday present one could hope for, see the
slides here.
Ohsawa develops a Hamiltonian formulation of the dynamics of the “shape" of
N point vortices on the plane and the sphere. If N = 3, it is the dynamics of the
shape of the triangle formed by three point vortices, regardless of the position
and orientation of the triangle on the plane/sphere.
For the planar case, reducing the basic equations of point vortex dynamics by the
special Euclidean group SE(2) yields a Lie-Poisson equation for relative config-
urations of the vortices. The shape dynamics is periodic in certain cases. The
approach can be extended to the spherical case by first lifting the dynamics from
the two-sphere to C 2 and then performing reductions by symmetries.
Figure 6.3: (green) Level set of quadratic Casimir C2 ellipsoid. (orange) Level set of
Hamiltonian h. The intersection is the unique periodic orbit of the symmetry-reduced
N = 3 vortex system. See also the corresponding figure for 3 vortices on the sphere
on p. 46 of the slide presentation. T. Ohsawa
I think Ohsawa discovery is that the system has a previously un-noted quartic
Casimir, whose invariance reduces the dimension of the symmetry reduced phase
space by one degree of freedom (dof). The implies that the symmetry-reduced
dynamics in the N = 3 case is 1 dof, i.e., integrable, see figure 6.3. In addition,
if the sum of the vortex circulations is 0, the N = 4 case is integrable. This fact
is not yet explained - my intuition is that the zero total circulation implies extra
rotational symmetry. For more vortices I expect the usual Hamiltonian mixed
phase space.
Parenthetically, statement that there is quartic order Casimir that is invariant un-
der the symmetry group can probably be written as a syzygy constraint on the
invariant polynomial basis in (Hilbert’s) theory of invariant polynomial bases.
In the Katz-Saporta and Efrati [15] example there is no quartic Casimir, so one
ends with a generic chaotic system. Ohsawa’s geometric technique works be-
cause of the simple symplectic structure on the point-vortex problem (there is no
‘momentum’), whereas Katz-Saporta and Efrati problem is a standard classical-
mechanical one on the cotangent bundle of a configuration space, with momen-
tum there. Ohsawa believes that one can apply the techniques developed by
Richard Montgomery to this setting as well. (Montgomery’s paper motivated
him to work on the point-vortex problem).
There are also examples in cardiac (!) dynamics where one must reduce 2D
Euclidean symmetry first, with similar outcome to yours, but no moment maps,
as such PDEs have no variational formulation (that I am aware of). Googling
“Barkley model" might do the trick. I do not think there is a variational (La-
grangina) formulation.
But that is the whole point - any flow with a symmetry has to have the symmetry
quotiented out. It’s easier to understand this for flows which are not symplectic
- in that case, every continuous symmetry parameter reduces the dimension of
the symmetry-reduced state space by one. The Hamiltonin case is a pain (or
bliss, if you love moment maps) because every continuous symmetry reduces
the dimension of the phase space by one degree of freedom (ie, by 2). Also
variational problems obey Noether’s theorem, our (dissipative) problems usually
do not. If I understand this right...
Kimberly Short Foulkes PhD thesis [7] Drift and Meander of Spiral Waves (2009)
might be an user friendly introduction for students that need to understand Eu-
clidean symmetry? Covers refs. [2, 8]. Page 14 solves a differential equation
with SE(2) symmetry. Appendix 8.9 discusses symmetries, and gives the con-
dition for equivariance.
(a) (b)
Figure 6.4: (a) Lorenz attractor plotted in [x̂, ŷ, z], the doubled-polar angle coordinates
(6.13), with points related by π-rotation in the [x, y] plane identified. Stable eigen-
vectors of E0 : e(3) and e(2) , along the z axis (6.12). Unstable manifold orbit W u (E0 )
(green) is a continuation of the unstable e(1) of E0 . (b) Blow-up of the region near E1 :
The unstable eigenplane of E1 defined by Re e(2) and Im e(2) , the stable eigenvector
e(3) . The descent of the E0 unstable manifold (green) defines the innermost edge of
the strange attractor. As it is clear from (a), it also defines its outermost edge. (E.
Siminos)
Predrag Heilman and Strichartz [13] Homotopies of Eigenfunctions and the Spectrum
of the Laplacian on the Sierpinski Carpet, arXiv:0908.2942, is not an obvious
read for us, but they compute a spectrum on a square domain, and we might
have to be mindful of it: “ Since all of our domains are invariant under the
D4 symmetry group, we can simplify the eigenfunction computations by reduc-
ing to a fundamental domain. On this domain we impose appropriate boundary
conditions according to the rep-resentation type. For the 1-dimensional repre-
sentation, we restrict to the sector 0 ≤ θ ≤ π/4 . Recall that the functions will
extend evenly when reflected about θ = 0 in the 1++ and 1– cases, and oddly
in the 1-+ and 1+- cases. Note that performing an even extension across a ray is
equivalent to imposing Neumann boundary conditions on that ray. Similarly, the
odd extension is equivalent to Dirichlet conditions. For the 2-dimensional rep-
resentation our fundamental domain is the sector 0 ≤ θ ≤ π/2 , and we impose
Neumann boundary conditions on the ray θ = 0 and Dirichlet conditions on the
ray θ = π/2. Note that our fundamental domains are simply connected. ”
This seems to be saying that one gets the 2-dimensional representation by dou-
bling the fundamental domain and mixing boundary conditions. Do you under-
stand that?
Boris Here is my present understanding of the fundamental domains issue: If you want
simple boundary conditions like Dirichlet or Neumann you stick to 1d represen-
tations only. They connect eigenfunction to itself at the fundamental domain
boundaries – otherwise you would need to connect pair of functions (would be
Predrag Your argument is in the spirit of Harter’s class operators construction (see
week 5) of higher-dimensional representations by using particular chains of sub-
groups, but I am not able to visualize how that larger fundamental domain (of
the lower-order subgroup) folds back into the small fundamental domain of the
whole group. By the time the dust settles, I have the symmetry factorization of
the determinants that we need, but I do not have a gut feeling for the boundary
conditions that you do, when it comes to higher-dimensional irreps.
Commentary
Remark 6.1. Lorenz equation. The Lorenz equation (6.4) is the most celebrated early illus-
tration of “deterministic chaos” [17] (but not the first - that honor goes to Dame Cartwright [3]
in 1945. Amusingly, Denisov and Ponomarev [6] argue that Ben F. Laposky might have been
the first to observe chaotic attractors as early as 1953, which, strictly speaking falls after 1945,
even in Russia). Lorenz’s 1963 paper, which can be found in reprint collections refs. [5, 12],
is a pleasure to read, and it is still one of the best introductions to the physics motivating such
models (read more about Lorenz here). The equations, a set of ODEs in R3 , exhibit strange
attractors. W. Tucker [24–26] has proven rigorously (via interval arithmetic) that the Lorenz
attractor is strange for the original parameters (no stable orbits) and that it has a long stable
periodic orbit for slightly different parameters. In contrast to the hyperbolic strange attractors
such as the weakly perturbed cat map [4], the Lorenz attractor is structurally unstable. Frøy-
land [9] has a nice brief discussion of Lorenz flow. Frøyland and Alfsen [10] plot many periodic
and heteroclinic orbits of the Lorenz flow; some of the symmetric ones are included in ref. [9].
Guckenheimer-Williams [11] and Afraimovich-Bykov-Shilnikov [1] offer an in-depth discus-
sion of the Lorenz equation. The most detailed study of the Lorenz equation was undertaken by
Sparrow [22]. For a geophysics derivation, see Rothman course notes [20]. For a physical inter-
pretation of ρ as “Rayleigh number,” see Jackson [14] and Seydel [21]. The Lorenz truncation
to 3 modes, however, is so drastic that the model bears no relation to the geophysical hydro-
dynamics problem that motivated it. Just for fun, as Lorentz was such a lovable weatherman,
in 1972 Willem Malkus constructed [18], by a feat of reverse engineering, a physical system,
a “water wheel”, popularized by Strogatz [23], that is described by Lorentz equations. You
can see it simulated on wolfram.com, and tested experimentally at http://www.ace.gatech.edu.
There is no deep physics in this lovely game, it is but a cute distraction. For detailed pictures of
Lorenz invariant manifolds consult Vol II of Jackson [14] and “Realtime visualization of invari-
ant manifolds” by Ronzan. The Lorenz attractor is a very thin fractal – as we shall see, stable
manifold thickness is of the order 10−4 – whose fractal structure has been accurately resolved
by D. Viswanath [27, 28]. If you wonder what analytic function theory has to say about Lorenz,
check ref. [29]. Modular flows are your thing? E. Ghys and J. Leys have a beautiful tale for you.
Refs. [16, 19] might also be of interest.
References
[1] V. S. Afraimovich, B. B. Bykov, and L. P. Shilnikov, “On the origin and structure
of the Lorenz attractor”, Dokl. Akad. Nauk SSSR 234, In Russian, 336–339
(1977).
[2] I. V. Biktasheva, D. Barkley, V. N. Biktashev, G. V. Bordyugov, and A. J. Foulkes,
“Computation of the response functions of spiral waves in active media”, Phys.
Rev. E 79, 056702 (2009).
[3] M. L. Cartwright and J. E. Littlewood, “On non-linear differential equations of
the second order”, J. London Math. Soc. 20, 180–189 (1945).
[4] S. C. Creagh, “Quantum zeta function for perturbed cat maps”, Chaos 5, 477–
493 (1995).
[5] P. Cvitanović, Universality in Chaos, 2nd ed. (Adam Hilger, Bristol, 1989).
[6] S. Denisov and A. V. Ponomarev, “Oscillons: An encounter with dynamical
chaos in 1953?”, Chaos 21, 023123 (2011).
[7] A. J. Foulkes, Drift and Meander of Spiral Waves, PhD thesis (Univ. of Liver-
pool, Liverpool, UK, 2020).
[8] A. J. Foulkes and V. N. Biktashev, “Riding a spiral wave: numerical simulation
of spiral waves in a comoving frame of reference”, Phys. Rev. E 81, 046702
(2010).
[9] J. Frøyland, Introduction to Chaos and Coherence (Taylor & Francis, Bristol,
1992).
[10] J. Frøyland and K. H. Alfsen, “Lyapunov-exponent spectra for the Lorenz model”,
Phys. Rev. A 29, 2928 (1984).
[11] J. Guckenheimer and R. Williams, “Structural stability of the Lorenz attractor”,
Publ. Math. IHES 50, 55–72 (1979).
[12] B.-L. Hao, Chaos II (World Scientific, Singapore, 1990).
[13] S. M. Heilman and R. S. Strichartz, “Homotopies of eigenfunctions and the spec-
trum of the Laplacian on the Sierpinski carpet”, Fractals 18, 1–34 (2010).
[14] E. A. Jackson, Perspectives of Nonlinear Dynamics, Vol. 1 (Cambridge Univ.
Press, Cambridge UK, 1989).
[15] O. Katz-Saporta and E. Efrati, “Self-driven fractional rotational diffusion of the
harmonic three-mass system”, Phys. Rev. Lett. 122, 024102 (2019).
[16] J. B. Laughlin and P. C. Martin, “Transition to turbulence of a statically stressed
fluid”, Phys. Rev. Lett. 33, 1189 (1974).
[17] E. N. Lorenz, “Deterministic nonperiodic flow”, J. Atmos. Sci. 20, 130–141
(1963).
Exercises
6.1. 3-disk symbolic dynamics. As periodic trajectories will turn out to be our main tool to
breach deep into the realm of chaos, it pays to start familiarizing oneself with them now
by sketching and counting the few shortest prime cycles. Show that the 3-disk pinball has
3 · 2n−1 itineraries of length n. List periodic orbits of lengths 2, 3, 4, 5, · · · . Verify that
the shortest 3-disk prime cycles are 12, 13, 23, 123, 132, 1213, 1232, 1323, 12123, · · · .
Try to sketch them. (continued in exercise 6.3)
A comment about exercise 6.1, exercise 6.2, and exercise 6.3: If parts of these problems
seem trivial - they are. The intention is just to check that you understand what these
symbolic dynamics codings are - the main message is that the really smart coding (fun-
damental domain) is 1-to-1 given by the group theory operations that map a point in the
fundamental domain to where it is in the full state space. This observation you might not
find deep, but it is - instead of absolute labels, in presence of symmetries one only needs to
keep track of relative motions, from domain to domain, does not matter which domain in
absolute coordinates - that is what group actions do. And thus the word ‘relative’ creeps
into this exposition.
6.2. Reduction of 3-disk symbolic dynamics to binary. (continued from exercise 6.1)
(a) Verify that the 3-disk cycles
{1 2, 1 3, 2 3}, {1 2 3, 1 3 2}, {12 13 + 2 perms.},
{121 232 313 + 5 perms.}, {121 323+ 2 perms.}, · · · ,
correspond to the fundamental domain cycles 0, 1, 01, 001, 011, · · · respectively.
(b) Check the reduction for short cycles in table 6.1 by drawing them both in the full
3-disk system and in the fundamental domain, as in figure 6.3.
(c) Optional: Can you see how the group elements listed in table 6.1 relate irreducible
segments to the fundamental domain periodic orbits?
(continued in exercise 6.3)
6.3. 3-disk fundamental domain cycles. Try to sketch 0, 1, 01, 001, 011, · · · . in the
fundamental domain, and interpret the symbols {0, 1} by relating them to topologically
distinct types of collisions. Compare with table 6.1. Then try to sketch the location of
periodic points in the Poincaré section of the billiard flow. The point of this exercise is
that while in the configuration space longer cycles look like a hopeless jumble, in the
Poincaré section they are clearly and logically ordered. The Poincaré section is always to
be preferred to projections of a flow onto the configuration space coordinates, or any other
subset of state space coordinates which does not respect the topological organization of
the flow.
6.4. C2 -equivariance of Lorenz system. Verify that the vector field in Lorenz equations
(6.4)
ẋ σ(y − x)
ẋ = v(x) = ẏ = ρx − y − xz (6.14)
ż xy − bz
is equivariant under the action of cyclic group C2 = {e, r} acting on R3 by a π rotation
about the z axis,
r(x, y, z) = (−x, −y, z) ,
as claimed in example 6.1.
6.5. Proto-Lorenz system. Here we quotient out the C2 symmetry by constructing an
explicit “intensity” representation of the desymmetrized Lorenz flow.
2. Show that this is the (Lorenz)/C2 quotient map for the Lorenz flow, i.e., that it
identifies points related by the π rotation (6.6).
3. Show that (6.15) is invertible. Where does the inverse not exist?
4. Compute the equilibria of proto-Lorenz and their stabilities. Compare with the
equilibria of the Lorenz flow.
5. Plot the strange attractor both in the original form (6.4) and in the proto-Lorenz
form (6.16)
45
40
35
30
W
25
20
15
10
5
0 200 400 600 800
V
for the Lorenz parameter values σ = 10, b = 8/3, ρ = 28. Topologically, does it
resemble more the Lorenz, or the Rössler attractor, or neither? (plot by J. Halcrow)
6. Show that a periodic orbit of the proto-Lorenz is either a periodic orbit or a relative
periodic orbit of the Lorenz flow.
7. Show that if a periodic orbit of the proto-Lorenz is also periodic orbit of the Lorenz
flow, their Floquet multipliers are the same. How do the Floquet multipliers of
relative periodic orbits of the Lorenz flow relate to the Floquet multipliers of the
proto-Lorenz?
8. Show that the coordinate change (6.15) is the same as rewriting
r
ṙ = (−σ − 1 + (σ + ρ − z) sin 2θ
2
+(1 − σ) cos 2θ)
1
θ̇ = (−σ + ρ − z + (σ − 1) sin 2θ
2
+(σ + ρ − z) cos 2θ)
r2
ż = −bz + sin 2θ . (6.17)
2
in variables
(u, v) = (r2 cos 2θ, r2 sin 2θ) ,
i.e., squaring a complex number z = x + iy, z 2 = u + iv.
Homework HW7
Total of 16 points = 100 % score. Bonus points accumulate, can help you later if you
miss a few problems.
83
GROUP THEORY - WEEK 7. DISCRETE FOURIER REPRESENTATION
(extra) Rant: Symmetrize you must. Karl Schwarzschild found his exact so-
lution in 1915, a month after the publication of Einstein’s theory of general
relativity, while serving on a World War I front. (3 min)
A meta truth; We live in The Matrix; Fourier transformation is just a matrix (10
min)
(extra) Rocket science needs complex numbers; Why Fourier? Digital image
processing! (8 min)
Exercises
7.1. Am I a group? Show that multiplication table
e a b c d f
e e a b c d f
a a e d b f c
b b d e f c a
c c b f e a d
d d f c a e b
f f c a d b e
describes a group. Or does it? (Hint: check whether this table satisfies the group axioms.)
7.2. Product of two groups. Let G1 and G2 be two finite groups. The elements of the
product set G = G1 × G2 are defined as pairs (g1 , g2 ), g1 ∈ G1 g2 ∈ G2 .
(a) Show that G is a group with the multiplication operation (g1 , g2 ) · (g10 , g20 ) =
(g1 g10 , g2 g20 ).
Let D1 be an irreducible representation of G1 and let D2 be an irreducible representation
of G2 . For each g = (g1 , g2 ) ∈ G define D(g) = D1 (g1 ) × D2 (g2 )
(b) Show that D = D1 × D2 is an irreducible representation of G. What are the
characters of D?
7.3. Laplacian is a non-local operator.
While the Laplacian is a simple tri-diagonal difference operator, its inverse (the “free”
propagator of statistical mechanics and quantum field theory) is a messier object. A way
to compute is to start expanding propagator as a power series in the Laplacian
∞
1 1 X 1
= 2 ∆n . (7.1)
m2 1−∆ m n=0 m2n
Space groups
Homework HW8
Bonus points
Exercise 8.3 Tight binding model 8 points
Total of 10 points = 100 % score. Extra points accumulate, can help you later if you
miss a few problems.
87
GROUP THEORY - WEEK 8. SPACE GROUPS
Gutkin lecture notes Lect. 7 Applications III. Energy Band Structure, Sect. 7.2
Lattice symmetries.
Liang and Cvitanović [23] A chaotic lattice field theory in one dimension (2022).
Also good reads: Dresselhaus et al. [11] chapter 9. Space Groups in Real Space
(click here), and Cornwell [9] chapter 7. Crystallographic Space Groups (click here).
Walt De Heer learned this stuff from Herzberg [15] Molecular Spectra and Molec-
ular Structure. Condensed matter people like Kittel [21] Introduction to Solid State
Physics, but I am not a fan, because simple group theoretical facts are there presented
as condensed matter phenomena.
Quinn and Yi [25] Solid State Physics: Principles and Modern Applications intro-
duction to space groups looks compact and sensible. Band structure of graphene.
Martin Mourigal found the Presqu’île Giens, May 2009 Contribution of Symmetries
in Condensed Matter Summer School very useful. Villain [30] Symmetry and group
theory throughout physics gives a readable overview. The overheads are here, many of
them are of potential interest (Mourigal recommended).
Canals and Schober [8] Introduction to group theory. It is very concise and precise,
a bastard child of Bourbaki and Hamermesh [13]. Space groups show up only once, on
p. 24: “By working with the cosets we have effectively factored out the translational
part of the problem.”
Ballou [1] An introduction to the linear representations of finite groups appears
rather formal (and very erudite).
Grenier, B. and Ballou [12] Crystallography: Symmetry groups and group repre-
sentations.
The word crystal stems from Greek ‘krustallas’ and means “solidified by the cold.”
Schober [28] Symmetry characterization of electrons and lattice excitations gives
an eminently readable discussion of space groups.
Rodríguez-Carvajal and Bourée [26] Symmetry and magnetic structures
Schweizer [29] Conjugation and co-representation analysis of magnetic structures
deals with black, white and gray groups that Martin tries not to deal with, so all Mouri-
gal groups are gray.
Villain discusses graphene in the Appendix A of Symmetry and group theory through-
out physics [30].
powerful way to think about these problems is Poincaré’s qualitative theory of solu-
tions of differential equations : analyse the geometry of their flows in their state space.
I know for a fact (from a study of cat maps and spatiotemporal cat maps - see links to
talks in ChaosBook.org/overheads/spatiotemporal; the papers are slowly being written
up) that in that case the translational eigenfunctions are hyperbolic sinhes and coshes,
rather than the sines and cosines we are used to as Cn eigenfunctions. For finite dis-
crete symmetries you saw that irreps were fine for linear problems, like coupled arrays
of springs, but symmetry reduction for a nonlinear problem like Lorenz equations re-
quired quite different techniques. For space group symmetries the analogous nonlinear
problems seem still quite unexplored.
century, but the field is undergoing a revival, as the study of topological insulators
requires diving deeper into crystallography than simply looking up the tables.
The translation group T , the set of translations t that put the crystallographic struc-
ture in coincidence with itself, constitutes the lattice. T is a normal subgroup of G. It
defines the Bravais lattice. Translations are of the form
t = tn = n1 a1 + n2 a2 + n3 a3 , nj ∈ Z .
The basis vectors aj span the unit cell. There are 6 simple (or primitive) unit cells
that contain a single point, specified by the lengths of the unit translations a, b, c and
pairwise angles α, β, γ between them. The most symmetric among them is the cubic
cell, with a = b = c and α = β = γ = 90o .
The lattice unit cell (primitive cell) is always a generating region (a tile that tiles the
entire space). The smallest generating region –the fundamental domain– is a minimal
region that generates the whole pattern through its images under all symmetries. At
each lattice point the identical group of “atoms” constitutes the motif. The lattice and
the motif completely characterize the crystal.
The cosets by translation subgroup T (the set all translations) form the factor (AKA
quotient) group G/T , isomorphic to the point group g (rotations). All irreducible rep-
resentations of a space group G can be constructed from irreducible representations of
g and T . This step, however, is tricky, as, due to the non-commutativity of translations
and rotations, the quotient group G/T is not a normal subgroup of the space group G.
The quantum-mechanical calculations are executed by approximating the infinite
crystal by a triply-periodic one, and going go to the reciprocal space by deploying CNj
discrete Fourier transforms. This implements the G/T quotienting by translations and
reduces the calculation to a finite Brilluoin zone. That is the content of the ‘Bloch
theorem’ of condensed matter physics. Further work is then required to reduce the
calculations to the point group irreps.
Point symmetry operations leave at least one point fixed. They are (a) inversion
through a point, (b) rotation around an axis, (c) roto-inversion around an axis and
through a point and (d) reflection through a mirror plane. The rotations have to be
compatible with the translation symmetry: in 3 spatial dimensions they can only be of
orders 1, 2, 3, 4, or 6. They can be proper (det = +1) or improper (det = −1).
The spectroscopists’ Schoenflies notation labels point groups as: cyclic Cn , dihe-
dral Cn0 , tetrahedral T and octahedral O rotation point groups, of order n = 1, 2, 3, 4, 6,
0
respectively. The superscript refers to either v (parallel mirror plane) or h (perpen-
dicular mirror plane). The crystallographer’s preferred classification is, however, the
international crystallographic (Hermann-Mauguin) notation.
(a) (b)
Figure 8.1: The shaded (or yellow) area indicates a fundamental domain, i.e., the
smallest part of the pattern whose repeats tile the entire plane. (a) For the most sym-
metric 2D square lattice, with point group p4mm, the fundamental domain is indicated
by the shaded triangle ΓΛRSX∆Γ which constitutes 1/8 of the Brillouin zone, and
contains the basic wave vectors and the high symmetry points (Fig. 10.2 of Dressel-
haus et al. [11]). (b) For the 2D square lattice with the glide and reflect point group
p4g the fundamental domain is indicated by the yellow triangle (Figure drawn by M.
von Gagern).
forming a normal subgroup. Rotations can only have order 1,2,3,4 or 6 (crystallo-
graphic restriction).
For wallpaper groups the Hermann-Mauguin notation begins with either p or c, for
a primitive cell or a face-centred cell. This is followed by a digit, n, indicating the
highest order of rotational symmetry: 1-fold (none), 2-fold, 3-fold, 4-fold, or 6-fold.
The first, resp. second of the next two symbols indicates the symmetry relative to one
translation axis of the pattern, referred to as the main, resp. second one. The symbols
are either m, g, or 1, for mirror, glide reflection, or none.
Section 9.3 Two-Dimensional Space Groups of Dresselhaus et al. [11] discusses
the most symmetric of the wallpaper groups, the tiling of a plane by squares, which in
the international crystallographic notation is denoted by #11, with point group p4mm.
We work out this space group in exercise 8.2. The largest invariant subgroup of C4v is
C4 . In that case, the space group is p4, or #10. Prefix p indicates that the unit cell is
primitive (not centered). This is a ‘simple’, or symmorphic group, which makes calcu-
lations easier. There is, however, the third, non-symmorphic two-dimensional square
space group p4g or #12 (p4gm), see Table B.10 of ref. [11]. If someone can explain
its ‘Biblical’ diagram to me, I would be grateful. The wiki explanation, reproduced
here as figure 8.1 (b), is the best one that I have found so far, but I’m still scratching
my head:) The Bravais lattice ‘unit cell’ is a square in all three cases. In the crystallo-
graphic literature the ChaosBook’s ‘fundamental domain’ makes an appearance only in
the reciprocal lattice, as the Brilloun zone depicted for p4mm in figure 8.1 (a). How-
ever, the ‘wallpaper groups’ wiki does call ‘fundamental domain’ the smallest part of
the configuration pattern that, when repeated, tiles the entire plane.
The quantum-mechanical calculations are carried out in the reciprocal space, in our
case with the full Γ point, k = 0, wave vector symmetry (see Table 10.1 of ref. [11]),
uij (x) = 1
2 (∂i uj + ∂j ui ) .
Kardar picks the triangular lattice, as its elastic energy is isotropic (i.e., invariant under
lattice rotations, see Landau and Lifshitz [22]). In terms of the Lamé coefficients λ and
µ,
Z
1
βH = d2 x (2µ uij uij + λ uii ujj )
2
Z
1
= − d2 x ui [2µ δij + (µ + λ) ∂i ∂j ] uj . (8.2)
2
(here we have assumed either infinite or doubly periodic lattice, so no boundary terms
from integration by parts), with the equations of motion something like (FIX!)
(Note that Kardar keeps time continuous, but discretizes space. In numerical compu-
tations time is discretized as well.) The symmetry of a square lattice permits an addi-
tional term proportional to ∂x2 u2x + ∂y2 u2y . In general, the number of independent elastic
constants depends on the dimensionality and rotational symmetry of the lattice in ques-
tion. In two dimensions, square lattices have three independent elastic constants, and
triangular lattices are “elastically isotropic” (i.e., elastic properties are independent of
direction and thus have only two [22]).
The Goldstone modes associated with the broken (PC: why “broken”?) transla-
tional symmetry are phonons, the normal modes of vibrations. Eq. (8.3) supports two
types of lattice normal modes, transverse and longitudinal.
The order parameter describing broken translational symmetry is
decreases at finite temperatures, and its correlations decay as hρG (x)ρ∗G (0)i . This is
the order parameter ChaosBook and Gaspard use in deriving formulas for deterministic
diffusion. Kardar computes this in Fourier space by approximating G · q with its
angular average G2 q 2 /2, ignoring the rotationally symmetry-breaking term cos q · x,
and getting only the asymptotics of the correlations right (the decay is algebraic).
The translational correlations are measured in diffraction experiments. The scat-
tering amplitude is the Fourier transform of ρG , and the scattered intensity at a wave-
vector q is proportional to the structure factor. At zero temperature, the structure factor
is a set of delta-functions (Bragg peaks) at the reciprocal lattice vectors.
The orientational order parameter that characterizes the broken rotational symmetry
of the crystal can be defined as
Ψ(x) = e6iθ(x) ,
where θ(x) is the angle between local lattice bonds and a reference axis. The factor of
6 accounts for the equivalence of the 6 possible D3 orientations of the triangular lattice.
(Kardar says the appropriate choice for a square lattice is exp(4iθ(x)) - shouldn’t the
factor be 8, the order of C4v ?) The order parameter has unit magnitude at T = 0, and
is expected to decrease due to fluctuations at finite temperature. The distortion u(x)
leads to a change in bond angle given by
θ(x) = − 12 (∂x uy − ∂y ux ) .
(This seems to be dimensionally wrong? For detailed calculations, see the above Kar-
dar lecture notes.)
(a) (b)
Figure 8.2: (a) The special points and the lines of symmetry in the first Brillouin zone of a
simple cubic lattice define its fundamental domain, the tetrahedron ΓXM R. (b) Just not to get
any ideas that this is easy: the fundamental domain for the first Brillouin zone of a bcc lattice.
(From Joshi [20].)
strange brew: the space is non-compact (homogeneity), while rotations are com-
pact (isotropy). That leads to the conceptually awkward situation of mixing a
group of additions (translations) with a group of multiplications (rotations). To
work with such group we first translate objects to the origin and then rotate them
with the respect to the origin. That’s not nice, because by translation invariance
any point is as good as any other, there is no preferred origin. There is no reason
why one should translate first, rotate second. What one needs is a formalism that
implements translations and rotations on the same footing.
If I understand Hestenes [16] right (also David Finkelstein and perhaps Holger
Beck Nielsen have told me things in this spirit) a way to accomplish that is to
replace the flat translational directions by a compact manifold where translations
and rotations are non-commuting multiplicative group operations.
A part of the Hestenes program is redoing crystallography. I have read Hestenes [17]
paper (but not the Hestenes and Holt [18] follow up). It looks very interesting,
but I will spare you from my comments here, as I do not know how to make this
formalism work for our purposes (character; explicit computations), so I should
not waste your time on that. If you do have a look at his, or at Coxeter [10]
discussion of planar tilings, please do report back to me.
The energy of a free electron is proportional to the square of its momentum, but
not so in a crystal. As this is reminiscent of massless elementary particles like
photons and neutrino’s, it has been renamed since ‘Dirac cones’, but Dirac has
nothing whatsoever to do with that. To learn more, talk to people from the Claire
Berger and Walt De Heer’s group [2] - I have extracted above history of graphene
from De Heer’s notes (the “con man” is my own angle on what went down with
this particular Nobel prize).
References
[1] R. Ballou, “An introduction to the linear representations of finite groups”, EPJ
Web Conf. 22, 00005 (2012).
[2] C. Berger, Z. Song, T. Li, X. Li, A. Y. Ogbazghi, R. Feng, Z. Dai, A. N. Marchenkov,
E. H. Conrad, P. N. First, and W. A. De Heer, “Ultrathin epitaxial graphite: 2D
electron gas properties and a route toward graphene-based nanoelectronics”, J.
Phys. Chem. B 108, 19912–19916 (2004).
[3] H. Boehm, R. Setton, and E. Stumpp, “Nomenclature and terminology of graphite
intercalation compounds”, Carbon 24, 241–245 (1986).
[4] H. P. Boehm., A. Clauss, G. O. Fischer, and U. Hofmann, “Dünnste Kohlenstoff-
Folien”, Z. Naturf. B 17, 150–153 (1962).
[5] L. P. Bouckaert, R. Smoluchowski, and E. P. Wigner, “Theory of Brillouin zones
and symmetry properties of wave functions in crystals”, Phys. Rev. 50, 58–67
(1936).
[6] I. B. Božović, “Irreducible representations of the symmetry groups of polymer
molecules. III. Consequences of time-reversal symmetry”, J. Phys. A 14, 1825
(1981).
[7] I. B. Božović, M. Vujičić, and F. Herbut, “Irreducible representations of the
symmetry groups of polymer molecules. I”, J. Phys. A 11, 2133 (1978).
[8] B. Canals and H. Schober, “Introduction to group theory”, EPJ Web Conf. 22,
00004 (2012).
[9] J. F. Cornwell, Group Theory in Physics: An Introduction (Academic, New York,
1997).
[10] H. S. M. Coxeter, Introduction to Geometry, 2nd ed. (Wiley, New York, 1989).
[11] M. S. Dresselhaus, G. Dresselhaus, and A. Jorio, Group Theory: Application to
the Physics of Condensed Matter (Springer, New York, 2007).
[12] B. Grenier and R. Ballou, “Crystallography: Symmetry groups and group repre-
sentations”, EPJ Web Conf. 22, 00006 (2012).
[13] M. Hamermesh, Group Theory and Its Application to Physical Problems (Dover,
New York, 1962).
[14] C. Herring, “Effect of time-reversal symmetry on energy bands of crystals”,
Phys. Rev. 52, 361–365 (1937).
[15] G. Herzberg, Molecular Spectra and Molecular Structure (Van Nostrand, Prin-
centon NJ, 1950).
[16] D. Hestenes, Space-time Algebra, 2nd ed. (Springer, 1966).
[17] D. Hestenes, “Point groups and space groups in geometric algebra”, in Appli-
cations of Geometric Algebra in Computer Science and Engineering, edited by
L. Dorst, C. Doran, and J. Lasenby (Birkhäuser, Boston, MA, 2002), pp. 3–34.
[18] D. Hestenes and J. W. Holt, “Crystallographic space groups in geometric alge-
bra”, J. Math. Phys. 48, 023514 (2007).
[19] D. F. Johnston, “Group theory in solid state physics”, Rep. Prog. Phys. 23, 66
(1960).
[20] A. W. Joshi, Elements of Group Theory for Physicists (New Age International,
New Delhi, India, 1997).
[21] C. Kittel, Introduction to Solid State Physics, 8th ed. (Wiley, 2004).
[22] L. D. Landau and E. M. Lifshitz, Theory of Elasticity, 3rd ed. (Pergamon Press,
Oxford, 1970).
[23] H. Liang and P. Cvitanović, “A chaotic lattice field theory in one dimension”, J.
Phys. A 55, 304002 (2022).
[24] J. W. Mintmire, B. I. Dunlap, and C. T. White, “Are Fullerene tubules metallic?”,
Phys. Rev. Lett. 68, 631–634 (1992).
[25] J. J. Quinn and K. . Yi, Solid State Physics: Principles and Modern Applications
(Springer, Berlin, 2009).
[26] J. Rodríguez-Carvajal and F. Bourée, “Symmetry and magnetic structures”, EPJ
Web Conf. 22, 00010 (2012).
[27] C. Schafhaeutl, “LXXXVI. On the combinations of carbon with silicon and iron,
and other metals, forming the different species of cast iron, steel, and malleable
iron”, Philos. Mag. Ser. 3 16, 570–590 (1840).
[28] H. Schober, “Symmetry characterization of electrons and lattice excitations”,
EPJ Web Conf. 22, 00012 (2012).
[29] J. Schweizer, “Conjugation and co-representation analysis of magnetic struc-
tures”, EPJ Web Conf. 22, 00011 (2012).
[30] J. Villain, “Symmetry and group theory throughout physics”, EPJ Web Conf. 22,
00002 (2012).
[31] P. R. Wallace, “The band theory of graphite”, Phys. Rev. 71, 622–634 (1947).
Exercises
8.1. Space group.
(a) Show that for any space group, the translations by vectors from Bravais lattice form
a normal subgroup.
(b) Can rotations of the lattice at a fixed point constitute a normal subgroup of a space
group?
(B. Gutkin)
8.2. Band structure of a square lattice. A charged particle (without spin) moves in a po-
tential created by an infinite square lattice of atoms, see figure 8.3.
(a) What are the symmetry groups of the Bravais and reciprocal lattices?
(b) Plot the 1st Brillouin zone. What is its symmetry? What is the corresponding
fundamental domain?
Let k be quasi-momentum and En (k) the energy of the nth band.
(c) At which points of the Brillouin zone is the group G(k) (the group which leaves
vector k invariant) nontrivial? What is it?
(d) What is the symmetry of En (k) as a function of k? At which points of the Brillouin
zone is the group velocity ∇En (k) equal 0?
(e) At which points of the Brillouin zone neighboring bands (generically) stick to each
other? How many bands can stick? Explain from the group theory prospective.
(f) Assume now that the lattice is slightly squeezed along one of the axis. What will
be the new symmetry of the system and its 1st Brillouin zone? Will the sticking
between bands be lifted or persiss?
(B. Gutkin)
8.3. Tight binding model. Verify your solution of exercise 8.2 within the 2-state tight bind-
ing model. Assume that particle can hop either from corner to corner of the square lattice
with coefficient t1 or from corner to the middle of the square with coefficient t2 (and vice
versa).
(a) Show the obtained energy bands Ei (k) as both contour- and 3-dimensional plots.
(b) Compare with the results from exercise 8.2.
Continuous groups
Homework HW9
Bonus points
Exercise 9.5 A fluttering flame front 4 points
Exercise 9.6 O(2) fundamental domain for a PDE (difficult) 10 points
101
GROUP THEORY - WEEK 9. CONTINUOUS GROUPS
lecture playlist
These lectures are about the basic ideas of how one goes from finite groups to the
continuous ones. We have worked one example out in week 2, the discrete Fourier
transform of example 2.6 Projection operators for cyclic group CN . The cyclic group
CN is generated by the powers of the rotation by 2π/N , and in the N → ∞ limit one
only needs to understand
P the algebra of T` , generators of infinitesimal transformations,
D(θ) = 1 + i ` θ` T` . Applied to functions, they turn out to be partial derivatives.
◦ Lie groups, sect. 9.3: Definition of a Lie group; Cyclic group CN → contin-
uous SO(2) plane rotations; Infinitesimal transformations; SO(2) generator of
rotations.
• What’s the payback? While for you the geometrically intuitive representation is
the set of rotation [2×2] matrices, group theory says no! They split into pairs of 1-
dimensional irreps, and the basic building blocks of our 2-dimensional rotations
on our kitchen table (forget quantum mechanics!) are the U(1) [1×1] complex
unit vector phase rotations.
Reading: C. K. Wong Group Theory notes, Chap 6 1D continuous groups,
Sects. 6.1-6.3 Irreps of SO(2). In particular, note that while geometrically
intuitive representation is the set of rotation [2×2] matrices, they split into
pairs of 1-dimensional irreps.
Reading: C. K. Wong Group Theory notes, Chap 6 1D continuous groups,
Sect. 6.6 completes discussion of Fourier analysis as continuum limit of
cyclic groups Cn , compares SO(2), O(2), discrete translations group, and
continuous translations group.
Chen, Ping and Wang [2] Group Representation Theory for Physicists, Sect
5.2 Definition of a Lie group, with examples (click here).
AWH Chapter 17 Group Theory, Sect. 17.7 Continuous groups (click here).
◦ Sect. 9.4 Character orthogonality theorem
Infinitesimal symmetries: Lie derivative (8 min)
Tell no Lie to plumbers (39 sec)
It’s a matter of no small pride for a card-carrying dirt physics theorist to claim
full and total ignorance of group theory (XX min)
Why did we move from orthogonal group O(n) to special orthogonal group
SO(n)? (3:32 min)
Why SU(n) rather than U(n)? (6:30 min)
Why is SU(n) dimension n2 − 1? (56 sec)
What are these “characters”? And why is there a Journal of Linear Algebra,
today? Inconclusive blah blah. (12 min)
Rant 1 - Is beauty symmetry? The first piece of art found in China is a perfect
disk carved out of jade. All of Bach is symmetries. (9 min)
Rant 2 - students find letter A beautifully symmetric, but Predrag finds zero ‘O’
the most beautiful grade. (1 min)
Rant 3 - SO(3) & SU(2) preview and a long rant - listen to it at your own risk.
Roger Penrose thoughts on quantum spacetime and quantum brain. Are laws of
physics time invariant? Waiting for dark energy to go away. Arrow of time. (17
min)
Rant 4 - SO(3) & SU(2) preview and a long rant - listen to it at your own risk.
Get this: math uses 2d complex vectors (spinors) to build our real 3d space. And
all we see - starlight, graphene, greenhouse effect, helioseismography, gravita-
tional wave detectors - it is all irreps! (12 min)
Rant 5 - Help me, I’m bullied by a mathematician. (3 min)
Rant 6 - you can always count on Prof. Z. (1/2 min)
6. How does a vector space come about from a set of traces, each of which I normally think
of as just a number, like the determinant? And finally,
7. How can we use our knowledge of classes/character to find irreducible representations,
since that seems to be an important goal in examining a group.
8. Exercise 4.2 (c) says to find the characters for this representation, which seems to imply
that character depends on representation. But I would’ve thought that character, which is
a trace of a matrix, is invariant under any similarity transform, which is how you get from
a reducible representation to an irreducible representation.
9. Do the multiplicities of irreducible representations correspond to the multiplicity of char-
acters (i.e. the number of elements in each class)? If so, why? (Or if not, why not?)
10. The same thing for classes, correct? Classes shouldn’t depend on representation b/c they
can be thought of as corresponding to a physical operation (e.g. transposition or cyclic
permutation), something which is independent of basis.
.
A Great framing for a discussion, thanks! I’ll probably reedit this post several times, every-
body’s input is very welcome. Items numbered as in above:
(2) My favorite step-by-step, pedagogical exposition are the chapters 2 Representation The-
ory and Basic Theorems and 3 Character of a Representation of Dresselhaus et al. [6].
There is too much material for our course, but if you want to understand it once for all
times, it’s worth your time.
(3) Correct.
(4) Correct. Note, however, that while every matrix representation has a trace, and thus a
character, you want to decompose this character into the sum of irrep characters, as it is
obvious after the block diagonalization has been attained.
(5) The unitary diagonalization matrix, whose entries are characters, takes character-weighted
sums of classes in order to project them onto irreps, just like what the Fourier represen-
tation does. The result (as we know from projection operators analysis), are mutually
orthogonal sub-spaces.
(6) Whenever you do not understand something about finite groups, ask yourself - how does
it work for finite lattice Fourier representation?
There the vector space comes via a unitary transformation from the configuration coor-
dinates (where each group element is represented by a full matrix) to the diagonalized,
irreducible subspaces coordinates (Fourier modes).
The unitary F matrix is full of ω ij , ie, characters of the cyclic group Cn . That’s where
the characters come from.
Now mess up C3 by adding a reflection. Dihedral group D3 , the group of rotations and re-
flections, has more symmetry constrains, it cannot have 6 irreps, as reflection invariance
mixes together the two senses of rotation. Now there are 3 classes, ie, kinds of things
the group does: nothing, flip, rotate. The unitary transformation that diagonalizes group
element matrices is now morally a smaller unitary [3 × 3] matrix from ‘classes’ in config-
uration space to ‘irreps’ in the diagonalized representation, where some sub-spaces must
have dimension higher than one.
The surprise, for me, is that the entries in the unitary diagonalization matrix can still be
written as traces of irreps, ie, characters. For me it is a calculation, a beautiful example
of mathematics leading us somewhere where our intuition falls short. If you find a good
intuitive explanation somewhere, please let us all know.
(7) That’s automatic, now. Each irrep has the projection operator associated with it; we
construct it as a sub-product of factors in Hamilton-Cayley formula. Now we know we
can write it -just as we did with the Fourier representation- as sum over all class group
actions, each weighted by a the irrep’s character.
(8) Characters are elements of the unitary matrix with one index running over classes, the
other over irreps. So you expect character to differ from representation to representation;
very clear from D3 character table. As always, you already know that from the Fourier
representation example.
(9) They do not. Dresselhaus et al. [6] has the answer - enter it here once you understand it.
(10) Correct.
These thoughts are spread over chapters of my book Group Theory - Birdtracks,
Lie’s, and Exceptional Groups [5] that you can steal from my website, but the book
itself is too sophisticated for this course.
All that is needed in practice is knowledge of the characteristic equation for the
invariant matrix M (see sect. 9.5.1). The characteristic equation is usually a simple
consequence of the algebraic relations satisfied by the primitive invariants, and the
eigenvalues λi are easily determined. The λi ’ s determine the projection operators
Pi , which in turn contain all relevant spectroscopic information: the rep dimension is
given by tr Pi , and the casimirs, 6-j’s, crossing matrices, and recoupling coefficients)
are traces of various combinations of Pi ’s. All these numbers are combinatoric; they
can often be interpreted as the number of different colorings of a graph, the number of
singlets, and so on.
The key results are the construction of projection operators from invariant matri-
ces, the Clebsch-Gordan coefficients rep of projection operators (9.26), the invariance
conditions (10.41) and the Lie algebra relations (10.45).
The basic idea is simple: a hermitian matrix can be diagonalized. If this matrix
is an invariant matrix, it decomposes the reps of the group into direct sums of lower-
dimensional reps. Most of computations to follow implement the spectral decomposi-
tion
M = λ1 P1 + λ2 P2 + · · · + λr Pr ,
which associates with each distinct root λi of invariant matrix M a projection operator
(9.5):
Y M − λj 1
Pi = .
λi − λj
j6=i
The exposition given here in sect. 9.5.1 is taken from refs. [3, 4]. Who wrote this down
first I do not know, but I like Harter’s exposition [9–11] best.
λ1 ... 0
.. 0 0
.
0 . . . λ1
λ2 0 ... 0
†
CMC =
0 λ2 .
(9.3)
0 .. .. .. 0
. . .
0 ... λ2
λ3 ...
0 0
.. ..
. .
λ1 − λ2
λ1 − λ2
λ1 − λ2
0
. ..
,
0
λ3 − λ2
λ3 − λ2
..
.
and so on, so the product over all factors (M − λ2 1)(M − λ3 1) . . . , with exception of
the (M − λ1 1) factor, has nonzero entries only in the subspace associated with λ1 :
1 0 0
0 1 0 0
0 0 1
0
Y Y
†
C (M − λj 1)C = (λ1 − λj ) .
j6=1 j6=1
0
0
0
..
.
In this way, we can associate with each distinct root λi a projection operator Pi ,
Y M − λj 1
Pi = , (9.5)
λi − λj
j6=i
which acts as identity on the ith subspace, and zero elsewhere. For example, the pro-
jection operator onto the λ1 subspace is
1
..
.
1
†
P1 = C 0 C .
(9.6)
0
..
.
0
di = tr Pi . (9.9)
It follows from the characteristic equation (9.4) and the form of the projection operator
(9.5) that λi is the eigenvalue of M on Pi subspace:
Hence, any matrix polynomial f (M) takes the scalar value f (λi ) on the Pi subspace
This, of course, is the reason why one wants to work with irreducible reps: they reduce
matrices and “operators” to pure numbers.
[M1 , M2 ] = 0 , (9.12)
Pi Pj = Pj Pi . (9.13)
(i) (i)
X
[M2 , M1 ] = λj [M2 , Pj ] = 0 . (9.14)
j
(i)
Now the characteristic equation for M2 (if nontrivial) can be used to decompose Vi
subspace.
In the diagonalized rep (9.6), the matrix G has a block diagonal form:
G1 0 0
CGC † = 0 G2 0 ,
X
G= C i Gi Ci . (9.17)
.. i
0 0 .
of the two terms in the diagonal representation of a projection operator (9.6). This
matrix has nonzero entries only in the dλ rows of subspace Vλ . We collect them in a
[dλ × d] rectangular matrix (Cλ )α
σ , α = 1, 2, . . . d, σ = 1, 2, . . . dλ :
The index α in (Cλ )ασ stands for all tensor indices associated with the d = n
p+q
-
p q
dimensional tensor space V ⊗ V̄ . In the birdtrack notation these indices are explicit:
b1
λ
(Cλ )σ , abpq...b1
= . (9.21)
...
...
...a2 a1
aq
(C λ )11 (C λ )d1λ
...
.. .. d
Cλ = . .
(C λ )ddλ
| {z }
dλ
b1
b2
a a ...aq σ λ σ
...
(C λ )b11...b
..
2
p
, = . (9.23)
..
aq
The tensors transform as tensor products of the defining rep. In general, tensors
transform as tensor products of various reps, with indices running over the correspond-
ing rep dimensions:
a1 = 1, 2, . . . , d1
a2 = 1, 2, . . . , d2
..
xaap+1 ...ap+q
1 a2 ...ap
where . (9.24)
ap+q = 1, 2, . . . , dp+q .
The action of the transformation g on the index ak is given by the [dk × dk ] matrix rep
Gk .
Clebsches are notoriously index overpopulated, as they require a rep label and a
tensor index for each rep in the tensor product. Diagrammatic notation alleviates this
index plague in either of two ways:
By the definition of clebsches (9.6), the λ rep projection operator can be written out in
terms of Clebsch-Gordan matrices C λ Cλ :
C λ Cλ = Pλ , (no sum on i)
a a2 ...ap α a a2 ...dp dq ...d1
(C λ )b11...bq
, (Cλ )α , cpd...c
q ...d1
2 c1
= (Pλ )b11...bq
, cp ...c2 c1 (9.26)
λ
= Pλ .
...
...
...
...
Cλ C µ = δλµ 1 ,
a a ...ap
(Cλ )β , b11...b
2
q
(C µ )apbq...a
...b1
2 a1
,α = δβα δλµ
λ μ λ μ
= . (9.27)
...
Here 1 is the [dλ × dλ ] unit matrix, and Cλ ’s are multiplied as [dλ × d] rectangular
matrices.
The completeness relation (9.8)
X
C λ Cλ = 1, ([d × d] unit matrix) ,
λ
a a2 ...ap α d
X
(C λ )b11...bq
, (Cλ )α , dcpq ...c
...d1
2 c1
= δca11 δca22 . . . δbqq
λ
X λ
= (9.28)
...
...
...
λ
C λ Pµ = δλµ C λ ,
Pλ C µ = δλµ C µ , (no sum on λ, µ) , (9.29)
in the defining space defines the little group of transformations in the remaining di-
rections), but the tensorial technology in this context is underdeveloped compared to
the canonical methods. And, as Stedman [13] rightly points out, if you need to check
your calculations against the existing literature, keeping track of phase conventions is
a necessity.
References
[1] N. B. Budanur and P. Cvitanović, “Unstable manifolds of relative periodic orbits
in the symmetry-reduced state space of the Kuramoto-Sivashinsky system”, J.
Stat. Phys. 167, 636–655 (2015).
[2] J.-Q. Chen, J. Ping, and F. Wang, Group Representation Theory for Physicists
(World Scientific, Singapore, 1989).
[3] P. Cvitanović, “Group theory for Feynman diagrams in non-Abelian gauge the-
ories”, Phys. Rev. D 14, 1536–1553 (1976).
[4] P. Cvitanović, Classical and exceptional Lie algebras as invariance algebras, Ox-
ford Univ. preprint 40/77, unpublished., 1977.
[5] P. Cvitanović, Group Theory: Birdtracks, Lie’s and Exceptional Groups (Prince-
ton Univ. Press, Princeton NJ, 2008).
[6] M. S. Dresselhaus, G. Dresselhaus, and A. Jorio, Group Theory: Application to
the Physics of Condensed Matter (Springer, New York, 2007).
[7] L. Frappat, P. Sorba, and A. Sciarrino, Dictionary on Lie algebras and superal-
gebras (Academic Press, 2000).
[8] M. Gell-Mann, The Eightfold Way: A Theory of Strong Interaction Symmetry
(CalTech, 1961).
[9] W. G. Harter, “Algebraic theory of ray representations of finite groups”, J. Math.
Phys. 10, 739–752 (1969).
[10] W. G. Harter, Principles of Symmetry, Dynamics, and Spectroscopy (Wiley, New
York, 1993).
[11] W. G. Harter and N. dos Santos, “Double-group theory on the half-shell and the
two-level system. I. Rotation and half-integral spin states”, Amer. J. Phys. 46,
251–263 (1978).
[12] R. Slansky, “Group Theory for Unified Model Building”, Phys. Rep. 79, 1–128
(1981).
[13] G. E. Stedman, Diagram Techniques in Group Theory (Cambridge U. Press,
Cambridge, 1990).
Exercises
9.1. Irreps of SO(2). Matrix
0 −i
T = (9.30)
i 0
is the generator of rotations in a plane.
(a) Use the method of projection operators to show that for rotations in the kth Fourier
mode plane, the irreducible 1D subspaces orthonormal basis vectors are
1 (k) (k)
e(±k) = √ ±e1 − i e2 .
2
In this equation t ≥ 0 is the time and x is the spatial coordinate. The subscripts x
and t denote partial derivatives with respect to x and t: ut = ∂u/d∂, uxxxx stands
for the 4th spatial derivative of u = u(x, t) at position x and time t. Consider the
form of equations under coordinate shifts x → x + ` and reflection x → −x. What
is the symmetry group of (9.31)?
(b) Expand u(x, t) in terms of its SO(2) irreducible components (hint: Fourier expan-
sion) and rewrite (9.31) as a set of linear ODEs for the expansion coefficients. What
are the eigenvalues of the time evolution operator? What is their degeneracy?
(c) Expand u(x, t) in terms of its O(2) irreducible components (hint: Fourier expan-
sion) and rewrite (9.31) as a set of linear ODEs. What are the eigenvalues of the
time evolution operator? What is their degeneracy?
(d) Interpret u = u(x, t) as a ‘flame front velocity’ and add a quadratic nonlinearity to
(9.31),
ut + 12 (u2 )x + uxx + νuxxxx = 0 , x ∈ [0, L] . (9.32)
This nonlinear equation is known as the Kuramoto-Sivashinsky equation, a baby
cousin of Navier-Stokes. What is the symmetry group of (9.32)?
(e) Expand u(x, t) in terms of its O(2) irreducible components (see exercise 9.3) and
rewrite (9.32) as an infinite tower of coupled nonlinear ODEs.
(f) What are the degeneracies of the spectrum of the eigenvalues of the time evolution
operator?
9.6. O(2) fundamental domain for Kuramoto-Sivashinsky equation. You have C2
discrete symmetry generated by flip σ, which tiles the space by two tiles.
• Is there a subspace invariant under this C2 ? What form does the tower of ODEs
take in this subspace?
• How would you restrict the flow (the integration of the tower of coupled ODEs) to
a fundamental domain?
This problem is indeed hard, a research level problem, at least for me and the grad students
in our group. Unlike the beautiful full-reducibility, character-orthogonality representation
theory of linear problems, in nonlinear problems symmetry reduction currently seems to
require lots of clever steps and choices of particular coordinates, and we am not at all sure
that our solution is the optimal one. Somebody looking at the problem with a fresh eye
might hit upon a solution much simpler than ours. Has happened before :)
Burak Budanur’s solution is written up in Budanur and Cvitanović [1] Unstable mani-
folds of relative periodic orbits in the symmetry-reduced state space of the Kuramoto-
Sivashinsky system sect. 3.2 O(2) symmetry reduction, eq. (17) (get it here).
9.7. Lie algebra from invariance. Derive the Lie algebra commutator and the Jacobi
identity as particular examples of the invariance condition, using both index and birdtracks
notations. The invariant tensors in question are “the laws of motion,” i.e., the generators
of infinitesimal group transformations in the defining and the adjoint representations.
Homework HW10
Bonus points
Exercise 10.5 Integrate the two-modes system 4 point
Exercise 10.7 Invariant subspace of the two-modes system 1 point
Exercise 10.8 Slicing the two-modes system 1 point
Exercise 10.9 The symmetry reduced two-modes flow (difficult) 6 points
119
GROUP THEORY - WEEK 10. LIE GROUPS, ALGEBRAS
There is way too much material in this week’s notes. Watch the main sequence of video
clips, that and recommended reading should suffice. The rest is optional.
• Bridging the step from discrete to continuous compact groups: invariant integra-
tion measures, characters, character orthonormality and completeness relations:
Rotations in 3 dimensions (30 min)
Lie algebra (21 min)
Birdtracks (6 min)
◦ Sect. 10.7 Lie groups for pedestrians is advanced material, safely ignored,
here only to whet your appetite for things not done in 19th century. It
is a very condensed extract of chapters 3 Invariants and reducibility and
4 Diagrammatic notation from Group Theory - Birdtracks, Lie’s, and Ex-
ceptional Groups [11]. I am usually reluctant to use birdtrack notations in
front of graduate students indoctrinated by their professors in the 1890’s
tensor notation, but now I’m emboldened by the very enjoyable article on
The new language of mathematics by Dan Silver [22].
◦ Ditto for sect. 10.10 Birdtracks - updated history. Your professor’s notation
is as convenient for actual calculations as -let’s say- long division using
roman numerals. So leave them wallowing in their early progressive rock
of 1968, King Crimsons of their youth. You chill to beats younger than
Windows 98, to grime, to trap, to hardvapour, to birdtracks.
◦ Go to week 16 to learn more.
• OK, I see that formally SU(2) ' SO(3), but who ordered ”spin?”
Rotations in 2 complex dimensions (42 min)
◦ Read sect. 10.3 SU(2) Pauli matrices
For overall clarity and pleasure of reading, I like Schwichtenberg [21]
(click here) discussion best. If you read anything for this week’s lectures,
read Schwichtenberg sects. 3.4 to 3.6.
This is difficult material, so it is OK if you do not get it this time around. None
of this will be on the final - the main point is that once you face a nonlinear
problem, nothing is easy - not even rotations on a circle.
ŷ2
y2
1.0 0.2
0.0
0.0 −0.2
−1.0
−1.5 −1.5 −1.5
0
0.75
x2 0
1.5 1.5
0
x1 x̂2 0.75 0 1.5 x̂ 1
(a) (b)
Figure 10.1: Two-modes flow before (a) and after (b) symmetry reduction by first
Fourier mode slice. Here a long trajectory (red and blue) starting on the unstable man-
ifold of the TW1 (red), until it falls on to the strange attractor (blue) and the shortest
relative periodic orbit 1 (magenta). Note that the relative equilibrium becomes an equi-
librium, and the relative periodic orbit becomes a periodic orbit after the symmetry
reduction.
The model has an unreasonably high number of parameters. After some experi-
mentation we fix or set to zero various parameters, and in the numerical examples that
follow, we settle for parameters set to
µ1 = −2.8 , µ2 = 1 , e1 = 0 , e2 = 1 ,
a1 = −1 , a2 = −2.66 , b1 = 0 , b2 = 0 , c1 = −7.75 , c2 = 1 , (10.3)
unless explicitly stated otherwise. For these parameter values the system exhibits
chaotic behavior. Experiment! If you find a more interesting behavior for some other
parameter values, please let us know. The simplified system of equations can now be
written as a 3-parameter {µ1 , c1 , a2 } two-modes system,
ż1 = µ1 z1 − z1 |z1 |2 + c1 z 1 z2
ż2 = (1 − i) z2 + a2 z2 |z1 |2 + z12 . (10.4)
In order to numerically integrate and visualize the flow, we recast the equations in real
variables by substitution z1 = x1 + i y1 , z2 = x2 + i y2 . The two-modes system (10.2)
is now a set of four coupled ODEs
exercise 10.5
ẋ1 = (µ1 − r2 ) x1 + c1 (x1 x2 + y1 y2 ) , r2 = x21 + y12
ẏ1 = (µ1 − r2 ) y1 + c1 (x1 y2 − x2 y1 )
ẋ2 = x2 + y2 + x21 − y12 + a2 x2 r2
ẏ2 = −x2 + y2 + 2 x1 y1 + a2 y2 r2 . (10.5)
Try integrating (10.5) with random initial conditions, for long times, times much
beyond which the initial transients have died out. What is wrong with this picture?
Figure 10.4 (a) is a mess. As we show here, the attractor is built up by a nice ‘stretch
& fold’ action, hidden from the view by the continuous symmetry induced drifts. That
exercise 10.6
is fixed by ‘quotienting’ model’s SO(2) symmetry, and reducing the dynamics to a
exercise 10.7
3-dimensional symmetry-reduced state space, figure 10.4 (b).
exercise 10.8
where the unit vector n determines the plane and the direction of the rotation by angle
ϕ. Here L1 , L2 , L3 are the generators of rotations along x, y, z axes respectively,
0 0 0 0 0 1 0 −1 0
L1 = i 0 0 1 , L2 = i 0 0 0 , L3 = i 1 0 0 , (10.8)
0 −1 0 −1 0 0 0 0 0
with eiφn2 ·L and e−iθn2 ·L mapping the vector n1 to n3 and back, so that the rotation
around axis n1 by angle θ is mapped to a rotation around axis n3 by the same θ. The
conjugacy classes of SO(3) thus consist of rotations by the same angle about all distinct
rotation axes, and are thus labelled the angle θ. As the conjugacy class depends only on
exercise 10.2
θ, the characters can only be a function of θ. For the 3-dimensional special orthogonal
representation, the character is
χ = 2 cos(θ) + 1 . (10.11)
sin(j + 1/2)θ
χ(j) (θ) = (10.12)
sin(θ/2)
To check that these characters are orthogonal to each other, one needs to define
the group integration over a parametrization of the SO(3) group manifold. A group
element is parametrized by the rotation axis n and the rotation angle θ ∈ (−π, π] ,
with n a unit vector which ranges over all points on the surface of a unit ball. Note
however, that a π rotation is the same as a −π rotation (n and −n point along the
same direction), and the n parametrization of SO(3) is thus a 2-dimensional surface of
a unit-radius ball with the opposite points identified.
The Haar measure for SO(3) requires a bit of work, here we just note that after the
integration over the solid angle (characters do not depend on it), the Haar measure is
dθ dθ
dg = dµ(θ) = (1 − cos(θ)) = sin2 (θ/2) . (10.13)
2π π
With this measure the characters are orthogonal, and the character orthogonality the-
exercise 10.3
orems follow, of the same form as for the finite groups, but with the group averages
replaced by the continuous, parameter dependant group integrals
Z
1 X
→ dg .
|G| G
g∈G
The good news is that, as explained in ChaosBook.org Chap. Relativity for cyclists
(and in Group Theory - Birdtracks, Lie’s, and Exceptional Groups [11]), one never
needs to actually explicitly construct a group manifold parametrizations and the corre-
sponding Haar measure.
(a) vector addition ”+” is defined in V such that V is an abelian group under addi-
tion, with identity element 0;
(b) the set is closed with respect to scalar multiplication and vector addition
a(x + y) = ax + ay , a, b ∈ F , x, y ∈ V
(a + b)x = ax + bx
a(bx) = (ab)x
1x = x, 0x = 0. (10.14)
Here the field F is either R, the field of reals numbers, or C, the field of complex
numbers. Given a subset V0 ⊂ V , the set of all linear combinations of elements of V0 ,
or the span of V0 , is also a vector space.
Dual space, dual basis. Under a general linear transformationPg ∈ GL(n, F), the
row of basis vectors transforms by right multiplication as e(j) = k (g−1 )j k e(k) , and
the column of xa ’s transforms by left multiplication as x0 = gx. Under left multiplica-
tion the column (row transposed) of basis vectors e(k) transforms as e(j) = (g† )j k e(k) ,
where the dual rep g† = (g−1 )> is the transpose of the inverse of g. This observation
motivates introduction of a dual representation space V̄ , the space on which GL(n, F)
acts via the dual rep g† .
Definition. If V is a vector representation space, then the dual space V̄ is the set of all
linear forms on V over the field F.
If {e(1) , · · · , e(d) } is a basis of V , then V̄ is spanned by the dual basis {e(1) , · · · , e(d) },
the set of d linear forms e(k) such that
where δjk is the Kronecker symbol, δjk = 1 if j = k, and zero otherwise. The compo-
nents of dual representation space vectors ȳ ∈ V̄ will here be distinguished by upper
indices
(y 1 , y 2 , . . . , y n ) . (10.15)
They transform under GL(n, F) as
y 0a = (g† )a b y b . (10.16)
(a) the set is closed with respect to multiplication T · T → T , so that for any two
elements tα , tβ ∈ T , the product tα · tβ also belongs to T :
r−1
X
tα · tβ = ταβ γ tγ , ταβ γ ∈ C ; (10.18)
γ=0
(tα + tβ ) · tγ = tα · tγ + tβ · tγ
tα · (tβ + tγ ) = tα · tβ + tα · tγ .
The set of numbers ταβ γ are called the structure constants. They form a matrix rep of
the algebra,
(tα )β γ ≡ ταβ γ , (10.19)
(tα · tβ ) · tγ = tα · (tβ · tγ ) ,
the algebra is associative. Typical examples of products are the matrix product
(tα · tβ )ca = (tα )ba (tβ )cb − (tα )bc (tβ )ab , tα ∈ V ⊗ V̄ (10.21)
Definition: A Lie group is a topological group G such that (i) G has the structure of
a smooth differential manifold, and (ii) the composition map G × G → G : (g, h) →
gh−1 is smooth, i.e., C∞ differentiable.
Do not be mystified by this definition. Mathematicians also have to make a living.
The compact Lie groups that we will deploy here are a generalization of the theory of
SO(2) ' U(1) rotations, i.e., Fourier analysis. By a ‘smooth differential manifold’
one means objects like the circle of angles that parameterize continuous rotations in a
plane, figure 10.2, or the manifold swept by the three Euler angles that parameterize
SO(3) rotations.
By ‘compact’ one means that these parameters run over finite ranges, as opposed
to parameters in hyperbolic geometries, such as Minkowsky SO(3, 1). The groups we
focus on here are compact by default, as their representations are linear, finite-dimen-
sional matrix subgroups of the unitary matrix group U(d).
Example 1. Circle group. A circle with a direction, figure 10.2, is invariant under rota-
tion by any angle θ ∈ [0, 2π), and the group multiplication corresponds to composition
of two rotations θ1 + θ2 mod 2π. The natural representation of the group action
is by a complex numbers of absolute value 1, i.e., the exponential eiθ . The composi-
tion rule is then the complex multiplication eiθ2 eiθ1 = ei(θ1 +θ2 ) . The circle group is
a continuous group, with infinite number of elements, parametrized by the continuous
parameter θ ∈ [0, 2π). It can be thought of as the n → ∞ limit of the cyclic group Cn .
Note that the circle divided into n segments is compact, in distinction to the infinite
lattice of integers Z, whose limit is a line (noncompact, of infinite length).
Figure 10.2: Circle group S 1 = SO(2), the symmetry group of a circle with directed
rotations, is a compact group, as its natural parametrization is either the angle φ ∈
[0, 2π), or the perimeter x ∈ [0, L).
g
x’ x
Manifold traced out by action of G
for all possible group elements g
(a) (b)
Figure 10.3: (a) Lie algebra fields {t1 , · · · , tN } span the tangent space of the group
orbit Mx at state space point x, see (10.26) (figure from WikiMedia.org). (b) A global
group transformation g : x → x0 can be pieced together from a series of infinitesimal
steps along a continuous trajectory connecting the two points. The group orbit of state
space point x ∈ Rd is the N -dimensional manifold of all actions of the elements of
group G on x.
norm, hTa x|xi + hx|Ta xi = 0 , hence the Lie algebra generators Ta are hermitian for,
Ta† = Ta . (10.25)
The flow field at the state space point x induced by the action of the group is given by
the set of N tangent fields
ta (x)i = (Ta )ij xj , (10.26)
which span the d-dimensional group tangent space at state space point x, parametrized
by δφ.
For continuous groups the Lie algebra, i.e., the algebra spanned by the set of N gen-
erators Ta of infinitesimal transformations, takes the role that the |G| group elements
play in the theory of discrete groups (see figure 10.3).
10.7.1 Invariants
One constructs the irreps of finite groups by identifying matrices that commute with
all group elements, and using their eigenvalues to decompose arbitrary representation
of the group into a unique sum of irreps. The same strategy works for the compact Lie
groups, (10.30), and is indeed the key idea that distinguishes the invariance groups clas-
sification developed in Group Theory - Birdtracks, Lie’s, and Exceptional Groups [11]
from the 19th century Cartan-Killing classification of Lie algebras.
q = Gq . (10.27)
Ma b = Ga c Gb d Mc d (10.29)
is an invariant matrix. Multiplying with Gb e and using the unitary, we find that the
invariant matrices commute with all transformations g ∈ G:
[G, M] = 0 . (10.30)
Definition. An invariance group G is the set of all linear transformations (10.28) that
preserve the primitive invariant relations (and, by extension, all invariant relations)
Unitarity guarantees that all contractions of primitive invariant tensors, and hence all
composed tensors h ∈ H, are also invariant under action of G. As we assume unitary
G, it follows that the list of primitives must always include the Kronecker delta.
then G is the full unitary group U (n) (invariance group of the complex norm |x|2 =
xb xa δba ), whose elements satisfy
G† G = 1 . (10.33)
Example 3. If we wish the z-direction to be invariant in our 3-dimensional space,
q = (0, 0, 1) is an invariant vector (10.27), and the invariance group is O(2), the group
of all rotations in the x-y plane.
allowed transformation; indeed, one of the main objectives of group theory is to define
the class of allowed transformations.
This subspace is called the adjoint space, and its special role warrants introduction
of special notation. We shall refer to this vector space by letter A, in distinction to
the defining space V . We shall denote its dimension by N , label its tensor indices by
i, j, k . . ., denote the corresponding Kronecker delta by a thin, straight line,
δij = , i, j = 1, 2, . . . , N , (10.36)
11
00 11
00
00
11 0011 00111100 11
00
00
11 1100 00
11
00
11 T 1100 = 11001100 + 00111100 − , (10.41)
00
11 1100
00
11 1100 00
11 00
11
00
11 1100 00
11 00
11
with a relative minus sign between lines flowing in opposite directions. The reader will
recognize this as the Leibnitz rule.
The invariance conditions take a particularly suggestive form in the birdtrack nota-
tion. Equation (10.39) amounts to the insertion of a generator into all external legs of
the diagram corresponding to the invariant tensor q:
0= + −
+ − . (10.42)
The insertions on the lines going into the diagram carry a minus sign relative to the
insertions on the outgoing lines.
As the simplest example of computation of the generators of infinitesimal transfor-
mations acting on spaces other than the defining space, consider the adjoint rep. Where
does the ugly word “adjoint” come from in this context is not obvious, but remember
it this way: this is the one distinguished representation, which is intrinsic to the Lie
algebra, with the explicit matrix elements (Ti )jk of the adjoint rep given by the fully
antisymmetric structure constants iCijk of the algebra (i.e., its multiplication table un-
der the commutator product). It’s the continuous groups analogue of the multiplication
table, or the regular representation for the finite groups. The factor i ensures their re-
ality (in the case of hermitian generators Ti ), and we keep track of the overall signs by
always reading indices counterclockwise around a vertex:
− iCijk = (10.43)
0110
1010
=− 11001010 . (10.44)
As all other invariant tensors, the generators must satisfy the invariance conditions
(10.42):
0=− + − .
Redrawing this a little and replacing the adjoint rep generators (10.43) by the structure
constants, we find that the generators obey the Lie algebra commutation relation
i j
− = (10.45)
Ti Tj − Tj Ti = iCijk Tk . (10.46)
0= + + .
− = . (10.47)
This is the Lie algebra commutator for the adjoint rep generators, known as the Jacobi
relation for the structure constants
Hence, the Jacobi relation is also an invariance statement, this time the statement that
the structure constants are invariant tensors.
opening into general relativity and set him on the path to his 1965 discovery celebrated
by this year’s prize.
A half of the 2020 Nobel Prize in Physics was awarded jointly to Reinhard Genzel
and Andrea Ghez for the discovery of –in Ghez’s words- ”The Monster at the heart of
the Milky Way," a black hole whose existence had been hypothesized since the early
1970s. In order to visually observe an object that famously does not emit any light,
precise measurements of stars moving in the black hole’s gravitational field had to be
carried out. The independent work of Genzel and Ghez mapping the positions of these
stars over many years has led to the clearest evidence yet that the center of our Milky
Way galaxy contains “The Monster”, that possibly every galaxy contains a black hole,
and that the environment near it looks nothing like what was expected.
Nobel Lecture: Roger Penrose, Nobel Prize in Physics 2020 (34 min)
Nobel Prize in Physics 2020 (56 min)
Abstract
Penrose slides for Predrag’s 1/2 of the presentation
2020 Nobel Prizes in Chemistry and Physics, Explained
Roger Penrose gets Nobel Prize. How David Ritz Finkelstein and Roger Penrose
met, and exchanged their lives’ paths.
Negative dimensions (6 min)
Andrea Ghez: ”The Monster at the Heart of our Galaxy"
Veritasium: ”The Infinite Pattern That Never Repeats"
France, Klein decided that it would be safer to return to Germany; Lie also decided to
go home to Norway. However (in a move that I think questions his geometric abilities),
Lie decided that to go from Paris to Norway, he would walk to Italy (and then presum-
ably take a ship to Norway). The trip did not go as Lie had planned. On the way, Lie
ran into some trouble–first some rain, and he had a habit of taking off his clothes and
putting them in his backpack when he walked in the rain (so he was walking to Italy
in the nude). Second, he ran into the French military (quite possibly while walking
in the nude) and they discovered in his sack (in addition to his hopefully dry clothing)
letters written to Klein in German containing the words ’lines’ and ’spheres’ (which the
French interpreted as meaning ’infantry’ and ’artillery’). Lie was arrested as a (insane)
German spy. However, due to intervention by Gaston Darboux, he was released four
weeks later and returned to Norway to finish his doctoral dissertation.”
algorithm is inefficient - the expression balloon quickly, the Young projection operators
soon become unwieldy and impractical, if not impossible to implement.
The Alcock-Zeilinger algorithm, based on the simplification rules of ref. [3], leads
to explicitly Hermitian and drastically more compact expressions for the projection op-
erators than the Keppeler-Sjödahl algorithm [14]. Alcock-Zeilinger fully supersedes
Cvitanović’s formulation, and any future full exposition of reduction of SU(N ) ten-
sor products into irreducible representations should be based on the Alcock-Zeilinger
algorithm.
References
[1] G. ’t Hooft, “A planar diagram theory for strong interactions”, Nucl. Phys. B 72,
461–473 (1974).
[2] A. C. Aitken, Determinants & Matrices (Oliver & Boyd, Edinburgh, 1939).
[3] J. Alcock-Zeilinger and H. Weigert, “Compact Hermitian Young projection op-
erators”, J. Math. Phys. 58, 051702 (2017).
[4] G. B. Arfken, H. J. Weber, and F. E. Harris, Mathematical Methods for Physi-
cists: A Comprehensive Guide, 7th ed. (Academic, New York, 2013).
[5] R. Brauer, “On algebras which are connected with the semisimple continuous
groups”, Ann. Math. 38, 857 (1937).
[6] G. P. Canning, “Diagrammatic group theory in quark models”, Phys. Rev. D 18,
395–410 (1978).
[7] J.-Q. Chen, J. Ping, and F. Wang, Group Representation Theory for Physicists
(World Scientific, Singapore, 1989).
[8] J. F. Cornwell, Group Theory in Physics: An Introduction (Academic, New York,
1997).
[9] P. Cvitanović, “Group theory for Feynman diagrams in non-Abelian gauge the-
ories”, Phys. Rev. D 14, 1536–1553 (1976).
[10] P. Cvitanović, Field Theory, Notes prepared by E. Gyldenkerne (Nordita, Copen-
hagen, 1983).
[11] P. Cvitanović, Group Theory: Birdtracks, Lie’s and Exceptional Groups (Prince-
ton Univ. Press, Princeton NJ, 2008).
[12] P. Cvitanović and A. D. Kennedy, “Spinors in negative dimensions”, Phys. Scr.
26, 5–14 (1982).
[13] H. Elvang, P. Cvitanović, and A. Kennedy, “Diagrammatic Young projection
operators for U(n)”, J. Math. Phys. 46, 043501 (2005).
[14] S. Keppeler and M. Sjödahl, “Hermitian Young operators”, J. Math. Phys. 55,
021702 (2014).
[15] J. E. Mandula, Diagrammatic techniques in group theory, Univ. of Southampton,
Notes taken by S. N. Coulson and A. J. G. Hey, 1981.
Exercises
10.1. Conjugacy classes of SO(3): Show that all SO(3) rotations (10.7) by the same angle
θ around any rotation axis n are conjugate to each other:
Check this for infinitesimal φ, and argue that from that it follows that it is also true for
finite φ. Hint: use the Lie algebra commutators (10.9).
10.2. The character of SO(3) 3-dimensional representation: Show that for the 3-dimen-
sional special orthogonal representation (10.7), the character is
χ = 2 cos(θ) + 1 . (10.50)
Hint: evaluate the character explicitly for Rx (θ), Ry (θ) and Rz (θ), then explain what is
the intuitive meaning of ‘class’ for rotations.
10.3. The orthonormality of SO(3) characters: Verify that given the Haar measure (10.13),
the characters (10.12) are orthogonal:
Z
0
hχ(j)|χ(j )i = dg χ(j) (g −1 ) χ(j ) (g) = δjj 0 .
0
(10.51)
G
10.4. U(1) equivariance of two-modes system for finite angles: Show that the vector field
in two-modes system (10.2) is equivariant under (10.22), the unitary group U(1) acting
on R4 ∼
= C2 as the k = 1 and 2 modes:
10.5. Integrate the two-modes system: Integrate (10.5) and plot a long trajectory of two-
modes in the 4d state space, (x1 , y1 , y2 ) projection, as in figure 10.4 (a). To save you time
(typing in (10.5) is tedious), we have prepared for you python code, and online graded
problem set here. If you do this exercise, please get started early, in order to make sure
that the autograder is working, and forward to us the grades that you receive from the
autograder.
10.6. SO(2) or harmonic oscillator slice: Construct a moving frame slice for action of
SO(2) on R2
(x, y) 7→ (x cos θ − y sin θ, x sin θ + y cos θ)
by, for instance, the positive y axis: x = 0, y > 0. Write out explicitly the group
transformation that brings any point back to the slice. What invariant is preserved by this
construction?
10.7. Invariant subspace of the two-modes system: Show that (0, 0, x2 , y2 ) is a flow invari-
ant subspace of the two-modes system (10.5), i.e., show that a trajectory with the initial
point within this subspace remains within it forever.
10.8. Slicing the two-modes system: Choose the simplest slice template point that fixes the
1. Fourier mode,
x̂0 = (1, 0, 0, 0) . (10.53)
(a) Show for the two-modes system (10.5), that the velocity within the slice, and the
phase velocity along the group orbit are
(b) Determine the chart border (the locus of point where the group tangent is either not
transverse to the slice or vanishes).
(c) What is its dimension?
(d) What is its relation to the invariant subspace of exercise 10.7?
(e) Can a symmetry-reduced trajectory cross the chart border?
10.9. The symmetry reduced two-modes flow: Pick an initial point x̂(0) that satisfies
the slice condition for the template choice (10.53) and integrate (10.54) & (10.55). Plot
the three dimensional slice hyperplane spanned by (x1 , x2 , y2 ) to visualize the symmetry
reduced dynamics. Does it look like figure 10.4 (b)?
Homework HW11
Bonus points
Exercise 11.2 Real and pseudo-real representations of SO(3) 4 points
Exercise 11.3 Total spin of N particles 5 points
Total of 4 points = 100 % score. Extra points accumulate, can help you later if you
miss a few problems.
143
GROUP THEORY - WEEK 11. SU(2) AND SO(3)
lecture playlist
Gutkin notes, Lect. 9 SU(2), SO(3) and their representations, Sects. 1-3.2.
14.1 Recap: Irreps of SO(2) are not what you would have expected (24:37 min)
(optional) Anthony Zee [3] Group Theory in a Nutshell for Physicists: Cartan
classification of Lie algebras and Dynkin diagrams (5 lecture course)
Henriette Roux Why is it that the Fourier transformation works? The presence of a
discrete but infinite translational symmetry in a system calls for its use of it to
diagonalize the matrix and thus make calculations easier, but exactly why is the
Fourier transform able to do this?
Henriette Roux How is this Fourier transform as we have studied in the space/point
groups section related to that which we have derived from the projection opera-
tors?
Henriette Roux As an extension of the Fourier transform, are there any equivalent
of Fourier transforms for rotations or other infinite but discrete symmetries as
well? So for example, if there is a system with a discrete but infinite rotational
symmetry, is there a “rotational” transform where the representing matrix is di-
agonalized? Are there whole classes of such transformations?
Henriette Roux You say that position and momentum are “dual” to each other, and
so is the real space and reciprocal space (I guess it’s the same thing as position
and momentum but just for argument sake). The commonality between these are
the fact that they can be Fourier transformed from one space to another. Does
this mean that unitary operations, eiHt , suggest a Fourier transform from the
“energy” or “frequency” space to “time” space as well?
Henriette Roux This seems very closely related to Noether’s theorem as well, is there
a way to explain this similarity?
Henriette Roux The special thing about Lie groups is that there exist analytic func-
tions which link g(a) and c = f (b, a) for g(c) = g(a)g(b). Does this need
for analytic functions come from the fact that to construct a group manifold, the
maps relating different “local” Euclidean spaces need to be C ∞ , or smooth? If
so, is there a reference we can refer to which explains how the Lie groups sat-
isfy all the other conditions of a manifold (establishing an open ball, building an
atlas and so on) as well? Just as an extension, how do you even study groups
which do not fall under the realm of a manifold? Don’t common functions like
differentials and integrations not apply in spaces outside a manifold?
Henriette Roux Why is that we Taylor expand the group in the first place? How is
this connected to the shift to left/right group operators?
The next few questions are about General Relativity, and how is what is covered
in this course applicable to GR:
Henriette Roux We keep to the first order in the expansion for g(θ) as we are con-
sidering the tangent space to the manifold. In the context of the GR, the tangent
space was defined as the space of directional derivatives at a point. In our case,
we are studying groups, which are not, in general, vectors (well I guess they can
be [1×1] vectors/matrices but that’s only specific irreps, so how do we understand
the concept of tangent space as you have define it?
- Or does it work out since Lie groups are always Abelian and thus have an
infinite number of 1D irreps?
- What happens if we keep the expansion to the 2nd order? Does the mathematics
change in any way? Is there a good reason to ignore the 2nd and higher order
expansions, not just in the physics sense (keeping to largest order of significance)
but in the mathematical way of understanding things?
Question 11.1. Henriette Roux asks
Q Why is this complex 2-dimensional vector called a ‘spinor’?
A Historical, as Arfken, Weber & Harris [1] explain: “It turns out that half-integral angular mo-
mentum states are needed to describe the intrinsic angular momentum of the electron and many
other particles. Since these particles also have magnetic moments, an intuitive interpretation is
that their charge distributions are spinning about some axis; hence the term spin. It is now un-
derstood that the spin phenomena cannot be explained consistently by describing these particles
as ordinary charge distributions undergoing rotational motion, [...] ”
Schwichtenberg [2]: “[...] spinors have properties that usual vectors do not have. For in-
stance, the factor 1/2 in the exponent. This factor shows us that a spinor1 is after a rotation by
2π not the same, but gets a minus sign. This is a pretty crazy property, because all objects we
deal with in everyday life are exactly the same after a rotation by 360o = 2π.
gives the length of a vector. Consider a SU(2) transformation (11.1) of this matrix,
U † (σ · x)U . Taking the determinant, we find the same expression as before:
det U (σ · x)U † = det U det (σ · x) det U † = det (σ · x) . (11.4)
Just as SO(3), SU(2) preserves the lengths of vectors.
To make the correspondence between SO(3) and SU(2) more explicit, consider a
SU(2) transformation on a complex two-component spinor
α
ψ= (11.5)
β
related to x by
1 2 i
x= (β − α2 ), y = − (α2 + β 2 ), z = αβ (11.6)
2 2
Check that a SU(2) transformation of ψ is equivalent to a SO(3) transformation on x.
From this equivalence, one sees that a SU(2) transformation has three real parameters
that correspond to the three rotation angles of SO(3). If we label the ”angles” for the
SU(2) transformation by α, β, and γ, we observe, for a ”rotation” about x̂
cos α/2 i sin α/2
Ux (α) = , (11.7)
i sin α/2 cos α/2
for a ”rotation” about ŷ,
cos β/2 sin β/2
Uy (β) = , (11.8)
− sin β/2 cos β/2
and for ”rotation” about ẑ,
eiγ/2
0
Uz (γ) = . (11.9)
0 e−iγ/2
Compare these three matrices to the corresponding SO(3) rotation matrices:
1 0 0 cos φ 0 sin φ
Rx (ζ) = 0 cos ζ sin ζ , Ry (φ) = 0 1 0
0 − sin ζ cos ζ − sin φ 0 cos φ
cos θ sin θ 0
Rz (θ) = − sin θ cos θ 0 (11.10)
0 0 1
They’re equivalent! Result: Half the rotation angle generated by SU(2) corresponds
to a rotation generated by SO(3).
What does this mean? At this point, probably best to switch to Schwichtenberg [2]
(click here) who explains clearly that SU(2) is a simply-connected group, and thus the
”mother" or covering group, or the double cover of SO(3). This means there is a two-
to-one map from SU(2) to SO(3); an SU(2) turn by 4π corresponds to an SO(3) turn
by 2π. So, the building blocks of your 3-dimensional world are not 3-dimensional real
vectors, but the 2-dimensional complex spinors! Quantum mechanics chose electrons
to be spin 1/2, and there is nothing Fox Channel can do about it.
References
[1] G. B. Arfken, H. J. Weber, and F. E. Harris, Mathematical Methods for Physi-
cists: A Comprehensive Guide, 7th ed. (Academic, New York, 2013).
[2] J. Schwichtenberg, Physics from Symmetry (Springer, Berlin, 2015).
[3] A. Zee, Group Theory in a Nutshell for Physicists (Princeton Univ. Press, Prince-
ton NJ, 2016).
Exercises
11.1. The characters of SO(3) representations: Show that for an irrep labeled by j, the
character of a conjugacy class labeled by θ
sin(j + 1/2)θ
χ(j) (θ) = (11.11)
sin(θ/2)
can be obtained by taking the trace of Rzj (θ). Verify that for j = 1 this character is the
three dimensional special orthogonal representation character (10.11).
11.2. Real and pseudo-real representations of SO(3). Recall (Gutkin notes, sect. 4.5
Representation Theory II, Sect. 5 5. Three types of representations) that there are exist
three types of representation which can be distinguished by the indicator (4.6):
Z +1 real
dµ(g)χ(`) (g 2 ) = 0 complex . (11.12)
G
−1 pseudo-real
Hint: The characters and Haar measure (10.13) of SO(3) are given by
sin ` + 12 θ
dθ
χ(`) (g) = , dµ(g) = sin2 (θ/2) (11.13)
sin 21 θ
π
(a) What are possible values s for the total spin of the system?
(b) Determine dimension of the subspace of V0 with 0 total spin. In other words: how
many times trivial representation enters into product:
(B. Gutkin)
11.5. Quadrupole transitions.
a) Write Q1 = xy, Q2 = zy, Q3 = x2 − y 2 and Q4 = 2z 2 − x2 − y 2 as components of
spherical tensor of rank 2. Hint: use spherical harmonics Ylm (θ, ϕ).
b) The last quantity Q4 is known as quadrupole moment. What are the selection rules for
transitions induced by Q4 in a system with SO(3) symmetry? In other words, for which
m, l and k, j the transition rates:
are non-zero?
c) By using Wigner-Eckart theorem write down the relationship between |hm l|Q4 |j ki|2
and |hm l|Q1 |j ki|2 in terms of Clebsch-Gordan coefficients.
(B. Gutkin)
Homework HW12
Bonus points
Exercise 12.5 The unbearable lightness of SO(4) Lie algebra 15 points
151
GROUP THEORY - WEEK 12. LORENTZ GROUP; SPIN
Lecture 23
2021-07-08 Predrag
SO(4) = SU(2) ⊗ SU(2); Lorentz group
Lecture 15 (Unedited) SU(2) irreps. SO(4) = SU(2) × SU(2). More importantly:
Minkowski metric, Lorentz group SO(1, 3) irreps are also labeled by pairs of
SU(2) × SU(2) irrep labels. (2:29:20 h)
For Lorentz group, read Schwichtenberg [3] Sect. 3.7 (click here).
For SO(n) see also birdtracks.eu Chapt. 10 Orthogonal groups, pp. 121-123.
For SO(4) = SU(2) ⊗ SU(2) see also birdtracks.eu sect. 20.3.1 SO(4) or Cartan
A1 + A1 algebra.
Peter Voit [5] Quantum Theory, Groups and Representations (2017) has a nice
calculation of spherical harmonics as SO(3) eigenvectors in polar coordinates.
Used as problem set #6 for Quantum Mechanics I Georgia Tech PHYS-6105,
November 17, 2022, should make it a section in these notes. S
John Wood’s (click here) notes and exercise 12.5 The unbearable lightness of
SO(4) Lie algebra. The challenge: achieve some elegance in deriving the SO(4)
commutator relations.
Ilya Kuprov [2] ‘What exactly is spin?’ (40 min). Starts out with very enter-
taining bits of physics history. Wlad Sobol writes: “Kuprov derives the Dirac
equation from Wigner symmetry theorem, from a product of two Casimir opera-
tors of the Poincaré group." Revisit sect. 10.8.1 Quaternionic speculations for a
different point of view.
xµ yµ = xµ ηµν y ν = x0 y 0 − x1 y 1 − x2 y 2 − x3 y 3
x 7→ x0 = ΛxΛ† (12.11)
and Λ explicitly acts as a Lorentz transformation (12.1), with Λµν = 12 tr (σ̄ µ Λσν Λ† ) .
The mapping is two-to-one, as two matrices ±Λ ∈ SL(2, C) generate the same Lorentz
transformation ΛxΛ† = (−Λ)x(−Λ)† . This Λ belong to the proper orthochronous
Lorentz group SO+ (1, 3), and it can be shown that SL(2, C) is simply connected and
is the double universal cover of the SO+ (1, 3).
Consider the fully antisymmetric Levi-Civita tensor ε = −ε−1 = −εT in two
dimensions
0 1
ε = iσ2 = . (12.13)
−1 0
This defines a symplectic (i.e., skew-symmetric) bilinear form hu, vi = −hv, ui on two
spinors u and v, elements of the two-dimensional complex vector (or spinor) space C2
1 1
u v
u= , v= , (12.14)
u2 v2
equipped with the symplectic form
hu, vi = u1 v 2 − u2 v 1 = uT εv . (12.15)
This symplectic form is SL(2, C)-invariant
hu, vi = uT εv = hΛu, Λvi = uT ΛT εΛv , (12.16)
so one can interpret the group acting on spinors as SL(2, C) ∼
= Sp(2, C) , the complex
symplectic group in two dimensions
Sp(2, C) = {Λ ∈ GL(2, C) | ΛT εΛ = ε} . (12.17)
Summary. The group of Lorentz transformations of spinors is the group SL(2,C)
of [2×2] complex matrices with determinant 1, i.e., the invariant tensor is the 2-index
Levi-Civita εAB . The SL(2,C) matrices are parametrized by three complex dimensions
and therefore six real ones (the matrices have four complex numbers and one complex
constraint on the determinant). This matches the 6 dimensions of the group manifold
associated with the Lorentz group SO(1, 3).
Andrew M. Steane writes “A spinor is the most basic mathematical object that can
be Lorentz-transformed.” His An introduction to spinors, arXiv:1312.3824, might help
you develop intuition about spinors.
Andrzej Trautman tracks the origin of spinors to Euclid, and General Relativity to
Clifford. He includes a letter from Hades saying, inter alia, “Unfortunately, it appears
that there is now in your world a race of vampires, called referees, who clamp down
mercilessly upon mathematicians unless they know the right passwords.”
n(n − 1)
N = tr PA = = . (12.18)
2
Young tableaux × = • + +
n(n−1) (n+2)(n−1)
Dimensions n2 = 1 + 2
+ 2
111
000
1 000
111
000
111 1
Projectors = n
+ 000
111 + −
000
111
000
111
n
Table 12.1: SO(n) Clebsch-Gordan series for V ⊗V , worked out in detail in Group The-
ory – Birdtracks, Lie’s, and Exceptional Groups, birdtracks.eu Chapt. 10 Orthogonal
groups.
The SO(2, 1) method can be extended to solve relativistic central force problems
(one of my students did his Ph.D. thesis on this 20 years ago).
Q: Is the geometry associated with these algebraic structures, as applied to central
force problems, explored?
References
[1] M. Bander and C. Itzykson, “Group theory and the hydrogen atom (I)”, Rev.
Mod. Phys. 38, 330–345 (1966).
[2] I. Kuprov, Spin: From Basic Symmetries to Quantum Optimal Control (Springer,
New York, 2022).
[3] J. Schwichtenberg, Physics from Symmetry (Springer, Berlin, 2015).
[4] E. Wigner, “On unitary representations of the inhomogeneous Lorentz group”,
Ann. Math. 40, 149–204 (1939).
[5] P. Woit, Quantum Theory, Groups and Representations (Springer, New York,
2017).
Exercises
12.1. Lie algebra of SO(4) and SU(2) ⊗ SU(2). One particle Hamiltonian with a central
potential has in general SO(3) symmetry group. It turns out, however, that for Coulomb
potential the symmetry group is actually larger - SO(4), rather than SO(3). This explains
why the energy level degeneracies in the hydrogen atom are anomalously large. So SO(4)
and its representations are of a special importance in atomic physics.
(a) Show that the Lie algebra so(4) of the group SO(4) is generated by real antisym-
metric 4 × 4 matrices.
(b) What is the dimension of so(4)?
[Mab , Mcd ] = Mad δbc + Mbc δad − Mac δbd − Mbd δac . (12.20)
(d) Show that Lie algebras of the groups SO(4) and SU(2) × SU(2) are isomorphic.
Path:
and
1 1
Ak = (Jk + Kk ) and Bk = (Jk − Kk ) .
2 2
(d.ii) Show that A and B satisfy the same commutation relations as two copies of
su(2).
(e) How does one construct irreps of so(4) out of irreps of su(2)?
(f) Are groups SO(4) and SU(2) ⊗ SU(2) isomorphic to each other?
(B. Gutkin)
12.2. SO(n) Clebsch-Gordan series for V ⊗V .
(a) Show that the product of two n-dimensional reps of SO(n) decomposes into three
irreps:
000
111
1 000
111 1
= + 000
111 + − . (12.21)
000
111
n 000
111 n
000
111
x 7→ x0 = ΛxΛ† . (12.25)
(a) Check that x0 is Hermitian, and the Minkowski scalar product (12.23) is preserved.
(b) Show that Λ explicitly acts as a Lorentz transformation x0µ = Λµν xν .
(c) Show that the mapping from a Λ ∈ SL(2, C) to the Lorentz transformation in
SO(1, 3) is two-to-one.
hu, vi = uT εv = u1 v 2 − u2 v 1 . (12.27)
12.5. The unbearable lightness of SO(4) Lie algebra. Download John Wood’s (click here)
notes. The challenge: achieve some elegance in deriving the SO(4) commutator bracket
relations, for example reduce the number of steps in the calculation by 30% or 50%.
The prize: a case of beer, details to be negotiated with John.
The challenges start on p. 9-8, following eq. (9.21), i.e., “(i)”, “(iv)”, and “(v)”. For
instance, on p. 9-11 John indicates all of the cancellations. These suggest that his solution
is “calculating zero” unnecessarily. One could take linear combinations of the operators
that possess these commutator bracket relations; but the combinations do not seem a priori
warranted on the basis of the dynamics of the problem.
(J. Wood)
Homework HW13
Bonus points
Exercise 13.3 SU(3) symmetry in 3D harmonic oscillator 5 points
Total of 10 points = 100 % score. Extra points accumulate, can help you later if you
miss a few problems.
161
GROUP THEORY - WEEK 13. SIMPLE LIE ALGEBRAS; SU(3)
In week 4 we learned that for finite groups there is one very special matrix rep, the
regular representation constructed from the group multiplication table, that is intrinsic
to the abstract group itself, and whose reduction yields all irreps of a given group.
In week 9 we saw that for continuous groups we need to study the Lie algebra
of the finite number of generators Tj , rather than the infinity of group elements g =
exp(iφ · T ). Here the finite group multiplication table is replaced by the ‘Lie product’,
i.e., the table of Lie commutators’ fully antisymmetric structure constants iCijk .
So far we have chosen the hermitian basis Tj . But non-hermitian bases are also
OK, as we know from raising / lowering operators of SU(2) of quantum mechanics
irrep constructions.
(a) (b)
Figure 13.1: (a) The meson (pseudoscalars) octet. (b) The quark triplet, the anti-quark
triplet and the gluon octet. (Wikipedia).
Exercises
13.1. Root system of simple Lie algebras.
a) Determine dimensions of Lie algebras so(N ), su(N ) and dimensions of their Cartan
subalgebras. What is the number of the positive roots for these Lie algebras?
b) Show that N × N diagonal matrices Hi with zero traces and uper/lower corner N × N
(a,b)
matrices E (a,b) with the elements Ei,j = δia δib provide Cartan-Weyl basis of su(N ).
To put it differently, show that E (a,b) are eigenstates for adjoint representation of Hi ’s.
(B. Gutkin)
13.2. Meson octet. In Gutkin lecture notes, Lect. 11 Strong interactions: flavor SU(3), the
meson octet, figure 13.1 (a)
π0 η
π+ K+
√ + √
2 6
Φ = π− π0
−√ + √η6 K0
2
K− K0 −√2η
6
0
π +
π 0 K+
√ 0 0 1 0 0
−2 0 0 η
= π π + 0 0 K + √ 0 1 0 (13.1)
,
−√ 0
2
K− K0 0 6 0 0 −2
0 0 0
is interpreted as arising from the adjoint representation of SU(3), i.e., the traceless part of
the quark-antiquark 3 ⊗ 3 = 1 ⊕ 8 outer product (see figure 13.1 (b)),
uu ud us
du dd ds . (13.2)
su sd ss
where we have replaced in (13.1) the constituent q ⊗ q combinations by the names of the
elementary particles they build.
Given the quark quantum numbers
Q I I3 Y B
u 2/3 1/2 1/2 1/3 1/3
d -1/3 1/2 -1/2 1/3 1/3
s -1/3 0 0 -2/3 1/3
Flavor SU(3)
Homework HW14
Bonus points
Exercise 14.1 Gell-Mann–Okubo mass formula 8 points
Exercise 15.3 Young tableaux for SU(3) 3 points
Exercise 15.4 Irrep projection operators for unitary groups 5 points
All points are bonus points. Extra points accumulate, can help you if you had missed a
few problems.
165
GROUP THEORY - WEEK 14. FLAVOR SU(3)
Figure 14.1: A lattice gauge theory calculation of the light QCD spectrum. Horizontal
lines and bands are the experimental values with their decay widths. The π, K and and
Ξ have no error bars because they are used to set the light and strange quark masses
and the overall scale respectively. From Scholarpedia.
The Gell-Mann-Okubo mass sum rules [1, 2, 4] are an easy consequence of the
approximate SU(3) flavor symmetry. Determination of quark masses is much harder -
they are parameters of the standard model, determined by optimizing the spectrum of
particle masses obtained by lattice QCD calculations as compared to the experimental
baryon and meson masses. The best determination of the mass spectrum as of 2012
is given in figure 14.1. Up, down quarks are about 3 and 6 MeV, respectively, with
strange quark mass about 100 MeV, all with large error brackets. As of 2021, I have
not found an update to figure 14.1, but the latest on the subject can probably be traced
in Georg von Hippel’s latticeqcd.blogspot.com.
References
[1] M. Gell-Mann, The Eightfold Way: A Theory of Strong Interaction Symmetry
(CalTech, 1961).
[2] M. Gell-Mann, “Symmetries of baryons and mesons”, Phys. Rev. 125, 1067–
1084 (1962).
[3] R. Murgan and A. Zender, “Energy eigenvalues of the three-dimensional quan-
tum harmonic oscillator from SU(3) cubic Casimir operator”, Eur. J. Phys. 40,
015405 (2018).
[4] S. Okubo, “Note on unitary symmetry in strong interactions”, Progr. Theor.
Phys. 27, 949–966 (1962).
Exercises
14.1. Gell-Mann–Okubo mass formula. The mass symmetry-breaking interaction for
an isospin multiplet is proportional to the 3rd component of the isospin operator, I3 .
Similarly, the symmetry-breaking interaction of SU(3) for the meson octet is given by the
8th component of the octet operator Y = λ8 . Derive the GMO mass formula for mesons
4 m2K − m2π
m2η = . (14.1)
3
by eliminating the parameter for the strength of this interaction, as in Gutkin lecture notes,
Lect. 11 Strong interactions: flavor SU(3) .
Homework HW15
Bonus points
Exercise 15.1 Representations of SU(3) 5 points
Exercise 15.2 Young tableaux for S5 3 points
Exercise 15.3 Young tableaux for SU(3) 3 points
Exercise 15.4 Irrep projection operators for unitary groups 5 points
169
GROUP THEORY - WEEK 15. MANY PARTICLE SYSTEMS. YOUNG
TABLEAUX
Young tableaux for SU(3) and SU(n) have not yet been covered in the lectures, but
you can easily learn them yourself, from, for example, Gutkin notes, Lect. 12 Young
tableaux. Boris Gutkin is a professor, beyond learning new stuff, so he follows old
fashioned references such as Fulton and Harris [6]. The resulting simple recipe with 0
explanation can be found, for example, here: Young diagrams by C.G. Wohl.
A modern exposition is given in Group Theory – Birdtracks, Lie’s, and Exceptional
Groups, birdtracks.eu Chapt. 9 Unitary groups. Currently I am a fan of the Alcock-
Zeilinger algorithm [1–3], based on the simplification rules of ref. [2], which leads to
explicitly Hermitian and compact expressions for the projection operators.
Probably best to read Alcock-Zeilinger course The Special Unitary Group, Bird-
tracks, and Applications in QCD notes [4]. Alcock-Zeilinger fully supersedes Cvi-
tanović’s formulation, and any future full exposition of birdtracks reduction of SU(N )
tensor products into irreducible representations should be based on the Alcock-Zeilin-
ger algorithm.
Probably best to read Alcock-Zeilinger course The Special Unitary Group, Bird-
tracks, and Applications in QCD notes [4]. You want to study these in detail if your
research leads you to study of multiparticle states.
References
[1] J. Alcock-Zeilinger and H. Weigert, “Transition operators”, J. Math. Phys. 58,
051702 (2016).
[2] J. Alcock-Zeilinger and H. Weigert, “Compact Hermitian Young projection op-
erators”, J. Math. Phys. 58, 051702 (2017).
[3] J. Alcock-Zeilinger and H. Weigert, “Simplification rules for birdtrack opera-
tors”, J. Math. Phys. 58, 051701 (2017).
[4] J. M. Alcock-Zeilinger, The Special Unitary Group, Birdtracks, and Applica-
tions in QCD, tech. rep. (Univ. Tübingen, 2018).
[5] P. Cvitanović, Group Theory: Birdtracks, Lie’s and Exceptional Groups (Prince-
ton Univ. Press, Princeton NJ, 2008).
[6] W. Fulton and J. Harris, Representation Theory (Springer, New York, 1991).
Exercises
15.1. Representations of SU(3). Any irrep of SU(3) can be labeled D(p, q) by its highest
weight λ = pλ1 + qλ2 , where λ1,2 are the two fundamental weights.
(a) Find all irreps D(p, q) of SU(3) with the dimensions less then 20 (see lecture notes
for the dimensions of D(p, q)).
(b) Draw the lattice Λ generated by λ1,2 and mark there all the weights v (i.e., lattice
nodes) which belong to irrep. D(3, 0). Is D(3, 0) a real irrep?
(c) Consider product (reducible) representation 3 ⊗ 3, where 3 = D(1, 0) is the fun-
damental irrep. Mark all the weights v on Λ which belong to 3 ⊗ 3. Using this find
out decomposition of 3 ⊗ 3 into irreps:
3 ⊗ 3 = ⊕ 4, =?, 4 =?
15.3. Young tableaux for SU(3). Solve exercise 15.1 (c,d) by using Young tableaux.
Remark: If Young tableaux for SU(3) are not covered in the lectures, learn them your-
self from, for example, birdtracks.eu Group Theory Birdtracks, Lie’s, and Exceptional
Groups. The resulting simple recipe with 0 explanation can be found, for example, here
C.G. Wohl.
(B. Gutkin)
15.4. Irrep projection operators for unitary groups. Derive projection operators and
dimensions for irreps of the Kronecker product of the defining and the adjoint reps of
SU(n) listed in Group Theory Birdtracks, Lie’s, and Exceptional Groups, birdtracks.eu
table 9.3. (Ignore “indices," we have not defined them here.)
Homework HW16
Bonus points
Exercise 16.1 Gravity tensors, part (a) 2 points
Exercise 16.1 Gravity tensors, part (b) 4 points
Exercise 16.1 Gravity tensors, part (c) 1 point
Exercise 16.1 Gravity tensors, part (d) 2 points
Exercise 16.1 Gravity tensors, part (e) 3 points
Exercise 16.1 Gravity tensors, part (f) 4 points
Exercise 16.1 Gravity tensors, part (g) 3 points
Exercise 16.1 Gravity tensors, part (h) 6 points
Exercise 16.1 Gravity tensors, part (i) 4 points
Exercise 16.1 Gravity tensors, part (j) 10 points
173
GROUP THEORY - WEEK 16. WIGNER 3- AND 6-J COEFFICIENTS
The webbook for cyclists (medium level): Tracks, Lie’s, and Exceptional Magic.
Most of the webbook at a cyclist pace, in 50 overheads.
The final result is invariant and highly elegant: any group-theoretical invariant
quantity can be expressed in terms of Wigner 3- and 6-j coefficients:
Sect. 5.1 Couplings and recouplings,
Sect. 5.2 Wigner 3n-j coefficients, and
Sect. 5.3 Wigner-Eckart theorem.
The rest is just bedside reading, nothing technical:
Sect. 4.8 Irrelevancy of clebsches and
Sect. 4.9 A brief history of birdtracks.
Woit writes here about the “The Stormy Onset of Group Theory in the New Quan-
tum Mechanics,” citing Bonolis [2] From the rise of the group concept to the
stormy onset of group theory in the New Quantum Mechanics. A saga of the
invariant characterization of physical objects, events and theories.
Chayut [3] From the periphery: the genesis of Eugene P. Wigner’s application of
group theory to quantum mechanics traces the origins of Wigner’s application
of group theory to quantum physics to his early work as a chemical engineer,
in chemistry and crystallography. “In the early 1920s, crystallography was the
only discipline in which symmetry groups were routinely used. Wigner’s early
training in chemistry exposed him to conceptual tools which were absent from
the pedagogy available to physicists for many years to come. This both enabled
and pushed him to apply the group theoretic approach to quantum physics. It
took many years for the approach first introduced by Wigner in the 1920s – and
whose reception by the physicists was initially problematical – to assume the
pivotal place it now holds.” Another historical exposition is given by Scholz [6]
Introducing groups into quantum theory (1926–1930).
So what is group theory good for? By identifying the symmetries, one can apply
group theory to determine good quantum numbers which describe a physical state (i.e.,
the irreps). Group theory then says that many matrix elements vanish, or shows how
are they related to others. While group theory does not determine the actual value of a
matrix element of interest, it vastly simplifies its calculation.
The old fashioned atomic physics, fixated on SO(3) / SU(2), is too explicit, with too
many bras and kets, too many square roots, too many deliriously complicated Clebsch-
Gordan coefficients that you do not need, and way too many labels, way too explicit for
you to notice that all of these are eventually summed over, resulting in a final answer
much simpler than any of the intermediate steps.
I wrote my book [4] Group Theory - Birdtracks, Lie’s, and Exceptional Groups to
teach you how to compute everything you need to compute, without ever writing down
a single explicit matrix element, or a single Clebsch-Gordan coefficient. There are two
versions. There is a particle-physics / Feynman diagrams version that is index free,
graphical and easy to use (at least for the low-dimensional irreps). The key insights
are already in Wigner’s book [8]: the content of symmetry is a set of invariant numbers
that he calls 3n-j’s. Then there are various mathematical flavors (Weyl group on Cartan
lattice, etc.), elegant, but perhaps too elegant to be computationally practical.
References
[1] S. L. Adler, J. Lieberman, and Y. J. Ng, “Regularization of the stress-energy ten-
sor for vector and scalar particles propagating in a general background metric”,
Ann. Phys. 106, 279–321 (1977).
[2] L. Bonolis, “From the Rise of the Group Concept to the Stormy Onset of Group
Theory in the New Quantum Mechanics. A saga of the invariant characterization
of physical objects, events and theories”, Rivista Nuovo Cim. 27, 1–110 (2005).
[3] M. Chayut, “From the periphery: the genesis of Eugene P. Wigner’s application
of group theory to quantum mechanics”, Found. Chem. 3, 55–78 (2001).
[4] P. Cvitanović, Group Theory: Birdtracks, Lie’s and Exceptional Groups (Prince-
ton Univ. Press, Princeton NJ, 2008).
[5] R. Penrose, “Applications of negative dimensional tensors”, in Combinatorial
mathematics and its applications, edited by D. J. A. Welsh (Academic, New
York, 1971), pp. 221–244.
[6] E. Scholz, “Introducing groups into quantum theory (1926-1930)”, Hist. Math.
33, 440–490 (2006).
[7] S. Weinberg, Gravitation and Cosmology: Principles and Applications of the
General Theory of Relativity (Wiley, New York, 1972).
[8] E. P. Wigner, Group Theory and Its Application to the Quantum Mechanics of
Atomic Spectra (Academic, New York, 1931).
Exercises
16.1. Gravity tensors. In this problem we will apply diagrammatic methods (“birdtracks”)
to construct and count the numbers of independent components of the “irreducible rank-
four gravity curvature tensors.” However, any notation that works for you is OK, as long
as you obtain the same irreps and their dimensions. The goal of this exercise (longish, as
much of it is the recapitulation of the material covered in the book) is to give you basic
understanding for how Young tableaux work for groups other than U(n). We start with
(a) The Riemann-Christoffel curvature tensor of general relativity has the following
symmetries (see, for example, Weinberg [7] or the Riemann curvature tensor wiki):
Rαβγδ = −Rβαγδ (16.1)
Rαβγδ = Rγδαβ (16.2)
Rαβγδ + Rβγαδ + Rγαβδ = 0. (16.3)
Introducing a birdtrack notation for the Riemann tensor
α
β
Rαβγδ =
γ R , (16.4)
δ
check that we can state the above symmetries as
Rαβγδ = −Rβαγδ
R = R , (16.5)
Rαβγδ = Rγδαβ
R = R , (16.6)
R + R + R =0. (16.7)
(c) Show that the third condition (16.7) says that R has no components in
the irrep:
R + R + R =3 R = 0 . (16.9)
Hence, the symmetries of the Riemann tensor are summarized by the irrep
projection operator [5]:
α α´
0 0 0 0 4β β´
(PR )αβγδ ,δ γ β α = γ´ (16.10)
3γ
δ δ´
(e) Compute the number of independent components of the Riemann tensor Rαβγδ by
taking the trace of the irrep projection operator:
n2 (n2 − 1)
dR = tr PR = . (16.12)
12
The U(n) 2-index tensors can be decomposed into a sum of their symmetric and antisym-
metric parts. Specializing to the subgroup SO(n), the rule is to lower all indices on all
tensors, and the symmetrization projection operator is written as
0 0
Sµν,ρσ = gρρ0 gσσ0 Sµν ,ρ σ
1
= (gµσ gνρ + gµρ gνσ )
2
From now on, we drop all arrows and g µν ’s and write the decomposition into symmetric
and antisymmetric parts as
= +
1 1
gµσ gνρ = (gµσ gνρ + gµρ gνσ ) + (gµσ gνρ − gµρ gνσ ) . (16.17)
2 2
T2 = = nT (16.19)
yields the trace and the traceless part projection operators. As T is symmetric, ST = T,
only the symmetric subspace is reduced by this invariant.
(f) Show that SO(n) 2-index tensors decompose into three irreps:
traceless symmetric:
1 1 1
(P2 )µν,ρσ = (gµσ gνρ + gµρ gνσ ) − gµν gρσ = − ,
2 n n
(16.20)
1 1
singlet: (P1 )µν,ρσ = gµν gρσ = , (16.21)
n n
1
antisymmetric: (P3 )µν,ρσ = (gµσ gνρ − gµρ gνσ ) = .
(16.22)
2
What are the dimensions of the three irreps?
(g) In the same spirit, the U(n) irrep is decomposed by the SO(n) intermediate
2-index state invariant matrix
Q= . (16.23)
Show that the intermediate 2-index subspace splits into three irreducible reps by
(16.20) – (16.22):
( )
1 1
Q = + − +
n n
= Q0 + QS + QA . (16.24)
P R QA = 0 . (16.25)
(Hint: The Riemann tensor is symmetric under the interchange of index pairs.)
(h) Fix the normalization of the remaining two projection operators by computing
Q2S , Q20 :
2
P0 = , (16.26)
n(n − 1)
( )
4 1
PS = − (16.27)
n−2 n
U(n) → SO(n)
→ + + ◦
PR = PW + PS + P0
n2 (n2 −1) (n+2)(n+1)n(n−3) (n+2)(n−1)
12
= 12
+ 2
+ 1
(16.28)
The projection operator for the SO(n) traceless irrep is:
PW = PR − PS − P0
4 4 2
PW = − + (16.29).
3 n−2 (n − 1)(n − 2)
(i) The above three projection operators project out the standard, SO(n)-irreducible
general relativity tensors:
Curvature scalar:
R=− R = Rµνµν (16.30)
1 1
Rµν − gµν R = − R +n R (16.31)
n
Weyl tensor:
The numbers of independent components of these tensors are given by the dimen-
sions of corresponding irreducible subspaces in (16.28).
What is the lowest dimension in which the Ricci tensor contributes? the Weyl tensor
contributes? Show that in 2, respectively 3 dimensions, we have
(j) The last example of this exercise is an application of birdtracks to general relativ-
ity index manipulations. The object is to find the characteristic equation for the
Riemann tensor in four dimensions.
The antisymmetrization tensor Aa1 a2 ... ,bp ...b2 b1 has nonvanishing components, only
if all lower (or upper) indices differ from each other. If the defining dimension is
smaller than the number of indices, the tensor A has no nonvanishing components:
1
2
=0 if p > n . (16.34)
...
...
p
This identity implies that for p > n, not all combinations of p Kronecker deltas are
linearly independent. A typical relation is the p = n + 1 case
0= R , (16.36)
R R R R
0 = 2 −4
−4
R R + 2R
R (16.37)
R +1 R R
( )
R2 R
− −2 .
2 2
This identity has been used by Adler et al., eq. (E2) in ref. [1].
183
GROUP THEORY - WEEK 17. AN OVERVIEW, AND THE EPILOGUE
week 4 Characters
Schur’s Lemma. Unitary matrices can be diagonalized, and from that follows
the Wonderful Orthogonality Theorem for Characters (coordinate independent,
intrinsic numbers), and the full reducibility of any representation of any finite
group.
week 5 Classes
The algebra of central or ‘all-commuting’ class operators, connects the reduction
in terms of characters to the projection operators of week 1. The key idea:
Define a group by what objects (primitive invariant tensors) it leaves invariant.
week 6 Fundamental domain
Dynamical systems application: the Lorenz flow, its C2 symmetry and its desym-
metrization: if the system is nonlinear, its symmetry reduction is not easy.
week 7 Discrete Fourier representation
So far, everything was finite and compact. Next: two distinct ways of going
infinite: (a) discrete translations, exemplified by deterministic diffusion and spa-
ce groups of week 8, and (b) continuous Lie groups, exemplified by rotations of
week 9.
week 8 Space groups
Translation group, Bravais lattice, wallpaper groups, reciprocal lattice, Brilluoin
zone.
week 9 Continuous groups
Lie groups. Matrix representations. Invariant tensors. Lie algebra. Adjoint
representation, Jacobi relation. Birdtracks.
Irreps of SO(2) and O(2) Clebsch-Gordan series (i.e., reduction of their prod-
ucts).
week 10 Lie groups, algebras; O(2) symmetry sliced
(a) Group integrals. SO(3) character orthogonality.
(b) Continuous symmetry reduction for a nonlinear system is much harder than
discrete symmetry reduction of week 7. “Slicing” is a research level topic.
week 11 SU(2) and SO(3)
SU(2) ' SO(3) correspondence leads to the next rude awakening; our 3-dim-
ensional Euclidean space is not fundamental! All irreps of SO(3) are built from
2-dimensional complex vectors, or 1/2 spins.
week 12 Lorentz group; spin
(a) We now loose compactness: even though the SO(1, 3) Lorentz invariance
group of the Minkowski space symmetries is not compact, its Lie algebra still
closes, as for the compact SO(4).
(b) SO(4) ' SU(2) ⊗ SU(2) correspondence leads to the Minkowski 4-dimen-
sional space not being fundamental either - all irreps of the Lorentz group are
built from combinations of 2-dimensional complex vectors, or spinors.
(c) together with general relativity, this leads to replacement of the Minkowski
continuum by a 4-dimensional spacetime (or quantum) foam, a candidate theory
of quantum gravity.
week 13 Simple Lie algebras; SU(3)
The next profound shift:
So far all our group notions were based on tangible, spatial intuition: permuta-
tions, reflections, rotations. But now Lie groups take on a life of their own.
(a) The SO(3) theory of angular momenta generalizes to Killing-Cartan lattices,
and a fully abstract enumeration of all possible semi-simple compact Lie groups.
(b) SU(2) is promoted to an internal isospin symmetry, decoupled from our Eu-
clidean spatial intuition. Modern particle physics is born, with larger and larger
internal symmetry groups, tacked onto higher and higher dimensional continuum
spacetimes.
week 14 Flavor SU(3)
Gell-Mann–Okubo formula. The next triumph of particle physics is yet another
departure; observed baryons and mesons are built up from quarks, particles by
assumption unobservable in isolation.
week 15 Young tableaux
We have come full circle now: as a much simpler alternative to the Cartan-
Killing construction, irreps of the finite symmetric group Sn classify the irreps
of the continuous SU(n) symmetry multi-particle states.
week 16 Wigner 3- and 6-j coefficients
The goal of group theory is to predict measurable numbers, numbers independent
of any particular choice of coordinate. The full reducibility says that any such
number is built from 3- and 6-j coefficients: they are the total content of group
theory.
If you are reading this in preparation for a final exam, think of it as an opportunity to
rethink the key ideas of this branch of mathematics, take with you the few essential
insights that may serve you well in your career later on.
improvements For this specific course, there are so many resources. It covers the best
parts of many books most probably, but from the student’s perspective every time
they open the middle of a new book they can find themselves looking at the very
first page at that book by following previous equations. Also the mathematical
aspect of the course can screen the physics side of the picture. As an experimen-
talist I got lost in the second half of the course. I found it easier to follow when
Dresselhaus or Tinkham book were being followed. Tinkham’s book was hard
for me at some point as well but these two books were easier to read for me than
other resources.
The course could be dramatically improved with a stronger emphasis on funda-
mentals of group structure and constructive examples throughout, accompanied
by a much more deliberate focus on a smaller set of topics. As an example of
a fundamental topic which I think has been inadequately covered: I am still un-
sure what distinguishes irreducible representations from any representation of a
group on an intuitive level. It is, as far as I can tell, something I should simply
know by now, but has not been adequately explained to me.
comments This course has been bewildering for the latter half of the semester. On
several occasions I have co-opted class time to have concepts from homework
and previous lectures explained carefully, which has been consistently unhelp-
ful. Overall, this course is by far the most stressful, confusing, and frustrating
course I have taken at Georgia Tech and on more than one occasion has given
me nightmares.
Instructor improvements The most needed improvement in the lectures is a stronger
focus on examples (especially those reinforcing fundamental concepts) which
reflect the Cvitanović expectations for the rest of the course.
to be CONTINUED