Skip to content

Marginals and Correlations

Marginals are reduced descriptions of subsystems. Correlations are features of the joint state that are not fixed by those marginal descriptions alone.

For a bipartite state ρAB\rho_{AB}, the marginal states are

ρA=Tr⁡BρAB,ρB=Tr⁡AρAB.\rho_A = \operatorname{Tr}_B\rho_{AB}, \qquad \rho_B = \operatorname{Tr}_A\rho_{AB}.

They determine all local measurement statistics on AA and BB separately. They generally do not determine the joint statistics of measurements on AA and BB together.

Helpful background. Reduced density operators explains local states; local measurement statistics distinguishes paired outcomes from separate marginals.

The word “marginal” is borrowed from probability theory. A joint probability distribution p(a,b)p(a,b) has marginals

pA(a)=∑bp(a,b),pB(b)=∑ap(a,b).p_A(a) = \sum_b p(a,b), \qquad p_B(b) = \sum_a p(a,b).

In quantum mechanics, the density-operator analogue is the reduced state:

ρA=Tr⁡BρAB,ρB=Tr⁡AρAB.\rho_A = \operatorname{Tr}_B\rho_{AB}, \qquad \rho_B = \operatorname{Tr}_A\rho_{AB}.

For local measurement effects {Ea}\{E_a\} on AA and {Fb}\{F_b\} on BB, the joint probabilities are

p(a,b)=Tr⁡AB[ρAB(Ea⊗Fb)].p(a,b) = \operatorname{Tr}_{AB} \bigl[ \rho_{AB}(E_a\otimes F_b) \bigr].

The marginals are

pA(a)=Tr⁡A(ρAEa),pB(b)=Tr⁡B(ρBFb).p_A(a) = \operatorname{Tr}_A(\rho_A E_a), \qquad p_B(b) = \operatorname{Tr}_B(\rho_B F_b).

Thus ρA\rho_A and ρB\rho_B are enough for local statistics, while ρAB\rho_{AB} is needed for joint statistics.

The same pair of marginals can come from many different joint states. For example,

∣Φ+⟩=12(∣00⟩+∣11⟩)\lvert\Phi^+\rangle = \frac{1}{\sqrt2} \bigl( \lvert00\rangle+\lvert11\rangle \bigr)

and

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

both satisfy

ρA=ρB=12I.\rho_A = \rho_B = \frac12 I.

They nevertheless have different joint states and different correlations. The first is pure and entangled. The second is mixed, separable, and classically correlated.

This is the quantum version of a familiar classical fact: knowing the marginal distributions of two random variables does not determine their joint distribution.

A bipartite density operator is a product state if

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

For bounded local observables (or unbounded observables with the required moments), every product-observable expectation in a product state factorizes:

⟨MA⊗NB⟩=Tr⁡AB[(ρA⊗ρB)(MA⊗NB)]=Tr⁡A(ρAMA)Tr⁡B(ρBNB)=⟨MA⟩⟨NB⟩.\begin{aligned} \langle M_A\otimes N_B\rangle &= \operatorname{Tr}_{AB} \bigl[ (\rho_A\otimes\rho_B)(M_A\otimes N_B) \bigr]\\ &= \operatorname{Tr}_A(\rho_A M_A) \operatorname{Tr}_B(\rho_B N_B)\\ &= \langle M_A\rangle\langle N_B\rangle. \end{aligned}

Equivalently, all connected correlations vanish:

⟨MA⊗NB⟩−⟨MA⊗IB⟩⟨IA⊗NB⟩=0\langle M_A\otimes N_B\rangle - \langle M_A\otimes I_B\rangle \langle I_A\otimes N_B\rangle = 0

for all local observables MAM_A and NBN_B.

Conversely, in finite dimensions, if this factorization holds for all local operators MAM_A and NBN_B, then ρAB=ρA⊗ρB\rho_{AB}=\rho_A\otimes\rho_B. Product operators span the operator space on HA⊗HB\mathcal H_A\otimes\mathcal H_B, so all matrix elements of ρAB\rho_{AB} are fixed by those product expectations.

There is a useful exception to the usual nonuniqueness of joint states: if either marginal is pure, the joint state must be a product. For example, ρA=∣0⟩⟨0∣\rho_A=|0\rangle\langle0| forces ρAB=∣0⟩⟨0∣⊗ρB\rho_{AB}=|0\rangle\langle0|\otimes\rho_B, even if ρB\rho_B is mixed. A pure subsystem cannot remain correlated with another subsystem. Exercise 4 proves this using positivity.

A correlation function is an expectation value involving observables on multiple subsystems. For two subsystems,

CMN=⟨MA⊗NB⟩=Tr⁡AB[ρAB(MA⊗NB)].C_{MN} = \langle M_A\otimes N_B\rangle = \operatorname{Tr}_{AB} \bigl[ \rho_{AB}(M_A\otimes N_B) \bigr].

The connected correlation is

CMNconn=⟨MA⊗NB⟩−⟨MA⟩⟨NB⟩.C^{\mathrm{conn}}_{MN} = \langle M_A\otimes N_B\rangle - \langle M_A\rangle \langle N_B\rangle.

Connected correlations remove the part already explained by the one-subsystem averages. They are the quantum analogue of covariance.

Nonzero connected correlation means the joint state is not a product state. It does not, by itself, mean the state is entangled. Separable mixed states can have nonzero connected correlations.

A connected correlation removes local bias

Section titled “A connected correlation removes local bias”

Consider the separable family

ρp=p∣00⟩⟨00∣+(1−p)∣11⟩⟨11∣,0≤p≤1.\rho_p=p|00\rangle\langle00|+(1-p)|11\rangle\langle11|, \qquad 0\le p\le1.

Every run gives matching ZZ outcomes, so ⟨Z⊗Z⟩=1\langle Z\otimes Z\rangle=1 for every pp. But the local mean is 2p−12p-1, and therefore

CZZconn=1−(2p−1)2=4p(1−p).C^{\mathrm{conn}}_{ZZ}=1-(2p-1)^2=4p(1-p).

At p=1/2p=1/2 the matching outcomes are individually unpredictable and the connected correlation is 11. At p=0p=0 or 11, the outcomes are already fixed locally; the state is product and the connected correlation vanishes despite the raw correlation remaining 11.

Thus a large raw product expectation can reflect local bias. A nonzero connected value establishes correlation, but the explicit separable preparation above shows that it does not by itself establish entanglement.

For two qubits, a useful summary is the Pauli correlation matrix

Tij=Tr⁡[ρAB(σi⊗σj)],i,j∈{x,y,z}.T_{ij} = \operatorname{Tr} \bigl[ \rho_{AB}(\sigma_i\otimes\sigma_j) \bigr], \qquad i,j\in\{x,y,z\}.

The local Bloch vectors are

ri=Tr⁡[ρAB(σi⊗I)],sj=Tr⁡[ρAB(I⊗σj)].r_i = \operatorname{Tr} \bigl[ \rho_{AB}(\sigma_i\otimes I) \bigr], \qquad s_j = \operatorname{Tr} \bigl[ \rho_{AB}(I\otimes\sigma_j) \bigr].

The connected Pauli correlation matrix is

Tijconn=Tij−risj.T^{\mathrm{conn}}_{ij} = T_{ij}-r_i s_j.

These quantities reconstruct the full two-qubit state:

ρAB=14[I⊗I+∑iriσi⊗I+∑jsjI⊗σj+∑i,jTijσi⊗σj].\rho_{AB}=\frac14\left[I\otimes I +\sum_i r_i\sigma_i\otimes I +\sum_j s_jI\otimes\sigma_j +\sum_{i,j}T_{ij}\sigma_i\otimes\sigma_j\right].

The formula follows from orthogonality of the Pauli operator basis under the trace inner product. The marginals supply six real coefficients, and the nine entries of TT supply the remaining information. Arbitrary coefficient choices need not produce a positive density operator.

For states with r=s=0\mathbf r=\mathbf s=0, three examples are:

Joint state Correlation matrix TT
Uncorrelated I4/4I_4/4 00
Correlated mixture ρcc\rho_{\mathrm{cc}} diag⁡(0,0,1)\operatorname{diag}(0,0,1)
Bell projector ∣Φ+⟩⟨Φ+∣\lvert\Phi^+\rangle\langle\Phi^+\rvert diag⁡(1,−1,1)\operatorname{diag}(1,-1,1)

All three have identical marginals. Their correlation matrices differ, but no single entry is an entanglement measure. The two-correlation bound shows how a combination of entries can certify entanglement.

In finite dimensions, the quantum mutual information is

I(A:B)=S(ρA)+S(ρB)−S(ρAB),I(A:B) = S(\rho_A)+S(\rho_B)-S(\rho_{AB}),

where

S(ρ)=−Tr⁡(ρlog⁡2ρ)S(\rho) = -\operatorname{Tr}(\rho\log_2\rho)

is the von Neumann entropy in bits. In infinite dimensions, individual entropies may diverge; use the relative-entropy definition of mutual information rather than subtracting infinities.

The mutual information measures total correlation, both classical and quantum. It satisfies

I(A:B)≥0,I(A:B)\ge0,

and

I(A:B)=0if and only ifρAB=ρA⊗ρB.I(A:B)=0 \qquad \text{if and only if} \qquad \rho_{AB}=\rho_A\otimes\rho_B.

It is not the same as entanglement entropy for mixed states. For a pure bipartite state, it equals twice the subsystem entropy:

I(A:B)=2S(ρA)=2S(ρB).I(A:B) = 2S(\rho_A) = 2S(\rho_B).

For a classically correlated separable state, it can be nonzero even though there is no entanglement.

The reduced-state entropy itself is developed in Subsystem Entropy, including examples where the same local entropy comes from entanglement, classical correlation, or uncorrelated local noise. The full entropy-based total-correlation measure is developed in Mutual Information.

In many-body physics, one often studies local or few-body observables at sites, modes, or spatial regions. A two-point function has the schematic form

⟨OiOj⟩,\langle O_i O_j\rangle,

and the connected two-point function is

⟨OiOj⟩c=⟨OiOj⟩−⟨Oi⟩⟨Oj⟩.\langle O_i O_j\rangle_{\mathrm c} = \langle O_i O_j\rangle - \langle O_i\rangle\langle O_j\rangle.

This is the same idea as CMNconnC^{\mathrm{conn}}_{MN}, expressed for observables labeled by position, site, mode, or field operator. Higher nn-point functions and their connected parts organize more detailed information about the state.

In QFT and statistical mechanics, correlation functions often become the primary objects: they diagnose phases, response, clustering, critical behavior, and particle content. The finite-dimensional reduced-state language here is a useful starting point. Connected correlation functions explains the ordering and clustering assumptions; continuum field correlations also require their distributional or regularization prescriptions.

Perfectly matching fixed outcomes need not be correlated. If both outcomes are already certain separately, their connected correlation is zero.

One zero covariance does not establish a product state. The complete product criterion tests a spanning set of local operators on both subsystems.

Correlations do not imply a causal influence. A joint preparation can correlate outcomes while every unconditioned local marginal obeys no-signalling.

  • Alexander Altland and Ben Simons, Condensed Matter Field Theory, 2nd ed., Cambridge University Press, 2010, correlation functions.
  • Michael A. Nielsen and Isaac L. Chuang, Quantum Computation and Quantum Information, 10th anniversary ed., Cambridge University Press, 2010, chapters 2 and 11.
  • John Watrous, The Theory of Quantum Information, Cambridge University Press, 2018, chapters 2 and 5. Author-hosted text.
  1. Show that product density operators have factorized product-observable expectations.
Solution

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

⟨MA⊗NB⟩=Tr⁡AB[(ρA⊗ρB)(MA⊗NB)]=Tr⁡A(ρAMA)Tr⁡B(ρBNB)=⟨MA⟩⟨NB⟩.\begin{aligned} \langle M_A\otimes N_B\rangle &= \operatorname{Tr}_{AB} \bigl[ (\rho_A\otimes\rho_B)(M_A\otimes N_B) \bigr]\\ &= \operatorname{Tr}_A(\rho_A M_A) \operatorname{Tr}_B(\rho_B N_B)\\ &= \langle M_A\rangle\langle N_B\rangle. \end{aligned}
  1. For the biased mixture ρp\rho_p above with p=0.9p=0.9, compute the raw and connected ZZ correlations.
Solution

The raw correlation is 11, and both local means are 2p−1=0.82p-1=0.8. Therefore the connected value is 1−0.82=0.361-0.8^2=0.36. The decrease from the unbiased value 11 reflects that much of the agreement is already explained by the strong local bias.

  1. Reconstruct the state with r=s=0\mathbf r=\mathbf s=0 and T=diag⁡(0,0,c)T=\operatorname{diag}(0,0,c). Determine the physical range of cc and whether the state is separable.
Solution ρ=14(I⊗I+cZ⊗Z)=14diag⁡(1+c,1−c,1−c,1+c).\rho=\frac14(I\otimes I+cZ\otimes Z) =\frac14\operatorname{diag}(1+c,1-c,1-c,1+c).

Positivity requires −1≤c≤1-1\le c\le1. In this range the diagonal entries are probabilities of four computational-basis product states, so every state in the family is separable. A nonzero connected ZZZZ correlation alone cannot certify entanglement.

  1. Show that a pure marginal ρA=∣a⟩⟨a∣\rho_A=|a\rangle\langle a| forces ρAB=∣a⟩⟨a∣⊗ρB\rho_{AB}=|a\rangle\langle a|\otimes\rho_B.
Solution

Set Q=(IA−∣a⟩⟨a∣)⊗IBQ=(I_A-|a\rangle\langle a|)\otimes I_B. Then Tr⁡(QρAB)=0\operatorname{Tr}(Q\rho_{AB})=0. Positivity implies QρAB1/2=0Q\rho_{AB}^{1/2}=0, since its squared Hilbert–Schmidt norm equals that trace. Hence the support of ρAB\rho_{AB} lies inside span⁡{∣a⟩}⊗HB\operatorname{span}\{|a\rangle\}\otimes\mathcal H_B. The state must therefore be ∣a⟩⟨a∣⊗τB|a\rangle\langle a|\otimes\tau_B; tracing AA identifies τB=ρB\tau_B=\rho_B.

  1. Why does vanishing connected correlation for one observable pair not prove that a state is product?
Solution

Product structure requires factorization for all local observables, not just one chosen pair. A state may have zero covariance for MAM_A and NBN_B while still being correlated for other observables. To conclude ρAB=ρA⊗ρB\rho_{AB}=\rho_A\otimes\rho_B in finite dimensions, product-observable expectations must factorize for an operator basis on each subsystem, or equivalently for all local operators.