Rank, Nullity, and the Fundamental Theorem of Linear Maps

A linear map T:V→WT:V\to W does only two things to the space it acts on: it crushes some directions to nothing, and it carries the rest through faithfully. The crushed directions form the kernel ker⁡T={v:Tv=0}\ker T=\{v: Tv=0\}; the directions that survive fill out the image im⁡T={Tv:v∈V}\operatorname{im}T=\{Tv: v\in V\}. Both are honest subspaces — the kernel inside the domain VV, the image inside the codomain WW — and that alone is a consequence of linearity. The dimension of the kernel is the nullity, the dimension of the image is the rank, and the whole subject of this article is a single, exact accounting identity between them. The Fundamental Theorem of Linear Maps — the rank-nullity theorem — says these two numbers are not free to roam: their sum is pinned at dim⁡V\dim V. The domain has a fixed budget of dimensions, and every one is spent exactly once, either annihilated into the kernel or transmitted into the image. The statement is a conservation law, and like the conservation laws of physics it converts a hard structural question into simple bookkeeping. It needs only that VV be finite-dimensional; WW may be as large as you like, even infinite-dimensional. That one identity does an astonishing amount of work. It decides when a linear system has solutions and how many; it proves that a square operator is injective exactly when it is surjective exactly when it is invertible; it measures the degeneracy of an energy level as the nullity of H−E 1H-E\,\mathbb{1}. It is the first structural theorem of linear algebra that a working physicist uses every day, almost always without naming it.

Key takeaways
  • A linear map T:V→WT:V\to W splits the domain into a part it kills (the kernel, dimension = nullity) and a part it transmits (the image, dimension = rank).
  • The Fundamental Theorem of Linear Maps is the conservation law dim⁡V=rank⁡T+null⁡T\dim V=\operatorname{rank}T+\operatorname{null}T — every domain dimension is spent exactly once.
  • The proof extends a basis of the kernel to a basis of VV and shows the images of the added vectors form a basis of im⁡T\operatorname{im}T; only dim⁡V<∞\dim V<\infty is needed.
  • In matrix form rank⁡A+dim⁡(null⁡A)=n\operatorname{rank}A+\dim(\operatorname{null}A)=n for A∈Fm×nA\in\mathbb{F}^{m\times n}, and remarkably column rank always equals row rank.
  • Immediate physics payoffs: a square operator is injective iff surjective iff invertible; a system Ax=bA\mathbf{x}=\mathbf{b} is solvable iff rank⁡[A ∣ b]=rank⁡A\operatorname{rank}[A\,|\,\mathbf{b}]=\operatorname{rank}A; the geometric multiplicity of an eigenvalue λ\lambda is n−rank⁡(A−λ1)n-\operatorname{rank}(A-\lambda\mathbb{1}).
  • Worked cleanly: a 3×43\times4 matrix of rank 22 has nullity 22 (so 2+2=42+2=4), and the defective matrix with eigenvalue 22 of algebraic multiplicity 22 has geometric multiplicity only 11.
  • Key idea: Every dimension of the domain is accounted for exactly once: dim⁡V=rank⁡T+null⁡T\dim V=\operatorname{rank}T+\operatorname{null}T. The map either kills a direction (kernel) or transmits it to the image — never both, never neither.

Where this lives today

Rank and nullity have quietly become the load-bearing numbers of computational science. Every time an engineer runs a singular value decomposition, the count of singular values above the noise floor is a numerical rank, and the near-null space it exposes is what principal component analysis, image compression, and recommender systems all exploit — a data matrix of nominal rank 10001000 that is "really" rank 2020 is the premise of low-rank approximation and compressed sensing. In control engineering, Kalman's controllability and observability tests are literally rank conditions: a system is controllable exactly when a certain matrix has full rank, i.e. nullity zero. In first-principles materials physics, phonon calculations enforce the acoustic sum rule precisely so the dynamical matrix has the correct three-dimensional null space of rigid translations. And in inverse problems — X-ray tomography, seismic imaging, gravitational-wave parameter estimation — a rank-deficient forward operator is the signature of information you did not measure, its nullity naming exactly the parameter combinations the data can never constrain. The theorem that began as bookkeeping for Victorian determinants is now the diagnostic that tells a scientist which questions the data can answer.

The physical picture

Think of TT as a pipeline that swallows an nn-dimensional input and delivers whatever comes out the far end. Some input directions are squeezed flat and emerge as zero — that information is gone, and those directions make up the kernel. Every direction that is not flattened has to appear at the output as a genuinely new, independent direction of the image. There is no third fate, so the input dimension is conserved: flattened plus transmitted equals the total. The first figure draws this ledger with real numbers — a map from a 55-dimensional domain whose kernel is 22-dimensional must have a 33-dimensional image, because 5=3+25=3+2 and there is nowhere else for the dimensions to go.

Figure
A linear map from a 5-dimensional domain V to W. A 2-dimensional kernel (red) is crushed to the single zero vector, while the complementary 3 dimensions (blue) are carried faithfully onto a 3-dimensional image. The ledger closes exactly: dim V = rank T + nullity T, that is 5 = 3 + 2.
Every domain dimension is spent once: crushed into the kernel, or carried into the image domain V dim = 5 transmitted part dim = 3 kernel nullity = 2 codomain W image rank = 3 0 T dim V = rank T + nullity T 5 = 3 + 2

The cleanest concrete instance is a projection. Drop every vector of R3\mathbb{R}^3 straight down onto the horizontal plane: the shadow is its image, and the vertical direction you project along is squashed to zero. The image is the whole plane, so the rank is 22; the kernel is the single vertical line, so the nullity is 11; and 2+1=32+1=3, the domain dimension. The second figure makes the trade-off vivid: you cannot enlarge the image without shrinking the kernel, because they are two halves of one fixed budget.

Figure
Orthogonal projection of ℝ³ onto the horizontal plane. A vector v drops to its shadow P v in the plane (the image, rank 2), while the vertical direction it is projected along is sent to zero (the kernel, nullity 1). The budget closes as 2 + 1 = 3.
Projection onto a plane: image = plane (rank 2), kernel = the vertical axis (nullity 1) image rank = 2 kernel nullity = 1 v P v 0 rank + nullity = 2 + 1 = 3 = dim ℝ³

Choosing a basis adapted to TT makes the theorem visible at a glance: take basis vectors inside the kernel, then add just enough more to fill out VV. The kernel vectors map to zero and the added ones map to an independent set spanning the image, so the added vectors number dim⁡V\dim V minus the nullity and their images number the rank. This is also why solving a square system is all-or-nothing: if no direction is crushed, none is missing from the image either, so existence and uniqueness arrive together.

The core equations

Everything begins with two sets attached to the map: the vectors it sends to zero, and the vectors it can produce. Linearity guarantees both are subspaces, so it makes sense to speak of their dimensions.

ker⁡T={v∈V:Tv=0},im⁡T={Tv:v∈V}(1)\ker T=\{v\in V: Tv=0\},\qquad \operatorname{im}T=\{Tv: v\in V\} \qquad (1)

where: TT — a linear map V→WV\to W obeying T(αu+βv)=α Tu+β TvT(\alpha u+\beta v)=\alpha\,Tu+\beta\,Tv; ker⁡T\ker T — the kernel or null space, a subspace of the domain VV; im⁡T\operatorname{im}T — the image or range, a subspace of the codomain WW. Reading: the kernel collects what is destroyed, the image what is reachable. Both are subspaces by linearity — if Tu=Tv=0Tu=Tv=0 then T(αu+βv)=0T(\alpha u+\beta v)=0, and images are closed the same way. The two numbers that matter are the dimensions of these subspaces. They are intrinsic to the map, unchanged by any choice of basis or by any similarity transform, which is exactly why they carry physical meaning.

rank⁡T=dim⁡(im⁡T),null⁡T=dim⁡(ker⁡T)(2)\operatorname{rank}T=\dim(\operatorname{im}T),\qquad \operatorname{null}T=\dim(\ker T) \qquad (2)

where: rank⁡T\operatorname{rank}T — the number of independent output directions the map produces; null⁡T\operatorname{null}T — the number of independent input directions it destroys. Reading: rank measures how much gets through, nullity how much is lost. The theorem itself ties them together. It is a single equation, and it is exact for any finite-dimensional domain regardless of how large the codomain is.

dim⁡V=rank⁡T+null⁡T(3)\dim V=\operatorname{rank}T+\operatorname{null}T \qquad (3)

where: dim⁡V\dim V — the dimension of the finite-dimensional domain. Reading: domain dimension = transmitted + annihilated. This is the conservation law; the derivation below turns "spend each dimension once" into a rigorous basis count. Fixing bases turns TT into a matrix AA, and the theorem becomes a statement about columns. The number of columns is the domain dimension; the pivots count the rank; the leftover free columns count the nullity.

rank⁡A+dim⁡(null⁡A)=n,A∈Fm×n(4)\operatorname{rank}A+\dim(\operatorname{null}A)=n,\qquad A\in\mathbb{F}^{m\times n} \qquad (4)

where: AA — the m×nm\times n matrix of T:Fn→FmT:\mathbb{F}^n\to\mathbb{F}^m; nn — number of columns, equal to dim⁡V\dim V; rank⁡A\operatorname{rank}A — number of pivot columns, equal to dim⁡(column space)\dim(\text{column space}). Reading: pivots + free columns = total columns, so dim⁡(null⁡A)=n−rank⁡A\dim(\operatorname{null}A)=n-\operatorname{rank}A — each non-pivot column is one free parameter of Ax=0A\mathbf{x}=\mathbf{0}. A companion fact makes the word "rank" unambiguous: the number of independent columns equals the number of independent rows. Row reduction produces the same pivot count whether you read it down the columns or across the rows.

rank⁡A=rank⁡AT(column rank=row rank)(5)\operatorname{rank}A=\operatorname{rank}A^{\mathsf T}\quad(\text{column rank}=\text{row rank}) \qquad (5)

where: ATA^{\mathsf T} — the transpose, whose column space is the row space of AA. Reading: the rr pivots of the echelon form count both independent rows and independent columns, so the two ranks coincide — a nontrivial theorem, not a definition. Finally, rank decides solvability. Appending the right-hand side b\mathbf{b} as an extra column must not create a new independent direction, or there is no solution at all.

Ax=b is solvable  ⟺  b∈im⁡A  ⟺  rank⁡[A ∣ b]=rank⁡A(6)A\mathbf{x}=\mathbf{b}\ \text{is solvable}\iff \mathbf{b}\in\operatorname{im}A\iff \operatorname{rank}[A\,|\,\mathbf{b}]=\operatorname{rank}A \qquad (6)

where: [A ∣ b][A\,|\,\mathbf{b}] — the augmented matrix. Reading: consistency is a rank condition, and when it holds the solution set is an affine copy of ker⁡A\ker A of dimension n−rank⁡An-\operatorname{rank}A — one point when the nullity is zero, a whole flat when it is positive.

Derivation AProof of the fundamental theorem of linear maps

The goal is to prove dim⁡V=rank⁡T+null⁡T\dim V=\operatorname{rank}T+\operatorname{null}T from scratch, using nothing but the definition of a basis and the linearity of TT. The strategy is the standard and beautiful one: start with a basis of the kernel, extend it to a basis of the whole domain, and then show that the images of exactly the added vectors form a basis of the image. Counting the vectors in each basis then delivers the identity with no arithmetic to spare. Throughout, let T:V→WT:V\to W with dim⁡V=n\dim V=n finite and null⁡T=k\operatorname{null}T=k.

Step 1. Choose a basis u1,…,uku_1,\dots,u_k of the kernel, so that by definition null⁡T=k\operatorname{null}T=k; if the kernel is trivial we take k=0k=0 and start from the empty set. Why: the kernel is where the map does its destroying, and pinning down a basis for it lets us separate the destroyed directions cleanly from everything else.

Step 2. Extend that independent set to a basis {u1,…,uk, v1,…,vn−k}\{u_1,\dots,u_k,\,v_1,\dots,v_{n-k}\} of VV; the basis-extension theorem guarantees this is always possible, and the number of added vectors is forced to be n−kn-k because a basis of an nn-dimensional space has exactly nn elements. Why: the added vjv_j are the candidate carriers of the image — they are precisely the directions the map does not kill.

Step 3. Show the images Tv1,…,Tvn−kTv_1,\dots,Tv_{n-k} span im⁡T\operatorname{im}T. Any output is TvTv for some v=∑iaiui+∑jbjvjv=\sum_i a_i u_i+\sum_j b_j v_j, and since Tui=0Tu_i=0 the kernel part drops out, leaving Tv=∑jbj TvjTv=\sum_j b_j\,Tv_j. Why: the kernel basis vectors contribute nothing to any output, so whatever the image contains is already built from the images of the added vectors alone.

Step 4. Show those images are linearly independent. Suppose ∑jcj Tvj=0\sum_j c_j\,Tv_j=0; then T ⁣(∑jcjvj)=0T\!\left(\sum_j c_j v_j\right)=0, so ∑jcjvj∈ker⁡T\sum_j c_j v_j\in\ker T and can be written ∑idiui\sum_i d_i u_i. Rearranged, ∑jcjvj−∑idiui=0\sum_j c_j v_j-\sum_i d_i u_i=0, a linear relation among the basis of VV, which forces every coefficient — in particular every cjc_j — to vanish. Why: independence is what upgrades "spanning set" to "basis", and it rests entirely on the ui,vju_i,v_j together being a basis of VV.

Step 5. Conclude. The TvjTv_j form a basis of the image, so rank⁡T=dim⁡(im⁡T)=n−k\operatorname{rank}T=\dim(\operatorname{im}T)=n-k, and adding back the nullity gives rank⁡T+null⁡T=(n−k)+k=n=dim⁡V\operatorname{rank}T+\operatorname{null}T=(n-k)+k=n=\dim V. Why: this is the whole theorem, and it used only dim⁡V<∞\dim V<\infty — the codomain was never assumed finite-dimensional, which is why it survives into the operator theory of quantum mechanics.

Go deeper · Derivation B · 1 minInjective equals surjective equals invertible for a square operator
Derivation BInjective equals surjective equals invertible for a square operator

The goal is to prove one of the most-used facts in applied mathematics: when a linear operator maps a finite-dimensional space to itself, being injective, being surjective, and being invertible are three names for the same condition, each equivalent to the nullity being zero. This equivalence fails badly in infinite dimensions, so the proof leans entirely on the rank-nullity theorem, which is what makes finiteness decisive. Let T:V→VT:V\to V be linear with dim⁡V=n<∞\dim V=n<\infty.

Step 1. Suppose TT is injective, i.e. null⁡T=0\operatorname{null}T=0. Rank-nullity gives rank⁡T=n−0=n\operatorname{rank}T=n-0=n, so im⁡T\operatorname{im}T is an nn-dimensional subspace of the nn-dimensional space VV, hence all of VV: TT is surjective. Why: a subspace whose dimension equals that of the whole space cannot be a proper subspace, so full rank forces the image to fill the codomain.

Step 2. Conversely suppose TT is surjective, i.e. rank⁡T=n\operatorname{rank}T=n. Rank-nullity gives null⁡T=n−n=0\operatorname{null}T=n-n=0, so the kernel is trivial and TT is injective. Why: the same conservation law run backwards — filling the image completely leaves no dimensions over to be destroyed.

Step 3. Each property implies the other, so together they make TT a bijection, and a bijective linear map has a linear inverse T−1T^{-1}. In matrix terms this is exactly the condition det⁡A≠0\det A\neq0. Why: existence and uniqueness of solutions to Tx=bT\mathbf{x}=\mathbf{b} are governed by surjectivity and injectivity respectively, and here they stand or fall together.

Step 4. Note where finiteness is essential: the right-shift operator S(x1,x2,… )=(0,x1,x2,… )S(x_1,x_2,\dots)=(0,x_1,x_2,\dots) on ℓ2\ell^2 is injective but not surjective, and the left shift is surjective but not injective. Why: with infinitely many dimensions the "budget" is inexhaustible, so a map can gain or lose room without contradiction — the equivalence is a purely finite-dimensional gift of rank-nullity.

Go deeper · Derivation C · 1 minCounting the solutions of a linear system
Derivation CCounting the solutions of a linear system

The goal is to extract, from rank-nullity alone, the complete decision tree for the equation Ax=bA\mathbf{x}=\mathbf{b}: when a solution exists, when it is unique, and when it is not, how large the family of solutions is. The answer is entirely dimensional — existence is a rank condition on the augmented matrix, and the size of the solution set is the nullity — so a physicist can read off the structure of a linear problem without ever computing a solution. Let A∈Fm×nA\in\mathbb{F}^{m\times n} and write r=rank⁡Ar=\operatorname{rank}A.

Step 1. Existence. A solution exists if and only if b\mathbf{b} lies in the image of AA, equivalently rank⁡[A ∣ b]=rank⁡A=r\operatorname{rank}[A\,|\,\mathbf{b}]=\operatorname{rank}A=r; if appending b\mathbf{b} raises the rank, the system is inconsistent and has no solution. Why: AxA\mathbf{x} is by definition a combination of the columns of AA, so it can equal b\mathbf{b} only when b\mathbf{b} adds no new column direction.

Step 2. Structure. When consistent, the general solution is x=xp+xh\mathbf{x}=\mathbf{x}_p+\mathbf{x}_h, one particular solution xp\mathbf{x}_p with Axp=bA\mathbf{x}_p=\mathbf{b} plus any homogeneous solution xh\mathbf{x}_h with Axh=0A\mathbf{x}_h=\mathbf{0}. Why: the difference of any two solutions lies in the kernel, so the solution set is a single translate of ker⁡A\ker A — an affine flat, not a subspace unless b=0\mathbf{b}=\mathbf{0}.

Step 3. Size. By rank-nullity the homogeneous solutions fill ker⁡A\ker A, of dimension n−rn-r, so the solution set is an affine space of dimension n−rn-r. Why: the nullity directly counts the free parameters, so it is literally the number of independent knobs left after the equation is imposed.

Step 4. Special cases. If r=nr=n (full column rank) the kernel is trivial and there is at most one solution; if in addition r=m=nr=m=n the matrix is square and invertible, and there is exactly one solution x=A−1b\mathbf{x}=A^{-1}\mathbf{b} for every b\mathbf{b}. Why: these are the invertible and overdetermined-consistent limits, and they are exactly the cases the square-operator equivalence above singles out.

Limiting cases

The single identity specializes into sharply different regimes depending on the shape of the map and how much of its domain survives:

CaseResultMeaning
Injective mapker⁡T={0}  ⟺  null⁡T=0  ⟺  rank⁡T=dim⁡V\ker T=\{0\}\iff \operatorname{null}T=0\iff \operatorname{rank}T=\dim VInjectivity is exactly a trivial kernel; the rank then equals the full domain dimension, which requires dim⁡W≥dim⁡V\dim W\geq\dim V.
Surjective mapim⁡T=W  ⟺  rank⁡T=dim⁡W\operatorname{im}T=W\iff \operatorname{rank}T=\dim WOnto means the image fills the codomain; combined with rank-nullity this forces null⁡T=dim⁡V−dim⁡W≥0\operatorname{null}T=\dim V-\dim W\geq0, so dim⁡V≥dim⁡W\dim V\geq\dim W.
Square operator V→VV\to Vnull⁡T=0  ⟺  rank⁡T=dim⁡V  ⟺  T invertible\operatorname{null}T=0\iff \operatorname{rank}T=\dim V\iff T\text{ invertible}When dim⁡V=dim⁡W\dim V=\dim W, injective, surjective, and bijective coincide — the workhorse behind invertibility tests and nonzero determinants.
Rank deficiency / degeneracyrank⁡T=r<dim⁡V⇒null⁡T=dim⁡V−r>0\operatorname{rank}T=r<\dim V\Rightarrow \operatorname{null}T=\dim V-r>0A nontrivial kernel means non-unique solutions and, for A−λ1A-\lambda\mathbb{1}, a degenerate (multi-dimensional) eigenspace.
Zero maprank⁡T=0, null⁡T=dim⁡V\operatorname{rank}T=0,\ \operatorname{null}T=\dim VThe entire domain is crushed to {0}\{0\}; the kernel is everything — the extreme opposite of an isomorphism.

Worked examples

1 · Rank, nullity, and the null space of a matrix (fully worked)

Find the rank, the nullity, and an explicit basis of the null space of

A=(123424680111),A∈R3×4.A=\begin{pmatrix} 1 & 2 & 3 & 4\\ 2 & 4 & 6 & 8\\ 0 & 1 & 1 & 1 \end{pmatrix},\qquad A\in\mathbb{R}^{3\times4}.

Step 1. Row reduce. The second row is exactly twice the first, so R2→R2−2R1=(0,0,0,0)R_2\to R_2-2R_1=(0,0,0,0) vanishes. Swapping the zero row to the bottom leaves two nonzero rows, (1,2,3,4)(1,2,3,4) with a pivot in column 1 and (0,1,1,1)(0,1,1,1) with a pivot in column 2. There are two pivots, so rank⁡A=2\operatorname{rank}A=2.

Step 2. Apply the theorem. The domain is R4\mathbb{R}^4, so n=4n=4 and rank-nullity gives dim⁡(null⁡A)=n−rank⁡A=4−2=2\dim(\operatorname{null}A)=n-\operatorname{rank}A=4-2=\mathbf{2}. Columns 3 and 4 are the non-pivot (free) columns, matching the two free parameters we are about to find.

Step 3. Solve Ax=0A\mathbf{x}=\mathbf{0}. With free variables x3=sx_3=s and x4=tx_4=t, the pivot row (0,1,1,1)(0,1,1,1) gives x2=−s−tx_2=-s-t, and the pivot row (1,2,3,4)(1,2,3,4) gives x1=−2x2−3x3−4x4=−s−2tx_1=-2x_2-3x_3-4x_4=-s-2t. Hence

x=s(−1−110)+t(−2−101).\mathbf{x}=s\begin{pmatrix}-1\\-1\\1\\0\end{pmatrix}+t\begin{pmatrix}-2\\-1\\0\\1\end{pmatrix}.

Substituting back confirms A(−1,−1,1,0)T=0A(-1,-1,1,0)^{\mathsf T}=\mathbf{0} and A(−2,−1,0,1)T=0A(-2,-1,0,1)^{\mathsf T}=\mathbf{0}, so these two vectors are a basis of the null space, of dimension 22.

Step 4. Cross-check. Column rank =2=2 (columns 1 and 2 are independent), row rank =2=2 (rows 1 and 3 are independent, row 2 dependent); they agree, illustrating equation (5). The ledger closes: rank⁡A+dim⁡(null⁡A)=2+2=4=n\operatorname{rank}A+\dim(\operatorname{null}A)=2+2=4=n. Rank 2, nullity 2, and 2 + 2 = 4.

2 · Geometric multiplicity and a defective matrix (fully worked)

For B=(210020003)B=\begin{pmatrix}2&1&0\\0&2&0\\0&0&3\end{pmatrix} find the geometric multiplicity of the eigenvalue λ=2\lambda=2 and decide whether BB is diagonalizable.

Step 1. The characteristic polynomial is det⁡(B−λ1)=(2−λ)2(3−λ)\det(B-\lambda\mathbb{1})=(2-\lambda)^2(3-\lambda), so λ=2\lambda=2 has algebraic multiplicity 2 and λ=3\lambda=3 has algebraic multiplicity 1.

Step 2. Form B−21=(010000001)B-2\mathbb{1}=\begin{pmatrix}0&1&0\\0&0&0\\0&0&1\end{pmatrix}. Its nonzero rows (0,1,0)(0,1,0) and (0,0,1)(0,0,1) are independent, so rank⁡(B−21)=2\operatorname{rank}(B-2\mathbb{1})=2.

Step 3. The geometric multiplicity is the nullity of B−21B-2\mathbb{1}, and rank-nullity on the 3×33\times3 operator gives dim⁡ker⁡(B−21)=3−rank⁡(B−21)=3−2=1\dim\ker(B-2\mathbb{1})=3-\operatorname{rank}(B-2\mathbb{1})=3-2=\mathbf{1}. Solving (B−21)x=0(B-2\mathbb{1})\mathbf{x}=\mathbf{0} gives x2=0, x3=0, x1x_2=0,\ x_3=0,\ x_1 free, so the single eigenvector is (1,0,0)T(1,0,0)^{\mathsf T}.

Step 4. Compare. The geometric multiplicity 11 is strictly less than the algebraic multiplicity 22, so BB has too few eigenvectors to form a basis: it is defective and not diagonalizable, and its structure is captured instead by the Jordan canonical form. The whole diagnosis was a nullity computation. Geometric multiplicity 1 < algebraic multiplicity 2, so B is defective.

3 · Counting molecular vibrational modes

How many independent vibrational modes does a nonlinear water molecule (H₂O, N=3N=3 atoms) have, and how many for linear carbon dioxide (CO₂, N=3N=3)?

Step 1. A molecule of NN atoms has 3N3N Cartesian coordinates, so its configuration space is 3N3N-dimensional. For both molecules 3N=93N=9.

Step 2. The zero-frequency modes — rigid translations and rotations of the whole molecule — form the kernel of the mass-weighted dynamical (Hessian) matrix, since a rigid displacement costs no restoring energy. A nonlinear molecule has 33 translations and 33 rotations, nullity 66; a linear molecule has 33 translations and only 22 rotations, nullity 55.

Step 3. By rank-nullity the vibrational modes are the rank of the dynamical matrix, 3N3N minus the nullity: water gives 9−6=39-6=\mathbf{3} vibrational modes, and CO₂ gives 9−5=49-5=\mathbf{4}. H₂O: 3 modes; CO₂: 4 modes — the standard 3N − 6 and 3N − 5 counts.

Physics of the portrait

The rank-nullity theorem carries no single discoverer's name because it grew up alongside the language of matrices themselves. Its vocabulary, though, is almost entirely the invention of one flamboyant Victorian.

Origins · 1 minJames Joseph Sylvester and how this was discovered
JJ
James Joseph Sylvester
1814–1897 · matrix algebra and invariant theory; coined the words "matrix", "minor", and "nullity"

Sylvester was one of the great word-coiners of mathematics, and much of the language of this article is his. Barred for decades from a professorship in England because he was Jewish, he worked as an actuary and a lawyer before holding chairs at Virginia, Johns Hopkins — where he founded the American Journal of Mathematics in 1878 — and finally Oxford. In an 1850 paper he introduced the term "matrix" for a rectangular array from which determinants (his "minors") could be cut, and in his 1884 lectures he introduced "nullity" for the dimension of what we now call the kernel, along with his "law of nullity" bounding the nullity of a product. His collaborator Arthur Cayley built the abstract algebra of matrices around these ideas, while Georg Frobenius gave the rank its rigorous modern theory in 1878, establishing the equality of row and column rank. The clean statement in terms of abstract linear maps — dim⁡V=rank⁡T+null⁡T\dim V=\operatorname{rank}T+\operatorname{null}T — belongs to the twentieth-century vector-space viewpoint of Peano, Noether, and van der Waerden, but the two nouns in its name are Sylvester's gift.

History
1844
— Grassmann's Ausdehnungslehre founds the ideas of linear independence and dimension
1850
— Sylvester coins the term "matrix" for a rectangular array of numbers
1858
— Cayley's memoir builds the abstract algebra of matrices, including the inverse
1878
— Frobenius gives the rank of a matrix its rigorous theory and proves row rank equals column rank
1884
— Sylvester introduces "nullity" and states his law of nullity for products
1930
— van der Waerden's Moderne Algebra frames rank plus nullity for abstract linear maps

Going further

Texts. Axler, Linear Algebra Done Right, 3rd ed., Ch. 3 (the Fundamental Theorem of Linear Maps) · Stone & Goldbart, Mathematics for Physics, Ch. 1 (Vector Spaces and Linear Operators) · Arfken, Weber & Harris, Mathematical Methods for Physicists, 7th ed., Ch. 2–3 · Riley, Hobson & Bence, Mathematical Methods for Physics and Engineering, 3rd ed., Ch. 8 · Strang, Introduction to Linear Algebra, Ch. 3–4 (the four fundamental subspaces) · Hassani, Mathematical Physics, Ch. 2 · Nakahara, Geometry, Topology and Physics, §2 (linear maps and their invariants). Key results to keep in mind: the basis-extension proof yields rank⁡=n−k\operatorname{rank}=n-k using only dim⁡V<∞\dim V<\infty; the worked 3×43\times4 matrix has rank 22, nullity 22, and the stated null-space basis; the defective matrix BB has geometric multiplicity 11 versus algebraic multiplicity 22; geometric multiplicity =n−rank⁡(A−λ1)=n-\operatorname{rank}(A-\lambda\mathbb{1}); and the molecular mode counts are 3N−63N-6 (nonlinear) and 3N−53N-5 (linear). No nonstandard numerical values are used anywhere in the article.

Builds on. Vector Spaces and Subspaces · Linear Independence, Basis, and Dimension · Linear Operators and Matrix Representations. Leads to. Eigenvalues, Eigenvectors, and the Characteristic Equation · Diagonalization and the Spectral Decomposition · Jordan Canonical Form and Defective Matrices · Adjoint and Hermitian Operators.

Ask the tutor: Ask me to row-reduce any matrix you like and read off its rank, nullity, null-space basis, and exactly how many solutions Ax=bA\mathbf{x}=\mathbf{b} has.

Why does the theorem need the domain to be finite-dimensional but not the codomain?

The proof counts basis vectors of VV, splitting them into a kernel basis and the added vectors whose images span the range. That count only makes sense if VV has finitely many basis vectors to distribute — the whole argument is a finite partition of the domain's dimensions. The codomain WW never gets counted; it only has to be large enough to contain the image. This is why the theorem survives into quantum mechanics, where operators act on finite-dimensional spin spaces mapping into infinite-dimensional function spaces. In genuinely infinite-dimensional domains the identity fails, and one needs the Fredholm index instead.

What is the difference between geometric and algebraic multiplicity?

The algebraic multiplicity of an eigenvalue λ\lambda is how many times (λ−x)(\lambda-x) divides the characteristic polynomial; the geometric multiplicity is dim⁡ker⁡(A−λ1)\dim\ker(A-\lambda\mathbb{1}), the number of independent eigenvectors, computed by rank-nullity as n−rank⁡(A−λ1)n-\operatorname{rank}(A-\lambda\mathbb{1}). Geometric multiplicity is always at least 11 and never exceeds the algebraic multiplicity. When they are equal for every eigenvalue the matrix is diagonalizable; when the geometric falls short — as in Worked Example 2, where a double eigenvalue supplies only one eigenvector — the matrix is defective and requires the Jordan canonical form. In physics the geometric multiplicity is the degeneracy of an energy level, so this gap is precisely where the tidy spectral picture breaks down.

Is the rank of a matrix the same as the number of its nonzero rows?

Only after you row-reduce it. In a general matrix the rows may be linearly dependent even though none is literally zero — the second row of Worked Example 1, (2,4,6,8)(2,4,6,8), is nonzero yet contributes nothing to the rank because it is twice the first. The rank is the number of pivots in the row echelon form, equivalently the number of independent rows, equivalently the number of independent columns. Counting visibly nonzero rows in the original matrix overcounts whenever there is hidden dependence. The safe recipe is always the same: reduce to echelon form and count the pivots, which simultaneously gives you the rank and, by subtraction from the number of columns, the nullity.

How does rank-nullity relate to solving A x = b?

It supplies the complete answer to both existence and uniqueness. Existence is the rank condition of equation (6): a solution exists iff appending b\mathbf{b} leaves the rank unchanged, i.e. b\mathbf{b} is in the column space. Uniqueness is the nullity: the solution set, when nonempty, is a translate of the null space, so it has dimension n−rank⁡An-\operatorname{rank}A. Zero nullity means at most one solution; positive nullity means an infinite affine family. This is why full column rank guarantees uniqueness and full row rank guarantees existence for every b\mathbf{b}, and why a square matrix with both — nonzero determinant — gives exactly one solution always. The theorem turns "how many solutions?" into two dimension counts.

Why is column rank always equal to row rank?

It is a genuine theorem, not a definition, and it can feel surprising because rows and columns look so different. One clean argument: row-reduce AA to echelon form. Row operations never change the row space, so the number of nonzero rows (the pivots) is the row rank; but those same pivots sit in specific columns, and the corresponding columns of the original AA are independent and span its column space, making the pivot count the column rank too. So a single integer, the number of pivots, is both ranks at once — equivalently rank⁡A=rank⁡AT\operatorname{rank}A=\operatorname{rank}A^{\mathsf T}. The consequence is that "rank" is unambiguous, which is quietly essential to the whole theory.

Do the kernel and image ever have any overlap?

They live in different spaces — the kernel in the domain VV, the image in the codomain WW — so in general asking whether they overlap is a category error. The confusion arises for an operator T:V→VT:V\to V mapping a space to itself, where both are subspaces of the same VV and can genuinely intersect. A pure projection has ker⁡T\ker T and im⁡T\operatorname{im}T complementary, meeting only at 00; but a nilpotent shift (x,y,z)↦(y,z,0)(x,y,z)\mapsto(y,z,0) has vectors that are both in the image and in the kernel. So for self-maps, rank plus nullity still equals dim⁡V\dim V, yet the two subspaces need not be complementary — that holds only for special operators such as projections and, with an inner product, normal operators.

Before this

Where this leads