Skip to content

Projectors

A projector is a linear operator PP that is idempotent:

P2=P.P^2=P.

Applying it once separates a vector into a retained component and a discarded component; applying it again changes nothing. This simple algebraic condition organizes subspace decompositions, constrained approximation, spectral resolutions, and ideal quantum alternatives.

This page develops the linear algebra. The interpretation of projectors as yes–no quantum propositions belongs to Projectors in Core Formalism, and probabilities and conditional states belong to Projective Measurement.

Let P:V→VP:V\to V satisfy P2=PP^2=P. Its two distinguished subspaces are

Ran⁡P={Pv:v∈V},ker⁡P={v∈V:Pv=0}.\begin{aligned} \operatorname{Ran}P &=\{P v:v\in V\},\\ \ker P &=\{v\in V:Pv=0\}. \end{aligned}

The operator is the identity on its range. Indeed, if u=Pvu=Pv, then Pu=P2v=Pv=uPu=P^2v=Pv=u. Every v∈Vv\in V has the decomposition

v=Pv+(I−P)v,Pv∈Ran⁡P,(I−P)v∈ker⁡P.\begin{aligned} v &=Pv+(I-P)v,\\ Pv &\in\operatorname{Ran}P,\\ (I-P)v &\in\ker P. \end{aligned}

The intersection of these subspaces is trivial: if w=Pvw=Pv and Pw=0Pw=0, then w=Pw=0w=Pw=0. Consequently,

V=Ran⁡P⊕ker⁡P.V=\operatorname{Ran}P\oplus\ker P.

Conversely, a direct-sum decomposition V=M⊕NV=M\oplus N defines a unique projector onto MM along NN:

P(m+n)=m,m∈M,n∈N.P(m+n)=m, \qquad m\in M,\quad n\in N.

Thus a projector remembers two subspaces, not merely its range. Different choices of the complementary subspace NN give different projectors with the same range. See Direct Sums for the underlying decomposition.

If Pv=λvPv=\lambda v for a nonzero vector vv, idempotence gives

λ2v=P2v=Pv=λv.\lambda^2v=P^2v=Pv=\lambda v.

Therefore every eigenvalue is either 00 or 11. More strongly, the minimal polynomial of PP divides x(x−1)x(x-1). Because this polynomial has distinct roots, every projector in finite dimensions is diagonalizable, even when it is not Hermitian. In a basis adapted to the range–kernel decomposition,

[P]=(Ir000),r=dim⁡Ran⁡P.[P]= \begin{pmatrix} I_r & 0\\ 0 & 0 \end{pmatrix}, \qquad r=\dim\operatorname{Ran}P.

It follows that

tr⁡P=rank⁡P=r.\operatorname{tr}P =\operatorname{rank}P =r.

The equality between trace and rank is basis independent. It is also a useful consistency check for an exact finite-dimensional projector, although floating-point traces need not be exact integers.

On an inner-product space, an orthogonal projector satisfies

P2=P,P†=P.P^2=P, \qquad P^\dagger=P.

Self-adjointness forces the discarded and retained subspaces to be orthogonal. For u=Pxu=Px and w∈ker⁡Pw\in\ker P,

⟨w,u⟩=⟨w,Px⟩=⟨Pw,x⟩=0.\langle w,u\rangle =\langle w,Px\rangle =\langle Pw,x\rangle =0.

Hence

ker⁡P=(Ran⁡P)⊥,V=Ran⁡P⊕⊥ker⁡P.\begin{aligned} \ker P &=(\operatorname{Ran}P)^\perp,\\ V &=\operatorname{Ran}P \mathbin{\oplus_\perp}\ker P. \end{aligned}

This is an equivalence: an idempotent is orthogonal precisely when its range is perpendicular to its kernel. There is therefore exactly one orthogonal projector onto a specified subspace MM.

The complementary operator

P⊥=I−PP^\perp=I-P

is the orthogonal projector onto M⊥M^\perp. For every vv,

⟨v,Pv⟩=∥Pv∥2≥0,⟨v,(I−P)v⟩=∥(I−P)v∥2≥0.\begin{aligned} \langle v,Pv\rangle &=\lVert Pv\rVert^2\ge 0,\\ \langle v,(I-P)v\rangle &=\lVert(I-P)v\rVert^2\ge 0. \end{aligned}

In operator order this is written 0≤P≤I0\le P\le I. Orthogonal projection cannot increase a norm:

∥Pv∥≤∥v∥.\lVert Pv\rVert\le\lVert v\rVert.

If P≠0P\ne0, its operator norm is exactly ∥P∥op=1\lVert P\rVert_{\mathrm{op}}=1, because equality holds for every nonzero vector in its range.

Let {u1,…,ur}\{u_1,\ldots,u_r\} be an orthonormal basis for a subspace MM. The orthogonal projector onto MM is

PM=∑a=1r∣ua⟩⟨ua∣.P_M =\sum_{a=1}^{r} \lvert u_a\rangle\langle u_a\rvert.

This formula is independent of the orthonormal basis chosen inside MM. For a one-dimensional subspace spanned by a nonzero, not necessarily normalized vector uu,

Pu=∣u⟩⟨u∣⟨u∣u⟩.P_u =\frac{\lvert u\rangle\langle u\rvert} {\langle u\vert u\rangle}.

Omitting the denominator is a common source of a non-idempotent outer product.

There is also a useful matrix formula. Put linearly independent spanning vectors into the columns of an n×rn\times r matrix BB. The Gram matrix B†BB^\dagger B is invertible, and

PM=B(B†B)−1B†.P_M =B(B^\dagger B)^{-1}B^\dagger.

For any xx, the residual is orthogonal to every column of BB:

B†(x−PMx)=0.B^\dagger(x-P_Mx)=0.

The formula is exact mathematics, but explicitly forming (B†B)−1(B^\dagger B)^{-1} is usually a poor numerical algorithm. A stable QR factorization B=QRB=QR gives

PM=QQ†,P_M=QQ^\dagger,

where the columns of QQ are orthonormal. A singular-value decomposition is preferable when the spanning vectors may be nearly linearly dependent.

Orthogonal projection gives the nearest vector in a subspace. Let MM be a subspace and let PMP_M be its orthogonal projector. For any m∈Mm\in M,

v−m=(v−PMv)+(PMv−m).v-m =\bigl(v-P_Mv\bigr) +\bigl(P_Mv-m\bigr).

The two terms on the right are orthogonal: the first lies in M⊥M^\perp and the second lies in MM. The Pythagorean theorem therefore gives

∥v−m∥2=∥v−PMv∥2+∥PMv−m∥2.\begin{aligned} \lVert v-m\rVert^2 &=\lVert v-P_Mv\rVert^2\\ &\quad+\lVert P_Mv-m\rVert^2. \end{aligned}

The first term is independent of mm, and the second vanishes only for m=PMvm=P_Mv. Thus

PMv=arg min⁡m∈M ∥v−m∥.P_Mv =\underset{m\in M}{\operatorname{arg\,min}}\, \lVert v-m\rVert.

This variational characterization is often more useful than the equation P2=PP^2=P. It underlies least-squares methods, basis truncation, and many variational approximations in quantum mechanics.

Consider the line in C3\mathbb C^3 spanned by

u=(11i),u†u=3.u= \begin{pmatrix} 1\\ 1\\ i \end{pmatrix}, \qquad u^\dagger u=3.

The orthogonal projector onto this line is

Pu=13(11−i11−iii1).P_u =\frac13 \begin{pmatrix} 1 & 1 & -i\\ 1 & 1 & -i\\ i & i & 1 \end{pmatrix}.

The conjugation in the row u†u^\dagger is essential. The matrix is Hermitian, has trace one, and satisfies Pu2=PuP_u^2=P_u. Acting on the first coordinate vector gives

Pu(100)=13(11i).P_u \begin{pmatrix} 1\\0\\0 \end{pmatrix} =\frac13 \begin{pmatrix} 1\\1\\i \end{pmatrix}.

The residual is orthogonal to uu, as the best-approximation theorem requires.

Idempotence alone does not imply orthogonality. For any a∈Ca\in\mathbb C,

Pa=(1a00)P_a= \begin{pmatrix} 1 & a\\ 0 & 0 \end{pmatrix}

satisfies Pa2=PaP_a^2=P_a. Its range and kernel are

Ran⁡Pa=span⁡{(1,0)T},ker⁡Pa=span⁡{(−a,1)T}.\begin{aligned} \operatorname{Ran}P_a &=\operatorname{span}\{(1,0)^T\},\\ \ker P_a &=\operatorname{span}\{(-a,1)^T\}. \end{aligned}

Unless a=0a=0, these subspaces are not orthogonal and Pa†≠PaP_a^\dagger\ne P_a. The operator projects onto the first coordinate axis along a tilted direction. Its norm is

∥Pa∥op=1+∣a∣2,\lVert P_a\rVert_{\mathrm{op}} =\sqrt{1+\lvert a\rvert^2},

so an oblique projection can amplify vectors substantially. This is one reason nearly parallel complementary subspaces lead to ill-conditioned decompositions.

A finite family of mutually orthogonal projectors satisfies

PaPb=δabPa.P_aP_b=\delta_{ab}P_a.

If the ranges span the whole space, the family is a resolution of the identity:

∑aPa=I.\sum_a P_a=I.

Every vector then decomposes into orthogonal components,

∣ψ⟩=∑aPa∣ψ⟩,∥ψ∥2=∑a∥Paψ∥2.\lvert\psi\rangle =\sum_a P_a\lvert\psi\rangle, \qquad \lVert\psi\rVert^2 =\sum_a\lVert P_a\psi\rVert^2.

Rank-one projectors from an orthonormal basis are the simplest example, but the PaP_a may have ranks greater than one. Higher-rank projectors encode degenerate subspaces without selecting a preferred basis inside them. Their role in decomposing Hermitian operators is developed in Spectral Decomposition.

Let PP and QQ be orthogonal projectors. If they commute, then

R=PQ=QPR=PQ=QP

is an orthogonal projector. Its range is the intersection

Ran⁡R=Ran⁡P∩Ran⁡Q.\operatorname{Ran}R =\operatorname{Ran}P\cap\operatorname{Ran}Q.

Commutation is essential: (PQ)†=QP(PQ)^\dagger=QP, so the product can be orthogonal only when PQ=QPPQ=QP. For commuting projectors,

P+Q−PQP+Q-PQ

is the orthogonal projector onto Ran⁡P+Ran⁡Q\operatorname{Ran}P+\operatorname{Ran}Q. These formulas are the subspace analogs of intersection and union for compatible yes–no alternatives.

For a self-adjoint matrix with distinct eigenvalues,

A=∑aaPa,PaPb=δabPa,∑aPa=I.A=\sum_a aP_a, \qquad P_aP_b=\delta_{ab}P_a, \qquad \sum_aP_a=I.

Functional calculus then gives f(A)=∑af(a)Paf(A)=\sum_a f(a)P_a. When the distinct eigenvalues are known, one can recover a projector algebraically:

Pa=∏b≠aA−bIa−b.P_a=\prod_{b\ne a}\frac{A-bI}{a-b}.

For a state ρ\rho, a yes–no projector has probability p=Tr⁡(ρP)p=\operatorname{Tr}(\rho P) and variance

(ΔρP)2=p(1−p).(\Delta_\rho P)^2=p(1-p).

These Born and state-update rules are physical postulates layered on top of the projector algebra.

If two orthogonal projectors do not commute, PQPQ is generally neither Hermitian nor idempotent. The sandwich PQPPQP is positive and obeys

0≤PQP≤P,0\le PQP\le P,

but is not generally a projector. Sequential quantum measurements therefore cannot be treated as ordinary set intersection unless the relevant projectors commute.

Unitary conjugation preserves orthogonal projection:

P⟼UPU†,Ran⁡(UPU†)=U(Ran⁡P).P\longmapsto UPU^\dagger, \qquad \operatorname{Ran}(UPU^\dagger)=U(\operatorname{Ran}P).

A projective measurement is specified by an orthogonal resolution of the identity {Pa}\{P_a\}. The projector algebra ensures that the outcome subspaces are exclusive and complete. The physical postulates then assign outcome probabilities and conditional states to those subspaces.

Those postulates are not consequences of idempotence. See Projective Measurement for the Born probabilities, repeatability, degeneracy, and Lüders state update. Generalized measurement effects are positive operators and need not be projectors.

Projectors also have a distinct approximation-theory role. Complementary subspaces can separate retained and eliminated sectors of a Hamiltonian; see Projection Methods.

In a Hilbert space, every closed subspace MM has a unique bounded orthogonal projector PMP_M. Closedness matters: the nearest-point limit must remain in MM. A nonclosed proper subspace has no bounded orthogonal projector whose range is exactly that subspace.

For any bounded idempotent PP, both ker⁡P\ker P and Ran⁡P\operatorname{Ran}P are closed because

ker⁡P=P−1({0}),Ran⁡P=ker⁡(I−P).\begin{aligned} \ker P &=P^{-1}(\{0\}),\\ \operatorname{Ran}P &=\ker(I-P). \end{aligned}

Infinite resolutions of the identity require a convergence statement. For a countable orthogonal family, the partial sums converge strongly when

∑aPaψ=ψ\sum_a P_a\psi=\psi

for every vector ψ\psi; norm convergence of the operators is generally too strong. The spectral theorem extends this idea from sums to projection-valued measures.

For a computed matrix P~\widetilde P, inspect both

ϵid=∥P~2−P~∥,ϵH=∥P~†−P~∥.\epsilon_{\mathrm{id}} =\lVert\widetilde P^2-\widetilde P\rVert, \qquad \epsilon_{\mathrm{H}} =\lVert\widetilde P^\dagger-\widetilde P\rVert.

The first tests idempotence; the second distinguishes an approximate orthogonal projector from an approximate oblique one. Eigenvalues should cluster near 00 and 11, and the trace should be close to the intended rank. These checks must be interpreted relative to matrix size, norm convention, and floating-point precision.

When constructing a projector from spanning vectors:

  • orthonormalize with QR when the numerical rank is clear;
  • use an SVD and an explicit singular-value tolerance when rank is uncertain;
  • avoid explicit matrix inversion;
  • check the residual orthogonality B†(x−P~x)B^\dagger(x-\widetilde Px);
  • do not repair a poor basis by merely rounding eigenvalues.

See Matrix Diagonalization for related conditioning and residual checks.

  • Assuming P2=PP^2=P automatically implies P†=PP^\dagger=P.
  • Saying “the projector onto MM” without specifying the complementary direction unless orthogonal projection is intended.
  • Using ∣u⟩⟨u∣\lvert u\rangle\langle u\rvert for an unnormalized vector.
  • Forgetting complex conjugation when forming an outer product.
  • Applying the best-approximation theorem to an oblique projector.
  • Multiplying orthogonal projectors and assuming the product is a projector without checking commutation.
  • Treating every positive measurement effect as a projector.
  • Expecting an infinite resolution of identity to converge in operator norm.
  1. Let

    P=(1200).P= \begin{pmatrix} 1 & 2\\ 0 & 0 \end{pmatrix}.

    Verify idempotence, find the range and kernel, and show that PP is not an orthogonal projector. Compute its operator norm.

Solution

Direct multiplication gives

P2=(1200)=P.P^2= \begin{pmatrix} 1 & 2\\ 0 & 0 \end{pmatrix} =P.

The output P(x,y)T=(x+2y,0)TP(x,y)^T=(x+2y,0)^T lies on the first coordinate axis, so

Ran⁡P=span⁡{(1,0)T}.\operatorname{Ran}P =\operatorname{span}\{(1,0)^T\}.

The kernel condition x+2y=0x+2y=0 gives

ker⁡P=span⁡{(−2,1)T}.\ker P =\operatorname{span}\{(-2,1)^T\}.

These lines are not orthogonal, and

P†=(1020)≠P.P^\dagger= \begin{pmatrix} 1 & 0\\ 2 & 0 \end{pmatrix} \ne P.

Finally,

PP†=(5000),PP^\dagger= \begin{pmatrix} 5 & 0\\ 0 & 0 \end{pmatrix},

so the largest singular value is 5\sqrt5 and ∥P∥op=5\lVert P\rVert_{\mathrm{op}}=\sqrt5.

  1. Let M⊂R3M\subset\mathbb R^3 be spanned by b1=(1,0,1)Tb_1=(1,0,1)^T and b2=(0,1,1)Tb_2=(0,1,1)^T. Find the orthogonal projector onto MM without first orthonormalizing b1,b2b_1,b_2.
Solution

A normal vector to both spanning vectors is n=(1,1,−1)Tn=(1,1,-1)^T. Therefore the projector onto M=n⊥M=n^\perp is the complement of the rank-one projector onto nn:

PM=I−nnTnTn=13(2−11−121112).\begin{aligned} P_M &=I-\frac{nn^T}{n^Tn}\\ &=\frac13 \begin{pmatrix} 2 & -1 & 1\\ -1 & 2 & 1\\ 1 & 1 & 2 \end{pmatrix}. \end{aligned}

The matrix is symmetric, has trace two, fixes both b1b_1 and b2b_2, and annihilates nn. These facts also verify that its range is MM.

  1. Prove the best-approximation theorem directly from the normal equation PM†=PMP_M^\dagger=P_M: show that PMvP_Mv is the unique minimizer of ∥v−m∥\lVert v-m\rVert over m∈Mm\in M.
Solution

Because PMP_M fixes MM, for any m∈Mm\in M,

PM(PMv−m)=PMv−m.P_M(P_Mv-m)=P_Mv-m.

The residual obeys

PM(v−PMv)=0,P_M(v-P_Mv)=0,

so it belongs to ker⁡PM=M⊥\ker P_M=M^\perp. Hence v−PMvv-P_Mv is orthogonal to PMv−mP_Mv-m, and

∥v−m∥2=∥v−PMv∥2+∥PMv−m∥2.\begin{aligned} \lVert v-m\rVert^2 &=\lVert v-P_Mv\rVert^2\\ &\quad+\lVert P_Mv-m\rVert^2. \end{aligned}

The right side is minimized exactly when the second term vanishes, which occurs only for m=PMvm=P_Mv.

  1. Let PP and QQ be commuting orthogonal projectors. Prove that PQPQ is the orthogonal projector onto Ran⁡P∩Ran⁡Q\operatorname{Ran}P\cap\operatorname{Ran}Q. Then show that P+Q−PQP+Q-PQ is an orthogonal projector.
Solution

Commutation and idempotence imply

(PQ)2=P2Q2=PQ,(PQ)†=QP=PQ.\begin{aligned} (PQ)^2 &=P^2Q^2=PQ,\\ (PQ)^\dagger &=QP=PQ. \end{aligned}

Thus PQPQ is an orthogonal projector. If x=PQyx=PQy, then Px=xPx=x and Qx=xQx=x, so its range lies in the intersection. Conversely, if Px=Qx=xPx=Qx=x, then PQx=xPQx=x, so every vector in the intersection belongs to the range.

Set S=P+Q−PQS=P+Q-PQ. Since all factors commute and are Hermitian, S†=SS^\dagger=S. Expanding and using P2=PP^2=P, Q2=QQ^2=Q, and (PQ)2=PQ(PQ)^2=PQ gives

S2=P+Q−PQ=S.S^2=P+Q-PQ=S.

Therefore SS is an orthogonal projector. It fixes both Ran⁡P\operatorname{Ran}P and Ran⁡Q\operatorname{Ran}Q, and its image lies in their sum, so its range is Ran⁡P+Ran⁡Q\operatorname{Ran}P+\operatorname{Ran}Q.

  • S. Axler, Linear Algebra Done Right, 3rd ed., Springer, 2015.
  • R. A. Horn and C. R. Johnson, Matrix Analysis, 2nd ed., Cambridge University Press, 2013.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013.
  • M. Reed and B. Simon, Methods of Modern Mathematical Physics, Vol. I: Functional Analysis, revised and enlarged ed., Academic Press, 1980.
  • G. H. Golub and C. F. Van Loan, Matrix Computations, 4th ed., Johns Hopkins University Press, 2013.
  • P. R. Halmos, Introduction to Hilbert Space and the Theory of Spectral Multiplicity, 2nd ed., Chelsea, 1957.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.