Skip to content

A pure state of a composite system is a product state when it assigns a pure state to each subsystem:

∣Ψ⟩AB=∣ψ⟩A⊗∣ϕ⟩B.\lvert\Psi\rangle_{AB} = \lvert\psi\rangle_A \otimes \lvert\phi\rangle_B.

Product states describe independent pure preparations. They can contain local superpositions and other fully quantum behavior, but they contain no entanglement or correlation between AA and BB.

Helpful background. Tensor products explains how subsystem spaces combine; product bases supplies the coefficient array used to test factorization.

Let

HAB=HA⊗HB.\mathcal H_{AB} = \mathcal H_A\otimes\mathcal H_B.

A nonzero vector ∣Ψ⟩AB\lvert\Psi\rangle_{AB} is a product vector across the split A∣BA\mid B if there exist nonzero vectors

∣ψ⟩A∈HA,∣ϕ⟩B∈HB\lvert\psi\rangle_A\in\mathcal H_A, \qquad \lvert\phi\rangle_B\in\mathcal H_B

such that

∣Ψ⟩AB=∣ψ⟩A⊗∣ϕ⟩B.\lvert\Psi\rangle_{AB} = \lvert\psi\rangle_A \otimes \lvert\phi\rangle_B.

For a normalized joint state, the factors may be chosen normalized because

∥Ψ∥=∥ψ∥ ∥ϕ∥.\lVert\Psi\rVert = \lVert\psi\rVert\,\lVert\phi\rVert.

The normalized factors are unique up to compensating phases. If

∣Ψ⟩=∣ψ⟩⊗∣ϕ⟩=∣ψ′⟩⊗∣ϕ′⟩,\lvert\Psi\rangle = \lvert\psi\rangle\otimes\lvert\phi\rangle = \lvert\psi'\rangle\otimes\lvert\phi'\rangle,

then, for normalized factors,

∣ψ′⟩=eiα∣ψ⟩,∣ϕ′⟩=e−iα∣ϕ⟩\begin{aligned} \lvert\psi'\rangle &= e^{i\alpha}\lvert\psi\rangle,\\ \lvert\phi'\rangle &= e^{-i\alpha}\lvert\phi\rangle \end{aligned}

for some real α\alpha. An additional common phase changes only the representative of the physical ray.

A pure state that is not a product state across the chosen split is entangled.

Productness is not a property of a vector without further structure. It is defined relative to a factorization

H≃HA⊗HB.\mathcal H \simeq \mathcal H_A\otimes\mathcal H_B.

The same abstract Hilbert space can admit different identifications of subsystems. A vector may be product relative to one factorization and entangled relative to another.

Once the split is fixed, changing local bases does not change productness. If UAU_A and UBU_B are unitary, then

(UA⊗UB)(∣ψ⟩A⊗∣ϕ⟩B)=UA∣ψ⟩A⊗UB∣ϕ⟩B,(U_A\otimes U_B) ( \lvert\psi\rangle_A\otimes\lvert\phi\rangle_B ) = U_A\lvert\psi\rangle_A \otimes U_B\lvert\phi\rangle_B,

which is still a product state. Productness is therefore invariant under local unitary changes of basis.

A global unitary that does not factor as UA⊗UBU_A\otimes U_B can map product states to entangled states.

Choose product bases {∣i⟩A}\{\lvert i\rangle_A\} and {∣j⟩B}\{\lvert j\rangle_B\}. A general pure state is

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

If the local factors are

∣ψ⟩A=∑iai∣i⟩A,∣ϕ⟩B=∑jbj∣j⟩B,\begin{aligned} \lvert\psi\rangle_A &= \sum_i a_i\lvert i\rangle_A,\\ \lvert\phi\rangle_B &= \sum_j b_j\lvert j\rangle_B, \end{aligned}

then bilinearity gives

∣ψ⟩A⊗∣ϕ⟩B=∑i,jaibj∣i⟩A⊗∣j⟩B.\begin{aligned} \lvert\psi\rangle_A\otimes\lvert\phi\rangle_B &= \sum_{i,j} a_i b_j \lvert i\rangle_A\otimes\lvert j\rangle_B. \end{aligned}

Thus a product state has coefficients

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

Organize the coefficients into the matrix

C=(cij).C=(c_{ij}).

The state is product exactly when the nonzero matrix CC has rank one. In that case,

C=a bT,C=\boldsymbol a\,\boldsymbol b^{\mathsf T},

where a\boldsymbol a and b\boldsymbol b collect the local amplitudes.

The matrix rank does not depend on local basis choices. An active local unitary UA⊗UBU_A\otimes U_B, with the basis held fixed, changes the coefficient matrix as

C⟼UACUBT,C \longmapsto U_A C U_B^{\mathsf T},

Multiplication by invertible matrices preserves rank. For a passive basis change, the matrices appear with inverse/conjugate conventions instead; the rank is again unchanged.

For a normalized bipartite pure state, the following statements are equivalent:

  1. ∣Ψ⟩\lvert\Psi\rangle factors as ∣ψ⟩A⊗∣ϕ⟩B\lvert\psi\rangle_A\otimes\lvert\phi\rangle_B.
  2. Its coefficient matrix has rank one.
  3. Its Schmidt rank is one.
  4. Its reduced state ρA\rho_A has rank one.
  5. Its reduced state ρB\rho_B has rank one.
  6. Either reduced state is pure: Tr⁡(ρA2)=Tr⁡(ρB2)=1\operatorname{Tr}(\rho_A^2)=\operatorname{Tr}(\rho_B^2)=1.
  7. Every bounded product-observable expectation factorizes.

The coefficient test is usually fastest for small pure states. Reduced-state purity is convenient when density operators are already in use. Schmidt rank is the most structural criterion. The proofs involving reduced states and Schmidt decomposition are developed in reduced density operators and Schmidt decomposition.

These criteria apply to pure bipartite states. Mixed states require a different separability definition.

Write a general two-qubit pure state as

∣Ψ⟩=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.

Its coefficient matrix is

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

A nonzero 2×22\times2 matrix has rank one exactly when its determinant vanishes. Therefore

∣Ψ⟩ is product⟺c00c11−c01c10=0.\begin{aligned} \lvert\Psi\rangle \text{ is product} &\quad\Longleftrightarrow\\ c_{00}c_{11}-c_{01}c_{10} &=0. \end{aligned}

For example, a candidate product

(a∣0⟩+b∣1⟩)⊗(c∣0⟩+d∣1⟩)\left( a\lvert0\rangle+b\lvert1\rangle \right) \otimes \left( c\lvert0\rangle+d\lvert1\rangle \right)

expands to

ac∣00⟩+ad∣01⟩+bc∣10⟩+bd∣11⟩.\begin{aligned} ac\lvert00\rangle + ad\lvert01\rangle + bc\lvert10\rangle + bd\lvert11\rangle. \end{aligned}

The determinant condition follows:

(ac)(bd)−(ad)(bc)=0.(ac)(bd)-(ad)(bc)=0.

The test also reconstructs factors. If c00≠0c_{00}\neq0 and the determinant vanishes, then

∣Ψ⟩=(c00∣0⟩+c10∣1⟩)⊗(∣0⟩+c01c00∣1⟩).\lvert\Psi\rangle = \left( c_{00}\lvert0\rangle+c_{10}\lvert1\rangle \right) \otimes \left( \lvert0\rangle + \frac{c_{01}}{c_{00}}\lvert1\rangle \right).

Overall normalization can then be redistributed between the two factors. If c00=0c_{00}=0, choose any nonzero row or column as the pivot.

For a normalized product pure state,

ρAB=∣Ψ⟩⟨Ψ∣=ρA⊗ρB,\rho_{AB} = \lvert\Psi\rangle\langle\Psi\rvert = \rho_A\otimes\rho_B,

where

ρA=∣ψ⟩⟨ψ∣,ρB=∣ϕ⟩⟨ϕ∣.\begin{aligned} \rho_A &= \lvert\psi\rangle\langle\psi\rvert,\\ \rho_B &= \lvert\phi\rangle\langle\phi\rvert. \end{aligned}

Taking either partial trace returns the other factor:

Tr⁡B(ρA⊗ρB)=ρA,Tr⁡A(ρA⊗ρB)=ρB.\begin{aligned} \operatorname{Tr}_B(\rho_A\otimes\rho_B) &= \rho_A,\\ \operatorname{Tr}_A(\rho_A\otimes\rho_B) &= \rho_B. \end{aligned}

Each subsystem therefore has its own pure state vector. This is the defining contrast with an entangled pure state, whose reduced states are mixed.

The local rays determine the product ray without any phase convention. Choosing phases is needed only to write particular vector representatives; observable predictions depend on the phase-independent density operators.

Let EaE_a be a measurement effect on AA and FbF_b an effect on BB. In a product state,

p(a,b)=Tr⁡ ⁣[(ρA⊗ρB)(Ea⊗Fb)]=Tr⁡(ρAEa)Tr⁡(ρBFb)=pA(a)pB(b).\begin{aligned} p(a,b) &= \operatorname{Tr}\!\left[ (\rho_A\otimes\rho_B) (E_a\otimes F_b) \right]\\ &= \operatorname{Tr}(\rho_AE_a) \operatorname{Tr}(\rho_BF_b)\\ &= p_A(a)p_B(b). \end{aligned}

This factorization holds for every pair of local measurements. For bounded local observables AA and BB,

⟨A⊗B⟩=⟨A⟩A⟨B⟩B.\langle A\otimes B\rangle = \langle A\rangle_A \langle B\rangle_B.

For unbounded observables, the corresponding moments must exist. Every well-defined connected local correlation then vanishes:

⟨A⊗B⟩−⟨A⊗I⟩⟨I⊗B⟩=0.\begin{aligned} \langle A\otimes B\rangle - \langle A\otimes I\rangle \langle I\otimes B\rangle &=0. \end{aligned}

One uncorrelated measurement setting does not prove that a state is product. The criterion requires factorization for a tomographically complete set, or equivalently for all local observables.

The state

∣0⟩A⊗∣0⟩B+∣1⟩B2\lvert0\rangle_A \otimes \frac{ \lvert0\rangle_B+\lvert1\rangle_B }{\sqrt2}

is a product state. Subsystem BB is in a coherent superposition, while subsystem AA is sharp in the computational basis.

Expanding the tensor product gives

∣00⟩+∣01⟩2.\frac{ \lvert00\rangle+\lvert01\rangle }{\sqrt2}.

The presence of several product-basis terms is not an entanglement test. What matters is whether their coefficients factor into one set of amplitudes for AA and one for BB.

Similarly,

∣+⟩A⊗∣−⟩B\lvert+\rangle_A\otimes\lvert-\rangle_B

contains four computational-basis terms:

∣+⟩⊗∣−⟩=12(∣00⟩−∣01⟩+∣10⟩−∣11⟩),\lvert+\rangle\otimes\lvert-\rangle = \frac12\bigl( \lvert00\rangle-\lvert01\rangle + \lvert10\rangle-\lvert11\rangle \bigr),

yet it is manifestly product.

Local operators preserve product form whenever the output is nonzero:

(MA⊗NB)(∣ψ⟩A⊗∣ϕ⟩B)=MA∣ψ⟩A⊗NB∣ϕ⟩B.(M_A\otimes N_B) ( \lvert\psi\rangle_A\otimes\lvert\phi\rangle_B ) = M_A\lvert\psi\rangle_A \otimes N_B\lvert\phi\rangle_B.

In particular, local unitary evolution

UA(t)⊗UB(t)U_A(t)\otimes U_B(t)

cannot create entanglement from a product state. A coupling or other nonproduct operation is required to turn an initially product pure state into an entangled one.

The converse is also useful: local unitaries cannot remove entanglement from an entangled pure state. They change local bases while preserving Schmidt coefficients.

Product States Versus Product Measurements

Section titled “Product States Versus Product Measurements”

A product state is a property of a preparation:

∣ψ⟩A⊗∣ϕ⟩B.\lvert\psi\rangle_A\otimes\lvert\phi\rangle_B.

A product measurement is a property of a measurement effect:

Ea⊗Fb.E_a\otimes F_b.

These notions are independent:

  • product measurements can be performed on entangled states;
  • product states can be measured by global, nonproduct measurements;
  • product measurements on product states give factorized probabilities;
  • product measurements on nonproduct states may or may not reveal correlations in a chosen setting.

The apparatus structure does not determine whether the input state is product.

Pure Product, Mixed Product, and Separable

Section titled “Pure Product, Mixed Product, and Separable”

For density operators, a product state has the form

ρAB=ρA⊗ρB.\rho_{AB} = \rho_A\otimes\rho_B.

The factors may themselves be mixed. This is stronger than separability.

A separable state may be a correlated mixture of products:

ρAB=∑kpk ρA(k)⊗ρB(k).\rho_{AB} = \sum_k p_k\, \rho_A^{(k)} \otimes \rho_B^{(k)}.

For example,

ρcc=12∣00⟩⟨00∣+12∣11⟩⟨11∣\rho_{\mathrm{cc}} = \frac12\lvert00\rangle\langle00\rvert + \frac12\lvert11\rangle\langle11\rvert

is separable but not product. Its computational-basis outcomes are classically correlated.

Thus the pure-state dichotomy

productorentangled\text{product} \quad\text{or}\quad \text{entangled}

does not extend by replacing kets with arbitrary density operators. Mixed-state entanglement is defined by the absence of every separable decomposition. Continue to Separable Mixed States for the general definition and examples.

For NN labeled subsystems, a fully product pure state is

∣Ψ⟩=⨂r=1N∣ψr⟩.\lvert\Psi\rangle = \bigotimes_{r=1}^{N} \lvert\psi_r\rangle.

Multipartite states can have intermediate factorization structure. For example,

∣0⟩A⊗∣00⟩BC+∣11⟩BC2\lvert0\rangle_A \otimes \frac{ \lvert00\rangle_{BC}+\lvert11\rangle_{BC} }{\sqrt2}

is product across A∣BCA\mid BC but is not fully product because BB and CC are entangled.

Therefore a multipartite claim should specify the partition: fully product, product across a chosen bipartition, or entangled across that bipartition.

Product structure is not restricted to qubit labels. For position and spin,

Ψ(x)=ψ(x)χ\Psi(\mathbf x) = \psi(\mathbf x)\chi

is product across the position–spin split when the spinor χ\chi is independent of x\mathbf x. A position-dependent spinor generally does not factor this way.

For two distinguishable spinless particles on a line,

Ψ(x1,x2)=ψ(x1)ϕ(x2)\Psi(x_1,x_2) = \psi(x_1)\phi(x_2)

is product across the particle split. These examples emphasize that factorization is always relative to a specified tensor-product decomposition, not to a preferred notation or basis.

Equal basis probabilities can conceal entanglement

Section titled “Equal basis probabilities can conceal entanglement”

Compare the two normalized states

∣ψ⟩=∣00⟩+∣01⟩+∣10⟩+∣11⟩2,|\psi\rangle=\frac{|00\rangle+|01\rangle+|10\rangle+|11\rangle}{2}, ∣χ⟩=∣00⟩+∣01⟩+∣10⟩−∣11⟩2.|\chi\rangle=\frac{|00\rangle+|01\rangle+|10\rangle-|11\rangle}{2}.

Both give probability 1/41/4 to each computational-basis outcome. Yet the first coefficient matrix has determinant zero and factors as ∣+⟩∣+⟩|+\rangle|+\rangle, while the second has determinant −1/2-1/2 and is entangled. Its marginal is I/2I/2.

The difference is a relative phase, not the probabilities in this one basis. This is why an entanglement test must retain complex amplitudes or use enough measurement settings to determine the needed coherences.

Superposition alone is not entanglement. Several product-basis terms may combine into one product of local superpositions. Test their coefficient matrix, including phases.

A zero determinant is only a two-qubit shortcut. A larger square matrix may have zero determinant while having rank greater than one; a rectangular matrix has no determinant. Rank one is the general product criterion.

Mixed product and separable are different. A mixed product state has no correlations. A separable mixture can have correlations without entanglement.

  • Ingemar Bengtsson and Karol Życzkowski, Geometry of Quantum States, 2nd ed., Cambridge University Press, 2017, chapters on composite systems and entanglement.
  • Michael A. Nielsen and Isaac L. Chuang, Quantum Computation and Quantum Information, 10th anniversary ed., Cambridge University Press, 2010, section 2.5.
  • John Watrous, The Theory of Quantum Information, Cambridge University Press, 2018, chapters 1 and 6. Author-hosted text.
  1. Factor the normalized state (∣00⟩+i∣01⟩+2∣10⟩+2i∣11⟩)/10(|00\rangle+i|01\rangle+2|10\rangle+2i|11\rangle)/\sqrt{10} into normalized one-qubit states.
Solution

The second row of the coefficient matrix is twice the first. Factoring and distributing the norm gives

∣Ψ⟩=∣0⟩+2∣1⟩5⊗∣0⟩+i∣1⟩2.|\Psi\rangle=\frac{|0\rangle+2|1\rangle}{\sqrt5} \otimes\frac{|0\rangle+i|1\rangle}{\sqrt2}.

The tensor-product coefficients use an ordinary product of amplitudes, with no complex conjugation of the second factor.

  1. Use the determinant test to classify
∣Ψ⟩=∣00⟩+∣01⟩+∣10⟩+∣11⟩2.\lvert\Psi\rangle = \frac{ \lvert00\rangle+\lvert01\rangle +\lvert10\rangle+\lvert11\rangle }{2}.
Solution

The coefficient matrix is

C=12(1111).C= \frac12 \begin{pmatrix} 1&1\\ 1&1 \end{pmatrix}.

Its determinant vanishes, so it has rank one. Explicitly,

∣Ψ⟩=∣0⟩+∣1⟩2⊗∣0⟩+∣1⟩2.\lvert\Psi\rangle = \frac{ \lvert0\rangle+\lvert1\rangle }{\sqrt2} \otimes \frac{ \lvert0\rangle+\lvert1\rangle }{\sqrt2}.
  1. Determine whether
∣χ⟩=12(∣00⟩−∣01⟩+∣10⟩−∣11⟩)\lvert\chi\rangle = \frac12 \left( \lvert00\rangle-\lvert01\rangle +\lvert10\rangle-\lvert11\rangle \right)

is product, and find its local factors.

Solution

The coefficient matrix is

C=12(1−11−1).C= \frac12 \begin{pmatrix} 1&-1\\ 1&-1 \end{pmatrix}.

Its rows are identical, so it has rank one. Factoring gives

∣χ⟩=∣0⟩+∣1⟩2⊗∣0⟩−∣1⟩2=∣+⟩⊗∣−⟩.\begin{aligned} \lvert\chi\rangle &= \frac{ \lvert0\rangle+\lvert1\rangle }{\sqrt2} \otimes \frac{ \lvert0\rangle-\lvert1\rangle }{\sqrt2} \\ &= \lvert+\rangle\otimes\lvert-\rangle. \end{aligned}
  1. Let ρAB=ρA⊗ρB\rho_{AB}=\rho_A\otimes\rho_B. Prove that local measurement effects EaE_a and FbF_b have factorized joint probabilities.
Solution

The joint probability is

p(a,b)=Tr⁡ ⁣[(ρA⊗ρB)(Ea⊗Fb)]=Tr⁡ ⁣[(ρAEa)⊗(ρBFb)]=Tr⁡(ρAEa)Tr⁡(ρBFb)=pA(a)pB(b).\begin{aligned} p(a,b) &= \operatorname{Tr}\!\left[ (\rho_A\otimes\rho_B) (E_a\otimes F_b) \right]\\ &= \operatorname{Tr}\!\left[ (\rho_AE_a)\otimes(\rho_BF_b) \right]\\ &= \operatorname{Tr}(\rho_AE_a) \operatorname{Tr}(\rho_BF_b)\\ &= p_A(a)p_B(b). \end{aligned}

The third line uses the trace rule for tensor-product operators.

  1. Show that the Bell state
∣Φ+⟩=∣00⟩+∣11⟩2\lvert\Phi^+\rangle = \frac{ \lvert00\rangle+\lvert11\rangle }{\sqrt2}

is not product by computing the reduced state of subsystem AA.

Solution

The joint density operator is

ρAB=12(∣00⟩⟨00∣+∣00⟩⟨11∣+∣11⟩⟨00∣+∣11⟩⟨11∣).\rho_{AB} = \frac12\bigl( \lvert00\rangle\langle00\rvert + \lvert00\rangle\langle11\rvert + \lvert11\rangle\langle00\rvert + \lvert11\rangle\langle11\rvert \bigr).

Tracing out BB removes the cross terms because ⟨0∣1⟩=0\langle0\mid1\rangle=0:

ρA=12∣0⟩⟨0∣+12∣1⟩⟨1∣=I2.\rho_A = \frac12\lvert0\rangle\langle0\rvert + \frac12\lvert1\rangle\langle1\rvert = \frac{I}{2}.

Its purity is

Tr⁡(ρA2)=12<1.\operatorname{Tr}(\rho_A^2)=\frac12<1.

A bipartite pure state is product only if its reduced states are pure, so the Bell state is entangled.

  1. Explain why
ρcc=12∣00⟩⟨00∣+12∣11⟩⟨11∣\rho_{\mathrm{cc}} = \frac12\lvert00\rangle\langle00\rvert + \frac12\lvert11\rangle\langle11\rvert

is separable but not product.

Solution

It is separable because the displayed expression is a convex mixture of the two product states ∣00⟩\lvert00\rangle and ∣11⟩\lvert11\rangle.

Its reduced states are both I/2I/2. If it were the product of those marginals, then

ρA⊗ρB=I2⊗I2=I44,\rho_A\otimes\rho_B = \frac{I}{2}\otimes\frac{I}{2} = \frac{I_4}{4},

which assigns probability 1/41/4 to all four computational-basis outcomes. By contrast, ρcc\rho_{\mathrm{cc}} assigns probability 1/21/2 to 0000 and 1111 and zero to 0101 and 1010. It is correlated and therefore not product.

  1. Prove that a local unitary UA⊗UBU_A\otimes U_B maps every product pure state to another product pure state. Does the converse evolution create entanglement?
Solution

Apply the operator to the product:

(UA⊗UB)(∣ψ⟩A⊗∣ϕ⟩B)=UA∣ψ⟩A⊗UB∣ϕ⟩B.(U_A\otimes U_B) ( \lvert\psi\rangle_A\otimes\lvert\phi\rangle_B ) = U_A\lvert\psi\rangle_A \otimes U_B\lvert\phi\rangle_B.

The result is explicitly product. The inverse is also local:

(UA⊗UB)†=UA†⊗UB†.(U_A\otimes U_B)^\dagger = U_A^\dagger\otimes U_B^\dagger.

Therefore local unitary evolution neither creates entanglement from a product state nor removes entanglement from an entangled pure state.

  1. Consider
∣Ω⟩=∣0⟩A⊗∣00⟩BC+∣11⟩BC2.\lvert\Omega\rangle = \lvert0\rangle_A \otimes \frac{ \lvert00\rangle_{BC}+\lvert11\rangle_{BC} }{\sqrt2}.

Classify its product structure across A∣BCA\mid BC and across B∣ACB\mid AC.

Solution

Across A∣BCA\mid BC, the state is explicitly product:

∣Ω⟩=∣0⟩A⊗∣Φ+⟩BC.\lvert\Omega\rangle = \lvert0\rangle_A \otimes \lvert\Phi^+\rangle_{BC}.

Across B∣ACB\mid AC, write

∣Ω⟩=∣0⟩B∣00⟩AC+∣1⟩B∣01⟩AC2.\lvert\Omega\rangle = \frac{ \lvert0\rangle_B\lvert00\rangle_{AC} + \lvert1\rangle_B\lvert01\rangle_{AC} }{\sqrt2}.

The two ACAC states are orthogonal, so this expression has Schmidt rank two across B∣ACB\mid AC. It is entangled across that split.

The state is therefore product across one bipartition without being fully product.