Skip to content

Product Bases

A product basis is the basis of a composite Hilbert space obtained by pairing basis vectors of each subsystem. If HA\mathcal H_A has basis vectors labeled by ii and HB\mathcal H_B has basis vectors labeled by jj, then a product basis of HA⊗HB\mathcal H_A\otimes\mathcal H_B is

∣i⟩A⊗∣j⟩B.\lvert i\rangle_A\otimes\lvert j\rangle_B.

Product bases are the bookkeeping layer behind coordinate arrays, matrix representations, qubit registers, two-particle wavefunctions, oscillator products, and angular-momentum addition.

Let subsystem AA have an orthonormal basis

{∣i⟩A}i∈I.\{\lvert i\rangle_A\}_{i\in I}.

In a finite-dimensional example, II may be the set {0,1,…,m−1}\{0,1,\ldots,m-1\}. For a harmonic oscillator, II may be the nonnegative integers. For a continuous position basis, the label is a continuous variable and one must treat the basis in the generalized Dirac sense rather than as a countable Hilbert-space basis.

An arbitrary state of AA may be expanded as

∣ψ⟩A=∑iai∣i⟩A.\lvert\psi\rangle_A = \sum_i a_i\lvert i\rangle_A.

Similarly, let subsystem BB have basis

{∣j⟩B}j∈J,\{\lvert j\rangle_B\}_{j\in J},

so that

∣ϕ⟩B=∑jbj∣j⟩B.\lvert\phi\rangle_B = \sum_j b_j\lvert j\rangle_B.

The labels ii and jj do different jobs: ii indexes the AA basis, while jj indexes the BB basis. Reusing one symbol for both factors is a common source of avoidable mistakes.

The product basis consists of all pairs

∣i⟩A⊗∣j⟩B,i∈I, j∈J.\lvert i\rangle_A\otimes\lvert j\rangle_B, \qquad i\in I,\ j\in J.

When the subsystem order has been declared, this may be abbreviated as

∣ij⟩AB.\lvert ij\rangle_{AB}.

Orthonormality follows from the tensor-product inner product:

AB⟨ij∣kl⟩AB=A⟨i∣k⟩A B⟨j∣l⟩B=δikδjl.{}_{AB}\langle ij\vert kl\rangle_{AB} = {}_{A}\langle i\vert k\rangle_A\, {}_{B}\langle j\vert l\rangle_B = \delta_{ik}\delta_{jl}.

Completeness is expressed as

∑i,j∣ij⟩⟨ij∣=IAB\sum_{i,j} \lvert ij\rangle\langle ij\rvert = I_{AB}

in the finite or countable orthonormal case.

Coordinate Representation of Composite States

Section titled “Coordinate Representation of Composite States”

Every vector in the composite Hilbert space can be expanded as

∣Ψ⟩=∑i,jcij ∣i⟩A⊗∣j⟩B.\lvert\Psi\rangle = \sum_{i,j} c_{ij}\, \lvert i\rangle_A\otimes\lvert j\rangle_B.

The coefficients cijc_{ij} form a two-index array. Normalization is

∑i,j∣cij∣2=1.\sum_{i,j}|c_{ij}|^2=1.

A product state has coefficients that factor:

cij=aibj.c_{ij}=a_i b_j.

If no such factorization exists, the pure state is entangled across the A∣BA|B split. For finite-dimensional bipartite pure states, this is equivalent to the coefficient matrix having rank greater than one.

For many subsystems,

H=H1⊗H2⊗⋯⊗HN,\mathcal H = \mathcal H_1\otimes\mathcal H_2\otimes\cdots\otimes\mathcal H_N,

a product basis vector is written

∣i1i2⋯iN⟩.\lvert i_1 i_2\cdots i_N\rangle.

This compact label hides an ordered tensor product:

∣i1⟩1⊗∣i2⟩2⊗⋯⊗∣iN⟩N.\lvert i_1\rangle_1 \otimes \lvert i_2\rangle_2 \otimes\cdots\otimes \lvert i_N\rangle_N.

The order must be stated. Changing order changes the coordinate representation even when the abstract Hilbert spaces are isomorphic.

For two qubits ordered as A,BA,B, this volume uses

∣00⟩,∣01⟩,∣10⟩,∣11⟩.\lvert 00\rangle,\quad \lvert 01\rangle,\quad \lvert 10\rangle,\quad \lvert 11\rangle.

Thus

∣Ψ⟩=c00∣00⟩+c01∣01⟩+c10∣10⟩+c11∣11⟩\lvert\Psi\rangle = c_{00}\lvert00\rangle +c_{01}\lvert01\rangle +c_{10}\lvert10\rangle +c_{11}\lvert11\rangle

corresponds to the coordinate column

(c00c01c10c11).\begin{pmatrix} c_{00}\\ c_{01}\\ c_{10}\\ c_{11} \end{pmatrix}.

For three qubits ordered as A,B,CA,B,C, the analogous order is

∣000⟩, ∣001⟩, ∣010⟩, ∣011⟩, ∣100⟩, ∣101⟩, ∣110⟩, ∣111⟩.\lvert000\rangle,\ \lvert001\rangle,\ \lvert010\rangle,\ \lvert011\rangle, \ \lvert100\rangle,\ \lvert101\rangle,\ \lvert110\rangle,\ \lvert111\rangle.

Different software libraries may use different endian conventions. A calculation that moves between mathematical notation and code must state the convention explicitly.

Matrix Representation of Composite Operators

Section titled “Matrix Representation of Composite Operators”

A local operator on AA has matrix elements

(OA)ik=A⟨i∣OA∣k⟩A.(O_A)_{ik} = {}_{A}\langle i\vert O_A\vert k\rangle_A.

On the composite product basis,

AB⟨ij∣(OA⊗IB)∣kl⟩AB=(OA)ikδjl.{}_{AB}\langle ij\vert (O_A\otimes I_B) \vert kl\rangle_{AB} = (O_A)_{ik}\delta_{jl}.

For a product operator,

AB⟨ij∣(A⊗B)∣kl⟩AB=AikBjl.{}_{AB}\langle ij\vert (A\otimes B) \vert kl\rangle_{AB} = A_{ik}B_{jl}.

These formulas are the index form of the tensor-product action. They are also the safest way to check whether a matrix implementation has the correct basis order.

If two distinguishable harmonic oscillators have number bases

{∣n⟩1}n=0∞,{∣m⟩2}m=0∞,\{\lvert n\rangle_1\}_{n=0}^{\infty}, \qquad \{\lvert m\rangle_2\}_{m=0}^{\infty},

then the product basis is

∣n,m⟩=∣n⟩1⊗∣m⟩2.\lvert n,m\rangle = \lvert n\rangle_1\otimes\lvert m\rangle_2.

For uncoupled oscillators,

H0=ℏω1(a1†a1+12)+ℏω2(a2†a2+12),H_0 = \hbar\omega_1 \left(a_1^\dagger a_1+\frac12\right) +\hbar\omega_2 \left(a_2^\dagger a_2+\frac12\right),

the product states are energy eigenstates. A coupling term can mix states with different pairs of occupation numbers.

For a spin-1/21/2 particle,

H=L2(R3)⊗C2.\mathcal H = L^2(\mathbb R^3)\otimes\mathbb C^2.

A useful product basis is the generalized position-spin basis

∣x⟩⊗∣↑⟩,∣x⟩⊗∣↓⟩.\lvert\mathbf x\rangle\otimes\lvert\uparrow\rangle, \qquad \lvert\mathbf x\rangle\otimes\lvert\downarrow\rangle.

A state is then represented by two wavefunction components:

Ψ(x)=(ψ↑(x)ψ↓(x)).\Psi(\mathbf x) = \begin{pmatrix} \psi_\uparrow(\mathbf x)\\ \psi_\downarrow(\mathbf x) \end{pmatrix}.

This is a product-basis expansion, not a statement that the spin and position state must factor.

For two distinguishable spinless particles on a line, the position basis is formally

∣x1⟩⊗∣x2⟩,\lvert x_1\rangle\otimes\lvert x_2\rangle,

and the wavefunction is

Ψ(x1,x2)=⟨x1,x2∣Ψ⟩.\Psi(x_1,x_2) = \langle x_1,x_2\vert\Psi\rangle.

The resolution of identity is formal:

I=∫dx1 dx2 ∣x1,x2⟩⟨x1,x2∣.I = \int dx_1\,dx_2\, \lvert x_1,x_2\rangle \langle x_1,x_2\rvert.

As usual with position kets, this belongs to Dirac notation or a rigged-Hilbert-space treatment. For ordinary calculations, the resulting square-integrable wavefunctions live in L2(R2)L^2(\mathbb R^2).

  • Reversing subsystem order while keeping the same coordinate column.
  • Confusing compact labels such as ∣01⟩\lvert01\rangle and ∣10⟩\lvert10\rangle.
  • Treating the coefficient array cijc_{ij} as basis-independent.
  • Forgetting that continuous basis labels require integrals and distributional normalization.
  • Assuming a product-basis expansion is itself evidence of a product state.
  • Using a software tensor-order convention without matching it to the written basis order.
  • P. A. M. Dirac, The Principles of Quantum Mechanics, 4th ed., Oxford University Press, 1958.
  • 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.
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press, 2010.
  1. In the default two-qubit basis order, what coordinate column represents the state
12(∣01⟩−∣10⟩)?\frac{1}{\sqrt2} \bigl( \lvert01\rangle-\lvert10\rangle \bigr)?
Solution

The default order is

∣00⟩,∣01⟩,∣10⟩,∣11⟩.\lvert00\rangle,\quad \lvert01\rangle,\quad \lvert10\rangle,\quad \lvert11\rangle.

Therefore the coordinate column is

(01/2−1/20).\begin{pmatrix} 0\\ 1/\sqrt2\\ -1/\sqrt2\\ 0 \end{pmatrix}.
  1. Show that the coefficients cij=aibjc_{ij}=a_i b_j form a rank-one matrix.
Solution

The coefficient matrix is the outer product of the column vector aa and the row vector bTb^T:

C=abT.C=ab^T.

Every column of CC is proportional to aa, with proportionality coefficient bjb_j. Hence the column space is at most one-dimensional. If both factors are nonzero, the rank is one.

  1. For a two-oscillator basis vector ∣n,m⟩\lvert n,m\rangle, which oscillator is affected by a1†a1⊗I2a_1^\dagger a_1\otimes I_2?
Solution

It acts on the first oscillator. On the basis vector,

(a1†a1⊗I2)∣n,m⟩=n∣n,m⟩.(a_1^\dagger a_1\otimes I_2) \lvert n,m\rangle = n\lvert n,m\rangle.

The second occupation number mm is unchanged.