11.7  Isometries and Unitary Maps

Chapter 2 asked which linear maps are invertible. Now that we can measure, we can ask a sharper question: which linear maps leave the measurements alone? A map that never changes a length never changes an angle either, and the two facts turn out to be the same fact. This section identifies these maps in five equivalent ways, collects them into the orthogonal and unitary groups, and classifies them completely in dimensions two and three, where they are exactly the rotations and reflections that the words suggest. At the end we drop linearity from the hypothesis and get it back as a conclusion.

As throughout the chapter, F is R or C, inner products are linear in the first slot, and spaces are finite-dimensional unless stated otherwise.

11.7.1. Maps that do not change length

Take the plane with its usual inner product and ask which linear maps move the picture rigidly. Turning the plane about the origin clearly does not change any length. Flipping it across a line through the origin does not either. Scaling by 2 obviously does. Shearing, (x1,x2)↦(x1+x2,x2), also does: it fixes e1 but sends e2 to (1,1), of length 2. So “does not change length” is a genuine restriction, and it is worth a name.

An isometry is a linear map that every ruler agrees with: the image of a vector is exactly as long as the vector.

Definition 11.7.1 (Isometry)

Let V and W be inner product spaces over F. A linear map T∈L(V,W) is an isometry if ‖Tv‖=‖v‖for every v∈V.

In words: the condition is required for every vector, not merely for the vectors of a basis, and it compares a norm in W with a norm in V. The word comes from iso (equal) and metron (measure). Some authors reserve “isometry” for a bijective such map; we do not, and Definition 11.7.2 below names the bijective case separately.

Two consequences are immediate. An isometry is injective, since Tv=0 forces ‖v‖=‖Tv‖=0 and hence v=0 (Theorem 2.2.3). And if Tv=λv with v≠0, then ‖v‖=‖λv‖=|λ|‖v‖ by Theorem 11.1.2 (b), so |λ|=1. We record the second observation, which the classification below will use twice.

Proposition 11.7.1 (Eigenvalues of an Isometry Have Modulus 1)

Let T∈L(V) be an isometry of an inner product space over F, and let λ∈F be an eigenvalue of T. Then |λ|=1. In particular, if F=R then λ=1 or λ=−1.

Proof. Let v≠0 satisfy Tv=λv (Definition 9.2.1). Then ‖v‖=‖Tv‖=‖λv‖=|λ|‖v‖ by Theorem 11.1.2 (b), and ‖v‖≠0, so dividing gives |λ|=1. The only real numbers of modulus 1 are ±1.

Example (Isometries and Near-Isometries)

Decide which of the following linear maps are isometries, with the standard inner products.

  • (a) T:R2→R3, T(x1,x2)=(x1,x2,0).
  • (b) TP:Fn→Fn for a permutation matrix P (Definition 3.4.2).
  • (c) T:C→C, T(z)=αz, where |α|=1.
  • (d) T:R2→R2, T(x1,x2)=(x1+x2,x2).
  • (e) The orthogonal projection PU onto a subspace U with {0}≠U≠V.

Solution.

  • (a) Yes: ‖Tx‖2=x12+x22+0=‖x‖2. Note that T is not surjective, so an isometry need not be invertible.

  • (b) Yes. A permutation matrix has exactly one 1 in each row and column, so Px is the list x1,…,xn in a different order, and ‖Px‖2=∑i|xσ(i)|2=∑i|xi|2, the same sum reordered.

  • (c) Yes: ‖αz‖=|αz|=|α||z|=|z|.

  • (d) No: ‖Te2‖=‖(1,1)‖=2≠1=‖e2‖. One vector is enough to refute.

  • (e) No: pick v∈U⊥ with v≠0, which exists because U≠V (Theorem 11.3.2 (a)). Then PUv=0 (Theorem 11.3.3 (a)), so ‖PUv‖=0≠‖v‖. A projection destroys length by design.

11.7.2. Five ways to say the same thing

The definition mentions only norms, but norms and inner products determine each other. In one direction this is Definition 11.1.3. In the other it is the following identity, which reconstructs the inner product from lengths alone; it is the reason that preserving length automatically preserves angle.

Lemma 11.7.2 (Polarization Identities)

Let V be an inner product space over F and let u,v∈V.

  • (a) If F=R, then ⟨u,v⟩=14(‖u+v‖2−‖u−v‖2).
  • (b) If F=C, then ⟨u,v⟩=14∑k=03ik‖u+ikv‖2.

Idea. Expand every norm squared with ‖u+w‖2=‖u‖2+2Re⟨u,w⟩+‖w‖2 (Theorem 11.1.2 (c)). In the real case the two expansions differ only in the middle term, which survives doubled. In the complex case the four terms ‖u‖2+‖w‖2 are the same in all four summands and are killed by ∑kik=0; what is left is a sum in which the ⟨u,v⟩ parts reinforce and the ⟨u,v⟩― parts cancel.

Proof.

  • (a) By Theorem 11.1.2 (c), ‖u±v‖2=‖u‖2±2⟨u,v⟩+‖v‖2, since over R the inner product is its own real part. Subtracting, ‖u+v‖2−‖u−v‖2=4⟨u,v⟩.

  • (b) Write z=⟨u,v⟩. Since |ik|=1, Theorem 11.1.2 (b) and (c) give ‖ikv‖=‖v‖ and ‖u+ikv‖2=‖u‖2+‖v‖2+2Re⟨u,ikv⟩=‖u‖2+‖v‖2+2Re(i−kz), where ⟨u,ikv⟩=ik―z=i−kz by conjugate-linearity in the second slot. Multiply by ik and sum over k=0,1,2,3. The first two terms contribute (‖u‖2+‖v‖2)∑kik=0, because 1+i−1−i=0. For the remaining terms use 2Re(w)=w+w―: ∑k=03ik⋅2Re(i−kz)=∑k=03ik(i−kz+ikz―)=∑k=03z+z―∑k=03i2k=4z+0, since ∑ki2k=1−1+1−1=0. Dividing by 4 gives the identity.

Now the main theorem. It says that length, angle, the adjoint, and orthonormal bases all detect the same class of maps, and it is the reason isometries are easy to work with: whichever of the five descriptions a problem hands you, you may switch to whichever one you want.

Theorem 11.7.3 (Characterizations of Isometries)

Let V and W be finite-dimensional inner product spaces over F, and let T∈L(V,W). The following are equivalent.

  • (a) T is an isometry: ‖Tv‖=‖v‖ for every v∈V.
  • (b) T preserves inner products: ⟨Tu,Tv⟩=⟨u,v⟩ for all u,v∈V.
  • (c) T∗T=idV.
  • (d) T sends every orthonormal list in V to an orthonormal list in W.
  • (e) T sends some orthonormal basis of V to an orthonormal list in W.
  • (f) For one (equivalently every) pair of orthonormal bases B of V and C of W, the matrix A=[T]BC∈Mm×n(F), where n=dim⁡V and m=dim⁡W, satisfies A∗A=In; that is, the n columns of A are orthonormal vectors of Fm.

If in addition dim⁡W=dim⁡V, then (a)–(f) are further equivalent to each of

  • (g1) T is invertible with T−1=T∗;
  • (g2) TT∗=idW.

Idea. Run the cycle (a) ⇒ (b) ⇒ (c) ⇒ (a). The first step is polarization: lengths are all we are given, and Lemma 11.7.2 converts them into an inner product. The second is the adjoint doing its job, ⟨Tu,Tv⟩=⟨u,T∗Tv⟩, plus nondegeneracy. The third is the same identity read backwards with u=v. Then (b) ⇒ (d) ⇒ (e) ⇒ (b) is a second, shorter cycle, and (c) ⇔ (f) is Theorem 11.6.3. The last two need a dimension count: T∗T=id makes T injective, and in equal dimensions injective means invertible.

Proof.

  • (a) ⇒ (b). Suppose ‖Tv‖=‖v‖ for all v. If F=R, then by Lemma 11.7.2 (a) applied in W and then in V, and by the linearity of T, ⟨Tu,Tv⟩=14(‖T(u+v)‖2−‖T(u−v)‖2)=14(‖u+v‖2−‖u−v‖2)=⟨u,v⟩. If F=C, the same two lines work with Lemma 11.7.2 (b), since T(u+ikv)=Tu+ikTv.

  • (b) ⇒ (c). For all u,v∈V, Definition 11.6.1 gives ⟨u,T∗Tv⟩=⟨Tu,Tv⟩=⟨u,v⟩, so ⟨u,T∗Tv−v⟩=0 for every u. Taking u=T∗Tv−v makes its norm zero, so T∗Tv=v for every v.

  • (c) ⇒ (a). ‖Tv‖2=⟨Tv,Tv⟩=⟨v,T∗Tv⟩=⟨v,v⟩=‖v‖2, and norms are nonnegative, so ‖Tv‖=‖v‖.

  • (d) ⇒ (d). Let (u1,…,uk) be orthonormal in V, so ⟨ui,uj⟩=δij. Then ⟨Tui,Tuj⟩=⟨ui,uj⟩=δij, which is exactly what it means for (Tu1,…,Tuk) to be orthonormal (Definition 11.2.1).

  • (e) ⇒ (e). V is finite-dimensional, so it has an orthonormal basis by Theorem 11.2.2, and (d) applies to it.

  • (f) ⇒ (b). Let B=(v1,…,vn) be an orthonormal basis of V with (Tv1,…,Tvn) orthonormal. Given u,v∈V, expand them in B as u=∑iaivi and v=∑jbjvj. Expanding both slots and using orthonormality twice, ⟨Tu,Tv⟩=∑i,jaibj―⟨Tvi,Tvj⟩=∑iaibi―=∑i,jaibj―⟨vi,vj⟩=⟨u,v⟩, since in both double sums only the terms with i=j survive.

  • (g) ⇔ (f). Fix orthonormal bases B and C and put A=[T]BC. By Theorem 11.6.3, [T∗]CB=A∗, and by Theorem 4.1.2, [T∗T]BB=A∗A. Since a map is determined by its matrix in fixed bases, T∗T=idV if and only if A∗A=In. The (i,j) entry of A∗A is the standard inner product of column j with column i, so this says the columns are orthonormal. As (c) does not mention the bases, the condition holds for one pair of orthonormal bases if and only if it holds for every pair.

Now assume dim⁡W=dim⁡V.

  • (c) ⇒ (g1). By (c) and Theorem 2.2.3, T is injective, so rankT=dim⁡V=dim⁡W by Theorem 2.3.1, and T is surjective. Hence T is invertible, and multiplying T∗T=idV on the right by T−1 gives T∗=T−1.

(g1) ⇒ (g2). TT∗=TT−1=idW.

(g2) ⇒ (c). TT∗=idW makes T surjective, hence invertible by the equal dimensions and Theorem 2.3.1, and multiplying on the left by T−1 gives T∗=T−1, so T∗T=idV. This proves the theorem.

Condition (f) is the one to use on a concrete matrix: check that the columns are orthonormal. Condition (c) is the one to use in a proof. Condition (e) is the cheapest to verify on an abstract map, since one basis suffices, and (d) is what you get out afterwards.

Warning. T∗T=id does not give TT∗=id when the dimensions differ. Let T:R2→R3 be T(x1,x2)=(x1,x2,0), with matrix A=(100100). Then A⊤A=I2, so T is an isometry; but AA⊤=(100010000)≠I3. Ask identity on what? The composite TT∗ is the orthogonal projection of R3 onto imT, the x1x2-plane, and it is the identity there and nowhere else. Only when T is surjective can both products be identities.

Quick check. Is (111−1) the matrix of an isometry of R2? If not, is some scalar multiple of it?

Solution. No: the first column has norm 2≠1, so the columns are not orthonormal and condition (f) of Theorem 11.7.3 fails. The columns are, however, orthogonal, so dividing by 2 normalizes both and 12(111−1) is an isometry.

11.7.3. Unitary and orthogonal maps

Conditions (g1) and (g2) describe the case that matters most, an isometry of a space with itself. It is the case with a name.

Definition 11.7.2 (Unitary and Orthogonal)

Let V be a finite-dimensional inner product space over F. An operator T∈L(V) is unitary if T∗T=TT∗=idV, equivalently if T is an invertible isometry. When F=R, such a T is also called orthogonal.

A matrix A∈Mn(C) with A∗A=In is a unitary matrix; a matrix A∈Mn(R) with A⊤A=In is an orthogonal matrix.

For a square matrix the one-sided condition A∗A=I already forces AA∗=I, by part (g2) of Theorem 11.7.3, or directly by Theorem 3.4.5. The warning above shows that squareness is not a formality.

Definition 11.7.3 (The Unitary and Orthogonal Groups)

For n≥1 set U(n):={A∈Mn(C):A∗A=In},O(n):={A∈Mn(R):A⊤A=In}, and SO(n):={A∈O(n):detA=1}, the special orthogonal group.

Proposition 11.7.4 (The Groups and Their Determinants)
  • (a) U(n) and O(n) are groups under matrix multiplication, and O(n)=U(n)∩Mn(R).
  • (b) If A∈U(n) then |detA|=1. If A∈O(n) then detA=1 or detA=−1.
  • (c) SO(n) is a subgroup of O(n).

Proof.

  • (a) I∗I=I, so I∈U(n). If A,B∈U(n) then (AB)∗(AB)=B∗A∗AB=B∗B=I, so AB∈U(n). Each A∈U(n) is invertible with A−1=A∗, and (A∗)∗A∗=AA∗=I, so A−1∈U(n). Associativity is inherited from matrix multiplication. The same three lines prove the claim for O(n), with A∗=A⊤ for a real matrix, which is also the last assertion.

  • (b) The determinant is the sum ∑σsgn(σ)aσ(1)1⋯aσ(n)n of Definition 7.4.1, and conjugation of complex numbers respects sums and products, so detA―=detA―. With Theorem 7.4.5 this gives detA∗=det(A―⊤)=detA―=detA―. Hence, by Theorem 7.5.4, 1=detIn=det(A∗A)=detA―detA=|detA|2. For real A the determinant is real, and the real numbers of modulus 1 are ±1.

  • (c) detI=1, and det is multiplicative (Theorem 7.5.4), so det(AB)=1 and det(A−1)=(detA)−1=1 whenever detA=detB=1.

Warning. |λ|=1 does not mean λ=±1, and detA=1 does not mean “rotation”. The matrix −I4∈SO(4) has determinant (−1)4=1, but it fixes no nonzero vector, so it is not a rotation about an axis in any sense. Likewise −I3∈O(3) has determinant −1 but is not a reflection in a plane: a reflection fixes a plane pointwise, and −I3 fixes only 0. The clean classification below is a feature of dimensions 2 and 3, not a general theorem.

11.7.4. The Frobenius norm and unitary factors

Section 1 made Mm×n(F) into an inner product space by declaring ⟨A,B⟩=tr(B∗A), the Frobenius inner product, and observing that it is nothing but the standard inner product of Fmn after the entries are listed in a single column. Its induced norm (Definition 11.1.3) was named there but never given a symbol. This is the place to give it one, because the single property that makes it useful is a statement about unitary factors, and unitary factors have only just arrived.

For A∈Mm×n(C), the Frobenius norm of A is ‖A‖F:=⟨A,A⟩=(∑i,j|aij|2)1/2. The subscript F is there because ‖⋅‖ with no subscript already means the norm of a vector, and a matrix will later be asked to be both. Two remarks on the formula. First, it treats A as a list of mn numbers and forgets that they are arranged in a rectangle: it is blind to the difference between (1000) and (0100). Second, it is computable with no eigenvalues, no row reduction, and no choices: add up mn squared moduli.

The trace form of the same quantity is the one that does the work.

Quick check. Why is ‖A‖F2=tr(A∗A)?

Solution. The (j,j) entry of A∗A is ∑iaij―aij=∑i|aij|2, the squared length of the j-th column. Summing over j sweeps every entry of A exactly once, so tr(A∗A)=∑i,j|aij|2. Equivalently, it is ⟨A,A⟩ read off the definition of the Frobenius inner product.

Now the property that makes this the right yardstick whenever orthonormal bases are in play: unitary factors are invisible to it.

Lemma 11.7.5 (Unitary Invariance of the Frobenius Norm)

Let A∈Mm×n(C), let U∈U(m) and let V∈U(n). Then ‖UA‖F=‖A‖F=‖AV‖F. In particular, if A is square and U is unitary then ‖U∗AU‖F=‖A‖F.

Proof. Using ‖B‖F2=tr(B∗B), then (UA)∗=A∗U∗ and U∗U=Im, ‖UA‖F2=tr(A∗U∗UA)=tr(A∗A)=‖A‖F2. For the right-hand factor, Theorem 0.9.6 (3) lets the V travel around the trace, and VV∗=In: ‖AV‖F2=tr(V∗A∗AV)=tr(A∗AVV∗)=‖A‖F2. Both norms are nonnegative reals, so equality of the squares gives equality of the norms. For the last statement, U∗ is unitary too, since (U∗)∗U∗=UU∗=I; apply the first equality to the left factor U∗ and the second to the right factor U.

The consequence to keep in mind is this. If T∈L(V) is an operator on a finite-dimensional complex inner product space and we compute ‖[T]BB‖F in an orthonormal basis B, the answer does not depend on which orthonormal basis we chose, because two such matrices differ by a unitary similarity (Theorem 4.2.3). So the Frobenius norm is a property of T, not of our coordinates — exactly like the trace and the determinant, and unlike, say, the largest entry.

11.7.5. Rotations and reflections of the plane

In dimension 2 the group O(2) can be written down completely, and the two determinant values split it into exactly the two families the words suggest. For θ∈R put Rθ:=(cos⁡θ−sin⁡θsin⁡θcos⁡θ),Mθ:=(cos⁡θsin⁡θsin⁡θ−cos⁡θ).

Diagram. Labels: cvec e_1, cvec e_2, R_theta cvec e_1, R_theta cvec e_2, the columns of R_theta are the rotated axes

The first matrix rotates the plane by θ: its columns are e1 and e2 turned through θ. The second is a reflection, and the picture says in which mirror.

Diagram. Labels: L, cvec v, M_theta cvec v, P_Lcvec v, theta /2, M_theta reflects in the line L at angle theta /2
Theorem 11.7.6 (Classification of the Plane Isometries)

Let A∈O(2).

  • (a) If detA=1, then A=Rθ for a unique θ∈[0,2π).
  • (b) If detA=−1, then A=Mθ for a unique θ∈[0,2π). In this case A2=I, the eigenvalues of A are 1 and −1, and A is the reflection in the line L=span((cos⁡θ2,sin⁡θ2)): it fixes L pointwise and negates L⊥.

Consequently O(2)={Rθ}∪{Mθ}, and SO(2)={Rθ:θ∈[0,2π)}.

Idea. Condition (f) of Theorem 11.7.3 says the two columns are orthonormal. A unit vector of R2 is (cos⁡θ,sin⁡θ) for exactly one θ∈[0,2π), and there are exactly two unit vectors orthogonal to it. So there are only two candidates once the first column is chosen, and the determinant tells them apart. The geometric description of Mθ is then a two-line verification: guess the fixed vector from the picture and check it.

Proof. Write A=(acbd)∈O(2). By Theorem 11.7.3 (f), the columns are orthonormal: a2+b2=1, c2+d2=1, and ac+bd=0. Since (a,b) is a unit vector of R2, there is a unique θ∈[0,2π) with a=cos⁡θ and b=sin⁡θ. The vectors orthogonal to (a,b) form a line, spanned by (−b,a), which is already a unit vector; so the unit vectors on that line are ±(−b,a), and (c,d)=±(−sin⁡θ,cos⁡θ). The two choices give (cos⁡θ−sin⁡θsin⁡θcos⁡θ)=Rθand(cos⁡θsin⁡θsin⁡θ−cos⁡θ)=Mθ, with determinants cos2⁡θ+sin2⁡θ=1 and −cos2⁡θ−sin2⁡θ=−1. This proves (a) and the first sentence of (b), and uniqueness of θ in both cases, since θ is already determined by the first column.

For the rest of (b), put u=(cos⁡θ2,sin⁡θ2) and w=(−sin⁡θ2,cos⁡θ2), two orthogonal unit vectors. Writing θ=θ2+θ2 and using the subtraction formulas for cosine and sine, Mθu=(cos⁡θcos⁡θ2+sin⁡θsin⁡θ2sin⁡θcos⁡θ2−cos⁡θsin⁡θ2)=(cos⁡θ2sin⁡θ2)=u, Mθw=(−cos⁡θsin⁡θ2+sin⁡θcos⁡θ2−sin⁡θsin⁡θ2−cos⁡θcos⁡θ2)=(sin⁡θ2−cos⁡θ2)=−w. So (u,w) is an orthonormal basis of eigenvectors with eigenvalues 1 and −1, and Mθ2 fixes both of them, hence Mθ2=I. Since L=span(u), we have L⊥=span(w), because w⊥u and dim⁡L⊥=1 (Theorem 11.3.2 (c)), so Mθ fixes L pointwise and negates L⊥, which is the description of a reflection. This proves the theorem.

So SO(2) is the circle of rotations, and the other half of O(2) is the circle of reflections. Multiplying two rotations adds the angles, RαRβ=Rα+β, so SO(2) is abelian; Exercise (C3: Two reflections make a rotation) shows what happens when two reflections are multiplied.

Example (Naming Two Plane Isometries)

Classify each of the following as a rotation or a reflection, and give the angle or the mirror line.

  • (a) A=12(1−331).
  • (b) B=15(344−3).

Solution.

  • (a) The columns 12(1,3) and 12(−3,1) have norm 121+3=1 and inner product 14(−3+3)=0, so A∈O(2). Its determinant is 14(1+3)=1, so A=Rθ with cos⁡θ=12 and sin⁡θ=32: rotation by θ=π/3.

  • (b) The columns 15(3,4) and 15(4,−3) are orthonormal, and detB=125(−9−16)=−1, so B=Mθ with cos⁡θ=35 and sin⁡θ=45. The mirror is L=span((cos⁡θ2,sin⁡θ2)), and rather than halve the angle we find the fixed line directly: Bx=x reads 3x1+4x2=5x1, that is x1=2x2. So L=span((2,1)). Check: B(2,1)=15(6+4,8−3)=(2,1), and B(−1,2)=15(−3+8,−4−6)=(1,−2)=−(−1,2), as a reflection must do to L⊥.

11.7.6. Rotations of space

In dimension 3 the answer is the one physics assumes: every orientation-preserving isometry of R3 is a rotation about a line. The proof is where Chapter 9 is needed, for one reason only: to produce an eigenvector.

Theorem 11.7.7 (Every Element of SO(3) Is a Rotation)

Let A∈SO(3). Then there is a unit vector u∈R3 with Au=u, and an orthonormal basis B=(u,w1,w2) of R3 in which [TA]BB=(1000cos⁡θ−sin⁡θ0sin⁡θcos⁡θ) for some θ. That is, A fixes the axis span(u) and rotates the plane span(u)⊥ by θ. Moreover trA=1+2cos⁡θ.

Idea. Three steps. ① Get a real eigenvalue: the characteristic polynomial is real of degree 3, and an odd-degree real polynomial has a real root. ② That eigenvalue is ±1, since TA is an isometry. ③ Split off the eigenline: its orthogonal complement is invariant, the restriction there lies in O(2), and Theorem 11.7.6 finishes. If the eigenvalue was −1 the restriction has determinant −1, so it is a reflection and supplies a +1 eigenvector inside the plane; feed that back into step ③.

Proof. Write T=TA, an isometry of R3 by Theorem 11.7.3.

Step 1: a real eigenvalue exists. By Theorem 7.10.2, pA∈R[x] is monic of degree 3, which is odd, so pA has a real root λ by Corollary 6.5.9, and λ is an eigenvalue of A by Theorem 9.2.1. By Proposition 11.7.1, λ=1 or λ=−1. Fix a unit eigenvector u0 for λ, possible after dividing by its norm.

Step 2: the orthogonal complement of an invariant line is invariant.

Claim. If U⊆R3 is a subspace with T(U)⊆U, then T(U⊥)⊆U⊥.

Proof. T|U:U→U is injective, because T is, so it is surjective by Theorem 2.5.7 applied on the finite-dimensional space U. Let w∈U⊥ and u∈U. Write u=Tu′ with u′∈U. Then ⟨Tw,u⟩=⟨Tw,Tu′⟩=⟨w,u′⟩=0, the middle step by Theorem 11.7.3 (b) and the last because u′∈U. So Tw∈U⊥.

Step 3: reduce to the plane. Put U=span(u0), which is T-invariant. By the claim, P:=U⊥ is T-invariant, and dim⁡P=2 by Theorem 11.3.2 (c). Choose an orthonormal basis (w1,w2) of P (Theorem 11.2.2); then B0=(u0,w1,w2) is an orthonormal basis of R3, and in it [T]B0B0=(λ000c11c120c21c22),C:=(c11c12c21c22), the first column because Tu0=λu0, and the two zeros in the first row because Tw1,Tw2∈P have no u0-component. This matrix is orthogonal, being the matrix of an isometry in orthonormal bases (Theorem 11.7.3 (f)), so its second and third columns are orthonormal, which says C⊤C=I2 and C∈O(2); and C is the matrix of T|P in the basis (w1,w2). Expanding the determinant along the first column, and using that similar matrices have equal determinants (Corollary 7.5.7, with Theorem 4.2.4 (a)), 1=detA=λdetC, so detC=λ.

Case λ=1. Then detC=1, so C=Rθ for some θ by Theorem 11.7.6 (a), and B=B0 and u=u0 are as required.

Case λ=−1. Then detC=−1, so by Theorem 11.7.6 (b) the restriction T|P is a reflection of the plane P and has 1 as an eigenvalue. Let u∈P be a unit eigenvector for 1, so Tu=u. Now rerun Steps 2 and 3 with u in place of u0. Because u was chosen with Tu=u, the eigenvalue this time is 1, so the rerun lands in the previous case and cannot return here.

Finally, A and [T]BB are matrices of the same operator, hence similar (Theorem 4.2.4 (a)), and similar matrices have equal traces (Theorem 4.2.5), so trA=1+2cos⁡θ. This proves the theorem.

The angle is determined by A up to sign, since reversing the order of w1,w2 replaces θ by −θ; the trace formula cos⁡θ=12(trA−1) gives |θ| with no computation at all. This settles O(3) as well: if detA=−1 then −A∈SO(3), so A is a rotation composed with the antipodal map.

Example (Axis and Angle of a Rotation)

Let A=13(2−21122−2−12). Verify that A∈SO(3) and find its axis and angle.

Solution. Orthogonality. The three columns are 13(2,1,−2), 13(−2,2,−1) and 13(1,2,2). Their squared norms are 19(4+1+4)=1, 19(4+4+1)=1 and 19(1+4+4)=1, and their pairwise inner products are 19(−4+2+2)=0, 19(2+2−4)=0 and 19(−2+4−2)=0. So the columns are orthonormal and A∈O(3) by Theorem 11.7.3 (f).

Determinant. Expanding along the first row, 27detA=2(4+2)−(−2)(2+4)+1(−1+4)=12+12+3=27, so detA=1 and A∈SO(3).

Axis. Solve Au=u, that is (3I−3A)u=0: (12−1−11−2211)u=0. Adding row one to row two gives 3u2−3u3=0, so u2=u3; then row one gives u1=−u2. The axis is span((−1,1,1)), and indeed A(−1,1,1)=13(−2−2+1,−1+2+2,2−1+2)=(−1,1,1).

Angle. trA=13(2+2+2)=2, so cos⁡θ=12(2−1)=12 and |θ|=π/3. The map is a rotation by 60∘ about the line through (−1,1,1).

Quick check. A matrix A∈SO(3) has trA=−1. What is its angle, and what is A2?

Solution. cos⁡θ=12(−1−1)=−1, so θ=π: a half-turn about the axis. Then A2 is the rotation by 2π about the same axis, which is I3. (In the basis of Theorem 11.7.7 the matrix is diag(1,−1,−1), whose square is I3.)

11.7.7. Rigid motions

We assumed linearity from the start. It is worth knowing that we did not have to: for maps of Rn that fix the origin, preserving distance forces linearity.

Theorem 11.7.8 (Distance-Preserving Maps Are Affine)

Let f:Rn→Rn be any function, not assumed linear, with ‖f(x)−f(y)‖=‖x−y‖for all x,y∈Rn.

  • (a) If f(0)=0, then f is linear and orthogonal.
  • (b) In general, f(x)=Ax+b for a unique A∈O(n) and a unique b∈Rn.

Idea. Setting y=0 turns the hypothesis into “preserves norms”, and expanding the squared distance then turns it into “preserves inner products” without any polarization at all. Once inner products are preserved, f sends the standard basis to an orthonormal basis, and Theorem 11.2.1 writes f(x) as an explicit linear expression in the coordinates of x. Part (b) is part (a) applied to x↦f(x)−f(0).

Proof.

  • (a) Assume f(0)=0. Taking y=0 gives ‖f(x)‖=‖x‖ for all x. Expanding both sides of the hypothesis with Theorem 11.1.2 (c), ‖f(x)‖2−2⟨f(x),f(y)⟩+‖f(y)‖2=‖x‖2−2⟨x,y⟩+‖y‖2, and the four outer terms cancel in pairs, leaving (∗)⟨f(x),f(y)⟩=⟨x,y⟩for all x,y. Put fi:=f(ei). By (∗), ⟨fi,fj⟩=⟨ei,ej⟩=δij, so (f1,…,fn) is an orthonormal list of n vectors in Rn, hence an orthonormal basis (it is independent by Theorem 11.1.4 and has n=dim⁡Rn elements). Now let x=(x1,…,xn) be arbitrary. By Theorem 11.2.1 applied to the basis (f1,…,fn), and then by (∗), f(x)=∑i=1n⟨f(x),fi⟩fi=∑i=1n⟨f(x),f(ei)⟩fi=∑i=1n⟨x,ei⟩fi=∑i=1nxifi. The right-hand side is linear in x, so f is linear; and it preserves norms, so it is an isometry, hence orthogonal by Definition 11.7.2, since it maps Rn to Rn.

  • (b) Put b:=f(0) and g(x):=f(x)−b. Then g preserves distances, since the translation cancels, and g(0)=0; by (a), g is orthogonal, say g(x)=Ax with A∈O(n), and f(x)=Ax+b. For uniqueness, evaluating at 0 recovers b, and then Ax=f(x)−b determines A. This proves the theorem.

So the rigid motions of Rn are exactly the maps “rotate or reflect, then translate”, and the linear ones are exactly the orthogonal maps. That is the sense in which O(n) is the symmetry group of Euclidean geometry, and it is why the same group reappears whenever a problem has a distance in it.

11.7.8. Exercises

A. Check your understanding

Exercise (A1)
  • (a) Define an isometry T∈L(V,W), and state two conditions equivalent to it.
  • (b) What must be checked about a square matrix A to decide whether it is unitary? State the check in terms of its columns.
  • (c) True or false: every isometry is invertible. Justify your answer.
  • (d) True or false: if detA=1 for A∈O(n), then A fixes a nonzero vector. Justify your answer.
  • (e) Give the two families of matrices making up O(2), and say how the determinant distinguishes them.
  • (f) An operator T on a complex inner product space is unitary. Why must every eigenvalue of T lie on the unit circle?

Solution.

  • (a) ‖Tv‖=‖v‖ for every v∈V (Definition 11.7.1). Equivalent: ⟨Tu,Tv⟩=⟨u,v⟩ for all u,v; and T∗T=idV (Theorem 11.7.3).

  • (b) That its columns are orthonormal in Fn, that is A∗A=In (Theorem 11.7.3 (f)). For a square matrix this already gives AA∗=In.

  • (c) False. T:R2→R3, T(x1,x2)=(x1,x2,0), is an isometry and is not surjective (Example (Isometries and Near-Isometries) (a)). An isometry is always injective, and it is invertible exactly when dim⁡W=dim⁡V.

  • (d) False for even n: −I4∈SO(4) has determinant 1 and fixes only 0. It is true for odd n, though not by the argument of Theorem 11.7.7, which uses the O(2) classification of a 2-dimensional complement. The general reason is a count: the nonreal eigenvalues of a real A come in conjugate pairs λ,λ― with λλ―=|λ|2=1 for an isometry, so they contribute 1 to the determinant and an even number to the count. With n odd, an odd number of eigenvalues are real, each ±1, and their product is detA=1; an odd number of factors ±1 with product 1 cannot all be −1, so 1 is an eigenvalue and some nonzero vector is fixed.

  • (e) The rotations Rθ, with determinant 1, and the reflections Mθ, with determinant −1 (Theorem 11.7.6).

  • (f) A unitary operator is an isometry, so Proposition 11.7.1 applies: |λ|=1.

B. Practice

Exercise (B1: Determine which are unitary)

Determine which of the following matrices are orthogonal or unitary. Justify your answer; for those that are not, name the failing condition.

  • (a) 13(2−21122−2−12) over R.
  • (b) 12(1ii1) over C.
  • (c) (1101) over R.
  • (d) 15(3443) over R.

Solution.

  • (a) Orthogonal, as computed in Example (Axis and Angle of a Rotation): the three columns are orthonormal.

  • (b) Unitary. The columns are 12(1,i) and 12(i,1), of norm 121+1=1, and their inner product is 12(1⋅i―+i⋅1―)=12(−i+i)=0. So A∗A=I2. (It is unitary but not orthogonal, since it is not real.)

  • (c) Not orthogonal. The second column (1,1) has norm 2≠1, so the columns are not orthonormal. Equivalently, the map stretches e2.

  • (d) Not orthogonal. Both columns are unit vectors, 15(3,4) and 15(4,3), but their inner product is 125(12+12)=2425≠0. Unit columns are not enough; they must also be orthogonal. (Its determinant is 125(9−16)=−725, neither 1 nor −1, which confirms it by Proposition 11.7.4.)

Exercise (B2: Name the plane isometry)

Each of the following lies in O(2). Decide whether it is a rotation or a reflection, and give the angle of rotation or the mirror line.

  • (a) (0110).
  • (b) 12(11−11).
  • (c) 113(51212−5).

Solution.

  • (a) Determinant −1, so a reflection. Ax=x reads x2=x1, so the mirror is L=span((1,1)), the line x2=x1. Check: A(1,−1)=(−1,1)=−(1,−1).

  • (b) Determinant 12(1+1)=1, so a rotation Rθ. Comparing with Rθ, cos⁡θ=12 and sin⁡θ=−12, so θ=−π/4: a clockwise half of a right angle.

  • (c) Determinant 1169(−25−144)=−1, so a reflection. Ax=x reads 5x1+12x2=13x1, that is 3x2=2x1, so L=span((3,2)). Check: A(3,2)=113(15+24,36−10)=(3,2).

Exercise (B3: Permutations are orthogonal)

Let σ∈Sn and let Pσ∈Mn(R) be its permutation matrix (Definition 3.4.2), whose j-th column is eσ(j).

  • (a) Prove that Pσ∈O(n).
  • (b) Deduce that detPσ=±1, and identify the sign.
  • (c) Show that Pσ−1=Pσ⊤=Pσ−1.

Solution.

  • (a) Column j of Pσ is eσ(j), by Definition 3.4.2. So each column is a unit vector, and for i≠j the columns eσ(i) and eσ(j) are distinct standard basis vectors, because σ is injective, hence orthogonal. By Theorem 11.7.3 (f), Pσ∈O(n).

  • (b) By Proposition 11.7.4 (b), detPσ=±1. The sign is sgn(σ): permuting the columns of In by σ multiplies the determinant by sgn(σ) (Theorem 7.2.1), and detIn=1.

  • (c) Pσ−1=Pσ⊤ holds because Pσ is orthogonal, by (a). For the second equality, (Pσ)ij=1 exactly when i=σ(j), so (Pσ⊤)ij=(Pσ)ji=1 exactly when j=σ(i), that is when i=σ−1(j); and that is the defining condition for (Pσ−1)ij=1.

C. Going deeper

Exercise (C1: Determinant as a homomorphism)
  • (a) Prove that det:O(n)→{1,−1} is a surjective group homomorphism with kernel SO(n), and deduce that exactly half of O(n) has determinant −1, in the sense that O(n) is the disjoint union of SO(n) and one coset A0SO(n).
  • (b) Prove that det:U(n)→{z∈C:|z|=1} is a surjective group homomorphism.
  • (c) Describe O(1) and U(1) completely.

Solution.

  • (a) det is multiplicative (Theorem 7.5.4) and takes only the values ±1 on O(n) (Proposition 11.7.4 (b)), so it is a homomorphism into {1,−1}. It is surjective: detIn=1, and A0:=diag(−1,1,…,1) is orthogonal, since its columns are ±ei, with detA0=−1. Its kernel is {A:detA=1}=SO(n) by definition. If detA=−1 then det(A0−1A)=1, so A∈A0SO(n); and no matrix has both determinants, so the union is disjoint.

  • (b) Multiplicativity again, and |detA|=1 on U(n) by Proposition 11.7.4 (b). For surjectivity, given z with |z|=1, the matrix diag(z,1,…,1) has orthonormal columns, hence lies in U(n), and its determinant is z (Theorem 7.4.6).

  • (c) A 1×1 matrix (a) satisfies a―a=1 exactly when |a|=1. So U(1) is the unit circle of C and O(1)={(1),(−1)}. Consistently, SO(1)={(1)} is trivial.

Exercise (C2: Linearity is a hypothesis, not a consequence)
  • (a) Let c:C→C be c(z)=z―, where C carries ⟨z,w⟩=zw―. Show that ‖c(z)‖=‖z‖ for every z, that c is not linear over C, and that c does not preserve the inner product. Which step of Theorem 11.7.3 fails?
  • (b) Let V be a complex inner product space and let T∈L(V) satisfy ⟨Tv,v⟩=0 for every v∈V. Prove that T=0.
  • (c) Show that (b) fails over R by exhibiting a nonzero T∈L(R2) with ⟨Tv,v⟩=0 for all v.

Hint for (b): apply the hypothesis to v=u+w and to v=u+iw.

Solution.

  • (a) ‖c(z)‖=|z―|=|z|=‖z‖. But c(iz)=iz―=−iz―=−ic(z)≠ic(z) unless z=0, so c is not C-linear. And ⟨c(1),c(i)⟩=⟨1,−i⟩=1⋅−i―=i, while ⟨1,i⟩=−i. The proof of (a) ⇒ (b) in Theorem 11.7.3 uses T(u+ikv)=Tu+ikTv, which is exactly what c does not satisfy.

  • (b) Let u,w∈V. Applying the hypothesis to u+w and expanding both slots, (1)0=⟨T(u+w),u+w⟩=⟨Tu,u⟩+⟨Tu,w⟩+⟨Tw,u⟩+⟨Tw,w⟩=⟨Tu,w⟩+⟨Tw,u⟩. Applying it to u+iw, and using ⟨Tu,iw⟩=−i⟨Tu,w⟩, ⟨T(iw),u⟩=i⟨Tw,u⟩ and ⟨T(iw),iw⟩=ii―⟨Tw,w⟩=0, 0=−i⟨Tu,w⟩+i⟨Tw,u⟩. Multiplying by i gives (2)⟨Tu,w⟩−⟨Tw,u⟩=0. Adding (1) and (2), 2⟨Tu,w⟩=0, so ⟨Tu,w⟩=0 for all u,w. Taking w=Tu gives ‖Tu‖2=0, so Tu=0 for every u, that is T=0.

  • (c) Let T be the rotation Rπ/2, that is T(x1,x2)=(−x2,x1). Then ⟨Tx,x⟩=−x2x1+x1x2=0 for every x, and T≠0. The complex proof breaks down because the second equation is unavailable: there is no i to substitute.

Exercise (C3: Two reflections make a rotation)

Let α,β∈R and let Mα,Mβ∈O(2) be the reflections of Theorem 11.7.6, in the lines Lα and Lβ at angles α/2 and β/2.

  • (a) Prove that MαMβ=Rα−β.
  • (b) Deduce that the composite of the reflections in two lines through the origin is the rotation by twice the angle from the second line to the first, and that every rotation of the plane is a product of two reflections.
  • (c) Deduce that O(2) is not abelian, by computing MβMα.

Solution.

  • (a) Multiplying out, MαMβ=(cos⁡αsin⁡αsin⁡α−cos⁡α)(cos⁡βsin⁡βsin⁡β−cos⁡β)=(c−ssc), where c=cos⁡αcos⁡β+sin⁡αsin⁡β=cos⁡(α−β) and s=sin⁡αcos⁡β−cos⁡αsin⁡β=sin⁡(α−β), by the subtraction formulas. The top-right entry is cos⁡αsin⁡β−sin⁡αcos⁡β=−sin⁡(α−β), as displayed. So MαMβ=Rα−β.

  • (b) The angle from Lβ to Lα is α2−β2, and (a) says the composite is the rotation by α−β, twice that. Conversely, given θ, choose α=θ and β=0; then Rθ=MθM0, a product of two reflections.

  • (c) By (a) with the roles exchanged, MβMα=Rβ−α=Rα−β−1. These two agree only when Rα−β2=I, that is when α−β is a multiple of π. Taking α=π/2 and β=0 gives Mπ/2M0=Rπ/2≠R−π/2=M0Mπ/2, so O(2) is not abelian.

Found a typo or a mistake? Report it, or read this page’s source

The text of this book is licensed under CC BY-NC-SA 4.0.