Skip to content

Multipartite Systems

A multipartite system is a composite system with more than two named subsystems. The Hilbert space has the form

H1⋯n=H1⊗H2⊗⋯⊗Hn.\mathcal H_{1\cdots n} = \mathcal H_1\otimes\mathcal H_2 \otimes\cdots\otimes \mathcal H_n.

This looks like a straightforward extension of the bipartite tensor product, but the conceptual structure changes sharply. There is no longer only one split A∣BA\vert B. There are many ways to group subsystems into blocks, many reduced density operators, many possible notions of locality, and many inequivalent patterns of entanglement.

The purpose of this page is to set the notation and the mental model. Concrete multipartite examples such as GHZ States and W States use this language repeatedly.

For nn distinguishable subsystems, the composite Hilbert space is the tensor product of the subsystem Hilbert spaces:

H=⨂j=1nHj.\mathcal H = \bigotimes_{j=1}^{n}\mathcal H_j.

The order of factors is part of the convention. If

∣ψj⟩∈Hj,\lvert\psi_j\rangle\in\mathcal H_j,

then a fully product vector is

∣ψ1⟩⊗∣ψ2⟩⊗⋯⊗∣ψn⟩,\lvert\psi_1\rangle\otimes \lvert\psi_2\rangle\otimes\cdots\otimes \lvert\psi_n\rangle,

often abbreviated as

∣ψ1ψ2⋯ψn⟩\lvert\psi_1\psi_2\cdots\psi_n\rangle

when the subsystem ordering is clear.

For nn qubits, a standard product basis is

∣x1x2⋯xn⟩,xj∈{0,1}.\lvert x_1x_2\cdots x_n\rangle, \qquad x_j\in\{0,1\}.

For example,

∣0101⟩=∣0⟩1⊗∣1⟩2⊗∣0⟩3⊗∣1⟩4.\lvert0101\rangle = \lvert0\rangle_1\otimes \lvert1\rangle_2\otimes \lvert0\rangle_3\otimes \lvert1\rangle_4.

A general pure state is a linear combination

∣Ψ⟩=∑x1,…,xncx1⋯xn∣x1⋯xn⟩\lvert\Psi\rangle = \sum_{x_1,\ldots,x_n} c_{x_1\cdots x_n} \lvert x_1\cdots x_n\rangle

for qubits, or the analogous expansion in product bases for higher-dimensional subsystems. The coefficients form a multi-index tensor rather than a matrix. This is the first reason multipartite entanglement is more complex than bipartite entanglement: there is no single matrix singular-value decomposition that simultaneously solves all subsystem splits.

Let

N={1,2,…,n}N=\{1,2,\ldots,n\}

be the set of subsystem labels. For a subset S⊆NS\subseteq N, define

HS=⨂j∈SHj,Sˉ=N∖S.\mathcal H_S = \bigotimes_{j\in S}\mathcal H_j, \qquad \bar S=N\setminus S.

With a fixed ordering convention, the full Hilbert space can be regarded as

HN≅HS⊗HSˉ.\mathcal H_N \cong \mathcal H_S\otimes\mathcal H_{\bar S}.

An operator supported on SS has the form

OS⊗ISˉ,O_S\otimes I_{\bar S},

where OSO_S acts on HS\mathcal H_S. The identity on all other factors is often left implicit after the support has been stated.

This support language is essential in many-body physics and quantum information. A one-site observable, a two-spin correlation, a block entropy, and a stabilizer generator are all operators with specified supports.

A partition of the subsystem labels is a collection of nonempty disjoint blocks

Π={S1,S2,…,Sm}\Pi = \{S_1,S_2,\ldots,S_m\}

whose union is all of NN:

Sa∩Sb=∅(a≠b),⋃α=1mSα=N.S_a\cap S_b=\varnothing \quad(a\ne b), \qquad \bigcup_{\alpha=1}^{m}S_\alpha=N.

The partition says which groups of subsystems are being treated as the parties. For example, for three subsystems A,B,CA,B,C, the partition

{A,B,C}\{A,B,C\}

treats all three as separate parties, while

{A,BC}\{A,BC\}

treats BB and CC as one composite block.

A pure state is product with respect to the partition Π\Pi if it can be written as

∣Ψ⟩=⨂α=1m∣ψα⟩Sα.\lvert\Psi\rangle = \bigotimes_{\alpha=1}^{m} \lvert\psi_\alpha\rangle_{S_\alpha}.

The strongest product condition is full productness:

∣Ψ⟩=∣ψ1⟩⊗∣ψ2⟩⊗⋯⊗∣ψn⟩,\lvert\Psi\rangle = \lvert\psi_1\rangle\otimes \lvert\psi_2\rangle\otimes\cdots\otimes \lvert\psi_n\rangle,

which is product with respect to the singleton partition

{{1},{2},…,{n}}.\{\{1\},\{2\},\ldots,\{n\}\}.

Productness with respect to a coarse partition is weaker. The state

∣0⟩A⊗12(∣00⟩BC+∣11⟩BC)\lvert0\rangle_A\otimes \frac{1}{\sqrt2} \bigl( \lvert00\rangle_{BC} + \lvert11\rangle_{BC} \bigr)

is product across A∣BCA\vert BC, but it is not fully product because BB and CC are entangled with each other.

A bipartition is a partition into two blocks:

S∣Sˉ.S\vert\bar S.

Every nonempty proper subset S⊂NS\subset N determines a bipartition, but S∣SˉS\vert\bar S and Sˉ∣S\bar S\vert S describe the same split. Therefore the number of nontrivial bipartitions of nn labeled subsystems is

2n−1−1.2^{n-1}-1.

For n=3n=3, the three nontrivial bipartitions are

A∣BC,B∣AC,C∣AB.A\vert BC, \qquad B\vert AC, \qquad C\vert AB.

For a pure state, each bipartition reduces to an ordinary bipartite problem. One may write a Schmidt decomposition across the grouped split:

∣Ψ⟩=∑rsr(S)∣r⟩S⊗∣r⟩Sˉ,sr(S)≥0.\lvert\Psi\rangle = \sum_r s_r^{(S)} \lvert r\rangle_S \otimes \lvert r\rangle_{\bar S}, \qquad s_r^{(S)}\ge0.

Thus every bipartition has its own Schmidt rank, Schmidt coefficients, and entanglement entropy:

SS(Ψ)=−Tr⁡(ρSlog⁡ρS),ρS=Tr⁡Sˉ∣Ψ⟩⟨Ψ∣.S_S(\Psi) = - \operatorname{Tr}(\rho_S\log\rho_S), \qquad \rho_S=\operatorname{Tr}_{\bar S} \lvert\Psi\rangle\langle\Psi\rvert.

These bipartite cuts are extremely useful, but they do not fully classify multipartite entanglement. A list of one-vs-rest entropies can miss distinctions between different kinds of tripartite and many-party correlations.

For a density operator ρN\rho_N on the full system, the reduced state on a subset SS is

ρS=Tr⁡SˉρN.\rho_S = \operatorname{Tr}_{\bar S}\rho_N.

It gives all expectation values of observables supported on SS:

⟨OS⟩ρ=Tr⁡S(ρSOS).\langle O_S\rangle_\rho = \operatorname{Tr}_S(\rho_S O_S).

For three systems,

ρAB=Tr⁡CρABC,ρA=Tr⁡BCρABC.\rho_{AB} = \operatorname{Tr}_C\rho_{ABC}, \qquad \rho_A = \operatorname{Tr}_{BC}\rho_{ABC}.

The reductions are consistent under further tracing:

ρA=Tr⁡BρAB=Tr⁡BCρABC.\rho_A = \operatorname{Tr}_B\rho_{AB} = \operatorname{Tr}_{BC}\rho_{ABC}.

This consistency is useful but not complete information. A collection of one-body reduced states usually does not determine the global state. Even collections of low-body marginals can miss global phase coherence or topological constraints.

The GHZ example is the standard warning. For

∣GHZ3(ϕ)⟩=12(∣000⟩+eiϕ∣111⟩),\lvert\mathrm{GHZ}_3(\phi)\rangle = \frac{1}{\sqrt2} \bigl( \lvert000\rangle + e^{i\phi}\lvert111\rangle \bigr),

tracing out one qubit gives

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

independent of the phase ϕ\phi. The two-qubit marginal remembers perfect ZZ-basis correlation but forgets the three-body coherence.

For two subsystems, pure-state entanglement is controlled by the Schmidt coefficients. For three or more subsystems, no single list of coefficients plays the same universal role.

Several distinctions become necessary:

  • A state may be entangled across one bipartition but product across another.
  • A state may have pairwise entanglement without genuine all-party entanglement.
  • A state may have genuine multipartite entanglement even though all two-party marginals are separable.
  • Different multipartite entangled states may be inequivalent under local operations even when their bipartite Schmidt ranks look similar.
  • Mixed-state multipartite separability has a hierarchy rather than one yes-or-no test.

For example, the three-qubit state

∣0⟩A⊗∣Φ+⟩BC\lvert0\rangle_A\otimes \lvert\Phi^+\rangle_{BC}

is entangled across B∣ACB\vert AC and C∣ABC\vert AB, but product across A∣BCA\vert BC. It is not genuinely tripartite entangled.

By contrast, GHZ States have global coherence among all parties, while W States distribute a single excitation coherently over many parties. These are not merely two examples of the same bipartite phenomenon; they represent different multipartite patterns.

  • Treating a multipartite system as if it had only one natural bipartition.
  • Calling a state “entangled” without saying which partition or separability notion is being used.
  • Inferring full multipartite entanglement from one entangled pair.
  • Inferring absence of multipartite entanglement from separable two-party marginals.
  • Assuming all multipartite entanglement can be measured by one entropy.
  • Forgetting that grouping subsystems changes what counts as local.
  • D. M. Greenberger, M. A. Horne, and A. Zeilinger, “Going Beyond Bell’s Theorem,” in Bell’s Theorem, Quantum Theory, and Conceptions of the Universe, Kluwer, 1989.
  • W. Dur, G. Vidal, and J. I. Cirac, “Three Qubits Can Be Entangled in Two Inequivalent Ways,” Physical Review A 62, 062314, 2000.
  • V. Coffman, J. Kundu, and W. K. Wootters, “Distributed Entanglement,” Physical Review A 61, 052306, 2000.
  • O. Guhne and G. Toth, “Entanglement Detection,” Physics Reports 474, 1-75, 2009.
  • R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, “Quantum Entanglement,” Reviews of Modern Physics 81, 865-942, 2009.
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press, 2010.
  1. Count the nontrivial bipartitions of nn labeled subsystems.
Solution

Every nonempty proper subset S⊂NS\subset N determines a split S∣SˉS\vert\bar S. There are 2n−22^n-2 such subsets. Since S∣SˉS\vert\bar S and Sˉ∣S\bar S\vert S are the same bipartition, divide by two:

2n−22=2n−1−1.\frac{2^n-2}{2} = 2^{n-1}-1.
  1. Let
ρN=ρ1⊗ρ2⊗⋯⊗ρn.\rho_N = \rho_1\otimes\rho_2\otimes\cdots\otimes\rho_n.

Show that the reduced state on S⊆NS\subseteq N is

ρS=⨂j∈Sρj.\rho_S = \bigotimes_{j\in S}\rho_j.
Solution

Trace over each factor in Sˉ\bar S. Since every density operator has trace one,

Tr⁡Sˉ(⨂j=1nρj)=(⨂j∈Sρj)∏k∈SˉTr⁡ρk=⨂j∈Sρj.\operatorname{Tr}_{\bar S} \left( \bigotimes_{j=1}^{n}\rho_j \right) = \left( \bigotimes_{j\in S}\rho_j \right) \prod_{k\in\bar S}\operatorname{Tr}\rho_k = \bigotimes_{j\in S}\rho_j.
  1. Consider
∣Ψ⟩ABC=∣0⟩A⊗12(∣00⟩BC+∣11⟩BC).\lvert\Psi\rangle_{ABC} = \lvert0\rangle_A\otimes \frac{1}{\sqrt2} \bigl( \lvert00\rangle_{BC} + \lvert11\rangle_{BC} \bigr).

Across which of the three one-vs-two bipartitions is this state product?

Solution

It is product across A∣BCA\vert BC by inspection:

∣Ψ⟩ABC=∣0⟩A⊗∣Φ+⟩BC.\lvert\Psi\rangle_{ABC} = \lvert0\rangle_A\otimes \lvert\Phi^+\rangle_{BC}.

It is not product across B∣ACB\vert AC or C∣ABC\vert AB, because the Bell state entangles BB and CC. Grouping ACAC together does not remove the correlation between BB and CC, and similarly for grouping ABAB together.

  1. Trace out qubit CC from
∣GHZ3(ϕ)⟩=12(∣000⟩+eiϕ∣111⟩).\lvert\mathrm{GHZ}_3(\phi)\rangle = \frac{1}{\sqrt2} \bigl( \lvert000\rangle + e^{i\phi}\lvert111\rangle \bigr).

Why does the answer not depend on ϕ\phi?

Solution

The density operator contains diagonal terms and cross terms:

12∣000⟩⟨000∣+12∣111⟩⟨111∣+e−iϕ2∣000⟩⟨111∣+eiϕ2∣111⟩⟨000∣.\frac12 \lvert000\rangle\langle000\rvert + \frac12 \lvert111\rangle\langle111\rvert + \frac{e^{-i\phi}}{2} \lvert000\rangle\langle111\rvert + \frac{e^{i\phi}}{2} \lvert111\rangle\langle000\rvert.

When tracing over CC, the cross terms contain ⟨1∣0⟩C\langle1\vert0\rangle_C or ⟨0∣1⟩C\langle0\vert1\rangle_C, so they vanish. The result is

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

The phase ϕ\phi is stored in three-body coherence, not in this two-body marginal.

  1. Why is “all two-party marginals are separable” not enough to conclude that a pure three-party state has no entanglement?
Solution

The three-qubit GHZ state is the standard counterexample. Its two-party reduced states are

12∣00⟩⟨00∣+12∣11⟩⟨11∣,\frac12 \lvert00\rangle\langle00\rvert + \frac12 \lvert11\rangle\langle11\rvert,

which are separable. Nevertheless the full pure state

12(∣000⟩+∣111⟩)\frac{1}{\sqrt2} \bigl( \lvert000\rangle+\lvert111\rangle \bigr)

is entangled across every one-vs-two bipartition. Entanglement can live in genuinely multipartite coherence that disappears from pairwise marginals.