Skip to content

Spectral Decomposition

The spectral decomposition of a finite-dimensional self-adjoint observable AA is

A=∑a∈σ(A)aPa,I=∑a∈σ(A)Pa.A = \sum_{a\in\sigma(A)}aP_a, \qquad I = \sum_{a\in\sigma(A)}P_a.

Here σ(A)\sigma(A) is the set of distinct eigenvalues, and PaP_a is the orthogonal projector onto the full eigenspace with eigenvalue aa. The projectors satisfy

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

This one structure organizes several parts of quantum mechanics:

  • it is the basis-independent form of diagonalization;
  • it separates a state into mutually orthogonal outcome sectors;
  • it supplies the projectors used by the Born rule;
  • it makes powers, exponentials, and other functions of AA transparent;
  • it exposes how degeneracy and commuting observables produce block structure;
  • and it extends to continuous spectra through a projection-valued measure.

The page develops these consequences without proving the general infinite-dimensional spectral theorem. That theorem and its domain conditions belong to Spectral Theorem, Practical Version.

Let AA be self-adjoint on a finite-dimensional Hilbert space H\mathcal H. Its distinct eigenvalues are real, and the Hilbert space is an orthogonal direct sum of eigenspaces:

H=⨁a∈σ(A)Ea,\mathcal H = \mathop{\bigoplus}_{a\in\sigma(A)} \mathcal E_a,

where

Ea=Ker⁡(A−aI).\mathcal E_a = \operatorname{Ker}(A-aI).

Let PaP_a be the orthogonal projector onto Ea\mathcal E_a. Then

PaPb=δabPa,Pa†=Pa,∑aPa=I,APa=PaA=aPa.\begin{aligned} P_aP_b &= \delta_{ab}P_a,\\ P_a^\dagger &=P_a,\\ \sum_aP_a &=I,\\ AP_a &=P_aA=aP_a. \end{aligned}

The operator is reconstructed by

A=∑aaPa.A=\sum_a aP_a.

The sum is over distinct spectral values. Degeneracy is contained in the rank of each PaP_a, not represented by repeating the value aa in this sum.

The finite-dimensional spectral theorem provides an orthonormal basis of eigenvectors. Group those eigenvectors by their distinct eigenvalues. Every state has a unique decomposition

∣ψ⟩=∑aPa∣ψ⟩.|\psi\rangle = \sum_aP_a|\psi\rangle.

Because Pa∣ψ⟩∈EaP_a|\psi\rangle\in\mathcal E_a,

APa∣ψ⟩=aPa∣ψ⟩.A P_a|\psi\rangle = aP_a|\psi\rangle.

Therefore

A∣ψ⟩=A∑aPa∣ψ⟩=∑aaPa∣ψ⟩.\begin{aligned} A|\psi\rangle &= A\sum_aP_a|\psi\rangle\\ &= \sum_a aP_a|\psi\rangle. \end{aligned}

Since this equality holds for every ∣ψ⟩|\psi\rangle,

A=∑aaPa.A=\sum_a aP_a.

This argument shows what the decomposition means operationally: resolve the state into eigenspace components, multiply each component by its spectral value, and add the results.

A state resolved into eigenspace components that are scaled by their eigenvalues and recombined

The identity resolves ∣ψ⟩|\psi\rangle into orthogonal components Pai∣ψ⟩P_{a_i}|\psi\rangle. On each eigenspace, AA acts as multiplication by the corresponding value aia_i.

Spectral Decomposition versus Diagonalization

Section titled “Spectral Decomposition versus Diagonalization”

Diagonalization is often written as

A=UDU†,A=UDU^\dagger,

where UU is unitary and DD is diagonal. This is a matrix description after an orthonormal eigenbasis has been chosen. Spectral decomposition instead writes the basis-independent operator identity

A=∑aaPa.A=\sum_a aP_a.

The distinction matters when an eigenvalue is degenerate:

  • the diagonal matrix repeats the value once for each basis vector in its eigenspace;
  • the spectral sum lists that value once and assigns it one higher-rank projector;
  • a different orthonormal basis inside the degenerate eigenspace changes UU but leaves PaP_a unchanged.

If DD has entries a1,…,ada_1,\ldots,a_d and ∣un⟩|u_n\rangle are the columns of UU, one may write

A=∑n=1dan∣un⟩⟨un∣.A = \sum_{n=1}^d a_n|u_n\rangle\langle u_n|.

When values repeat, this sum is over basis vectors rather than distinct outcomes. Grouping equal values recovers the spectral-projector form.

If each eigenvalue aa has a one-dimensional eigenspace and ∣a⟩|a\rangle is a normalized eigenvector, then

Pa=∣a⟩⟨a∣.P_a = |a\rangle\langle a|.

The decomposition becomes

A=∑aa∣a⟩⟨a∣,A = \sum_a a|a\rangle\langle a|,

with completeness relation

∑a∣a⟩⟨a∣=I.\sum_a|a\rangle\langle a|=I.

Acting on a state

∣ψ⟩=∑aca∣a⟩|\psi\rangle = \sum_a c_a|a\rangle

gives

A∣ψ⟩=∑aaca∣a⟩.A|\psi\rangle = \sum_a a c_a|a\rangle.

This is the familiar textbook form. It is safe only when the sum is understood to include a complete eigenbasis and degeneracies are handled correctly.

Suppose the eigenspace Ea\mathcal E_a has dimension gag_a. Choose any orthonormal basis

{∣a,λ⟩:λ=1,…,ga}.\left\lbrace |a,\lambda\rangle: \lambda=1,\ldots,g_a \right\rbrace.

The eigenspace projector is

Pa=∑λ=1ga∣a,λ⟩⟨a,λ∣.P_a = \sum_{\lambda=1}^{g_a} |a,\lambda\rangle \langle a,\lambda|.

The spectral decomposition remains

A=∑aaPa.A=\sum_a aP_a.

The index aa labels a distinct outcome; λ\lambda labels basis vectors inside that outcome subspace. A unitary change of basis among the ∣a,λ⟩|a,\lambda\rangle leaves PaP_a invariant.

This is why a degenerate sharp measurement is specified by the full PaP_a, not by an arbitrary rank-one projector inside Ea\mathcal E_a. The resulting state-update question is treated in Degenerate Measurements and Lüders Rule.

For a self-adjoint AA, the projector associated with a distinct eigenvalue is fixed by AA itself. It does not depend on a chosen eigenbasis.

In finite dimensions this can be seen through Lagrange interpolation. For each a∈σ(A)a\in\sigma(A), define

qa(x)=∏b∈σ(A)b≠ax−ba−b.q_a(x) = \prod_{\substack{b\in\sigma(A)\\b\neq a}} \frac{x-b}{a-b}.

This polynomial satisfies

qa(b)=δabq_a(b)=\delta_{ab}

on the spectrum. Applying it to AA gives

Pa=qa(A)=∏b∈σ(A)b≠aA−bIa−b.P_a = q_a(A) = \prod_{\substack{b\in\sigma(A)\\b\neq a}} \frac{A-bI}{a-b}.

On an eigenvector with value cc, this product returns the vector if c=ac=a and zero otherwise. It therefore projects onto the entire aa eigenspace, including all degeneracy.

If AA has exactly two distinct values aa and bb, then

Pa=A−bIa−b,Pb=A−aIb−a.P_a = \frac{A-bI}{a-b}, \qquad P_b = \frac{A-aI}{b-a}.

These formulas extract the projectors without first writing eigenvectors. They are especially useful for qubits, parity sectors, and other binary observables.

The phrase projector-valued decomposition emphasizes that the operator is not merely a list of eigenvalues. It is the pairing

a⟷Paa \longleftrightarrow P_a

between each distinct value and the subspace on which that value acts.

The values answer “what number is reported?” The projectors answer “which component of the state supports that outcome?” Both are needed to reconstruct AA.

Two observables can have the same set of numerical eigenvalues but different projectors. For example, σz\sigma_z and σx\sigma_x both have values ±1\pm1, but their eigenspaces differ. Their statistics in a given state therefore differ.

Conversely, assigning different distinct numerical labels to the same family of projectors changes the operator but keeps the same sharp partition of Hilbert space. This observation becomes precise under functions of operators and coarse graining below.

For a normalized state, define

∣ψa⟩=Pa∣ψ⟩.|\psi_a\rangle = P_a|\psi\rangle.

The sector components satisfy

⟨ψa∣ψb⟩=0(a≠b),\langle\psi_a|\psi_b\rangle=0 \qquad (a\neq b),

and

∣ψ⟩=∑a∣ψa⟩.|\psi\rangle = \sum_a|\psi_a\rangle.

Their squared norms obey

∑a∥ψa∥2=1.\sum_a\|\psi_a\|^2=1.

When ∥ψa∥>0\|\psi_a\|>0, one may define the normalized sector vector

∣ψ^a⟩=Pa∣ψ⟩∥Paψ∥.|\widehat\psi_a\rangle = \frac{P_a|\psi\rangle}{\|P_a\psi\|}.

Then

∣ψ⟩=∑a: p(a)>0p(a)∣ψ^a⟩,|\psi\rangle = \sum_{a:\,p(a)>0} \sqrt{p(a)} |\widehat\psi_a\rangle,

where the phase of each sector is inherited from Pa∣ψ⟩P_a|\psi\rangle. This is an orthogonal decomposition of the original state, not yet a claim that a measurement has occurred.

For a normalized pure state, the Born probability of outcome aa is

p(a)=⟨ψ∣Pa∣ψ⟩=∥Paψ∥2.p(a) = \langle\psi|P_a|\psi\rangle = \|P_a\psi\|^2.

Completeness guarantees normalization:

∑ap(a)=∑a⟨ψ∣Pa∣ψ⟩=⟨ψ∣I∣ψ⟩=1.\begin{aligned} \sum_a p(a) &= \sum_a\langle\psi|P_a|\psi\rangle\\ &= \langle\psi|I|\psi\rangle\\ &=1. \end{aligned}

For a density operator ρ\rho,

p(a)=Tr⁡(ρPa).p(a) = \operatorname{Tr}(\rho P_a).

Spectral decomposition therefore separates the observable into the precise operators used to compute its outcome statistics. The canonical calculation workflow is Born Rule for Discrete Spectra.

Orthogonality of the spectral projectors gives

An=∑aanPaA^n = \sum_a a^nP_a

for every nonnegative integer nn. Hence

⟨An⟩ψ=∑aanp(a).\langle A^n\rangle_\psi = \sum_a a^n p(a).

In particular,

⟨A⟩ψ=∑aap(a)\langle A\rangle_\psi = \sum_a a p(a)

and

(ΔA)2=∑a(a−⟨A⟩ψ)2p(a).(\Delta A)^2 = \sum_a \left(a-\langle A\rangle_\psi\right)^2 p(a).

These are the ordinary moments and variance of the probability distribution induced by the state and the spectral projectors. Spectral decomposition does not replace the Born rule; it identifies the event operators to which the rule is applied.

If ff is defined on the finite spectrum of AA, then

f(A)=∑af(a)Pa.f(A) = \sum_a f(a)P_a.

The eigenspaces remain the same while each value aa is replaced by f(a)f(a). Examples include

A2=∑aa2Pa,e−iAt=∑ae−iatPa,∣A∣=∑a∣a∣Pa.\begin{aligned} A^2 &= \sum_a a^2P_a,\\ e^{-iAt} &= \sum_a e^{-iat}P_a,\\ |A| &= \sum_a |a|P_a. \end{aligned}

An inverse exists exactly when 0∉σ(A)0\notin\sigma(A), in which case

A−1=∑a1aPa.A^{-1} = \sum_a\frac{1}{a}P_a.

The full construction, including square roots, branch choices, and infinite-dimensional domains, is the canonical subject of Functions of Operators.

Suppose B=f(A)B=f(A). Its values are f(a)f(a), but distinct values of aa need not remain distinct.

If ff is one-to-one on σ(A)\sigma(A), then AA and BB have the same spectral projectors. They assign different numerical labels to the same sharp outcome subspaces.

If f(a)=f(b)f(a)=f(b) for two different values, the corresponding subspaces merge. For a value β\beta of BB, its projector is

Qβ=∑a: f(a)=βPa.Q_\beta = \sum_{a:\,f(a)=\beta}P_a.

The observable BB is then a coarse graining of AA: it cannot distinguish values that ff maps to the same label.

This distinction is physically useful. The operator A2A^2, for example, does not distinguish eigenvalues aa and −a-a even when AA does.

Let A=∑aaPaA=\sum_a aP_a have distinct values aa. For any operator BB in finite dimensions,

B=∑a,bPaBPb.B = \sum_{a,b}P_aBP_b.

The terms PaBPbP_aBP_b are blocks mapping Eb\mathcal E_b into Ea\mathcal E_a. The commutator is

[A,B]=∑a,b(a−b)PaBPb.[A,B] = \sum_{a,b} (a-b)P_aBP_b.

Therefore

[A,B]=0[A,B]=0

if and only if all off-diagonal spectral blocks vanish:

PaBPb=0(a≠b).P_aBP_b=0 \qquad (a\neq b).

Equivalently,

B=∑aPaBPa,B = \sum_aP_aBP_a,

so BB preserves every eigenspace of AA.

If AA is nondegenerate, each block is one-dimensional and a commuting self-adjoint BB is diagonal in the same eigenbasis. If AA is degenerate, BB can act nontrivially within each Ea\mathcal E_a. Diagonalizing those restrictions is the first step toward a Complete Set of Commuting Observables.

If a time-independent Hamiltonian has discrete spectral decomposition

H=∑EEPE,H=\sum_E E P_E,

then

U(t)=e−iHt/ℏ=∑Ee−iEt/ℏPE.U(t) = e^{-iHt/\hbar} = \sum_E e^{-iEt/\hbar}P_E.

The evolved state is

∣ψ(t)⟩=∑Ee−iEt/ℏPE∣ψ(0)⟩.|\psi(t)\rangle = \sum_E e^{-iEt/\hbar} P_E|\psi(0)\rangle.

Each energy sector acquires a phase. Every vector inside a degenerate energy eigenspace receives the same phase, so the state component within that eigenspace is not resolved further by HH alone.

The energy-sector weights are constant:

∥PEψ(t)∥2=∥PEψ(0)∥2.\|P_E\psi(t)\|^2 = \|P_E\psi(0)\|^2.

The dynamical interpretation is developed in Energy Eigenstates and Unitary Time Evolution.

Let ga=rank⁡Pag_a=\operatorname{rank}P_a be the degeneracy of aa in a finite-dimensional space. Spectral decomposition makes several basis-independent quantities immediate:

Tr⁡A=∑agaa,det⁡A=∏aaga,∥A∥op=max⁡a∣a∣.\begin{aligned} \operatorname{Tr}A &= \sum_a g_a a,\\ \det A &= \prod_a a^{g_a},\\ \|A\|_{\mathrm{op}} &= \max_a|a|. \end{aligned}

More generally,

Tr⁡f(A)=∑agaf(a).\operatorname{Tr}f(A) = \sum_a g_a f(a).

The determinant formula includes multiplicity. It vanishes when zero is an eigenvalue, exactly when AA is not invertible.

Consider

A=(2i−i2).A = \begin{pmatrix} 2&i\\ -i&2 \end{pmatrix}.

Its distinct eigenvalues are 11 and 33. Because there are only two values, the projectors can be extracted directly:

P1=3I−A2=12(1−ii1),P3=A−I2=12(1i−i1).\begin{aligned} P_1 &= \frac{3I-A}{2} = \frac12 \begin{pmatrix} 1&-i\\ i&1 \end{pmatrix},\\ P_3 &= \frac{A-I}{2} = \frac12 \begin{pmatrix} 1&i\\ -i&1 \end{pmatrix}. \end{aligned}

They satisfy

P1+P3=I,P1P3=0,P_1+P_3=I, \qquad P_1P_3=0,

and the reconstruction is

A=P1+3P3.A=P_1+3P_3.

For the state ∣0⟩=(1,0)T|0\rangle=(1,0)^{\mathsf T},

p(1)=⟨0∣P1∣0⟩=12,p(1) = \langle0|P_1|0\rangle = \frac12,

and similarly p(3)=1/2p(3)=1/2. Therefore

⟨A⟩=1⋅12+3⋅12=2,\langle A\rangle = 1\cdot\frac12 +3\cdot\frac12 =2,

which agrees with the direct matrix element ⟨0∣A∣0⟩=2\langle0|A|0\rangle=2.

For a unit vector n\boldsymbol n, define

σn=n⋅σ.\sigma_{\boldsymbol n} = \boldsymbol n\cdot\boldsymbol\sigma.

Its values are +1+1 and −1-1, with projectors

P+=I+σn2,P−=I−σn2.P_+ = \frac{I+\sigma_{\boldsymbol n}}{2}, \qquad P_- = \frac{I-\sigma_{\boldsymbol n}}{2}.

The spectral decomposition is

σn=(+1)P++(−1)P−=P+−P−.\sigma_{\boldsymbol n} = (+1)P_+ +(-1)P_- = P_+-P_-.

For spin angular momentum

Sn=ℏ2σn,S_{\boldsymbol n} = \frac{\hbar}{2} \sigma_{\boldsymbol n},

the projectors are unchanged while the values become ±ℏ/2\pm\hbar/2. The sharp alternatives are the same subspaces; only their numerical labels have been rescaled.

For a self-adjoint operator with continuous or mixed spectrum, a simple sum of normalizable eigenprojectors is generally unavailable. The spectral theorem assigns a projection-valued measure EAE_A. The decomposition becomes the spectral integral

A=∫Rλ dEA(λ),A = \int_{\mathbb R} \lambda\,dE_A(\lambda),

with identity resolution

I=EA(R)=∫RdEA(λ).I = E_A(\mathbb R) = \int_{\mathbb R}dE_A(\lambda).

For a Borel set Δ⊆R\Delta\subseteq\mathbb R, the projector EA(Δ)E_A(\Delta) selects the spectral component whose values lie in Δ\Delta. For a normalized pure state,

Pr⁡(A∈Δ)=⟨ψ∣EA(Δ)∣ψ⟩.\Pr(A\in\Delta) = \langle\psi|E_A(\Delta)|\psi\rangle.

At a discrete eigenvalue,

EA({a})=Pa.E_A(\{a\})=P_a.

At a purely continuous spectral point, the singleton projector is typically zero even though every interval around the point has nonzero spectral projector. This is why the formal notation

A∼∫a∣a⟩⟨a∣ daA \sim \int a|a\rangle\langle a|\,da

is a useful heuristic but not the rigorous starting point: generalized kets ∣a⟩|a\rangle need not be Hilbert-space vectors.

If AA is unbounded, the spectral integral for AA is defined only on vectors whose second spectral moment is finite. Writing

μψ(Δ)=⟨ψ∣EA(Δ)∣ψ⟩,\mu_\psi(\Delta) = \langle\psi|E_A(\Delta)|\psi\rangle,

the domain is

D(A)={∣ψ⟩:∫Rλ2 dμψ(λ)<∞}.\mathcal D(A) = \left\lbrace |\psi\rangle: \int_{\mathbb R} \lambda^2\,d\mu_\psi(\lambda) <\infty \right\rbrace.

The projectors EA(Δ)E_A(\Delta) themselves are bounded and everywhere defined. This distinction lets spectral probabilities remain meaningful even when the state does not have a finite expectation value of A2A^2.

For point, continuous, essential, and mixed spectra, see Discrete and Continuous Spectra. For the rigorous working theorem, continue to Spectral Theorem, Practical Version.

In finite dimensions, every self-adjoint operator admits an orthogonal spectral decomposition with real values. More generally, every finite-dimensional normal operator has an orthogonal spectral decomposition, though its values may be complex.

A general nonnormal operator need not have an orthonormal eigenbasis. It may have nonorthogonal eigenvectors, oblique spectral projectors, or Jordan blocks. Writing an arbitrary matrix as if it had the observable decomposition

A=∑aaPaA=\sum_a aP_a

with mutually orthogonal PaP_a is therefore unjustified.

For unbounded quantum observables, formal Hermiticity or symmetry on a test domain is not enough. The projection-valued spectral theorem requires a specified self-adjoint operator; see Hermitian vs Self-Adjoint Operators.

For a finite-dimensional self-adjoint observable:

  1. Find the distinct eigenvalues. Keep track of their degeneracies.
  2. Determine each eigenspace. Solve (A−aI)∣ψ⟩=0(A-aI)|\psi\rangle=0.
  3. Construct PaP_a. Sum rank-one projectors over an orthonormal basis of the eigenspace, or use the polynomial formula.
  4. Check orthogonality. Verify PaPb=0P_aP_b=0 for a≠ba\neq b.
  5. Check completeness. Verify ∑aPa=I\sum_aP_a=I.
  6. Reconstruct the operator. Verify A=∑aaPaA=\sum_a aP_a.
  7. Compute statistics. Use p(a)=⟨ψ∣Pa∣ψ⟩p(a)=\langle\psi|P_a|\psi\rangle or p(a)=Tr⁡(ρPa)p(a)=\operatorname{Tr}(\rho P_a).
  8. Apply functions spectrally. Replace aa by f(a)f(a) while keeping the projectors, then account for any merged values.

The completeness and reconstruction checks catch most missing-degeneracy and sign errors.

  • Summing over eigenvectors without tracking repeated values. The basis-independent spectral sum is over distinct eigenvalues and full eigenspace projectors.
  • Replacing a degenerate projector by one rank-one term. A sharp outcome aa corresponds to all of Ea\mathcal E_a.
  • Treating diagonalization as basis independent. The matrix DD and unitary UU depend on eigenbasis choices; the projectors PaP_a do not.
  • Applying ff to matrix entries. The rule is f(A)=∑af(a)Paf(A)=\sum_a f(a)P_a, not entrywise evaluation in an arbitrary basis.
  • Forgetting that a function can merge outcomes. If f(a)=f(b)f(a)=f(b), the spectral projector of f(A)f(A) is Pa+PbP_a+P_b for that combined value.
  • Using eigenvalues without projectors to predict statistics. Equal spectra do not imply equal observables; the eigenspaces matter.
  • Confusing state decomposition with measurement update. The identity ∣ψ⟩=∑aPa∣ψ⟩|\psi\rangle=\sum_aP_a|\psi\rangle is linear algebra, not a collapse postulate.
  • Writing a continuous spectral integral as an ordinary sum of normalizable eigenstates. Use the projection-valued measure for the exact statement.
  • Ignoring domains for unbounded operators. The spectral measure is everywhere defined, but AA and unbounded functions f(A)f(A) are not.
  • Applying the orthogonal formula to a nonnormal matrix. Nonnormal operators need not admit mutually orthogonal spectral projectors.

This page owns the physical reading of an observable as values paired with orthogonal outcome subspaces. Neighboring canonical pages carry the deeper pieces:

  • A finite-dimensional self-adjoint observable has the unique decomposition A=∑aaPaA=\sum_a aP_a over its distinct eigenvalues.
  • The projectors are mutually orthogonal and resolve the identity.
  • Nondegenerate values give rank-one projectors; degenerate values give higher-rank eigenspace projectors.
  • The formula is the basis-independent version of unitary diagonalization.
  • Finite spectral projectors can be recovered as interpolation polynomials in AA.
  • A state decomposes into orthogonal sectors Pa∣ψ⟩P_a|\psi\rangle, whose squared norms are Born probabilities.
  • Moments, functions, and time evolution follow by applying the corresponding scalar expression to each spectral value.
  • Commuting operators are block diagonal with respect to the spectral projectors.
  • A non-injective function of AA coarse-grains its spectral outcomes.
  • Continuous and mixed spectra replace the sum by A=∫λ dEA(λ)A=\int\lambda\,dE_A(\lambda).
  • The projection-valued theorem requires self-adjointness and, for unbounded operators, explicit domain control.
  • P. A. M. Dirac, The Principles of Quantum Mechanics, 4th ed., Oxford University Press, 1958.
  • J. von Neumann, Mathematical Foundations of Quantum Mechanics, Princeton University Press, 1955.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013.
  • M. Reed and B. Simon, Methods of Modern Mathematical Physics, Volume I: Functional Analysis, Academic Press, 1980.
  • S. Axler, Linear Algebra Done Right, 3rd ed., Springer, 2015.

Exercise 1: Verify a two-level decomposition

Section titled “Exercise 1: Verify a two-level decomposition”

For

A=(2i−i2),A = \begin{pmatrix} 2&i\\ -i&2 \end{pmatrix},

verify directly that the matrices

P1=12(1−ii1),P3=12(1i−i1)P_1 = \frac12 \begin{pmatrix} 1&-i\\ i&1 \end{pmatrix}, \qquad P_3 = \frac12 \begin{pmatrix} 1&i\\ -i&1 \end{pmatrix}

are orthogonal projectors and reconstruct AA.

Solution

Both matrices are self-adjoint. Direct multiplication gives

P12=P1,P32=P3,P1P3=0.P_1^2=P_1, \qquad P_3^2=P_3, \qquad P_1P_3=0.

Their sum is

P1+P3=I.P_1+P_3=I.

Finally,

P1+3P3=(2i−i2)=A.P_1+3P_3 = \begin{pmatrix} 2&i\\ -i&2 \end{pmatrix} =A.

Let AA have distinct eigenvalues −1-1, 00, and 22. Write each spectral projector as a polynomial in AA.

Solution

Using the interpolation formula,

P−1=A(A−2I)3,P0=−(A+I)(A−2I)2,P2=A(A+I)6.\begin{aligned} P_{-1} &= \frac{A(A-2I)}{3},\\ P_0 &= -\frac{(A+I)(A-2I)}{2},\\ P_2 &= \frac{A(A+I)}{6}. \end{aligned}

Each polynomial equals one on its target value and zero on the other two. Consequently P−1+P0+P2=IP_{-1}+P_0+P_2=I on the whole Hilbert space.

Exercise 3: Degenerate outcome probabilities

Section titled “Exercise 3: Degenerate outcome probabilities”

Let

A=a(∣1⟩⟨1∣+∣2⟩⟨2∣)+b∣3⟩⟨3∣,A = a \left( |1\rangle\langle1| +|2\rangle\langle2| \right) +b|3\rangle\langle3|,

where a≠ba\neq b, and let

∣ψ⟩=∣1⟩+i∣2⟩+2∣3⟩2.|\psi\rangle = \frac{|1\rangle+i|2\rangle+\sqrt2|3\rangle}{2}.

Find the probabilities of aa and bb and the expectation value of AA.

Solution

The spectral projectors are

Pa=∣1⟩⟨1∣+∣2⟩⟨2∣,Pb=∣3⟩⟨3∣.P_a = |1\rangle\langle1| +|2\rangle\langle2|, \qquad P_b = |3\rangle\langle3|.

Therefore

p(a)=14+14=12,p(b)=24=12.p(a) = \frac14+\frac14 = \frac12, \qquad p(b) = \frac{2}{4} = \frac12.

The expectation value is

⟨A⟩=a+b2.\langle A\rangle = \frac{a+b}{2}.

An observable has values −2-2, 00, and 33 with probabilities 1/41/4, 1/21/2, and 1/41/4. Compute its expectation value and variance.

Solution

The expectation value is

⟨A⟩=−2(14)+3(14)=14.\langle A\rangle = -2\left(\frac14\right) +3\left(\frac14\right) = \frac14.

The second moment is

⟨A2⟩=4(14)+9(14)=134.\langle A^2\rangle = 4\left(\frac14\right) +9\left(\frac14\right) = \frac{13}{4}.

Hence

(ΔA)2=134−(14)2=5116.(\Delta A)^2 = \frac{13}{4} -\left(\frac14\right)^2 = \frac{51}{16}.

Let

A=−2P−2+P1+2P2,A = -2P_{-2} +P_1 +2P_2,

where the three projectors are mutually orthogonal and complete. Find the spectral decomposition of A2A^2.

Solution

Applying f(a)=a2f(a)=a^2 gives

A2=4P−2+P1+4P2.A^2 = 4P_{-2} +P_1 +4P_2.

The value 44 occurs on two subspaces, so its spectral projector is their sum:

Q4=P−2+P2.Q_4 = P_{-2}+P_2.

The distinct-value decomposition is therefore

A2=1P1+4Q4.A^2 = 1P_1 +4Q_4.

Squaring has erased the distinction between the values −2-2 and 22.

Exercise 6: A commuting operator preserves eigenspaces

Section titled “Exercise 6: A commuting operator preserves eigenspaces”

Let A=∑aaPaA=\sum_a aP_a have distinct eigenvalues. Suppose [A,B]=0[A,B]=0. Show that

PaBPb=0(a≠b).P_aBP_b=0 \qquad (a\neq b).
Solution

Multiply the commutator by PaP_a on the left and PbP_b on the right:

0=Pa[A,B]Pb=(a−b)PaBPb.\begin{aligned} 0 &= P_a[A,B]P_b\\ &= (a-b)P_aBP_b. \end{aligned}

For a≠ba\neq b, the scalar a−ba-b is nonzero, so PaBPb=0P_aBP_b=0. Thus BB has no matrix blocks connecting different eigenspaces of AA.

Exercise 7: Time evolution by energy sectors

Section titled “Exercise 7: Time evolution by energy sectors”

Suppose

H=E0P0+E1P1,P0+P1=I.H=E_0P_0+E_1P_1, \qquad P_0+P_1=I.

Derive U(t)U(t) and show that the probabilities of the two energy values are constant in time.

Solution

The spectral rule gives

U(t)=e−iE0t/ℏP0+e−iE1t/ℏP1.U(t) = e^{-iE_0t/\hbar}P_0 +e^{-iE_1t/\hbar}P_1.

Because PjPk=δjkPjP_jP_k=\delta_{jk}P_j,

PjU(t)∣ψ(0)⟩=e−iEjt/ℏPj∣ψ(0)⟩.P_jU(t)|\psi(0)\rangle = e^{-iE_jt/\hbar} P_j|\psi(0)\rangle.

Taking the squared norm removes the phase:

∥Pjψ(t)∥2=∥Pjψ(0)∥2.\|P_j\psi(t)\|^2 = \|P_j\psi(0)\|^2.

Let XX be the position operator on L2(R)L^2(\mathbb R) and let RR be a measurable region. The spectral projector acts as

(EX(R)ψ)(x)=1R(x)ψ(x).\bigl(E_X(R)\psi\bigr)(x) = \mathbf1_R(x)\psi(x).

Verify the intersection rule

EX(R)EX(S)=EX(R∩S)E_X(R)E_X(S) = E_X(R\cap S)

and identify the probability assigned to RR.

Solution

Indicator functions obey

1R(x)1S(x)=1R∩S(x).\mathbf1_R(x)\mathbf1_S(x) = \mathbf1_{R\cap S}(x).

Therefore multiplication by the first two indicators in succession equals multiplication by the indicator of the intersection. The Born probability is

⟨ψ∣EX(R)∣ψ⟩=∫Rψ∗(x)1R(x)ψ(x) dx=∫R∣ψ(x)∣2 dx.\begin{aligned} \langle\psi|E_X(R)|\psi\rangle &= \int_{\mathbb R} \psi^*(x)\mathbf1_R(x)\psi(x)\,dx\\ &= \int_R|\psi(x)|^2\,dx. \end{aligned}