Problems in Matrix Calculs
Problems in Matrix Calculs
in
Matrix Calculus
by
Willi-Hans Steeb
International School for Scientific Computing
at
University of Johannesburg, South Africa
Yorick Hardy
Department of Mathematical Sciences
at
University of South Africa
Preface
by Willi-Hans Steeb
World Scientific Publishing, Singapore 2006
ISBN 981 256 916 2
http://www.worldscibooks.com/mathematics/6202.html
by Willi-Hans Steeb
World Scientific Publishing, Singapore 2007
ISBN 981-256-916-2
http://www.worldscibooks.com/physics/6515.html
v
The International School for Scientific Computing (ISSC) provides certifi-
cate courses for this subject. Please contact the author if you want to do
this course or other courses of the ISSC.
[email protected]
[email protected]
Home page of the author:
http://issc.uj.ac.za
vi
vii
Contents
Preface v
Notation xi
1 Basic Operations 1
2 Linear Equations 16
6 Decomposition of Matrices 86
7 Functions of Matrices 90
viii
16 Inequalities 160
21 Miscellaneous 176
Bibliography 201
Index 206
ix
xi
Notation
:= is defined as
∈ belongs to (a set)
∈
/ does not belong to (a set)
∩ intersection of sets
∪ union of sets
∅ empty set
N set of natural numbers
Z set of integers
Q set of rational numbers
R set of real numbers
R+ set of nonnegative real numbers
C set of complex numbers
Rn n-dimensional Euclidean space
space of column vectors with n real components
Cn n-dimensional complex linear space
space of column vectors with n complex components
H Hilbert
√ space
i −1
<z real part of the complex number z
=z imaginary part of the complex number z
|z| modulus of complex number z
|x + iy| = (x2 + y 2 )1/2 , x, y ∈ R
T ⊂S subset T of set S
S∩T the intersection of the sets S and T
S∪T the union of the sets S and T
f (S) image of set S under mapping f
f ◦g composition of two mappings (f ◦ g)(x) = f (g(x))
x column vector in Cn
xT transpose of x (row vector)
0 zero (column) vector
k.k norm
x · y ≡ x∗ y scalar product (inner product) in Cn
x×y vector product in R3
A, B, C m × n matrices
det(A) determinant of a square matrix A
tr(A) trace of a square matrix A
rank(A) rank of matrix A
AT transpose of matrix A
xii
A conjugate of matrix A
A∗ conjugate transpose of matrix A
A† conjugate transpose of matrix A
(notation used in physics)
A−1 inverse of square matrix A (if it exists)
In n × n unit matrix
I unit operator
0n n × n zero matrix
AB matrix product of m × n matrix A
and n × p matrix B
A•B Hadamard product (entry-wise product)
of m × n matrices A and B
[A, B] := AB − BA commutator for square matrices A and B
[A, B]+ := AB + BA anticommutator for square matrices A and B
A⊗B Kronecker product of matrices A and B
A⊕B Direct sum of matrices A and B
δjk Kronecker delta with δjk = 1 for j = k
and δjk = 0 for j 6= k
λ eigenvalue
real parameter
t time variable
Ĥ Hamilton operator
The Pauli spin matrices are used extensively in the book. They are given
by
0 1 0 −i 1 0
σx := , σy := , σz := .
1 0 i 0 0 −1
In some cases we will also use σ1 , σ2 and σ3 to denote σx , σy and σz .
Chapter 1
Basic Operations
a · (b × c) ≡ b · (c × a) ≡ c · (a × b)
1 1 1
b1 = a2 × a3 , b2 = a3 × a1 , b3 = a1 × a2
Va Va Va
1
2 Problems and Solutions
where × denotes the vector product. Are the vectors linearly independent?
T
Problem 3. Consider the normalized vector v0 = ( 1 0 0) in R3 .
Find three normalized vectors v1 , v2 , v3 such that
3
X 1
vj = 0, vjT vk = − (j 6= k) .
j=0
3
Problem 4. Let x ∈ R3 .
(i) Find all solutions of
1 x1 0
0 × x2 = 0 .
0 x3 1
(ii) Find all solutions of
√
0 x1 1/ 2
1 × x2 = 0 .
√
0 x3 1/ 2
Discuss.
calculate?
Problem 9. Let
x1 y1
x = x2 , y = y2
x3 y3
be two normalized vectors in R3 . Assume that xT y = 0, i.e. the vectors
are orthogonal. Is the vector x × y a unit vector again? Here × denotes
the vector product.
Problem 13. (i) Consider the Hilbert space M2 (R) of the 2 × 2 matrices
over R. Show that the matrices
1 0 1 1 1 1 1 1
, , ,
0 0 0 0 1 0 1 1
are linearly independent.
(ii) Use the Gram-Schmidt orthonormalization technqiue to find an or-
thonormal basis for M2 (R).
in matrix form.
0 1 0 0 1
1 0 1 0 1
A= , b= .
0 1 0 1 1
0 0 1 0 1
is normal?
Problem 31. Show that if hermitian matrices S and T are positive semi-
definite and commute (ST = T S), then their product ST is also positive
semi-definite. We have to show that
(ST u)∗ u ≥ 0
for all u ∈ Cn .
Problem 45. Can one find an orthogonal matrix over R such that
0 1 0 0
RT R= ?
0 0 1 0
Problem 46. Let 0n be the n×n zero matrix and In be the n×n identity
matrix. Find an invertible 2n × 2n matrix T such that
−1 0n In 0 −In
T T = .
−In 0n In 0
det g = 1, ηg ∗ η = g −1
where n = 1, 2, . . . and
H(0) = (1) .
Find H(1), H(2), and H(3).
is a projection matrix.
Problem 62. Consider the vector space of 2 × 2 matrices over R and the
matrices
1 1 0 1 0 0 1 0
A1 = , A2 = , A3 = , A4 = .
0 0 0 1 1 1 1 0
Are these matrices linearly independent? Which of these matrices are nor-
mal matrices?
Problem 65. (i) Consider the two-dimensional Euclidean space and let
e1 , e2 be the standard basis
1 0
e1 = , e2 = .
0 1
|xi − xj ||xk − x` |
rijk` := .
|xi − x` ||xk − xj |
A2 = −I2 , A∗ = −A .
Extend to 3 × 3 matrices.
Problem 73. Consider the normalized state u and the permutation ma-
trix P , respectively
1 0 1 0
1
u = √ 1 , P = 0 0 1 .
3 −1 1 0 0
0 1 0 0
0 0 1 0
P = .
0 0 0 1
1 0 0 0
Problem 75. Let 0 ≤ θ < π/4. Note that sec(x) := 1/ cos(x). Consider
the matrix
sec(2θ) −i tan(2θ)
A(θ) = .
i tan(2θ) sec(2θ)
Basic Operations 13
cos(x) 0 − sin(x) 0
0 cos(x) 0 − sin(x)
A(x) =
sin(x) 0 cos(x) 0
0 sin(x) 0 cos(x)
cosh(x) 0 sinh(x) 0
0 cosh(x) 0 sinh(x)
B(x) =
sinh(x) 0 cosh(x) 0
0 sinh(x) 0 cosh(x)
Show that
0 2
k a − b2 − c2 2ab 2ac k
0 1
` = 2ab −a2 + b2 − c2 2bc ` .
a2 + b2 + c2
m0 2ac 2bc −a2 − b2 + c2 m
the faces defined by two of these vectors, normalized to the area of the face,
is given by
1 1 1
n1 = v2 × v3 , n2 = v3 × v1 , n3 = v1 × v2 .
2 2 2
(ii) Show that
2 2 2
v1 = n2 × n3 , v2 = n3 × n1 , v3 = n1 × n2
3V 3V 3V
where V is the volume of the tetrahedron given by
r
1 2
V = (v1 × v2 ) · v3 = (n1 × n2 ) · n3 .
3! 9
S2 := { (x1 , x2 ) : 1 ≥ x1 ≥ x2 ≥ 0 } .
S3 := { (x1 , x2 , x3 ) : 1 ≥ x1 ≥ x2 ≥ x3 ≥ 0 } .
Linear Equations
and
1 1 1 1
A√ =√ . (2)
2 1 2 1
(ii) Do these matrices form a group under matrix multiplication?
x1 + 2x2 − 4x3 + x4 = 3
2x1 − 3x2 + x3 + 5x4 = −4
7x1 − 10x3 + 13x4 = 0 .
16
Linear Equations 17
Let
1 0 1
A = 0 1 0 .
1 0 1
Find out whether system (1) or system (2) has a solution.
Let
2 1 1
A = 1 2 1.
1 1 2
Find out whether system (1) or system (2) has a solution.
Let
1 0 1 1
A = 0 1 0, c = 1 .
1 0 1 1
Find out whether system (1) or system (2) has a solution.
[x] = { y ∈ V : y ∼ x } .
Y + CE + DX = 0n , AE + BX = 0n .
20 Problems and Solutions
Assume that A has an inverse. Eliminate the matrix E and solve the system
for Y
and
x1
W = : x1 ∈ R .
0
(i) Is
3 1 4 −3 3 4
∼ , ∼ , ∼ ?
0 0 1 1 0 1
Problem 18. For the three-body problem the following linear transfor-
mation plays a role
1
X(x1 , x2 , x3 ) = (x1 + x2 + x3 )
3
1
x(x1 , x2 , x3 ) = √ (x1 − x2 )
2
1
y(x1 , x2 , x3 ) = √ (x1 + x2 − 2x3 ) .
6
(i) Find the inverse transformation.
(ii) Introduce polar coordinates
1
x(r, φ) = r sin φ, y(r, φ) = r cos φ, r2 = ((x1 −x2 )2 +(x2 −x3 )2 +(x3 −x1 )2 ) .
3
Express (x1 − x2 ), (x2 − x3 ), (x3 − x1 ) using this coordinates.
Problem 19. Let α ∈ [0, 2π). Find all solutions of the linear equation
cos α sin α x1 b1
= .
sin α cos α x2 b2
and then solve the linear equation. Consider the cases ∆x = ∆y = 1/3 and
∆x = ∆y = 1/4.
n · (x − p) = 0
Problem 22. (i) The equation of a line in the Euclidean space R2 passing
through the points (x1 , y1 ) and (x2 , y2 ) is given by
Apply this equation to the points in R2 given by (x1 , y1 ) = (1, 1/2), (x2 , y2 ) =
(1/2, 1). Consider the unit square with the corner points (0, 0), (0, 1), (1, 0),
(1, 1) and the map
We can consider this as a 2 input AND-gate. Show that the line constructed
above classifies this map.
22 Problems and Solutions
trA = 0, A = A∗ , A2 = I2 .
A2 = tr(A)A .
23
24 Problems and Solutions
The first column of the matrices A and B agree, but the second column of
the two matrices differ. Is
det(A + B) = 2(det(A) + det(B))?
x0 y0 z0 1
1 x1 y1 z1 1
V = det .
6 x2 y2 z2 1
x3 y3 z3 1
Let
(x0 , y0 , z0 ) = (0, 0, 0), (x1 , y1 , z1 ) = (0, 0, 1),
tr(AB) = 0 .
where det[A]{1,...,k} denotes the k-th principal minor, that is, [A]{1,...,k} is
the submatrix by the first k rows and columns of A. The cost of computing
the determinant of a tridiagonal matrix using this recursion is linear in n,
while the cost is cubic for a general matrix. Apply this recursion relation
to calculate the determinant of the 4 × 4 matrix
0 1 0 0
1 1 2 0
A= .
0 2 2 3
0 0 3 3
1 1 0 ... 0 0
1 1 1 ... 0 0
0 1 1 ... 0 0
Mn =
... ... ... .. ....
. . .
0 0 0 ... 1 1
0 0 0 ... 1 1
with the initial conditions det M3 = det M4 = 1. Show that the solution is
√
2 3 nπ π
det(Mn ) = cos − (mod 2) .
3 3 6
Traces, Determinants and Hyperdeterminants 27
Problem 28. Consider the Hilbert space M4 (C) of all 4 × 4 matrices over
C with the scalar product hA, Bi := tr(AB ∗ ), where A, B ∈ M4 (C). The
γ-matrices are given by
0 0 0 −i 0 0 0 −1
0 0 −i 0 0 0 1 0
γ1 = , γ2 =
0 i 0 0 0 1 0 0
i 0 0 0 −1 0 0 0
0 0 −i 0 1 0 0 0
0 0 0 i 0 1 0 0
γ3 = , γ4 =
i 0 0 0 0 0 −1 0
0 −i 0 0 0 0 0 −1
Traces, Determinants and Hyperdeterminants 29
and
0 0 −1 0
0 0 0 −1
γ5 = γ1 γ2 γ3 γ4 = .
−1 0 0 0
0 −1 0 0
We define the 4 × 4 matrices
i
σjk := [γj , γk ], j<k
2
where j = 1, 2, 3, k = 2, 3, 4 and [ , ] denotes the commutator.
(i) Calculate σ12 , σ13 , σ14 , σ23 , σ24 , σ34 .
(ii) Do the 16 matrices
I4 , γ1 , γ2 , γ3 , γ4 , γ5 , γ5 γ1 , γ5 γ2 , γ5 γ3 , γ5 γ4 , σ12 , σ13 , σ14 , σ23 , σ24 , σ34
form a basis in the Hilbert space M4 (C)? If so is the basis orthogonal?
Problem 35. The n×n permutation matrices form a group under matrix
multiplications. Show that
det(In − P ) = 0
where S2n is the symmetric group and sgn(σ) is the signature of permuta-
tion σ. Consider the case with n = 2, i.e.
0 a12 a13 a14
−a 0 a23 a24
A = 12 .
−a13 −a23 0 a34
−a14 −a24 −a34 0
Calculate Pf(A).
Show that
1
det(A) = (−1)n(n−1)/2 (n + 1)nn−1 .
2
tr((V2 V1 )k )
1/2 1/2
can be written as tr(V k ), where V := V2 V1 V2 .
Find the determinant of M . Thus show that the inverse of M exists. Find
the inverse of M .
Is the matrix
1
P3 = I3 − B
1 + z1 z̄1 + z2 z̄2
a projection matrix?
tr(Aj ) = 0, for j = 1, 2, . . . , n .
√
Problem 52. Consider the golden mean number τ = ( 5 − 1)/2 and the
matrix √
√τ τ
F = .
τ −τ
Find tr(F ) and det(F ). Since det(F ) 6= 0 we have an inverse. Find F −1 .
e1 + e2 , e1 + e3 , e2 + e3 .
36 Problems and Solutions
(e1 + e2 e1 + e3 e2 + e3 ) .
Show that this matrix has an inverse and find the inverse.
Problem 62. (i) Let M be a 2×2 matrix over R. Assume that tr(M ) = 0.
Show that
M 2 = − det(M )I2 .
(ii) Show that
√
p
M sin( det(M ))
e = cos( M )I2 + p M.
det(M )
tr(R−1 A) 6= −1 .
38 Problems and Solutions
V 2 = tr(V )V .
1 0 0 ... ... 0 z1
0 1 0 ... ... 0 z2
0 0 1 ... ... 0 z3
. .. .. .. ..
.. ..
A = .. . . . . . . .
. .. .. .. .. .. ..
.. . . . . . .
0 0 0 ... ... 1 zn
z1 z2 z3 ... ... zn 1
e−iφ
0 0
U (φ) = 0 1 0 .
0 0 eiφ
a11 0 0 a12
0 b11 b12 0
A ? B := .
0 b21 b22 0
a21 0 0 a22
(i) Find the determinant and trace of A ? B. Express the result using tr(A),
tr(B), det(A), det(B).
(ii) Assume that the inverse of A and B exists. Is
(A ? B)−1 = A−1 ? B −1 ?
Consider
v1 = (1, 0, 0), v2 = (0, 1, 0), v3 = (0, 0, 1) .
Find the area of the triangle using both expressions. Discuss. The triangle
could be one of the faces of a tetrahedron.
Problem 81. A tetrahedron has four triangular faces. Given the coordi-
nates of the four vertices
x0 y0 z0 1
1 x y1 z1 1
V = det 1
3! x2 y2 z2 1
x3 y3 z3 1
√ √
(i) Given the four vertices (1, 0, − 2), (2, 0, 0), (0, 0, 0), (1, − 2, 0) find the
volume.
(ii) Derive an equation for surface area of a tretrahedron given by coordi-
nates. Apply it to the vertices given in (i).
Problem 84. Let 00 = 11 = 0, 01 = 1, 10 = −1, i.e. we consider the
2 × 2 matrix
0 1
= .
−1 0
Then the determinant of a 2 × 2 matrix A2 = (aij ) with i, j = 0, 1 can be
defined as
1 1 1 1
1 XXX X
det A2 := ij `m ai` ajm .
2 i=0 j=0 m=0
`=0
Thus
det A2 = a00 a11 − a01 a10 .
In analogy the hyperdeterminant of the 2 × 2 × 2 array A3 = (aijk ) with
i, j, k = 0, 1 is defined as
1 1 1 1 1 1
1 X X X X X X
DetA3 := − ii0 jj 0 kk0 mm0 nn0 pp0 aijk ai0 j 0 m anpk0 an0 p0 m0 .
2 0
ii =0 jj =0 kk =0 mm0 =0 nn0 =0 pp0 =0
0 0
Calculate DetA3 .
A = (ajk` ), j, k, ` = 0, 1
Assume that det(S) = 1, i.e. S ∈ SL(2, C). Show that Det(AS) = Det(A).
This is a typical problem to apply computer algebra. Write a SymbolicC++
program or Maxima program that solves the problem.
0 a1 0 0 0 a1 0 0
1 a1 0 a2 0 1 −a1 0 a2 0
A= , B= .
2 0 a2 0 a3 2i 0 −a2 0 a3
0 0 a3 0 0 0 −a3 0
det(A + B) = det(A) .
(M )jk = xj−1
k , j, k = 1, 2, 3
Eigenvalues and
Eigenvectors
0 0 0 1
0 0 0 1
A= .
0 0 0 1
1 1 1 0
A2 − AtrA + I2 det A = 0
43
44 Problems and Solutions
and therefore
(trA)2 = trA2 + 2 det A .
Hint. Apply the Cayley-Hamilton theorem.
Problem 5. Find all 2 × 2 matrices over R that admit only one eigen-
vector.
Can one find a condition on the parameter a so that A has only one eigen-
vector?
ej = V ∗ Aj V
A
(I2 − A)−1 = I2 + A + A2 + · · ·
0 1 0 0
0 1 0
0 1 0 0 1 0
, 0 0 1 , .
c1 c2 0 0 0 1
c1 c2 c3
c1 c2 c3 c4
Problem 14. (i) Find the eigenvalues and normalized eigenvectors of the
3 × 3 matrix
0 1 0
M = 1 0 1 .
0 1 0
(ii) Use the normalized eigenvectors to construct a 3 × 3 matrix R such that
RM R−1 is a diagonal matrix.
46 Problems and Solutions
0 1 0 0 1
0 1 0 1
0 1 1 1 0 1 0 0
1 0 1 0
A3 = 1 0 1, A4 = , A5 = 0 1 0 1 0 .
0 1 0 1
1 1 0 0 0 1 0 1
1 0 1 0
1 0 0 1 0
(ii) Extend the results from (i) to find the largest eigenvalue of the sym-
metric n × n matrix
0 1 0 ··· 0 0 1
1 0 1 ··· 0 0 0
0 1 0 ··· 0 0 0
An = .. .. .. . . . .. .. .. .
. . . . . .
0 0 0 ... 1 0 1
1 0 0 ... 0 1 0
−1 α 0 0
α −1/2 α 0
.
0 α 1/2 α
0 0 α 1
2k
4k
1 ..
≡ eiπ/n ,
√ . , 1 ≤ k ≤ n.
n
2(n−1)k
1
a 0 0 b
0 a 0 b
.
0 0 a b
b b b 0
a 0 0 b
0 a 0 b
.
0 0 −a b
b b b 0
Show that the two matrices have the same (normalized) eigenvectors. Find
the commutator [S, T ].
and
I2 02
γ0 = .
02 −I2
Find γ1 γ2 γ3 γ0 and tr(γ1 γ2 γ3 γ0 ).
1 c 0 0
c 2 2c 0
A= .
0 2c 3 c
0 0 c 4
λ1 + λ2 + λ3 + λ4 = 10
λ1 λ2 + λ1 λ3 + λ1 λ4 + λ2 λ3 + λ2 λ4 + λ3 λ4 = 35 − 6c2
λ1 λ2 λ3 + λ1 λ2 λ4 + λ1 λ3 λ4 + λ2 λ3 λ4 = 50 − 30c2
λ1 λ2 λ3 λ4 = 24 − 30c2 + c4
1 0 0 0 1
0 1 1 1 0
A = 0 1 1 1 0 .
0 1 1 1 0
1 0 0 0 1
using the trace and the determinant of the matrix and the information that
two eigenvalues are the same.
50 Problems and Solutions
0 0 1 0
0 0 0 1
X= .
b11 b12 0 0
b21 b22 0 0
If λ ∈ λ(A, B) and
Ax = λBx, x 6= 0
then x is referred to as an eigenvector of A − λB. Note that λ may be
finite, empty or infinite.
Let
0 0 0 1 1 0 0 1
0 0 1 0 1 0 1 1 0
A= , B=√ .
0 1 0 0 2 0 1 −1 0
1 0 0 0 1 0 0 −1
Find the eigenvalue of the pencil.
Eigenvalues and Eigenvectors 51
Problem 43. Let AT = (1/2 , 1/2)T . Find the eigenvalues of AAT and
AT A.
Problem 44. We know that a hermitian matrix has only real eigenvalues.
Can we conclude that a matrix with only real eigenvalues is hermitian?
Problem 45. Let A be an n×n matrix over C. Show that the eigenvalues
of A∗ A are nonnegative.
Cjk := x|j−k| , j, k = 1, . . . , N .
1 − x2
λ(φ) = .
1 − 2x cos φ + x2
Show that the eigenvalues of Tn are real and not of absolute value 1.
(ii) Show that the eigenvalues of the m×m matrix AAT are also eigenvalues
of AT A.
0 0 0 a14
0 0 a13
0 a12 1 0 0 a24
A2 = , A3 = 1 0 a23 , A4 = .
1 a22 0 1 0 a34
0 1 a33
0 0 1 a44
ln(A) .
e = U ∗ AU is a diagonal matrix.
transforms A into a diagonal matrix, i.e. A
Eigenvalues and Eigenvectors 55
0 1
0 0 1
0 0 0 1
H= ..
.
1
eiφ 0 0
D = diag(1, ω, ω 2 , . . . , ω 2j )
where ω := exp(i2π/(2j + 1). Is H unitary? Find ωDH − HD.
tr(A) = c1 , tr(A2 ) = c2 .
Ax = ρx, yT A = ρy, yT x = 1 .
Show that B has a positive spectral radius. Find the vectors x and y.
Problem 63. Show that the resolvent satisfies the so-called resolvent
equation
Rλ − Rµ = (λ − µ)Rλ Rµ .
√
Problem 64. Let τ = (1 + 5)/2 be the golden ratio. Consider the
modular matrix
1 1
M= .
1 0
Find the eigenvalues λ1 , λ2 and the corresponding normalized eigenvec-
tors u1 , u2 . Find the projection matrices Π1 and Π2 onto the associated
eigendirections.
A = B −1 DB
is not normal.
0 1
0 1
0 1
H= . ..
.
1
eiφ 0
G = diag(1, ω, ω 2 , · · · ω n−1 )
58 Problems and Solutions
Problem 76. (i) Find the eigenvalues and eigenvectors of the orthogonal
matrices
1 0 0 cos β sin β 0
R = 0 cos α sin α , S = − sin β cos β 0 .
0 − sin α cos α 0 0 1
Eigenvalues and Eigenvectors 59
Problem 81. Find all 2 × 2 matrices over the real numbers with only
one 1-dimensional eigenspace, i.e. all eigenvectors are linearly dependent.
60 Problems and Solutions
(A − λIn )j x = 0n
for some j ∈ {1, 2, . . .}, where In is the n × n identity matrix and 0n is the
n-dimensional zero vector. For j = 1 we find the eigenvectors. It follows
that x is a generalized eigenvector of A corresponding to λ if and only if
(A − λIn )n x = 0n .
and generally
√ √
1/ 2 0√ ... 0 0 0 ... 0√ 1/ 2
0 1/ 2 ... 0 0 0 ... 1/ 2 0
.. .. .. .. .. .. ..
. . ... . . . ... . .
√ √
0 0 ... 1/ 2 0 1/ 2 ... 0 0
An =
0 0 ... 0√ 1 0√ ... 0 0 .
0
0 ... 1/ 2 0 −1/ 2 ... 0 0
. .. .. .. .. .. ..
..
√. ... . . . ... .√ .
0
√ 1/ 2 ... 0 0 0 . . . −1/ 2 0√
1/ 2 0 ... 0 0 0 ... 0 −1/ 2
Eigenvalues and Eigenvectors 61
where c ∈ C.
plays a role for the chemical compounds ZnS and N aCl. Find the eigen-
values and eigenvectors of A. Then find the inverse of A. Find all x such
that Ax = x.
62 Problems and Solutions
Write a C++ program that calculates the right-hand side of the inequality
for a given matrix. Apply the complex class of STL. Apply it to the matrix
i 0 0 i
0 2i 2i 0
A= .
0 3i 3i 0
4i 0 0 4i
66 Problems and Solutions
Problem 106. Find the condition on a11 , a12 , b11 , b12 such that
a11 0 0 a12 1 1
0 b11 b12 0 1 1
= λ
0 b12 b11 0 1 1
a12 0 0 a11 1 1
i.e. we have an eigenvalue equation.
Problem 110. (i) Find the spectral decomposition of the normal matrix
0 1 0
A = 1 0 1 .
0 1 0
A2 ≡ AA = cA
(ii) Take the determinant of both sides of the equation. Discuss. Study the
cases that A is invertible and non-invertible.
(iii) Study the case
−z
e 1
A(z) = , z ∈ C.
1 ez
(iv) Study
(A ⊗ A)2 = c(A ⊗ A) .
(v) Let A be a 2 × 2 matrix and
a11 0 0 a12
0 a11 a12 0
A ? A := .
0 a21 a22 0
a21 0 0 a22
Problem 115. (i) Consider the Pauli spin matrices for describing a spin- 21
system
0 1 0 −i 1 0
σ1 = , σ2 = , σ3 = .
1 0 i 0 0 −1
n1 sx + n2 sy + n3 sz .
n1 sx + n2 sy + n3 sz .
0 0 0 a14
0 0 a13
1 0 a23 , 1 0 0 a24
.
0 1 0 a34
0 1 a33
0 0 1 a44
(iv) Apply the result from (ii) and (iii) to the unitary 2 × 2 matrix
cos θ sin θ
U (θ) = .
− sin θ cos θ
(v) Apply the result from (ii) and (iii) to the 2 × 2 unitary matrix
cos θ −eiφ sin θ
V (θ, φ) = .
e−iφ sin θ cos θ
(vi) Every hermitian matrix H can be written as H = iK, where K is a
skew-hermitian matrix. Find H for the examples given above.
Problem 125. Let a11 , a22 ∈ R and a12 ∈ C. Consider the hermitian
matrix
a11 a12
H=
ā12 a22
with the real eigenvalues λ1 and λ2 . What conditions are impost on the
matrix elements of H if λ1 = λ2 ?
Problem 126. (i) Consider the spin matrices for describing a spin- 21
system
1 0 1 1 0 −i 1 1 0
s1 = , s2 = , s3 = .
2 1 0 2 i 0 2 0 −1
and the spin matrices for describing a spin-1 system
0 1 0 0 −i 0 1 0 0
1 1
p1 = √ 1 0 1 , p2 = √ i 0 −i , p3 = 0 0 0 .
2 0 1 0 2 0 i 0 0 0 −1
72 Problems and Solutions
Ĥ
K̂ = = s1 ⊗ p1 ⊗ s1 + s2 ⊗ p2 ⊗ s2 + s3 ⊗ p3 ⊗ s3 .
~ω
Problem 127. s`(3, R) is the rank 2 Lie algebra with Cartan matrix
2 −1
C=
1 2
0 eiφ1 0 0
iφ2
e 0 0 0
A(φ1 , φ2 , φ3 , φ4 ) = iφ3 ,
0 0 0 e
0 0 eiφ4 0
0 0 eiφ5 0
iφ6
0 0 0 e
B(φ5 , φ6 , φ7 , φ8 ) = iφ7
e 0 0 0
0 eiφ8 0 0
and A(φ1 , φ2 , φ3 , φ4 )B(φ5 , φ6 , φ7 , φ8 ). Find the eigenvalues of these matri-
ces.
Problem 135. We know that any n × n unitary matrix has only eigenval-
ues λ with |λ| = 1. Assume that a given n × n matrix has only eigenvalues
with |λ| = 1. Can we conclude that the matrix is unitary?
(vvT )v = v(vT v) .
Discuss.
Problem 139. (i) Find the eigenvalues and eigenvectors of the matrix
a11 0 0 0
a21 a22 0 0
A= .
0 a32 a33 0
0 0 a43 a44
b11 0 0 0 0 0
b21 b22 0 0 0 0
0 b32 b33 0 0 0
B= .
0 0 b43 b44 0 0
0 0 0 b54 b55 0
0 0 0 0 b65 b66
These matrices are the so-called staircase matrices. Extend the results to
the n × n case.
0 0 0 0 0 1
1 0 0 0 0 0
0 1 0 0 0 0
A= .
0 0 1 0 0 0
0 0 0 1 0 0
0 0 0 0 1 0
Find the trace and determinant of these matrices. Show that for the matrix
A(α) the eigenvalues depend on α but the eigenvectors do not. Show that
for the matrix B(α) the eigenvalues do not depend on α but the eigenvectors
do.
(ii) Let α ∈ R. Consider the matrices
cosh α sinh α cosh α − sinh α
C(α) = , D(α) = .
sinh α cosh α sinh α − cosh α
Find the trace and determinant of these matrices. Show that for the matrix
C(α) the eigenvalues depend on α but the eigenvectors do not. Show that
for the matrix D(α) the eigenvalues do not depend on α but the eigenvectors
do.
76 Problems and Solutions
Av = λBv .
(ii) Use the spectral decomposition to find the matrices A1 and A2 such
that P1 = exp(A1 ), P2 = exp(A2 ).
Commutators and
Anticommutators
78
Commutators and Anticommutators 79
where µ ∈ R and µ 6= 0?
A1 ⊗A2 ⊗A3 +A2 ⊗A3 ⊗A1 +A3 ⊗A1 ⊗A2 −A1 ⊗A3 ⊗A2 −A2 ⊗A1 ⊗A3 −A3 ⊗A2 ⊗A1 .
[H, A] = 0, [H, B] = 0 .
Show that
[H ⊕ In + In ⊕ H, A ⊕ B] = 0
where ⊕ denotes the direct sum.
Problem 25. Show that any two 2 × 2 matrices which commute with the
matrix
0 1
−1 0
commute with each other.
∂/∂xn
c ⊕ A, A⊕c
Let
B1 0 B
e1 0
B= , B
e=
0 B2 0 B
e2
and
0 F1 0 Fe1
F = , Fe = .
F2 0 Fe2 0
Find the commutators [B.B],
e [B, F ] and the anticommutator [F, Fe]+ .
Commutators and Anticommutators 83
Problem 38. Find all 2 × 2 matrices over C such that the commutator
is an invertible diagonal matrix D, i.e. d11 6= 0 and d22 6= 0.
1 0 0 0
0 0 1 0
P = .
0 1 0 0
0 0 0 1
Decomposition of
Matrices
will be called a plane rotation through φ in the plane span (ep , er ). Let
Q = (qij )1≤i,j≤3 be a rotation matrix. Show that there exist angles φ ∈
[0, π), θ, ψ ∈ (−π, π] called the Euler angles of Q such that
86
Decomposition of Matrices 87
A∗ A = I k
where Ik is the k × k unit matrix. Find the n × n matrix AA∗ using the
singular value decomposition. Calculate tr(AA∗ ).
A∗ A = I k
0 1 0 0 0 0 1 0
−1 0 1 0 e = 0 0 0 1
Ω= , Ω .
0 −1 0 1 −1 0 0 0
0 0 −1 0 0 −1 0 0
Ω = P ΩQ
e ?
0 1 0 ... 0 0 0
−1 0 1 ... 0 0
0 −1 0 ... 0 0 0
Ω= ..
.
0 0 0 ... −1 0 1
0 0 0 ... 0 −1 0
Decomposition of Matrices 89
and
0n In
Ω
e=
−In 0n
where 0n is the n × n zero matrix and In is the n × n identity matrix. Can
one find 2n × 2n permutation matrices P , Q such that Ω = P ΩQ?e
AX = XB . (1)
V −1 AV = JA , W −1 BW = JB
JA Y = Y JB
where Y := V −1 XW .
Chapter 7
Functions of Matrices
cosh ζ 0 0 sinh ζ
0 cosh ζ sinh ζ 0
S(ζ) = .
0 sinh ζ cosh ζ 0
sinh ζ 0 0 cosh ζ
(i) Show that the matrix is invertible, i.e. find the determinant.
(ii) Calculate the inverse of S(ζ).
(iii) Calculate
d
A := S(ζ)
dζ ζ=0
Note that
∞
Z r
2 π (b2 −4ac)/(4a)
dqe−(aq +bq+c)
= e . (1)
−∞ a
90
Functions of Matrices 91
(ii) Let
0 1 0
B = 0 0 1
1 0 0
with B 3 = I3 . Calculate exp(B) using the result from (i). Calculate exp(B)
applying the Cayley-Hamilton theorem.
0 1 0 ... 0
0 0 1 ... 0
. . .. . . .
A= . . . ..
. . . .
0 0 0 ... 1
0 0 0 ... 0
if AB = BA.
94 Problems and Solutions
Find exp(A).
[A, B]+ = AB + BA = 0n .
Let a, b ∈ C. Calculate
eaA+bB .
[A, B] = AB − BA = 0n .
Let a, b ∈ C. Calculate
eaA+bB .
Find exp(tA).
n · σ := n1 σ1 + n2 σ2 + n3 σ3 .
Then
1 1
f (θn · σ) ≡ (f (θ) + f (−θ))I2 + (f (θ) − f (−θ))(n · σ) .
2 2
Apply this identity to f (x) = sin(x).
Calculate exp(M ).
Let
0 1
A= .
1 0
Calculate etA using this method.
eλ1 t − eλ2 t
[λ1 , λ2 ] :=
λ1 − λ2
[λ1 , . . . , λk ] − [λ2 , . . . , λk+1 ]
[λ1 , . . . , λk+1 ] := , k ≥ 2.
λ1 − λk+1
Let
0 1
A= .
1 0
Calculate etA using this method.
is closely related the the resolvent (λIn − A)−1 trough the formula
N1 λn−1 + N2 λn−2 + · · · + Nn N (λ)
(λIn − A)−1 = =
λn + a1 λn−1 + · · · + an p(λ)
where the adjugate matrix N (λ) is a polynomial in λ of degree n − 1 with
constant n × n coefficient matrices N1 , . . . , Nn . The Laplace transform of
the matrix exponential is the resolvent
Let
1 0 0 1
0 −1 1 0
A= .
0 1 −1 0
1 0 0 1
Find the Nk matrices and the coefficients ak and thus calculate the resol-
vent.
is closely related the the resolvent (λIn − A)−1 trough the formula
N1 λn−1 + N2 λn−2 + · · · + Nn N (λ)
(λIn − A)−1 = =
λn + a1 λn−1 + · · · + an p(λ)
98 Problems and Solutions
1
N1 = In , a1 = − tr(AN1 )
1
1
N2 = AN1 + a1 In , a2 = − tr(AN2 )
2
..
.
1
Nn = ANn−1 + an−1 In , an = − tr(ANn )
n
0 = ANn + an In .
Show that
p0 (λ)
tr(L(et A)) =
p(λ)
where p0 (λ) = dp(λ)/dλ.
Calculate f (A).
eA BeA .
Set f () = eA BeA , where is a real parameter. Then differentiate with
respect to . For = 1 we have eA BeA .
Let
2 −1 1 1
A= , B= .
−1 1 1 2
Calculate the left and right-hand side of the integral representation.
Problem 51. Consider the Pauli spin matrices σx , σy and σz . Can one
find an α ∈ R such that
exp(iασz )σx exp(−iασz ) = σy ?
102 Problems and Solutions
up to second order in .
Find exp(tB), where t ∈ R and thus solve the initial value problem of the
matrix differential equation
dA
= BA(t) .
dt
RA + A∗ R = −B .
[eA , eB ]
104 Problems and Solutions
(ii) Use the method given above to calculate exp(iK), where the hermitian
2 × 2 matrix K is given by
a b
K= , a, c ∈ R, b ∈ C .
b c
ezA Be−zA = B .
Obviously e−zA is the inverse of ezA . Show that, if this condition is satisfied,
one has [A, B] = 0n , where 0n is the n × n zero matrix. If ezA would be
unitary we have U BU ∗ = B.
Find the commutators [Σx , Σy ], [Σy , Σz ], [Σz , Σx ] and compare with the
commutators [σx , σy ], [σy , σz ], [σz , σx ].
106 Problems and Solutions
Problem 71. We know that for any n × n matrix A over C the matrix
exp(A) is invertible with the inverse exp(−A). What about cos(A) and
cosh(A)?
is closely related the the resolvent (λIn − A)−1 trough the formula
Let
1 0 0 1
0 −1 1 0
A= .
0 1 −1 0
1 0 0 1
Find the Nk matrices and the coefficients ak and thus calculate the resol-
vent.
w ∗ eA v
Find cos(A) and the inverse of this matrix. Find cos(B) and the inverse of
this matrix. Find the commutators [A, B] and [cos(A), cos(B)]. Discuss.
V = v0 I2 + v1 σ1 + v2 σ2 + v3 σ3
dR(t) T
H(t) = i~ R (t)
dt
and show it is hermitian.
110 Problems and Solutions
exp(c(A − I2 )) = A .
where
∞
X
C(zA, zB) = cj (A, B)z j .
j=1
c1 (A, B) = A + B
1
c2 (A, B) = [A, B]
2
1 1
c3 (A, B) = [A, [A, B]] − [B, [A, B]
12 12
1
c4 (A, B) = − [A, [B, [A, B] .
24
Linear Differential
Equations
du(t)
+ Au(t) = g(t), u(0) = u0 (1)
dt
where g(t) = (g1 (t), g2 (t), . . . , gn (t))T . The solution of the initial value
problem is
Z t
u(t) = e−tA u0 + e−(t−τ )A g(τ )dτ . (2)
0
(i) Discretize the system with the implicit Euler method with step size h.
(ii) Compare the two solutions of the two systems for the matrix
0 0 1
A = 0 1 0,
1 0 0
initial values u0 = (1 , 1 , 1)T with g(t) = (1, 0, 1)T and the step size
h = 0.1.
112
Linear Differential Equations 113
that K −1 (t) exists for all t. Assume that the time-evolution of L is given
by
L(t) = K(t)L(0)K −1 (t) .
(i) Show that L(t) satisfies the matrix differential equation
dL
= [L, B](t)
dt
where [ , ] denotes the commutator and
dK −1
B=− K (t) .
dt
(ii) Show that if L(t) is hermitian and K(t) is unitary, then the matrix B(t)
is skew-hermitian.
which does not depend on the initial time t0 . The matrix exp(T R) is called
the monodromy matrix of the periodic system (1). Calculate
tr exp(T R) .
Problem 4. Let
1 1
A= .
0 1
114 Problems and Solutions
Problem 5. Solve the initial value problem for the matrix differential
equation
dA
[B, A()] =
d
where A() and B are 2 × 2 matrices with
0 1
B= .
1 0
X(t) = exp(Ω(t))
d2 X dX
2
+a + bX = 0 .
dt dt
Linear Differential Equations 115
d2 u du
2
=A = A2 u .
dt dt
Thus we can write
du dv
= v = Au, = A2 u = Av
dt dt
or in matrix form
du/dt 0n In u(0)
=
dv/dt A2 0n v(0)
kM k = max kM xk
kxk=1
where kxk denotes the Euclidean norm. Show that if kU1 − U2 k ≤ then
kU1 v − U2 vk ≤ .
116
Norms and Scalar Products 117
Let
0 1 2
A = 3 4 5 .
6 7 8
Calculate ρ(A) and the right-hand side of the inequality.
for all x, y ∈ C.
Find the condition on t such that ρ(A(t)) < 1, where ρ(A(t)) denotes the
spectral radius of A(t).
U ∗ AU = diag(λ1 , . . . , λn ) + N
118 Problems and Solutions
where
p
X (p + q − j)!p!
Npq (A) := Aj
j=0
(p + q)!j!(p − j)!
q
X (p + q − j)!q!
Dpq (A) := (−A)j .
j=0
(p + q)!j!(q − j)!
Let
0 6 0 0
0 0 6 0
A= .
0 0 0 6
0 0 0 0
Calculate kR11 − eA k, where k · k denotes the 2-norm.
kAkHS = kAkO .
Calculate k[A, B]k, kABk, kBAk, kAk, kBk and pthus verify the inequality
for these matrices. The norm is given by kCk = tr(CC ∗ ).
Problem 17. Consider the Hilbert space Rn . The scalar product hx, yi
x, y ∈ Rn is given by
n
X
hx, yi := xT y = xj yj .
j=1
120 Problems and Solutions
p
Thus the norm is given by kxk = hx, xi. Show that
Au − λu = v .
(i) Show that for λ nonreal (i.e. it has an imaginary part) the vector v
cannot vanish unless u vanishes.
(ii) Show that for λ nonreal we have
1
k(A − λIn )−1 vk ≤ kvk .
|=λ|
Is B invertible?
(ii) Let A, B, C be invertible n × n matrices over R with
Is
kA−1 k ≥ kB −1 k ≥ kC −1 k ?
Norms and Scalar Products 121
Show that kM kF is the square root of the sum of the squares of the singular
values of M .
kM − AkF .
kΣ − A0 kF .
where kAxk denotes the Euclidean norm. Is ρ(A) ≤ kAk? Prove or dis-
prove.
Chapter 10
0 1 1 0 0 0 0
1 0 1 1 0 0 0
1 1 0 1 0 1 0
A = 0 1 1 0 0 1 0 .
0 0 0 0 0 1 0
0 0 1 1 1 0 1
0 0 0 0 0 1 0
123
124 Problems and Solutions
Hadamard Product
A•B=B •A
A • (B + C) = A • B + A • C
A • (λB) = λ(A • B) .
n
Y
det(A • B) ≥ det B ajj (1)
j=1
First show that A and B are positive definite and then calculate the left
and right-hand side of (1).
125
126 Problems and Solutions
kA • Bk ≤ kA∗ Ak kB ∗ Bk
Problem 4. If V and W are matrices of the same order, then their Schur
product V • W is defined by (entrywise multiplication)
If all entries of V are nonzero, then we say that X is Schur invertible and
define its Schur inverse, V (−) , by V (−) • V = J, where J is the matrix with
all 1’s.
XC M := CM, ∆C M := C • M, YC := M C T .
XA ∆B XA = ∆B XA ∆B .
We call this the braid relation. Give an example for a one-sided Jones pair.
(Aej ) • (Bek ), j, k = 1, . . . , n .
Hadamard Product 127
sT = ( 1 1 ... 1)
A • A−1
is invertible?
where n = 1, 2, . . . and
H(0) = (1) .
Find H(1), H(2) and H(3).
Unitary Matrices
U = In − 2uu∗ .
129
130 Problems and Solutions
A∗ U A = U
U1 = V U2 V ∗ ?
Problem 13. (i) Find the condition on the n × n matrix A over C such
that In + A is a unitary matrix.
(ii) Let B be an 2 × 2 matrix over C. Find all solutions of the equation
B + B ∗ + BB ∗ = 02 .
A + A−1 = I2 .
U = α1 σ1 + α2 σ2 + α3 σ3
is a unitary matrix.
unitary?
(ii) Is the 4 × 4 matrix
1 σ0 σ1
√
2 −iσ2 σ3
unitary?
unitary?
1
hU, V i := tr(U V ∗ ) .
n
Λ := { Σj C k Ω` : 0 ≤ j ≤ 3, 0 ≤ k ≤ 1, 0 ≤ ` ≤ 2 }
1 0 0 0 0 0 0 1
0 ω 0 0 1 0 0 0
σ= , Γ= .
0 0 ω2 0 0 1 0 0
3
0 0 0 ω 0 0 1 0
Unitary Matrices 135
Let c > 0. The four-state Potts quantum chain is defined by the Hamilton
operator
N
1 X
(σj + σj2 + σj3 ) + c(Γj Γ3j+1 + Γ2j Γ2j+1 + Γ3j Γj+1 )
Ĥ = − √
π c j=1
where N is the number of sites and one imposes cyclic boundary conditions
N + 1 ≡ 1. Let N = 2. Find the eigenvalues and eigenvectors of Ĥ.
Problem 37. Let U be a unitary and symmetric matrix. Show that there
exists a unitary and symmetric matrix V such that U = V 2 .
Problem 43. (i) What are the conditions on φ11 , φ12 , φ21 , φ22 ∈ R such
that iφ
1 e 11 eiφ12
U (φ11 , φ12 , φ21 , φ22 ) = √ iφ iφ
2 e 21 e 22
is a unitary matrix?
(ii) What are the condition on φ11 , φ12 , φ21 , φ22 ∈ R such that U (φ11 , φ12 , φ21 , φ22 )
is an element of SU (2)?
Chapter 13
Numerical Methods
Cx(k+1) = Rx(k) + b, k = 0, 1, 2, . . . .
Let
4 −1 0 4 0 0
A = −1 4 −1 , C = 0 4 0
0 −2 4 0 0 4
and
3 0
b = 2, x(0) = 0 .
2 0
The iteration converges if ρ(C −1 R) < 1, where ρ(C −1 R) denotes the spec-
tral radius of C −1 R. Show that ρ(C −1 R) < 1. Perform the iteration.
138
Numerical Methods 139
Use k = 2.
1
keA − fj,k (A)k ≤ kAkk+1 ekAk (2)
j k (k + 1)!
Let
0 1
A= .
0 0
Find f2,2 (A) and eA . Calculate the right-hand side of the inequality (2).
Numerical Methods 141
M xt
xt+1 = , t = 0, 1, . . .
kM xt k
1 0 0 1 1
1 0 1 1 0 0
B=√ and x0 = .
2 0 1 −1 0 0
1 0 0 −1 0
Chapter 14
Binary Matrices
142
Binary Matrices 143
is given by
det A = a11 a22 a33 + a12 a23 a31 + a13 a21 a32
−a13 a22 a31 − a11 a23 a32 − a12 a21 a33 .
det B = (b11 · b22 · b33 ) ⊕ (b12 · b23 · b31 ) ⊕ (b13 · b21 · b32 )
⊕(b13 · b22 · b31 ) ⊕ (b11 · b23 · b32 ) ⊕ (b12 · b21 · b33 ) .
Problem 3. The finite field GF (2) consists of the elements 0 and 1 (bits)
which satisfies the following addition (XOR-operation) and multiplication
(AND-operation) tables
⊕ 0 1 · 0 1
0 0 1 0 0 0
1 1 0 1 0 1
Ty = s
T −1 s = y .
For T we select the Hadamard matrix. The 2n ×2n Hadamard matrix H(n)
is recursively defined as
H(n − 1) H(n − 1)
H(n) = , n = 1, 2, . . .
H(n − 1) −H(n − 1)
Find the truth table, the vector y and then calculate, using H(3), the
spectral coefficients sj , (j = 0, 1, . . . , 7).
Problem 1. The Pauli matrix σx is not only hermitian, unitary and his
own inverse, but also a permutation matrix. Find all 2 × 2 matrices A such
that
σx−1 Aσx = A .
145
146 Problems and Solutions
Show that these matrices form a group, i.e. a subgroup of the 3 × 3 per-
mutation matrices.
Show that the Lie algebra s`(2, R) has no proper nontrivial ideals.
hA, Bi := tr(AB T ) .
hX, Y i = 0
Groups, Lie Groups and Lie Algebras 147
Let z ∈ C. Find
0 b−c −i(b + c)
a b
exp z , exp z c − b 0 2ia .
c −a
i(b + c) −2ia 0
(i) Map the standard basis of R4 into SU (2) and express these matrices
using the Pauli spin matrices and the 2 × 2 identity matrix.
(ii) Map the Bell basis
1 1 0 0
1 0 1 0 1 1 1 1
√ , √ , √ , √
2 0 0 2 1 2 −1
2
1 −1 0 0
148 Problems and Solutions
into SU (2) and express these matrices using the Pauli spin matrices and
the 2 × 2 identity matrix.
Problem 12. In the decomposition of the simple Lie algebra s`(3, R) one
finds the 3 × 3 matrices
a11 a12 0 0 0 b13
A = a21 a22 0 , B = 0 0 b23 .
0 0 −a11 − a22 b31 b32 0
Consider
0 0 0 1 1 0 0 0 0 0 0 0
0 0 1 0 0 1 0 0 0 0 0 0
X=
e , H=
e , Y =
e .
0 0 0 0 0 0 −1 0 0 1 0 0
0 0 0 0 0 0 0 −1 1 0 0 0
Find the commutators
[X,
e H],
e [X,
e Ye ], [Ye , H].
e
Problem 18. Consider the Lie group SL(n, C), i.e. the n × n matrices
over C with determinant 1. Can we find A, B ∈ SL(n, C) such that [A, B]
is an element of SL(n, C)?
Both are elements of the non-compact Lie group SL(2, C). Can one finds
α ∈ C such that the commutator [A(α), B] is again an element of SL(2, C)?
Problem 22. The simple Lie algebra s`(2, R) has a basis given by the
matrices
1 0 0 1 0 0
H= , E= , F =
0 −1 0 0 1 0
with the commutation relations
[H, E] = 2E, [H, F ] = −2F, [E, F ] = H .
The universal enveloping algebra U (s`(2, R)) of the Lie algebra s`(2, R) is
the associative algebra with generators H, E, F and the relations
HE − EH = 2E, HF − F H = −2F, EF − F E = H .
Find a basis of the Lie algebra s`(2, R) so that all matrices are invert-
ible. Find the inverse matrices of these matrices. Give the commutation
relations.
Problem 23. A Chevalley basis for the semisimple Lie algebra s`(3, R)
is given by
0 0 0 0 1 0 0 0 1
X1 = 0 0 1 , X2 = 0 0 0 , X3 = 0 0 0
0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0
Y1 = 0 0 0 , Y2 = 1 0 0 , Y3 = 0 0 0
0 1 0 0 0 0 1 0 0
0 0 0 1 0 0
H1 = 0 1 0 , H2 = 0 −1 0
0 0 −1 0 0 0
where Yj = XjT for j = 1, 2, 3. The Lie algebra has rank 2 owing to H1 , H2
and [H1 , H2 ] = 0. Another basis could be formed by looking at the linear
combinations
Uj = Xj + Yj , Vj = Xj − Vj .
Groups, Lie Groups and Lie Algebras 151
Using the expansion coefficients hjk one can associate to the hermitian
matrix H the function
N
X −1 N
X −1
h(q, q) = hjk e2πi(jp+kq)
j=0 k=1
Problem 29. (i) The standard basis for the vector space of the 2 × 2
matrices is given by
1 0 0 1 0 0 0 0
E11 = , E12 = , E21 = , E22 = .
0 0 0 0 1 0 0 1
We define the star composition of two 2 × 2 matrices as the 4 × 4 matrix
a11 0 0 a12
0 b11 b12 0
A ? B := .
0 b21 b22 0
a21 0 0 a22
Groups, Lie Groups and Lie Algebras 153
form a basis of the Lie algebra s`(2, R). Do the nine 4 × 4 matrices
X ? X, X ? Y, X ? H, Y ? X, Y ? Y, Y ? H, H ? X, H ? Y, H ? H
form a basis of the simple Lie algebra s`(2, R). Define the matrices
Discuss.
CAC T = A .
Note that C T = C −1 .
CAC T = A .
Note that C T = C −1 .
Show that
1 X
g
|G|
g∈G
is a projection matrix.
where ajk ∈ R, ajk > 0 for j, k = 1, 2 and a12 = a21 . We also have φjk ∈ R
for j, k = 1, 2 and impose φ12 = φ21 . What are the conditions on ajk and
φjk such that I2 + iA is a unitary matrix?
SAS −1 = A .
Discuss.
Let (3-sphere)
Show that SU (2) can be identified as a real manifold with the 3-sphere S 3 .
Find
z = xyx−1 y −1 .
Show that xz = zx and yz = zy.
(iii) The derived subgroup (or commutator subgroup) of a group G is the
subgroup [G, G] generated by the set of commutators of every pair of ele-
ments of G. Find [G, G] for the Heisenberg group.
(iv) Let
0 a c
A = 0 0 b
0 0 0
and a, b, c ∈ R. Find exp(A).
158 Problems and Solutions
(v) The Heisenberg group is a simple connected Lie group whose Lie algebra
consists of matrices
0 a c
L = 0 0 b .
0 0 0
Find the commutators [L, L0 ] and [L, L0 ], L0 ]], where [L, L0 ] := LL0 − L0 L.
Inequalities
(i) Let
1/2 −1/2 1/2 0
A= , B= .
−1/2 1/2 0 1/2
Calculate the left-hand side and the right-hand side of the inequality.
(ii) When is the inequality an equality?
treA+B ≤ tr(eA eB ) .
160
Inequalities 161
(i) Let
2 −1 1 1 4 2
A= , B= , C= .
−1 1 1 2 2 5
Braid Group
Let n ≥ 2. The braid group Bn of n strings has n−1 generators {σ1 , σ2 , . . . , σn−1 }
satisfying the relations
gives a particular braid. The length of the braid is the total number of
used letters, while the minimal irreducible length (referred sometimes as
the primitive word) is the shortest non-contractible length of a particular
braid which remains after applying all the group relations given above.
162
Braid Group 163
1 1
Bj−1 u = c1 e1 + · · · + (cj+1 − δ)ej + (cj − β)ej+1 + · · · + cn en .
γ α
Show that the linear operators Bj satisfy the braid condition
if
γβ + δ = αδ + β .
Problem 7. If V and W are matrices of the same order, then their Schur
product V • W is defined by (entrywise multiplication)
If all entries of V are nonzero, then we say that X is Schur invertible and
define its Schur inverse, V (−) , by V (−) • V = J, where J is the matrix with
all 1’s.
XC M := CM, ∆C M := C • M, YC := M C T .
XA ∆B XA = ∆B XA ∆B .
We call this the braid relation. Give an example for a one-sided Jones pair.
where L, R ∈ SL(2, Z). Thus L−1 , R−1 ∈ SL(2, Z). Show that
RL−1 R = L−1 RL−1 .
1 i 0 0 i 0 0 0
1 i 1 0 0 0 1 0 0
B4 = √ , B5 = .
2 0 0 1 −i 0 0 i 0
0 0 −i 1 0 0 0 1
Are the matrices unitary? Is (braid-like relation)
vec Operator
and
m X
n
∗
X
vec−1
m×n x = (ej,n ⊗ ei,m ) x ei,m ⊗ e∗j,n .
i=1 j=1
Show that
vec−1
m×n (vecm×n A) = A.
167
168 Problems and Solutions
A ⊗ B = vec−1
ms×nt (LA,s×t (vecs×t B))
where
n
X
LA,s×t := (In ⊗ It ⊗ A ⊗ Is ) ej,n ⊗ It ⊗ ej,n ⊗ Is .
j=1
Star Product
169
170 Problems and Solutions
a11 0 0 a12
0 b11 b12 0
A ? B := .
0 b21 b22 0
a21 0 0 a22
A?B =A⊗B?
a11 0 0 a12
0 b11 b12 0
A ? B := .
0 b21 b22 0
a21 0 0 a22
a11 a12 0 0
a21 a22 0 0 T
= P (A ? B)P .
0 0 b11 b12
0 0 b21 b22
a11 0 0 a12
0 b11 b12 0
A ? B := .
0 b21 b22 0
a21 0 0 a22
The 2 × 2 matrices
1 0 0 1
E= , N=
0 1 1 0
a11 0 0 a12
0 b11 b12 0
A ? B := .
0 b21 b22 0
a21 0 0 a22
a11 a12 0 0
a21 a22 0 0
P (A ? B)P T = .
0 0 b11 b12
0 0 b21 b22
b11 0 0 b12
0 a11 a12 0
A ? B := .
0 a21 a22 0
b21 0 0 b22
Is A ? B invertible?
a11 0 0 a12
0 b11 b12 0
A ? B := .
0 b21 b22 0
a21 0 0 a22
172 Problems and Solutions
Problem 12. (i) Given the eigenvalues of A and B. What can be said
about the eigenvalues of A ? B?
(ii) Can one find 4 × 4 permutation matrices P and Q such that
P (A ? B)Q = A ⊕ B ?
Here ⊕ denotes the direct sum
A1 02 02 A2
02 B1 B2 02
A ? B :=
02 B3 B4 02
A3 02 02 A4
Let √
√
1/ 2 0 1/ 2
M = 0√ 1 0√ .
1/ 2 0 −1/ 2
Find the eigenvalues of M and M ♦M .
Nonnormal Matrices
Problem 3. Let a > 0, b ≥ 0 and φ ∈ [0, π]. What are the conditions a,
b, φ such that
0 a
A(a, b, φ) =
eiφ b 0
is a normal matrix?
174
Nonnormal Matrices 175
Miscellaneous
c0 c1 c2 c3 0 0 0 0
c3 −c2 c1 −c0 0 0 0 0
0 0 c0 c1 c2 c3 0 0
0 0 c3 −c2 c1 −c0 0 0
D8 = .
0 0 0 0 c0 c1 c2 c3
0 0 0 0 c3 −c2 c1 −c0
c2 c3 0 0 0 0 c0 c1
c1 −c0 0 0 0 0 c3 −c2
176
Miscellaneous 177
Se = T −1 ST ?
Problem 4. Let
0 0 1 0 1 i 0 0
0 0 0 1 1 0 0 −i 1
J = , P =√ .
−1 0 0 0 2 i 1 0 0
0 −1 0 0 0 0 1 −i
p1 (x) = a0 + a1 x + · · · + an xn , p2 (x) = b0 + b1 x + · · · + bm xm
178 Problems and Solutions
p2 (x)
r(x) = .
p1 (x)
Show that AT = A−1 by showing that the column of the matrix are nor-
malized and pairwise orthonormal.
1
P0 (x) = 1, P1 (x) = x, P2 (x) = (3x2 − 1), . . . .
2
Calculate the infinite dimensional matrix A = (Ajk )
Z +1
dPk (x)
Ajk = Pj (x) dx
−1 dx
Problem 9. Consider the vector space R3 and the vector product ×. The
vector product is not associative. The associator of three vectors u, v, w
is defined by
ass(u × (v × w)) := (u × v) × w − u × (v × w) .
Calculate
d
tr(A + B)j =0
.
d
X = A + S + cIn
Problem 17. Find the partial differential equation given by the condition
0 ∂u/∂x1 ∂u/∂x2
det ∂u/∂x1 ∂ 2 u/∂x21 ∂ 2 u/∂x1 ∂x2 .
2
∂u/∂x2 ∂ u/∂x2 ∂x1 ∂ 2 u/∂x22
Problem 19. Consider vectors in the vector space R3 and the vector
product. Consider the mapping of the vectors into 3 × 3 skew-symmetric
matrices
a 0 c −b
b ↔ −c 0 a .
c b −a 0
Calculate
a1 a2
b1 × b2
c1 c2
and [M1 , M2 ], where
0 c1 −b1 0 c2 −b2
M1 = −c1 0 a1 , M2 = −c2 0 a2 .
b1 −a1 0 b2 −a2 0
Discuss.
Problem 20. Let u1 (t), u2 (t), u3 (t) ∈ R3 . Solve the initial value problem
of the nonlinear autonomous system of first order differential equations
Problem 21. Let u(t) ∈ R3 . Solve the initial value problem for the
differential equation
d2 u du
=u×
dt2 dt
where × denotes the vector product.
0 0 0 t
0 0 t
0 t 0 0 t 0
A1 = 1, A2 = , A3 = 0 1 0, A4 =
t r 0 t r 0
t 0 r
t 0 0 r
182 Problems and Solutions
0 0 0 0 t
0 0 0 t 0
A5 = 0 0 1 0 0 .
0 t 0 r 0
t 0 0 0 r
Thus the even dimensional matrix A2n has t along the skew-diagonal and
r along the lower main diagonal. Otherwise the entries are 0. The odd
dimensional matrix A2n+1 has t along the skew-diagonal except 1 at the
centre and r along the lower main diagonal. Otherwise the entries are 0.
Find the eigenvalues of these matrices.
Problem 25. Let Cn×N be the vector space of all n×N complex matrices.
Let Z ∈ Cn×N . Then Z ∗ ≡ Z̄ T , where T denotes transpose. One defines a
Gaussian measure µ on Cn×N by
1
dµ(Z) := exp(−tr(ZZ ∗ ))dZ
π nN
where dZ denotes the Lebesgue measure on Cn×N . The Fock space F(Cn×N )
consists of all entire functions on Cn×N which are square integrable with
respect to the Gaussian measure dµ(Z). With the scalar product
Z
hf |gi := f (Z)g(Z)dµ(Z), f, g ∈ F(Cn×N )
Cn×N
one has a Hilbert space. Show that this Hilbert space has a reproducing
kernel K. This means a continuous function K(Z, Z 0 ) : Cn×N × Cn×N → C
such that Z
f (Z) = K(Z, Z 0 )f (Z 0 )dµ(Z 0 )
Cn×N
Miscellaneous 183
Problem 26. The vector space of all n×n matrices over C form a Hilbert
space with the scalar product defined by
hA, Bi := tr(AB ∗ ) .
such a pair?
form a group under matrix multiplication. If not add the elements to find
a group. Here 02 is the 2 × 2 zero matrix.
Let √
√
1/ 2 0 1/ 2
M = 0√ 1 0√ .
1/ 2 0 −1/ 2
Find the eigenvalues of M and M ♦M .
(i) Show that they form an orthonormal basis in the Hilbert space of the
2 × 2 matrices with the scalar product hX, Y i = tr(XY ∗ ).
(ii) Find the multiplication table.
S 0n
Miscellaneous 185
C T C + CC T = I2 , C 2 = 02 .
and
in this form.
Consider the map for the special cases θ = 0, φ = 0 and θ = π/4, φ = π/4.
mink[A, B] − I2 k
where [ , ] denotes the commutator and for the norm k k consider the Frobe-
nius norm and max-norm.
lim dj = ∞ .
j→∞
p(x) = x2 − sx + d, s, d ∈ C .
A = S −1 BS .
Are the two matrices normal? Are the matrices similar? Are the matrices
X ⊗ Y and Y ⊗ X similar?
Show that R−1 = R∗ = R. Use these properties and tr(R) to find all the
eigenvalues of R. Find the eigenvectors.
(ii) Let
0 1 0 0 0 1
A1 = 1 0 0 , A2 = 0 0 1 .
0 0 1 1 1 0
Calculate RA1 R−1 and RA2 R−1 . Discuss.
Miscellaneous 189
Problem 51. (i) Find the conditions on the 2 × 2 matrices over C such
that
ABA = BAB .
Find solutions where AB 6= BA, i.e. [A, B] 6= 02 .
(ii) Find the conditions on the 2 × 2 matrices A and B such that
A⊗B⊗A=B⊗A⊗B.
Find solutions where AB 6= BA, i.e. [A, B] 6= 02 .
where α, β, γ are the three Euler angles with the range 0 ≤ α < 2π,
0 ≤ β ≤ π and 0 ≤ γ < 2π. Show that
−iα/2
cos(β/2)e−iγ/2 −e−iα/2 sin(β/2)eiγ/2
e
U (α, β, γ) = . (1)
e−iα/2 sin(β/2)e−iγ/2 eiα/2 cos(β/2)eiγ/2
hA, Bi := tr(AB ∗ ), A, B ∈ H .
Find the commutator [A(t), B(t)]. Discuss. What is the condition such
that [A(t), B(t)] = 02 .
projection matrices?
{ Dj P k : j, k = 0, 1, 2, . . . , n − 1 }
P D = ωDP, P j Dk = ω jk Dk P j .
X = ζP + ζ −1 P −1 + ηD + η −1 D−1
tr(AB) = 0 .
det(AB) = 0 .
Problem 69. How many 3×3 binary matrices can one form which contain
three 1’s? Write down these matrices. Which of them are invertible?
Find Rs Vs Rs .
(iii) Let s be 1/2, 3/2, . . . with n = 2s + 1 and the n × n matrix
1 1
.. .
. ..
1
1 1
Rs = √ .
2 1 −1
. . ..
. .
1 −1
Find Rs Vs Rs .
(i) Find the condition on λj , aj and µk , bk such that tr(AB ∗ ) = 0, i.e. the
two matrices are orthogonal to each other.
(ii) Find the condition on λj , aj and µk , bk such that [A, B] = 0n , i.e. the
the commutator of the matrices vanishes.
ezA Be−zA .
Show that
d(ln(det(A))) = tr(A−1 dA)
where d denotes the exterior derivative.
f ◦ π = π ◦ A, π : Rn → Tn := Rn /Zn .
#Fix(f ) := { x∗ ∈ Tn : f (x∗ ) = x∗ }
Let n = 2 and
2 1
A= .
1 1
Show that det(I2 − A) 6= 0 and find #Fix(f ).
Ĥ = ~ωσ 1 · σ 2 ≡ ~ω(σx ⊗ σx + σy ⊗ σy + σz ⊗ σz ) .
Show that this function also admits the fixed point (0, 0). Find the func-
tional matrix at this fixed point
∂h1 /∂x1 ∂h1 /∂x2
.
∂h2 /∂x1 ∂h2 /∂x2 (0,0)
An+1 = An Bn , Bn+1 = An .
S 2 = In , (T S)2 = In .
Show that
ST S −1 = T −1 , ST −1 S = T .
U 7→ U V W, V 7→ U V, W 7→ U .
Calculate dV (t)/dt and then find the commutator [dV (t)/dt, V (t)].
(ii) Let V (t) be a 2 × 2 matrix where all the entries are smooth functions
of t. Calculate dV (t)/dt and then find the conditions on the entries such
that [dV (t)/dt, V (t)] = 02 .
Problem 94. (i) Find all invertible 2 × 2 matrices over R such that
0 1 0 0
T T −1 = .
0 0 1 0
Armstrong M. A.
Groups and Symmetry
Springer (1988)
Bredon G. E.
Introduction to Compact Transformation Groups
Academic Press (1972)
Bronson R.
Matrix Operations
Schaum’s Outlines, McGraw-Hill (1989)
Bump D.
Lie Groups
Springer (2000)
Campoamor-Stursberg R.
“The structure of the invariants of perfect Lie algebras”, J. Phys. A: Math.
Gen. 6709–6723 (2003)
Carter R. W.
Simple Groups of Lie Type
John Wiley (1972)
201
202 Bibliography
DasGupta Ananda
American J. Phys. 64 1422–1427 (1996)
Dixmier J.
Enveloping Algebras
North-Holland (1974)
Englefield M. J.
Group Theory and the Coulomb Problem
Wiley-Interscience, New York (1972)
Fuhrmann, P. A.
A Polynomial Approach to Linear Algebra
Springer (1996)
Gallian J. A.
Contemporary Abstract Algebra, Sixth edition
Houghton Mifflin (2006)
Harville D. A.
Matrix Algebra from a Statistician’s Perspective
Springer (1997)
Grossman S. I.
Elementary Linear Algebra, Third Edition
Wadsworth Publishing, Belmont (1987)
Hall B. C.
Lie Groups, Lie Algebras, and Representations: an elementary introduction
Springer (2003)
Helgason S.
Groups and Geometric Analysis, Integral Geometry, Invariant Differential
Operators and Spherical Functions
Academic Press (1984)
Humphreys J. E.
Introduction to Lie Algebras and Representation Theory
Springer (1972)
Isham C. J.
Modern Diffferential Geometry
World Scientific (1989)
Johnson D. L.
Presentation of Groups
Cambridge University Press (1976)
Jones H. F.
Groups, Representations and Physics
Adam Hilger, Bristol (1990)
Lang S.
Linear Algebra
Addison-Wesley, Reading (1968)
Miller W.
Symmetry Groups and Their Applications
Academic Press, New York (1972)
Ohtsuki T.
Quantum Invariants
World Scientific, Singapore (2002)
Olver P. J.
Applications of Lie Groups to Differential Equations, second edition, Grad-
uate Text in Mathematics, vol. 107, Springer, New York (1993)
Steeb W.-H.
Matrix Calculus and Kronecker Product with Applications and C++ Pro-
grams
World Scientific Publishing, Singapore (1997)
Steeb W.-H.
Continuous Symmetries, Lie Algebras, Differential Equations and Com-
puter Algebra
World Scientific Publishing, Singapore (1996)
Steeb W.-H.
Hilbert Spaces, Wavelets, Generalized Functions and Quantum Mechanics
Kluwer Academic Publishers, Dordrecht (1998)
Steeb W.-H.
Problems and Solutions in Theoretical and Mathematical Physics,
Second Edition, Volume I: Introductory Level
Bibliography 205
Steeb W.-H.
Problems and Solutions in Theoretical and Mathematical Physics,
Second Edition, Volume II: Advanced Level
World Scientific Publishing, Singapore (2003)
Sternberg S.
Lie Algebras 2004
Varadarajan V. S.
Lie Groups, Lie Algebras and Their Representations
Springer (2004)
Watkins D. S.
Fundamentals of matrix computations
Wiley
Wawrzynczyk A.
Group Representations and Special Functions
D. Reidel (1984)
Whittaker E. T.
A treatise on the analytical dynamics of particles and rigid bodies
Wybourne B. G.
Classical Groups for Physicists
John Wiley, New York (1974)
Index
206
Index 207
Quotient space, 20
Ternary commutator, 81
Ternutator, 81
Tetrahedron, 14, 24, 40
Three-body problem, 20
Trace norm, 30
Triangle, 24