Unitary and Orthogonal Operators

Every inner-product space carries a geometry — lengths, angles, orthonormal frames — and one question organizes an enormous amount of physics and mathematics: which linear maps leave that geometry exactly intact? The answer is a single, beautifully rigid class. On a complex space they are the unitary operators, defined by U†U=IU^\dagger U = I; on a real space they are the orthogonal operators, defined by QTQ=IQ^\mathsf{T} Q = I. These are the linear analogues of rigid motions that fix the origin: rotations, reflections, and (in the complex case) phase rotations. Nothing is stretched, sheared, or squashed. What makes the subject so powerful is that one algebraic equation is secretly a whole dictionary of geometric statements. From U†U=IU^\dagger U = I alone it follows that UU preserves every inner product and every norm, sends orthonormal bases to orthonormal bases, has columns that form an orthonormal set, is invertible with U−1=U†U^{-1} = U^\dagger, has all its eigenvalues on the unit circle ∣λ∣=1|\lambda| = 1, and has ∣det⁡U∣=1|\det U| = 1. Because UU commutes with its own adjoint it is normal, so the spectral theorem hands us an orthonormal eigenbasis for free. These operators assemble into the continuous symmetry groups that scaffold modern physics — the unitary group U(n)U(n), the special unitary group SU(n)SU(n), the orthogonal group O(n)O(n), and the rotation group SO(n)SO(n). Unitarity is why quantum probabilities are conserved under time evolution, why observables diagonalize by a change of orthonormal frame, and why the stable workhorses of computation never amplify error. This article builds the linear-algebraic core; the group and Lie-algebra structure grows straight out of it.

Key takeaways
  • A unitary operator satisfies U†U=UU†=IU^\dagger U = U U^\dagger = I; an orthogonal operator satisfies QTQ=QQT=IQ^\mathsf{T} Q = Q Q^\mathsf{T} = I. They are the same idea over C\mathbb{C} and over R\mathbb{R}.
  • That one equation is equivalent to preserving the inner product, ⟨Ux,Uy⟩=⟨x,y⟩\langle Ux, Uy\rangle = \langle x,y\rangle — hence all lengths and angles. These are the rigid rotations, reflections, and phase rotations about the origin.
  • The columns of UU form an orthonormal basis, so a unitary matrix is the dictionary that translates between two orthonormal frames.
  • All eigenvalues lie on the unit circle, ∣λ∣=1|\lambda| = 1; the determinant obeys ∣det⁡U∣=1|\det U| = 1, with det⁡Q=±1\det Q = \pm 1 splitting rotations from improper motions.
  • Unitary operators are normal, so the spectral theorem gives an orthonormal eigenbasis U=∑keiθkPkU = \sum_k e^{i\theta_k}P_k with unit-modulus eigenvalues.
  • Key idea: one equation, U†U=IU^\dagger U = I (or QTQ=IQ^\mathsf{T} Q = I in the real case), is simultaneously equivalent to preservation of the inner product, of all norms and angles, of orthonormality of the columns, and to invertibility with U−1=U†U^{-1} = U^\dagger — so whenever a transformation must keep probabilities, lengths, or orthonormal frames intact, it must be unitary or orthogonal.
The core insight

A single equation — U†U=IU^\dagger U = I (or QTQ=IQ^\mathsf{T} Q = I in the real case) — is equivalent to preservation of the inner product, of all norms and angles, of orthonormality of the columns, and to invertibility with U−1=U†U^{-1} = U^\dagger. Whenever a transformation must keep probabilities, lengths, or orthonormal frames intact, it must be unitary or orthogonal.

Where this lives today

Unitarity is a daily working constraint across physics and computation. In quantum computing every gate is a unitary matrix — the Hadamard and SU(2)SU(2) gates of Examples 2 and 3 are literal hardware operations, and an algorithm is a product of unitaries steering a state around the unit sphere without losing a drop of probability; unitary synthesis asks how to decompose an arbitrary U∈SU(2n)U \in SU(2^n) into available gates. In metrology, SO(3)SO(3) governs the attitude estimation that keeps satellites and drones pointed, best done through unit quaternions, the double cover SU(2)→SO(3)SU(2) \to SO(3). In large-scale computation, orthogonal transforms are prized because their condition number is exactly one: Householder reflections and Givens rotations power the QR eigenvalue algorithm, and the discrete Fourier transform — unitary up to normalization — underlies the FFT. And in fundamental physics the unitarity of the S-matrix — that scattering conserve total probability — is what forced the prediction, and discovery, of the Higgs boson from the electroweak amplitudes.

The physical picture

Picture a rigid tripod of three mutually perpendicular unit arrows glued at the origin. An orthogonal operator on R3\mathbb{R}^3 is any way of re-orienting that tripod — spinning it, flipping it — without bending an arrow or changing an angle between two of them. That indestructible rigidity is precisely QTQ=IQ^\mathsf{T} Q = I: the columns of QQ are the images of the tripod's arrows, and the equation says those images are still orthonormal. A unitary operator does the same for complex "arrows", with one extra freedom — each complex coordinate may also pick up a phase eiθe^{i\theta}, an internal rotation invisible to real geometry but essential to quantum amplitudes. The first figure makes the two-dimensional version concrete with real numbers.

Figure
An orthogonal operator rotates an orthonormal frame. R at 30° sends ê₁ and ê₂ to the columns (0.866, 0.500) and (−0.500, 0.866); each image stays unit length and the two stay perpendicular, so QᵀQ = I and det Q = +1 — a proper rotation in SO(2)
QᵀQ = I : the image of an orthonormal frame is still orthonormal ê₁ ê₂ Qê₁ Qê₂ θ = 30° R(30°) = ( 0.866 −0.500 ) ( 0.500 0.866 ) Qê₁ · Qê₂ = 0 ‖Qê₁‖ = ‖Qê₂‖ = 1 det Q = +1 (proper)

The condition ⟨Ux,Uy⟩=⟨x,y⟩\langle Ux, Uy\rangle = \langle x, y\rangle says the operator moves the whole space the way a rotation does — no direction favored, no pair sheared, no volume inflated. This is why unitarity is the natural home of conservation: preserved lengths mean the squared norm of a quantum state is frozen, so total probability cannot leak away. The deepest fact is spectral: because UU commutes with U†U^\dagger it is diagonalizable in an orthonormal eigenbasis, and in each eigendirection it merely multiplies by a unit-modulus eiθe^{i\theta}. The whole operator is a bundle of independent rotations, one per eigenplane, every eigenvalue pinned to the unit circle. The second figure shows that circle.

Figure
Unitary spectrum lives on the unit circle. Every eigenvalue is e^{iθ} of modulus one; the rotation R at 30° carries the conjugate pair e^{±i30°} = 0.866 ± 0.500 i, while a real orthogonal map can also carry the real eigenvalues +1 and −1 that build reflections
The punchline: unitary ⇒ every eigenvalue sits on |λ| = 1 Re λ Im λ θ e^{iθ} e^{−iθ} +1 −1 rotation R(30°) λ = e^{±i30°} = 0.866 ± 0.500 i real spectrum +1 : fixed axis −1 : reflection det Q = (+1)(−1)^(#−1)

The core equations

The two definitions look almost identical, and that is the point: an orthogonal operator is the real specialization of a unitary one, the adjoint †\dagger (conjugate-transpose) collapsing to the transpose T\mathsf{T} for real entries. We state both, then unpack the geometric fact they encode.

U†U=UU†=I⟺U−1=U†(1)U^\dagger U = U U^\dagger = I \quad\Longleftrightarrow\quad U^{-1} = U^\dagger \qquad (1)

where: UU — operator on a complex inner-product space VV; U†U^\dagger — Hermitian adjoint, ⟨U†x,y⟩=⟨x,Uy⟩\langle U^\dagger x, y\rangle = \langle x, U y\rangle; II — identity. Reading: the adjoint is the inverse — undoing UU costs only a conjugate-transpose.

QTQ=QQT=I⟺Q−1=QT(2)Q^\mathsf{T} Q = Q Q^\mathsf{T} = I \quad\Longleftrightarrow\quad Q^{-1} = Q^\mathsf{T} \qquad (2)

where: QQ — operator on a real inner-product space; QTQ^\mathsf{T} — transpose, the real adjoint. Reading: identical to (1); everything below specializes by †→T\dagger \to \mathsf{T}. The whole geometric meaning is compressed into one line: UU preserves the inner product, hence length and angle. This is the equation to remember — the algebra U†U=IU^\dagger U = I and the geometry below are logically equivalent.

⟨Ux,Uy⟩=⟨x,y⟩⟹∥Ux∥=∥x∥(3)\langle Ux, Uy\rangle = \langle x, y\rangle \qquad \Longrightarrow \qquad \lVert Ux\rVert = \lVert x\rVert \qquad (3)

where: ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle — inner product on VV; ∥x∥=⟨x,x⟩\lVert x\rVert = \sqrt{\langle x, x\rangle} — induced norm. Reading: norm preservation is the case y=xy = x; via polarization the two are equivalent over C\mathbb{C}. In components (orthonormal basis), U†U=IU^\dagger U = I becomes a statement about the columns of UU: they form an orthonormal set — the most practical hand-check for unitarity.

∑i(U)ik‾ (U)iℓ=δkℓ(4)\sum_{i} \overline{(U)_{ik}}\,(U)_{i\ell} = \delta_{k\ell} \qquad (4)

where: (U)ik(U)_{ik} — matrix entry in an orthonormal basis; (⋅)‾\overline{(\cdot)} — complex conjugate; δkℓ\delta_{k\ell} — Kronecker delta. Reading: column kk dotted with column ℓ\ell gives δkℓ\delta_{k\ell} — columns and rows are orthonormal frames. Norm preservation constrains the spectrum: an eigenvector's length cannot change, so its eigenvalue cannot have modulus other than one.

Uv=λv  ⇒  ∣λ∣=1  ⇒  λ=eiθ(5)U v = \lambda v \;\Rightarrow\; |\lambda| = 1 \;\Rightarrow\; \lambda = e^{i\theta} \qquad (5)

where: λ\lambda — eigenvalue; v≠0v \neq 0 — eigenvector; θ\theta — real phase angle. Reading: a real orthogonal operator therefore admits only ±1\pm 1 as real eigenvalues, the rest in conjugate pairs e±iθe^{\pm i\theta}. The determinant inherits the same discipline, and separates the two connected pieces of the real orthogonal group.

∣det⁡U∣=1,det⁡Q=±1(6)|\det U| = 1, \qquad \det Q = \pm 1 \qquad (6)

where: det⁡\det — determinant. Reading: det⁡Q=+1\det Q = +1 is a proper rotation in SO(n)SO(n); det⁡Q=−1\det Q = -1 is an improper motion (rotation plus reflection), and the two never connect continuously. Because U†U=UU†U^\dagger U = U U^\dagger, UU is normal, and the spectral theorem gives an orthonormal eigenbasis: UU is a sum of phases times orthogonal projectors.

U=∑keiθk Pk,∑kPk=I,PjPk=δjkPk(7)U = \sum_{k} e^{i\theta_k}\, P_k, \qquad \sum_k P_k = I, \qquad P_j P_k = \delta_{jk} P_k \qquad (7)

where: PkP_k — orthogonal projector onto the kk-th eigenspace; eiθke^{i\theta_k} — unit-modulus eigenvalue. Reading: UU is a bundle of independent phase rotations. Two operational forms close the set. First, a unitary is the change between orthonormal frames, a geometry-preserving similarity.

A′=U†A U(8)A' = U^\dagger A\, U \qquad (8)

where: AA — operator in the original orthonormal basis; A′A' — the same operator in the new frame; UU — unitary whose columns are the new basis vectors. Reading: diagonalizing a Hermitian observable means choosing UU so A′A' is diagonal. Second, the canonical building block of every real orthogonal matrix is the 2×22\times 2 rotation block, alongside ±1\pm 1 entries.

R(θ)=(cos⁡θ−sin⁡θsin⁡θcos⁡θ),det⁡R=cos⁡2θ+sin⁡2θ=1(9)R(\theta) = \begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix}, \qquad \det R = \cos^2\theta + \sin^2\theta = 1 \qquad (9)

where: θ\theta — rotation angle; R(θ)∈SO(2)R(\theta) \in SO(2). Reading: stacked with ±1\pm 1s, this block reconstructs any rotation or reflection in any dimension.

Derivation AOne equation, four geometric facts

The goal is to prove that the algebraic definition U†U=IU^\dagger U = I is not merely sufficient for the geometry but equivalent to it, so that "preserves the inner product", "preserves every norm", and "U†U=IU^\dagger U = I" are three names for one property. The forward direction is a line of adjoint algebra; the reverse — recovering the full inner product from norms alone — is where the complex polarization identity does the real work.

Step 1. Using the defining property of the adjoint, ⟨Ux,Uy⟩=⟨x,U†Uy⟩\langle Ux, Uy\rangle = \langle x, U^\dagger U y\rangle, and if U†U=IU^\dagger U = I this collapses to ⟨x,y⟩\langle x, y\rangle. So the algebra preserves every inner product. Why: the adjoint is defined precisely so UU can be slid across the bracket, converting the operator equation into a statement about brackets.

Step 2. Specialize to norms by setting y=xy = x: ∥Ux∥2=⟨Ux,Ux⟩=⟨x,x⟩=∥x∥2\lVert Ux\rVert^2 = \langle Ux, Ux\rangle = \langle x, x\rangle = \lVert x\rVert^2. Preserving inner products preserves all lengths. Why: the norm is the inner product on the diagonal, so it carries no more information than the bracket it comes from.

Step 3. The subtle reverse: suppose only ∥Ux∥=∥x∥\lVert Ux\rVert = \lVert x\rVert for all xx. The complex polarization identity ⟨x,y⟩=14∑n=03in ∥x+iny∥2\langle x, y\rangle = \tfrac14\sum_{n=0}^{3} i^{n}\,\lVert x + i^{n} y\rVert^2 rebuilds every inner product from four norms, each of which UU preserves. Why: over C\mathbb{C} the norm secretly contains the full bracket — this fails on a merely normed space, which is why the theorem is special to inner-product spaces.

Step 4. Close the loop: inner-product preservation gives ⟨x,(U†U−I)y⟩=0\langle x, (U^\dagger U - I) y\rangle = 0 for all x,yx, y, and non-degeneracy forces U†U−I=0U^\dagger U - I = 0. The three characterizations are equivalent, and column orthonormality (4) and unit-modulus eigenvalues (5) follow as corollaries. Why: an operator that vanishes inside every bracket is the zero operator — this is what "non-degenerate" buys us.

Go deeper · Derivation B · 1 minEigenvalues on the unit circle and the spectral form
Derivation BEigenvalues on the unit circle and the spectral form

The aim is to show that a unitary operator not only has unit-modulus eigenvalues but automatically supplies orthogonal eigenvectors for distinct eigenvalues — exactly the raw material the spectral theorem needs to assemble an orthonormal eigenbasis. We derive the modulus bound from norm preservation, prove the orthogonality from inner-product preservation, and package the result as the projector decomposition of equation (7).

Step 1. Apply norm preservation to Uv=λvUv = \lambda v: ∥v∥=∥Uv∥=∣λ∣ ∥v∥\lVert v\rVert = \lVert Uv\rVert = |\lambda|\,\lVert v\rVert, so dividing by ∥v∥≠0\lVert v\rVert \neq 0 gives ∣λ∣=1|\lambda| = 1, hence λ=eiθ\lambda = e^{i\theta}. Why: an eigenvector is the one direction where UU acts by pure scaling, so if lengths are frozen the scale factor can only rotate.

Step 2. Take Uv=λvUv = \lambda v, Uw=μwUw = \mu w with λ≠μ\lambda \neq \mu. Preservation gives ⟨v,w⟩=⟨Uv,Uw⟩=λ‾ μ ⟨v,w⟩\langle v, w\rangle = \langle Uv, Uw\rangle = \overline{\lambda}\,\mu\,\langle v, w\rangle, the first factor conjugating because the bracket is antilinear in its first slot. Why: pushing both vectors through UU leaves the bracket unchanged yet multiplies out the eigenvalues, so the bracket equals a scalar multiple of itself.

Step 3. Use λ‾=1/λ\overline{\lambda} = 1/\lambda to get (1−μ/λ)⟨v,w⟩=0(1 - \mu/\lambda)\langle v, w\rangle = 0; since μ≠λ\mu \neq \lambda the prefactor is nonzero, forcing ⟨v,w⟩=0\langle v, w\rangle = 0. Eigenvectors of distinct eigenvalues are orthogonal. Why: unit modulus turns λ‾\overline{\lambda} into 1/λ1/\lambda, the trick that collapses the self-consistency equation to orthogonality.

Step 4. Within each eigenspace apply Gram–Schmidt; across eigenspaces the vectors are already orthogonal, so the space acquires an orthonormal eigenbasis and U=∑keiθkPkU = \sum_k e^{i\theta_k} P_k with ∑kPk=I\sum_k P_k = I. Why: orthogonality between eigenspaces plus orthonormality within them is a complete orthonormal eigenbasis — the projector form of the spectral theorem for normal operators.

Go deeper · Derivation C · 1 minReal orthogonal matrices are stacks of rotation blocks
Derivation CReal orthogonal matrices are stacks of rotation blocks

The goal is the canonical form of a real orthogonal matrix: every QQ is real-orthogonally similar to a block-diagonal matrix built from 2×22\times 2 rotations R(θ)R(\theta) together with ±1\pm 1 entries. This is the precise sense in which "orthogonal" means "rotation and reflection", and why a rigid rotation in any dimension organizes into independent rotation planes. The one obstruction — a real matrix may have genuinely complex eigenvalues — is resolved by folding each conjugate pair into a single real invariant plane.

Step 1. Fix the spectrum. Since QQ is orthogonal it is unitary (∣λ∣=1|\lambda| = 1), and since QQ is real, non-real eigenvalues come in conjugate pairs. Hence λ∈{+1,−1,e±iθ}\lambda \in \{+1, -1, e^{\pm i\theta}\}. Why: reality forces complex roots to appear with their conjugates, and the unit-circle bound removes everything except ±1\pm 1 and conjugate phase pairs.

Step 2. Turn a conjugate pair into a real plane. Write v=a+ibv = a + i b (a,ba, b real); expanding Qv=eiθvQv = e^{i\theta} v into real and imaginary parts gives Qa=cos⁡θ a−sin⁡θ bQa = \cos\theta\, a - \sin\theta\, b, Qb=sin⁡θ a+cos⁡θ bQb = \sin\theta\, a + \cos\theta\, b, so span{a,b}\mathrm{span}\{a, b\} is a QQ-invariant plane. Why: a real matrix cannot single out a complex direction, but it preserves the real plane spanned by its real and imaginary parts, on which it acts as an honest rotation.

Step 3. Normalize a,ba, b to an orthonormal pair; there QQ is the standard block R(θ)R(\theta). Stacking one block per conjugate pair and a ±1\pm 1 per real eigenvalue gives PTQP=R(θ1)⊕⋯⊕(±1)⊕⋯P^\mathsf{T} Q P = R(\theta_1) \oplus \cdots \oplus (\pm 1) \oplus \cdots. Why: an orthonormal basis inside each invariant plane makes the local action the textbook rotation matrix, and orthogonal blocks assemble into an orthogonal whole.

Step 4. Read off the determinant: each block gives det⁡R=+1\det R = +1, so det⁡Q=(−1)#(−1)\det Q = (-1)^{\#(-1)}, set by the number of −1-1 eigenvalues — even keeps QQ in SO(n)SO(n), odd makes it improper. Why: the sign of the determinant is the parity of reflections, the invariant separating the two disconnected components of O(n)O(n). *Note: for R(θ)R(\theta) the determinant is cos⁡2θ+sin⁡2θ=1\cos^2\theta + \sin^2\theta = 1, confirming the block form gives a proper rotation.*

Limiting cases

Pushing the definition to its edges shows what unitarity guarantees — and, in the last row, where finite-dimensional intuition breaks.

CaseResultMeaning
Real specializationU∈Rn×n, U†=UT⇒QTQ=IU \in \mathbb{R}^{n\times n},\ U^\dagger = U^\mathsf{T} \Rightarrow Q^\mathsf{T} Q = IA real unitary operator is an orthogonal operator; the Hermitian adjoint collapses to the transpose.
One dimensionU=eiθ, ∣U∣=1U = e^{i\theta},\ |U| = 1On C\mathbb{C} a unitary is a pure phase; on R\mathbb{R} an orthogonal operator is ±1\pm 1 — identity or reflection.
Determinant subgroupsU(n):det⁡U=eiϕ; SU(n):det⁡U=1; O(n):det⁡Q=±1; SO(n):det⁡Q=1U(n):\det U = e^{i\phi};\ SU(n):\det U = 1;\ O(n):\det Q = \pm 1;\ SO(n):\det Q = 1det⁡=+1\det = +1 selects proper rotations; O(n)O(n) splits into two disconnected sheets.
Isometry vs unitary (infinite dim.)S†S=I, SS†≠IS^\dagger S = I,\ S S^\dagger \neq IThe unilateral shift is norm-preserving but not surjective — an isometry that is not unitary.

Worked examples

1 · A rotation is orthogonal — check every claim with numbers (fully worked)

Take the plane rotation by θ=30°\theta = 30°, R=(cos⁡30°−sin⁡30°sin⁡30°cos⁡30°)R = \begin{pmatrix} \cos 30° & -\sin 30° \\ \sin 30° & \cos 30° \end{pmatrix}. Verify orthogonality, the determinant, and the eigenvalues.

Step 1. Insert cos⁡30°=0.8660\cos 30° = 0.8660, sin⁡30°=0.5000\sin 30° = 0.5000, so R=(0.8660−0.50000.50000.8660)R = \begin{pmatrix} 0.8660 & -0.5000 \\ 0.5000 & 0.8660 \end{pmatrix} with columns c1=(0.8660,0.5000)c_1 = (0.8660, 0.5000), c2=(−0.5000,0.8660)c_2 = (-0.5000, 0.8660).

Step 2. Check RTR=IR^\mathsf{T} R = I column by column: c1⋅c1=0.7500+0.2500=1.0000c_1 \cdot c_1 = 0.7500 + 0.2500 = \mathbf{1.0000}, identically c2⋅c2=1.0000c_2 \cdot c_2 = 1.0000, and the cross term c1⋅c2=−0.4330+0.4330=0c_1 \cdot c_2 = -0.4330 + 0.4330 = \mathbf{0}. Columns orthonormal, so RTR=IR^\mathsf{T} R = I — equation (4) made arithmetic.

Step 3. Determinant: det⁡R=0.7500+0.2500=+1\det R = 0.7500 + 0.2500 = \mathbf{+1}. Positive determinant means a proper rotation in SO(2)SO(2) — no reflection.

Step 4. Eigenvalues from λ2−(tr⁡R)λ+det⁡R=0\lambda^2 - (\operatorname{tr} R)\lambda + \det R = 0, i.e. λ2−2(0.8660)λ+1=0\lambda^2 - 2(0.8660)\lambda + 1 = 0. The discriminant is 3−4=−1<03 - 4 = -1 < 0, giving λ=0.8660±0.5000 i=e±i30°\lambda = 0.8660 \pm 0.5000\, i = \mathbf{e^{\pm i 30°}}, both on the unit circle (∣λ∣=1|\lambda| = 1) as equation (5) demands.

Step 5. A real rotation has no real eigenvector — nothing is left pointing the same way — which is why the eigenvalues are the complex pair e±iθe^{\pm i\theta} acting on the single rotation plane. RTR=IR^\mathsf{T} R = I, det⁡R=+1\det R = +1, eigenvalues e±i30°=0.866±0.500 ie^{\pm i30°} = 0.866 \pm 0.500\, i — a proper rotation of SO(2)SO(2).

2 · The Hadamard gate — orthogonal and unitary, but a reflection

The Hadamard gate is H=12(111−1)H = \tfrac{1}{\sqrt 2}\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}. Show it is both orthogonal and unitary, find its determinant and eigenvalues, and identify the motion.

Step 1. HH is real symmetric, so HT=H†=HH^\mathsf{T} = H^\dagger = H. Then HTH=H2=12(2002)=IH^\mathsf{T} H = H^2 = \tfrac12\begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix} = I, so HH is orthogonal, and being real with H†H=IH^\dagger H = I it is unitary too. It is also an involution: H2=IH^2 = I, so H=H−1H = H^{-1}.

Step 2. Determinant: det⁡H=12[(1)(−1)−(1)(1)]=−1\det H = \tfrac12\big[(1)(-1) - (1)(1)\big] = \mathbf{-1}. Negative determinant means HH is improper — not a pure rotation.

Step 3. Eigenvalues: from H2=IH^2 = I every λ2=1\lambda^2 = 1, so λ=±1\lambda = \pm 1; the trace tr⁡H=0\operatorname{tr} H = 0 fixes them as {+1,−1}\{+1, -1\}, consistent with det⁡H=−1\det H = -1. Both have ∣λ∣=1|\lambda| = 1, as required.

Step 4. Identify the motion. An improper 2×22\times 2 orthogonal matrix with eigenvalues +1,−1+1, -1 is a reflection (the +1+1 eigenvector is the mirror axis, the −1-1 is flipped). Matching HH to the reflection form (cos⁡2φsin⁡2φsin⁡2φ−cos⁡2φ)\begin{pmatrix} \cos 2\varphi & \sin 2\varphi \\ \sin 2\varphi & -\cos 2\varphi \end{pmatrix} gives cos⁡2φ=1/2\cos 2\varphi = 1/\sqrt 2, so φ=22.5°\varphi = 22.5°. HH is orthogonal and unitary with det⁡=−1\det = -1 and eigenvalues ±1\pm 1 — a reflection about the 22.5°22.5° axis, which is why it maps ∣0⟩,∣1⟩|0\rangle, |1\rangle to the two equal superpositions.

3 · A genuinely complex unitary — a qubit rotation in SU(2)SU(2)

Consider U=12(1ii1)U = \tfrac{1}{\sqrt 2}\begin{pmatrix} 1 & i \\ i & 1 \end{pmatrix}. Confirm it is unitary, find its determinant and eigenvalues, and connect it to a Hermitian generator.

Step 1. Adjoint (conjugate-transpose): U†=12(1−i−i1)U^\dagger = \tfrac{1}{\sqrt 2}\begin{pmatrix} 1 & -i \\ -i & 1 \end{pmatrix}. Then U†U=12(1+1i−i−i+i1+1)=IU^\dagger U = \tfrac12\begin{pmatrix} 1+1 & i - i \\ -i + i & 1 + 1 \end{pmatrix} = I, so UU is unitary — the conjugation is essential; the plain transpose would not give II.

Step 2. Determinant: det⁡U=12[(1)(1)−(i)(i)]=+1\det U = \tfrac12\big[(1)(1) - (i)(i)\big] = \mathbf{+1}. Unitary with det⁡=1\det = 1 places UU in SU(2)SU(2).

Step 3. Eigenvalues: trace 2\sqrt 2, det⁡=1\det = 1, so λ2−2 λ+1=0\lambda^2 - \sqrt 2\,\lambda + 1 = 0 gives λ=12(1±i)=e±iπ/4\lambda = \tfrac{1}{\sqrt2}(1 \pm i) = \mathbf{e^{\pm i\pi/4}}, on the unit circle.

Step 4. Connect to a generator. With X=(0110)X = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} (Hermitian), eiαX=cos⁡α I+isin⁡α Xe^{i\alpha X} = \cos\alpha\, I + i\sin\alpha\, X; at α=π/4\alpha = \pi/4 this is 12(I+iX)=U\tfrac{1}{\sqrt2}(I + iX) = U. So U=ei(π/4)XU = e^{i(\pi/4)X} — a Bloch-sphere rotation generated by XX, and det⁡U=ei(π/4)tr⁡X=1\det U = e^{i(\pi/4)\operatorname{tr} X} = 1 confirms SU(2)SU(2). UU is unitary with det⁡=1\det = 1 and eigenvalues e±iπ/4e^{\pm i\pi/4}; it equals ei(π/4)Xe^{i(\pi/4)X} — the recipe "unitary = exponential of ii times Hermitian".

Physics of the portrait

Unitary and orthogonal operators are how pure linear algebra became the grammar of physical symmetry — a thread from Cauchy's spectral theorem, through Hermite's forms, to the moment quantum mechanics made unitarity synonymous with conservation of probability. No one drew that synthesis more sharply than Hermann Weyl.

Origins · 1 minHermann Weyl and how this was discovered
HW
Hermann Weyl
1885–1955 · linear symmetry made the language of quantum mechanics

Hermann Weyl was the mathematician who taught physics to speak in groups of unitary and orthogonal operators. A student of Hilbert at Göttingen, he saw that the symmetries of nature are exactly the transformations that preserve an inner product. In his 1928 Gruppentheorie und Quantenmechanik he set out the modern view: quantum states live in a Hilbert space, observables are Hermitian operators, and every symmetry — rotation, reflection, phase change — is a unitary (or anti-unitary) operator leaving ⟨ψ∣ψ⟩\langle\psi|\psi\rangle invariant. Unitarity was no longer a technical condition on matrices; it was the mathematical form of "probability is conserved". His 1939 masterwork The Classical Groups then worked out the invariants and representations of U(n)U(n), SU(n)SU(n), O(n)O(n), and SO(n)SO(n), the very groups this article's operators generate, and his "gauge" idea — a local phase, unitary at every point of spacetime — became the foundation of the Standard Model. Weyl's conviction, that the content of a physical law is the group under which it is invariant, turned U†U=IU^\dagger U = I into a research program that has organized theoretical physics ever since.

History
1829
— Cauchy proves the eigenvalues of a real symmetric matrix are real, founding spectral theory
1855
— Hermite studies Hermitian forms and shows their eigenvalues are real
1897
— Hurwitz constructs the invariant integral on the orthogonal and unitary groups
1928
— Weyl publishes Gruppentheorie und Quantenmechanik, casting symmetries as unitary operators
1932
— Stone proves that one-parameter unitary groups are exponentials of self-adjoint generators

Going further

Texts. Arfken and Weber, Mathematical Methods for Physicists, 7th ed., ch. 3 (matrices and orthogonal/unitary transformations) · Hassani, Mathematical Physics, 2nd ed., chs. on operators and the spectral decomposition · Stone and Goldbart, Mathematics for Physics, ch. on linear operators and Hilbert space · Nakahara, Geometry, Topology and Physics, 2nd ed., §4 on Lie groups U(n)U(n), SU(n)SU(n), O(n)O(n), SO(n)SO(n) · Byron and Fuller, Mathematics of Classical and Quantum Physics, chs. 3–5 · Halmos, Finite-Dimensional Vector Spaces, §§ on unitary and orthogonal transformations. Key results to keep in mind: the complex polarization identity rebuilds ⟨x,y⟩\langle x, y\rangle from four norms; ∣λ∣=1|\lambda| = 1 follows from ∥Uv∥=∥v∥\lVert Uv\rVert = \lVert v\rVert; ∣det⁡U∣2=det⁡U‾ det⁡U=1|\det U|^2 = \overline{\det U}\,\det U = 1; R(30°)R(30°) has orthonormal columns, det⁡=+1\det = +1, and eigenvalues e±i30°=0.866±0.500 ie^{\pm i30°} = 0.866 \pm 0.500\,i; the Hadamard gate has det⁡=−1\det = -1 and eigenvalues ±1\pm 1; and U=12(1ii1)U = \tfrac{1}{\sqrt2}\big(\begin{smallmatrix} 1 & i \\ i & 1 \end{smallmatrix}\big) is in SU(2)SU(2) with eigenvalues e±iπ/4e^{\pm i\pi/4} and equals ei(π/4)Xe^{i(\pi/4)X}.

Prerequisites. Inner-Product Spaces and Norms · Adjoint and Hermitian Operators · Orthogonality and the Gram–Schmidt Process · Linear Operators and Matrix Representations. Leads to. Diagonalization and the Spectral Decomposition · Change of Basis and Similarity Transformations · Eigenvalues, Eigenvectors, and the Characteristic Equation · Determinants and Their Properties.

Ask the tutor: Want to see why demanding a unitary S-matrix predicted the Higgs, or how unit quaternions in SU(2)SU(2) beat rotation matrices for spacecraft attitude — just ask.

Is an orthogonal operator just a real unitary operator?

Exactly. The two definitions are the same equation, A†A=IA^\dagger A = I, read over different fields. For complex entries the adjoint †\dagger means conjugate-transpose and the operator is unitary; for real entries conjugation does nothing, so †\dagger collapses to the transpose T\mathsf{T} and the operator is orthogonal. Every orthogonal matrix is a real unitary matrix, and every theorem about unitaries specializes by deleting the conjugation. The one visible difference is spectral: a real orthogonal matrix has a real characteristic polynomial, so its non-real eigenvalues are forced into conjugate pairs e±iθe^{\pm i\theta}, whereas a complex unitary can place a single eigenvalue anywhere on the unit circle.

Why are the eigenvalues forced onto the unit circle?

Because a unitary cannot change any length, and an eigenvector is the one direction where it acts by pure scaling. If Uv=λvUv = \lambda v then ∥v∥=∣λ∣ ∥v∥\lVert v\rVert = |\lambda|\,\lVert v\rVert, so ∣λ∣=1|\lambda| = 1. The only thing a unitary may do along an eigendirection is rotate the complex amplitude by a phase eiθe^{i\theta} — never stretch or shrink it. This is also why unitary evolution conserves probability: the squared norm is a sum of ∣amplitude∣2|\text{amplitude}|^2 terms, and if every eigen-amplitude only rotates, the total is frozen. A matrix with an eigenvalue off the unit circle is, by this one fact, not unitary.

If det U has modulus 1, is U automatically unitary?

No, and this is a common trap. ∣det⁡U∣=1|\det U| = 1 is a necessary consequence of unitarity but nowhere near sufficient — the determinant is a single number, while unitarity is n2n^2 conditions on the columns. The clean counterexample is the shear (1101)\begin{pmatrix} 1 & 1 \\ 0 & 1 \end{pmatrix}: its determinant is 11, yet it stretches the vector (0,1)(0,1) into (1,1)(1,1), changing its length from 11 to 2\sqrt 2, and its columns (1,0)(1,0) and (1,1)(1,1) are not orthogonal. It preserves area but not lengths or angles, so it is neither orthogonal nor unitary. To certify unitarity you must check U†U=IU^\dagger U = I directly, or equivalently verify that the columns form an orthonormal set.

What does the spectral form U = Σ e^{iθₖ} Pₖ actually buy you?

It says a unitary operator, however complicated in coordinates, is really just a set of independent phase rotations. Because UU is normal, the spectral theorem gives an orthonormal eigenbasis; grouping eigenvectors by eigenvalue yields orthogonal projectors PkP_k, and UU acts as multiplication by eiθke^{i\theta_k} on each eigenspace. Practically, you can compute any function of UU by acting on the phases: Um=∑keimθkPkU^m = \sum_k e^{im\theta_k}P_k, log⁡U=∑kiθkPk\log U = \sum_k i\theta_k P_k. It is also the bridge to Stone's theorem, which writes a one-parameter unitary group as U(t)=e−iHtU(t) = e^{-iHt} with Hermitian HH — the skeleton of quantum time evolution.

Why does every real rotation break into 2×2 blocks?

Because a real rotation has no reason to fix a single direction unless the dimension forces it to. A genuine rotation in a plane sends every vector elsewhere, so it has no real eigenvector there; instead it has the conjugate pair e±iθe^{\pm i\theta}, and the real and imaginary parts of its eigenvector span a plane the rotation simply spins. The canonical form collects one 2×22\times 2 block per rotating plane, plus a +1+1 per fixed axis and a −1-1 per reflected axis. In odd dimensions one real eigenvalue is always left over, which is why every rotation of 3-space fixes an axis — Euler's rotation theorem, seen here as a counting fact about the blocks.

How is unitarity connected to Hermitian operators and observables?

Through the exponential map, the single most useful bridge in the subject. If HH is Hermitian (H†=HH^\dagger = H) then U=eiHU = e^{iH} is unitary, since U†=e−iH=U−1U^\dagger = e^{-iH} = U^{-1}. Physically, Hermitian operators are the observables (their eigenvalues are the measured values), and exponentiating ii times an observable produces the symmetry it generates: energy generates time translation via e−iHte^{-iHt}, momentum generates translations, angular momentum generates rotations. So Hermitian and unitary operators are two faces of one structure — the Lie algebra of generators and the Lie group of unitaries they exponentiate to — which is why diagonalizing an observable (equation 8) and evolving a state (Stone's theorem) use the same machinery.

Before this

Where this leads