Mathematics · Upper Division

Linear Algebra for Physics

Matrices, eigenvalues, and inner product spaces are the language of quantum mechanics, classical mechanics (normal modes), and every numerical simulation. This chapter develops the tools physicists use most.

PrerequisitesVectors(Ch.V)Calculus(Ch.22)BasicmatrixmultiplicationVectors (Ch. V) \cdot Calculus (Ch. 22) \cdot Basic matrix multiplication
Learning Goals
  • Identify vector spaces and inner product spaces, and apply the Cauchy-Schwarz inequality.
  • Find eigenvalues and eigenvectors of a matrix by solving the characteristic equation.
  • State the spectral theorem for Hermitian operators and explain its role in quantum mechanics.
  • Compute the singular value decomposition and use it for low-rank approximation and least-squares fitting.
  • Determine whether two Hermitian operators share an eigenbasis by checking if they commute.

LA.1 Vector Spaces and Inner Products

Definition LA.1Vector Space
A vector space V over a field F (usually ℝ or ℂ) is a set with addition and scalar multiplication satisfying 8 axioms (commutativity, associativity, distributivity, identity, inverses). Key examples in physics:ℝⁿ:column vectors, forces, positions.  ℂⁿ: quantum states (spinors, finite-dimensional Hilbert spaces).L2(a,bL^{2}(a,b: square-integrable functions — the Hilbert space of quantum wavefunctions.Polynomials, matrices, operators — all are vector spaces under appropriate operations.

An inner product ⟨u, v⟩ on a complex vector space satisfies: linearity in v, conjugate-linearity in u, ⟨v,v⟩ ≥ 0 with equality iff v = 0. For ℂⁿ: ⟨u, v⟩ = u†v = Σᵢ uᵢ* vᵢ. For L²: ⟨f, g⟩ = ∫ f*(x) g(x) dx.

u,vuv(CauchySchwarzinequality)⟨u, v⟩ \le \|u\| \|v\| \qquad (Cauchy-Schwarz inequality)(LA.1)

LA.2 Eigenvalues and Eigenvectors

For a linear operator A: V → V, a nonzero vector v with Av = λv is an eigenvector with eigenvalue λ. Eigenvalues satisfy:

det(A\lambdaI)=0(characteristicequation)det(A - \lambdaI) = 0 \qquad (characteristic equation)(LA.2)

For an n×n matrix, this is an nth-degree polynomial — giving n eigenvalues (counting multiplicity). The eigenvalues determine: stability of linear systems (Re(λ) < 0 → stable), principal stresses in mechanics, normal mode frequencies, and the spectrum of quantum observables.

Example LA.1Normal Modes of Two Coupled Oscillators

Twomassesmconnectedbysprings(ktowalls,κbetweenthem).FindnormalfrequenciesTwo masses m connected by springs (k to walls, \kappa between them). Find normal frequencies.

Equations of motion:mx¨1=kx1κ(x1x2)mx¨2=kx2κ(x2x1mẍ_{1} = -kx_{1} - \kappa(x_{1}-x_{2}) \qquad mẍ_{2} = -kx_{2} - \kappa(x_{2}-x_{1}
Matrix form:ẍ = (1/m)KxK=[[k+κ,κ],[κ,k+κ-(1/m)Kx \qquad K = [[k+\kappa, -\kappa],[-\kappa, k+\kappa
Eigenvalues of K:det(K\muI)=0(k+κμ)2κ2=0μ=kork+2κdet(K - \muI) = 0 \to (k+\kappa-\mu)^{2} - \kappa^{2} = 0 \to \mu = k or k+2\kappa
Normal frequencies:ω1=(k/m)(inphasemode),ω2=((k+2κ)/m)(outofphasemode\omega_{1} = \sqrt(k/m) (in-phase mode), \omega_{2} = \sqrt((k+2\kappa)/m) (out-of-phase mode
Eigenvectors:v1=(1,1)/2(bothmassesmovetogether),v2=(1,1)/2(massesmoveoppositelyv_{1} = (1,1)/\sqrt2 (both masses move together), v_{2} = (1,-1)/\sqrt2 (masses move oppositely

LA.3 Hermitian Operators and Spectral Theorem

Definition LA.2Hermitian (Self-Adjoint) Operators
An operator A is Hermitian (selfadjoint)ifA=A(Aequalsitsconjugatetranspose).Keypropertiesself-adjoint) if A = A^\dagger (A equals its conjugate transpose). Key properties1.Alleigenvaluesarereal(provedbyvAv=Avvλ=λ1. All eigenvalues are real (proved by ⟨v|A|v⟩ = ⟨Av|v⟩ \to \lambda = \lambda*2. Eigenvectors for different eigenvalues are orthogonal3.(Spectraltheorem)Ahasacompleteorthonormaleigenbasis:A=iλivivi3. (Spectral theorem) A has a complete orthonormal eigenbasis: A = \sum_{i} \lambda_{i} |v_{i}⟩⟨v_{i}|In quantum mechanics, observables are Hermitian operators — this guarantees real measured values and allows complete state expansion in the eigenstate basis.

Analogously, unitary operators U satisfy U†U = I (preserve inner products). All eigenvalues lie on the unit circle (|λ| = 1). Time evolution e^(−iHt/ℏ) is unitary when H is Hermitian — unitarity is quantum mechanical probability conservation.

LA.4 Singular Value Decomposition

Every m×n matrix A can be written as:

A=U\sumV(SVD)A = U\sumV^\dagger \qquad (SVD)(LA.3)

where U (m×m) and V (n×n) are unitary and Σ is diagonal with non-negative real entries σ₁ ≥ σ₂ ≥ ... ≥ 0 (singular values). The number of nonzero σᵢ is the rank of A. SVD is the foundation of:

Principal Component Analysis (PCA): the columns of V are principal axes — directions of maximum variance in data.

Low-rank approximation: A ≈ U_kΣ_kV_k† (keep only k largest singular values) — used in image compression, recommendation systems, quantum state compression.

Least squares: the pseudoinverse A⁺ = VΣ⁺U† solves the least-squares problem min‖Ax − b‖ when A is non-square.

Theorem LA.1Spectral Theorem for Commuting Observables
Two Hermitian operators A and B can be simultaneously diagonalized(shareacompleteorthonormaleigenbasis)ifandonlyiftheycommute:[A,B]=agonalized (share a complete orthonormal eigenbasis) if and only if they commute: [A, B] =les can be simultaneously measured(sharpvalues)onlyiftheycommute.Positionandmomentumdontcommute([x^,p^]=sured (sharp values) only if they commute. Position and momentum don't commute ([x̂, p̂] =arp.L2andLzcommutebothcanbesimultaneouslyspecifiedarp. L^{2} and Lz commute — both can be simultaneously specified.
Example LA.2Diagonalizing the Inertia Tensor

ArigidbodyhasinertiatensorI(symmetric3×3matrix).FindtheprincipalaxesA rigid body has inertia tensor I (symmetric 3\times3 matrix). Find the principal axes.

Problem:In an arbitraryframe,theinertiatensorIhasoffdiagonalelements(productsofinertia).L=arbitrary frame, the inertia tensor I has off-diagonal elements (products of inertia). L =sLnotparalleltoωingenerals L not parallel to \omega in general.
Spectral theorem:SinceIissymmetric(realHermitian),ithas3realeigenvaluesI1,I2,I3(principalmomSince I is symmetric (real Hermitian), it has 3 real eigenvalues I_{1}, I_{2}, I_{3} (principal moments) and orthogonal eigenvectors (principal axes).
In principal frame:I = diag(I_{1}, I_{2}, I_{3}). Euler's equations simplify: I_{1}\omegȧ_{1} = (I_{2}-I_{3})\omega_{2}\omega_{3}, etc.
Symmetry:Abodywith3foldorhighersymmetryaxishastwoequalprincipalmomentsaxisymmetricA body with 3-fold or higher symmetry axis has two equal principal moments \to axisymmetric top.SphericalsymmetryI1=I2=I3noprecessiontop. Spherical symmetry \to I_{1} = I_{2} = I_{3} \to no precession
Definition LA.3Common Traps
  • Basis vectors are choices: physical vectors and states are independent of representation.
  • Eigenvectors can be degenerate: a degenerate eigenspace needs a basis choice.
  • Hermitian operators have real eigenvalues: this is why they represent observables in quantum mechanics.
  • Matrix multiplication is not commutative: operator order can change the result.
Exercises — LA.1–LA.4 Linear Algebra for Physics
1.
FindthelargereigenvalueofthematrixA=[[4,1],[2,3Find the larger eigenvalue of the matrix A = [[4,1],[2,3.
Straightforward
2.
FindtherankofthematrixA=[[1,2,3],[2,4,6Find the rank of the matrix A = [[1,2,3],[2,4,6.
Straightforward
3.ApplyGramSchmidtorthogonalizationtothevectorsv1=(1,1,0),v2=(1,0,1),v3=(0,1,1)toApply Gram-Schmidt orthogonalization to the vectors v_{1}=(1,1,0), v_{2}=(1,0,1), v_{3}=(0,1,1) to produce an orthonormal basis.
Straightforward
4.Write the matrix representations of the harmonic oscillator Hamiltonian H and position operator x in the energy eigenstate basis {|0⟩, |1⟩, |2⟩, |3⟩}. What selection rule does this imply?
Intermediate
5.UsetheCayleyHamiltontheoremtofindarecursionforAnwhereA=[[1,1],[1,0]].HowdoUse the Cayley-Hamilton theorem to find a recursion for A^{n} where A = [[1,1],[1,0]]. How does this relate to Fibonacci numbers?
Intermediate
6.DiscussthepositionandmomentumoperatorsasoperatorsonL2(R).WhydotheynothavenoDiscuss the position and momentum operators as operators on L^{2}(ℝ). Why do they not have normalizable eigenstates? How does the continuous spectrum relate to the uncertainty principle?
Challenging
Key Takeaways
  • Innerproductu,v:CauchySchwarzgivesuvu,vInner product ⟨u,v⟩: Cauchy-Schwarz gives \|u\|\|v\| \ge |⟨u,v⟩|.
  • Eigenvaluesfromdet(A\lambdaI)=0;eigenvectorsspantheeigenspaceEigenvalues from det(A-\lambdaI)=0; eigenvectors span the eigenspace.
  • Hermitian:A=A.Eigenvaluesreal,eigenvectorsorthogonal,spectraltheoremholdsHermitian: A=A^\dagger. Eigenvalues real, eigenvectors orthogonal, spectral theorem holds.
  • Unitary:U\daggerU=I.Eigenvaluesonunitcircle.Timeevolutionisunitary(probabilityconservUnitary: U^\daggerU=I. Eigenvalues on unit circle. Time evolution is unitary (probability conserved).
  • SVD:A=U\sumV.FoundationofPCA,lowrankapproximation,andleastsquaresSVD: A = U\sumV^\dagger. Foundation of PCA, low-rank approximation, and least squares.
  • Commuting Hermitian operators share eigenbases — simultaneously measurable observables.