Skip to content

A GHZ state is a multipartite entangled state of the form

∣GHZn⟩=12(∣00⋯0⟩+∣11⋯1⟩)\lvert\mathrm{GHZ}_n\rangle = \frac{1}{\sqrt2} \bigl( \lvert00\cdots0\rangle + \lvert11\cdots1\rangle \bigr)

for n≥2n\ge2 qubits. The name comes from Greenberger, Horne, and Zeilinger.

For n=2n=2, this is the Bell state ∣Φ+⟩\lvert\Phi^+\rangle. For n≥3n\ge3, it is a canonical example of genuinely multipartite entanglement: the coherence is shared among all parties, but tracing out one party destroys the pure GHZ coherence among the rest.

Helpful background. Multipartite partitions explains the different cuts; conditional states distinguishes keeping a measurement outcome from discarding a subsystem.

For nn qubits with product basis

∣x1x2⋯xn⟩,xk∈{0,1},\lvert x_1x_2\cdots x_n\rangle, \qquad x_k\in\{0,1\},

the standard nn-qubit GHZ state is

∣GHZn+⟩=12(∣0⟩⊗n+∣1⟩⊗n).\lvert\mathrm{GHZ}_n^+\rangle = \frac{1}{\sqrt2} \bigl( \lvert0\rangle^{\otimes n} + \lvert1\rangle^{\otimes n} \bigr).

The minus-phase version is

∣GHZn−⟩=12(∣0⟩⊗n−∣1⟩⊗n).\lvert\mathrm{GHZ}_n^-\rangle = \frac{1}{\sqrt2} \bigl( \lvert0\rangle^{\otimes n} - \lvert1\rangle^{\otimes n} \bigr).

More generally, one may write

∣GHZn(ϕ)⟩=12(∣0⟩⊗n+eiϕ∣1⟩⊗n).\lvert\mathrm{GHZ}_n(\phi)\rangle = \frac{1}{\sqrt2} \bigl( \lvert0\rangle^{\otimes n} + e^{i\phi}\lvert1\rangle^{\otimes n} \bigr).

A local phase rotation on one qubit can change ϕ\phi, so the plus state is the standard representative unless the phase is being used as part of an interferometric or stabilizer calculation.

For any nontrivial bipartition of the nn qubits into a subset AA and its complement Aˉ\bar A, the GHZ state has Schmidt form

∣GHZn+⟩=12∣0A⟩∣0Aˉ⟩+12∣1A⟩∣1Aˉ⟩,\lvert\mathrm{GHZ}_n^+\rangle = \frac{1}{\sqrt2} \lvert0_A\rangle\lvert0_{\bar A}\rangle + \frac{1}{\sqrt2} \lvert1_A\rangle\lvert1_{\bar A}\rangle,

where

∣0A⟩=∣0⟩⊗∣A∣,∣1A⟩=∣1⟩⊗∣A∣.\lvert0_A\rangle = \lvert0\rangle^{\otimes |A|}, \qquad \lvert1_A\rangle = \lvert1\rangle^{\otimes |A|}.

Thus every nontrivial bipartition has Schmidt rank

SR⁡=2,\operatorname{SR}=2,

and entanglement entropy

SA=1S_A=1

bit.

This does not mean the state is just a collection of Bell pairs. GHZ entanglement is shared globally. For n≥3n\ge3, its two-qubit reduced states are not Bell states. Across a cut containing at least two qubits on both sides, the rank-two state is not maximally entangled on the full subsystem spaces; a larger Schmidt rank would be possible.

Let RR be a nonempty proper subset of the nn qubits. Tracing out at least one qubit gives

ρR=12∣0R⟩⟨0R∣+12∣1R⟩⟨1R∣,\rho_R = \frac12 \lvert0_R\rangle\langle0_R\rvert + \frac12 \lvert1_R\rangle\langle1_R\rvert,

where

∣0R⟩=∣0⟩⊗∣R∣,∣1R⟩=∣1⟩⊗∣R∣.\lvert0_R\rangle=\lvert0\rangle^{\otimes |R|}, \qquad \lvert1_R\rangle=\lvert1\rangle^{\otimes |R|}.

The off-diagonal terms disappear because the traced-out states

∣0Rˉ⟩and∣1Rˉ⟩\lvert0_{\bar R}\rangle \quad \text{and} \quad \lvert1_{\bar R}\rangle

are orthogonal.

For three qubits, tracing out qubit CC gives

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

This is a classically correlated separable state. It has perfect ZZ-basis correlation, but no two-qubit entanglement.

For ∣R∣=1|R|=1, this reduces to I/2I/2. The same formula holds for every phase ϕ\phi: all proper marginals are insensitive to that phase.

GHZ entanglement is fragile under loss of a subsystem. If one qubit is discarded, the remaining state is not a smaller pure GHZ state. It is the separable mixture

12∣00⋯0⟩⟨00⋯0∣+12∣11⋯1⟩⟨11⋯1∣.\frac12 \lvert00\cdots0\rangle\langle00\cdots0\rvert + \frac12 \lvert11\cdots1\rangle\langle11\cdots1\rvert.

This mixture can have strong classical correlations but no entanglement across a partition of the remaining individual qubits.

Do not confuse tracing out with measuring and conditioning. If the last qubit of ∣GHZ3+⟩\lvert\mathrm{GHZ}_3^+\rangle is measured in the XX basis and the outcome is kept, the remaining two qubits are conditionally projected into a Bell state:

∣+⟩C⟹12(∣00⟩AB+∣11⟩AB),\lvert+\rangle_C \quad \Longrightarrow \quad \frac{1}{\sqrt2} \bigl( \lvert00\rangle_{AB} + \lvert11\rangle_{AB} \bigr),

while

∣−⟩C⟹12(∣00⟩AB−∣11⟩AB).\lvert-\rangle_C \quad \Longrightarrow \quad \frac{1}{\sqrt2} \bigl( \lvert00\rangle_{AB} - \lvert11\rangle_{AB} \bigr).

Each outcome has probability 1/21/2. The nonselective average over both outcomes returns the separable mixture above. If the outcome is communicated, a local ZZ correction in the minus branch makes the final pair ∣Φ+⟩|\Phi^+\rangle in every branch. This uses existing three-party entanglement and a classical message; it does not create entanglement from a separable input.

GHZ states have perfect agreement in the computational basis. For any pair i,ji,j,

⟨ZiZj⟩GHZ=1.\langle Z_i Z_j\rangle_{\mathrm{GHZ}} = 1.

Each single-qubit ZZ expectation vanishes:

⟨Zi⟩GHZ=0.\langle Z_i\rangle_{\mathrm{GHZ}} = 0.

The global phase coherence appears in an nn-body XX correlation:

⟨X1X2⋯Xn⟩GHZ+=1,\langle X_1X_2\cdots X_n\rangle_{\mathrm{GHZ}^+} = 1,

and

⟨X1X2⋯Xn⟩GHZ−=−1.\langle X_1X_2\cdots X_n\rangle_{\mathrm{GHZ}^-} = -1.

For an arbitrary phase,

⟨X1⋯Xn⟩=cos⁡ϕ,⟨Y1X2⋯Xn⟩=sin⁡ϕ.\langle X_1\cdots X_n\rangle=\cos\phi, \qquad \langle Y_1X_2\cdots X_n\rangle=\sin\phi.

Together these correlations distinguish the relative phase that all proper marginals miss. For the plus state, a compact stabilizer description is:

X1X2⋯Xn∣GHZn+⟩=∣GHZn+⟩,X_1X_2\cdots X_n\lvert\mathrm{GHZ}_n^+\rangle = \lvert\mathrm{GHZ}_n^+\rangle,

and

Z1Zk∣GHZn+⟩=∣GHZn+⟩,k=2,…,n.Z_1Z_k\lvert\mathrm{GHZ}_n^+\rangle = \lvert\mathrm{GHZ}_n^+\rangle, \qquad k=2,\ldots,n.

The Z1ZkZ_1Z_k operators encode the fact that all computational-basis bits agree. The global X1X2⋯XnX_1X_2\cdots X_n operator encodes the coherence between the all-zero and all-one branches.

Three-party correlations and a GHZ contradiction

Section titled “Three-party correlations and a GHZ contradiction”

For ∣GHZ3+⟩|\mathrm{GHZ}_3^+\rangle, the exact Pauli correlations are

⟨XXX⟩=+1,⟨XYY⟩=⟨YXY⟩=⟨YYX⟩=−1.\langle XXX\rangle=+1,\qquad \langle XYY\rangle=\langle YXY\rangle=\langle YYX\rangle=-1.

A local model assigning a fixed ±1\pm1 value to each party’s XX and YY outcomes cannot satisfy all four. Multiplying the last three predicted products cancels each local YY value twice and requires the product of the three XX values to be −1-1, contradicting the first prediction. The quantum observables at a given site do not have simultaneous preassigned values. The ideal argument assumes the perfect correlations and locality/measurement-choice assumptions; experiments with finite noise use statistical inequalities. See entanglement and foundations.

GHZ states appear in quantum information as examples of multipartite resources:

These applications are context-dependent. A GHZ state is not automatically the best resource for every multipartite task. Its loss sensitivity is part of why W states and Graph States are studied separately.

Losing one qubit completely removes the coherence from the remaining marginal. By contrast, independent phase-flip noise E(ρ)=(1−q)ρ+qZρZ\mathcal E(\rho)=(1-q)\rho+qZ\rho Z on each qubit multiplies the full state’s off-diagonal term by (1−2q)n(1-2q)^n. The qubits remain present, but the X1⋯XnX_1\cdots X_n contrast decays rapidly with nn for fixed 0<q<1/20<q<1/2. Thus a metrological advantage for an ideal GHZ state does not by itself establish an advantage under a specified loss or noise model.

  1. Starting from ∣0⟩⊗n|0\rangle^{\otimes n}, apply a Hadamard gate to the first qubit and then controlled-NOT gates from it to every other qubit. Show that the result is ∣GHZn+⟩|\mathrm{GHZ}_n^+\rangle.
Solution

The Hadamard creates (∣0⟩+∣1⟩)∣0⋯0⟩/2(|0\rangle+|1\rangle)|0\cdots0\rangle/\sqrt2. Each controlled-NOT leaves the zero-control branch unchanged and flips one target in the one-control branch. After all n−1n-1 gates the two terms are ∣0⋯0⟩|0\cdots0\rangle and ∣1⋯1⟩|1\cdots1\rangle, with equal amplitudes. Their orthogonality confirms normalization.

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

The projector contains diagonal terms and cross terms:

12(∣000⟩⟨000∣+∣000⟩⟨111∣+∣111⟩⟨000∣+∣111⟩⟨111∣).\frac12 \Bigl( \lvert000\rangle\langle000\rvert + \lvert000\rangle\langle111\rvert + \lvert111\rangle\langle000\rvert + \lvert111\rangle\langle111\rvert \Bigr).

Tracing over CC removes the cross terms because ⟨0∣1⟩=0\langle0\vert1\rangle=0. Thus

ρAB=12∣00⟩⟨00∣+12∣11⟩⟨11∣.\rho_{AB} = \frac12 \lvert00\rangle\langle00\rvert + \frac12 \lvert11\rangle\langle11\rvert.
  1. For ∣GHZn(ϕ)⟩|\mathrm{GHZ}_n(\phi)\rangle, compute the parity distribution when all qubits are measured in XX. What happens at ϕ=π/2\phi=\pi/2?
Solution

The product PP of all ±1\pm1 outcomes has mean ⟨P⟩=cos⁡ϕ\langle P\rangle=\cos\phi. Thus Pr⁡(P=±1)=(1±cos⁡ϕ)/2\Pr(P=\pm1)=(1\pm\cos\phi)/2. At ϕ=π/2\phi=\pi/2, the two parities are equally probable, even though the state remains maximally entangled across every one-versus-rest cut. The Y1X2⋯XnY_1X_2\cdots X_n correlator is then one and reveals the phase.

  1. Find the Schmidt rank across any nontrivial bipartition of ∣GHZn+⟩\lvert\mathrm{GHZ}_n^+\rangle.
Solution

For a subset AA and complement Aˉ\bar A,

∣GHZn+⟩=12∣0A⟩∣0Aˉ⟩+12∣1A⟩∣1Aˉ⟩.\lvert\mathrm{GHZ}_n^+\rangle = \frac{1}{\sqrt2} \lvert0_A\rangle\lvert0_{\bar A}\rangle + \frac{1}{\sqrt2} \lvert1_A\rangle\lvert1_{\bar A}\rangle.

The two terms use orthonormal vectors on each side, so this is a Schmidt decomposition with two nonzero coefficients. The Schmidt rank is 22.

  1. Measure qubit CC of ∣GHZ3+⟩\lvert\mathrm{GHZ}_3^+\rangle in the XX basis. What state remains on ABAB after the ++ outcome?
Solution

Use

∣0⟩=12(∣+⟩+∣−⟩),∣1⟩=12(∣+⟩−∣−⟩).\lvert0\rangle = \frac{1}{\sqrt2} \bigl( \lvert+\rangle+\lvert-\rangle \bigr), \qquad \lvert1\rangle = \frac{1}{\sqrt2} \bigl( \lvert+\rangle-\lvert-\rangle \bigr).

Substituting for qubit CC and keeping the ∣+⟩C\lvert+\rangle_C component gives an unnormalized state proportional to

∣00⟩AB+∣11⟩AB.\lvert00\rangle_{AB} + \lvert11\rangle_{AB}.

After normalization, the conditional state is

12(∣00⟩AB+∣11⟩AB).\frac{1}{\sqrt2} \bigl( \lvert00\rangle_{AB} + \lvert11\rangle_{AB} \bigr).