This page collects frequently used formulas from composite systems, entanglement, identical particles, Fock space, and second-quantized notation. It is a quick reference, not a substitute for the explanatory pages linked below. For named example states, use Common Composite States . For choosing which entanglement test applies to a state class, use the Entanglement Diagnostic Table . For reproducible numerical checks, use Computational Notebooks .
Unless otherwise stated, logarithms in entropy formulas are natural logarithms.
For distinguishable subsystems A A A and B B B ,
H A B = H A ⊗ H B . \mathcal H_{AB}
=
\mathcal H_A\otimes\mathcal H_B. H A B = H A ⊗ H B .
Dimensions multiply:
dim ( H A ⊗ H B ) = dim H A dim H B . \dim(\mathcal H_A\otimes\mathcal H_B)
=
\dim\mathcal H_A\,\dim\mathcal H_B. dim ( H A ⊗ H B ) = dim H A dim H B .
If { ∣ i ⟩ A } \{\lvert i\rangle_A\} {∣ i ⟩ A } and { ∣ j ⟩ B } \{\lvert j\rangle_B\} {∣ j ⟩ B } are bases, then
{ ∣ i ⟩ A ⊗ ∣ j ⟩ B } i j \{\lvert i\rangle_A\otimes\lvert j\rangle_B\}_{ij} {∣ i ⟩ A ⊗ ∣ j ⟩ B } ij
is a product basis. A general pure state is
∣ Ψ ⟩ = ∑ i j C i j ∣ i ⟩ A ∣ j ⟩ B . \lvert\Psi\rangle
=
\sum_{ij}
C_{ij}\,
\lvert i\rangle_A\lvert j\rangle_B. ∣ Ψ ⟩ = ij ∑ C ij ∣ i ⟩ A ∣ j ⟩ B .
Local operators are embedded as
O A ↦ O A ⊗ I B , O B ↦ I A ⊗ O B . O_A
\mapsto
O_A\otimes I_B,
\qquad
O_B
\mapsto
I_A\otimes O_B. O A ↦ O A ⊗ I B , O B ↦ I A ⊗ O B .
For an uncoupled Hamiltonian,
H 0 = H A ⊗ I B + I A ⊗ H B . H_0
=
H_A\otimes I_B
+I_A\otimes H_B. H 0 = H A ⊗ I B + I A ⊗ H B .
Interactions add terms that are not local to one factor:
H = H 0 + V A B . H
=
H_0+V_{AB}. H = H 0 + V A B .
A pure product state has the form
∣ Ψ ⟩ = ∣ ψ ⟩ A ⊗ ∣ ϕ ⟩ B . \lvert\Psi\rangle
=
\lvert\psi\rangle_A\otimes\lvert\phi\rangle_B. ∣ Ψ ⟩ = ∣ ψ ⟩ A ⊗ ∣ ϕ ⟩ B .
In a product basis, this means
C i j = a i b j C_{ij}
=
a_i b_j C ij = a i b j
for some coefficient vectors a i a_i a i and b j b_j b j .
A bipartite mixed state is separable if it can be written as
ρ A B = ∑ r p r ρ A ( r ) ⊗ ρ B ( r ) , p r ≥ 0 , ∑ r p r = 1. \rho_{AB}
=
\sum_r
p_r\,
\rho_A^{(r)}\otimes\rho_B^{(r)},
\qquad
p_r\ge0,
\qquad
\sum_r p_r=1. ρ A B = r ∑ p r ρ A ( r ) ⊗ ρ B ( r ) , p r ≥ 0 , r ∑ p r = 1.
A pure bipartite state is entangled if it is not a product state. A mixed state is entangled if it is not separable.
For a pure state,
ρ = ∣ Ψ ⟩ ⟨ Ψ ∣ . \rho
=
\lvert\Psi\rangle\langle\Psi\rvert. ρ = ∣ Ψ ⟩ ⟨ Ψ ∣ .
The reduced density operators are
ρ A = Tr B ρ A B , ρ B = Tr A ρ A B . \rho_A
=
\operatorname{Tr}_B\rho_{AB},
\qquad
\rho_B
=
\operatorname{Tr}_A\rho_{AB}. ρ A = Tr B ρ A B , ρ B = Tr A ρ A B .
In a product basis,
( ρ A ) i i ′ = ∑ j ( ρ A B ) i j , i ′ j . (\rho_A)_{ii'}
=
\sum_j
(\rho_{AB})_{ij,i'j}. ( ρ A ) i i ′ = j ∑ ( ρ A B ) ij , i ′ j .
Equivalently,
Tr B ( ∣ i ⟩ ⟨ i ′ ∣ ⊗ ∣ j ⟩ ⟨ j ′ ∣ ) = ⟨ j ′ ∣ j ⟩ ∣ i ⟩ ⟨ i ′ ∣ . \operatorname{Tr}_B
\bigl(
\lvert i\rangle\langle i'\rvert
\otimes
\lvert j\rangle\langle j'\rvert
\bigr)
=
\langle j'\vert j\rangle\,
\lvert i\rangle\langle i'\rvert. Tr B ( ∣ i ⟩ ⟨ i ′ ∣ ⊗ ∣ j ⟩ ⟨ j ′ ∣ ) = ⟨ j ′ ∣ j ⟩ ∣ i ⟩ ⟨ i ′ ∣ .
Local expectation values are computed from the reduced state:
Tr A B [ ρ A B ( O A ⊗ I B ) ] = Tr A ( ρ A O A ) . \operatorname{Tr}_{AB}
\bigl[
\rho_{AB}
(O_A\otimes I_B)
\bigr]
=
\operatorname{Tr}_A(\rho_A O_A). Tr A B [ ρ A B ( O A ⊗ I B ) ] = Tr A ( ρ A O A ) .
For a product density operator,
Tr B ( ρ A ⊗ ρ B ) = ρ A Tr ρ B = ρ A . \operatorname{Tr}_B(\rho_A\otimes\rho_B)
=
\rho_A\,\operatorname{Tr}\rho_B
=
\rho_A. Tr B ( ρ A ⊗ ρ B ) = ρ A Tr ρ B = ρ A .
For a local observable M A M_A M A ,
Tr A B [ ρ A B ( M A ⊗ I B ) ] = Tr A ( ρ A M A ) . \operatorname{Tr}_{AB}
\bigl[
\rho_{AB}(M_A\otimes I_B)
\bigr]
=
\operatorname{Tr}_A(\rho_A M_A). Tr A B [ ρ A B ( M A ⊗ I B ) ] = Tr A ( ρ A M A ) .
For local POVM effects { E a } \{E_a\} { E a } on A A A ,
p ( a ) = Tr A B [ ρ A B ( E a ⊗ I B ) ] = Tr A ( ρ A E a ) . p(a)
=
\operatorname{Tr}_{AB}
\bigl[
\rho_{AB}(E_a\otimes I_B)
\bigr]
=
\operatorname{Tr}_A(\rho_A E_a). p ( a ) = Tr A B [ ρ A B ( E a ⊗ I B ) ] = Tr A ( ρ A E a ) .
For local measurements on both subsystems,
p ( a , b ) = Tr A B [ ρ A B ( E a ⊗ F b ) ] , p(a,b)
=
\operatorname{Tr}_{AB}
\bigl[
\rho_{AB}(E_a\otimes F_b)
\bigr], p ( a , b ) = Tr A B [ ρ A B ( E a ⊗ F b ) ] ,
and the marginal is
∑ b p ( a , b ) = Tr A ( ρ A E a ) . \sum_b p(a,b)
=
\operatorname{Tr}_A(\rho_A E_a). b ∑ p ( a , b ) = Tr A ( ρ A E a ) .
The marginal states are
ρ A = Tr B ρ A B , ρ B = Tr A ρ A B . \rho_A
=
\operatorname{Tr}_B\rho_{AB},
\qquad
\rho_B
=
\operatorname{Tr}_A\rho_{AB}. ρ A = Tr B ρ A B , ρ B = Tr A ρ A B .
For local observables M A M_A M A and N B N_B N B ,
C M N = ⟨ M A ⊗ N B ⟩ = Tr A B [ ρ A B ( M A ⊗ N B ) ] . C_{MN}
=
\langle M_A\otimes N_B\rangle
=
\operatorname{Tr}_{AB}
\bigl[
\rho_{AB}(M_A\otimes N_B)
\bigr]. C M N = ⟨ M A ⊗ N B ⟩ = Tr A B [ ρ A B ( M A ⊗ N B ) ] .
The connected correlation is
C M N c o n n = ⟨ M A ⊗ N B ⟩ − ⟨ M A ⟩ ⟨ N B ⟩ . C^{\mathrm{conn}}_{MN}
=
\langle M_A\otimes N_B\rangle
-
\langle M_A\rangle\langle N_B\rangle. C M N conn = ⟨ M A ⊗ N B ⟩ − ⟨ M A ⟩ ⟨ N B ⟩ .
For a product state,
ρ A B = ρ A ⊗ ρ B ⟹ C M N c o n n = 0 \rho_{AB}
=
\rho_A\otimes\rho_B
\quad
\Longrightarrow
\quad
C^{\mathrm{conn}}_{MN}=0 ρ A B = ρ A ⊗ ρ B ⟹ C M N conn = 0
for all local observables M A , N B M_A,N_B M A , N B .
The quantum mutual information is
I ( A : B ) = S ( ρ A ) + S ( ρ B ) − S ( ρ A B ) . I(A:B)
=
S(\rho_A)+S(\rho_B)-S(\rho_{AB}). I ( A : B ) = S ( ρ A ) + S ( ρ B ) − S ( ρ A B ) .
Equivalently,
I ( A : B ) = D ( ρ A B ∥ ρ A ⊗ ρ B ) . I(A:B)
=
D
\bigl(
\rho_{AB}
\Vert
\rho_A\otimes\rho_B
\bigr). I ( A : B ) = D ( ρ A B ∥ ρ A ⊗ ρ B ) .
For a projective measurement { Q b } \{Q_b\} { Q b } on subsystem B B B ,
p ( b ) = Tr A B [ ρ A B ( I A ⊗ Q b ) ] . p(b)
=
\operatorname{Tr}_{AB}
\bigl[
\rho_{AB}(I_A\otimes Q_b)
\bigr]. p ( b ) = Tr A B [ ρ A B ( I A ⊗ Q b ) ] .
The unnormalized conditional state of A A A is
ρ ~ A ∣ b = Tr B [ ( I A ⊗ Q b ) ρ A B ( I A ⊗ Q b ) ] . \widetilde\rho_{A\vert b}
=
\operatorname{Tr}_B
\bigl[
(I_A\otimes Q_b)\rho_{AB}(I_A\otimes Q_b)
\bigr]. ρ A ∣ b = Tr B [ ( I A ⊗ Q b ) ρ A B ( I A ⊗ Q b ) ] .
For p ( b ) > 0 p(b)>0 p ( b ) > 0 ,
ρ A ∣ b = ρ ~ A ∣ b p ( b ) . \rho_{A\vert b}
=
\frac{\widetilde\rho_{A\vert b}}{p(b)}. ρ A ∣ b = p ( b ) ρ A ∣ b .
The nonselective average returns the reduced state:
∑ b p ( b ) ρ A ∣ b = ρ A . \sum_b p(b)\rho_{A\vert b}
=
\rho_A. b ∑ p ( b ) ρ A ∣ b = ρ A .
If
ρ A = ∑ k = 1 r p k ∣ k ⟩ A ⟨ k ∣ A , \rho_A
=
\sum_{k=1}^{r}
p_k
\lvert k\rangle_A\langle k\rvert_A, ρ A = k = 1 ∑ r p k ∣ k ⟩ A ⟨ k ∣ A ,
then a canonical purification is
∣ Ψ ⟩ A R = ∑ k = 1 r p k ∣ k ⟩ A ∣ k ⟩ R . \lvert\Psi\rangle_{AR}
=
\sum_{k=1}^{r}
\sqrt{p_k}\,
\lvert k\rangle_A\lvert k\rangle_R. ∣ Ψ ⟩ A R = k = 1 ∑ r p k ∣ k ⟩ A ∣ k ⟩ R .
It satisfies
Tr R ( ∣ Ψ ⟩ ⟨ Ψ ∣ A R ) = ρ A . \operatorname{Tr}_R
\bigl(
\lvert\Psi\rangle\langle\Psi\rvert_{AR}
\bigr)
=
\rho_A. Tr R ( ∣ Ψ ⟩ ⟨ Ψ ∣ A R ) = ρ A .
The minimal purifying dimension is
dim H R = rank ρ A . \dim\mathcal H_R
=
\operatorname{rank}\rho_A. dim H R = rank ρ A .
Unitary transformations on the purifying system do not change ρ A \rho_A ρ A :
( I A ⊗ U R ) ∣ Ψ ⟩ A R purifies the same ρ A . (I_A\otimes U_R)\lvert\Psi\rangle_{AR}
\quad
\text{purifies the same }\rho_A. ( I A ⊗ U R ) ∣ Ψ ⟩ A R purifies the same ρ A .
The standard two-qubit Bell states are
∣ Φ ± ⟩ = 1 2 ( ∣ 00 ⟩ ± ∣ 11 ⟩ ) , \lvert\Phi^\pm\rangle
=
\frac{1}{\sqrt2}
\bigl(
\lvert00\rangle
\pm
\lvert11\rangle
\bigr), ∣ Φ ± ⟩ = 2 1 ( ∣ 00 ⟩ ± ∣ 11 ⟩ ) ,
and
∣ Ψ ± ⟩ = 1 2 ( ∣ 01 ⟩ ± ∣ 10 ⟩ ) . \lvert\Psi^\pm\rangle
=
\frac{1}{\sqrt2}
\bigl(
\lvert01\rangle
\pm
\lvert10\rangle
\bigr). ∣ Ψ ± ⟩ = 2 1 ( ∣ 01 ⟩ ± ∣ 10 ⟩ ) .
Each Bell state has maximally mixed one-qubit reductions:
ρ A = ρ B = 1 2 I . \rho_A
=
\rho_B
=
\frac12 I. ρ A = ρ B = 2 1 I .
For example,
Tr B ( ∣ Φ + ⟩ ⟨ Φ + ∣ ) = 1 2 I . \operatorname{Tr}_B
\bigl(
\lvert\Phi^+\rangle\langle\Phi^+\rvert
\bigr)
=
\frac12 I. Tr B ( ∣ Φ + ⟩ ⟨ Φ + ∣ ) = 2 1 I .
With ∣ 0 ⟩ = ∣ ↑ ⟩ \lvert0\rangle=\lvert\uparrow\rangle ∣ 0 ⟩ = ∣ ↑ ⟩ and ∣ 1 ⟩ = ∣ ↓ ⟩ \lvert1\rangle=\lvert\downarrow\rangle ∣ 1 ⟩ = ∣ ↓ ⟩ ,
∣ 1 , 1 ⟩ = ∣ ↑ ↑ ⟩ , ∣ 1 , 0 ⟩ = 1 2 ( ∣ ↑ ↓ ⟩ + ∣ ↓ ↑ ⟩ ) , ∣ 1 , − 1 ⟩ = ∣ ↓ ↓ ⟩ , ∣ 0 , 0 ⟩ = 1 2 ( ∣ ↑ ↓ ⟩ − ∣ ↓ ↑ ⟩ ) . \begin{aligned}
\lvert1,1\rangle
&=
\lvert\uparrow\uparrow\rangle,\\
\lvert1,0\rangle
&=
\frac{1}{\sqrt2}
\bigl(
\lvert\uparrow\downarrow\rangle
+
\lvert\downarrow\uparrow\rangle
\bigr),\\
\lvert1,-1\rangle
&=
\lvert\downarrow\downarrow\rangle,\\
\lvert0,0\rangle
&=
\frac{1}{\sqrt2}
\bigl(
\lvert\uparrow\downarrow\rangle
-
\lvert\downarrow\uparrow\rangle
\bigr).
\end{aligned} ∣ 1 , 1 ⟩ ∣ 1 , 0 ⟩ ∣ 1 , − 1 ⟩ ∣ 0 , 0 ⟩ = ∣ ↑↑ ⟩ , = 2 1 ( ∣ ↑↓ ⟩ + ∣ ↓↑ ⟩ ) , = ∣ ↓↓ ⟩ , = 2 1 ( ∣ ↑↓ ⟩ − ∣ ↓↑ ⟩ ) .
The triplet states are symmetric under exchange, while the singlet is antisymmetric:
P 12 ∣ 1 , m ⟩ = ∣ 1 , m ⟩ , P 12 ∣ 0 , 0 ⟩ = − ∣ 0 , 0 ⟩ . P_{12}\lvert1,m\rangle
=
\lvert1,m\rangle,
\qquad
P_{12}\lvert0,0\rangle
=
-\lvert0,0\rangle. P 12 ∣ 1 , m ⟩ = ∣ 1 , m ⟩ , P 12 ∣ 0 , 0 ⟩ = − ∣ 0 , 0 ⟩ .
The singlet correlations are isotropic:
⟨ ( σ ⋅ a ) ⊗ ( σ ⋅ b ) ⟩ s i n g l e t = − a ⋅ b . \left\langle
(\boldsymbol\sigma\cdot\mathbf a)
\otimes
(\boldsymbol\sigma\cdot\mathbf b)
\right\rangle_{\mathrm{singlet}}
=
-\mathbf a\cdot\mathbf b. ⟨ ( σ ⋅ a ) ⊗ ( σ ⋅ b ) ⟩ singlet = − a ⋅ b .
The standard n n n -qubit GHZ state is
∣ G H Z n + ⟩ = 1 2 ( ∣ 0 ⟩ ⊗ n + ∣ 1 ⟩ ⊗ n ) . \lvert\mathrm{GHZ}_n^+\rangle
=
\frac{1}{\sqrt2}
\bigl(
\lvert0\rangle^{\otimes n}
+
\lvert1\rangle^{\otimes n}
\bigr). ∣ GHZ n + ⟩ = 2 1 ( ∣ 0 ⟩ ⊗ n + ∣ 1 ⟩ ⊗ n ) .
Tracing out any nonempty proper complement of a subset R R R gives
ρ R = 1 2 ∣ 0 R ⟩ ⟨ 0 R ∣ + 1 2 ∣ 1 R ⟩ ⟨ 1 R ∣ . \rho_R
=
\frac12
\lvert0_R\rangle\langle0_R\rvert
+
\frac12
\lvert1_R\rangle\langle1_R\rvert. ρ R = 2 1 ∣ 0 R ⟩ ⟨ 0 R ∣ + 2 1 ∣ 1 R ⟩ ⟨ 1 R ∣ .
For any nontrivial bipartition A ∣ A ˉ A\vert\bar A A ∣ A ˉ ,
SR = 2 , S A = 1 bit . \operatorname{SR}=2,
\qquad
S_A=1
\quad
\text{bit}. SR = 2 , S A = 1 bit .
The standard n n n -qubit W state is
∣ W n ⟩ = 1 n ∑ k = 1 n ∣ 0 ⋯ 010 ⋯ 0 ⟩ , \lvert W_n\rangle
=
\frac{1}{\sqrt n}
\sum_{k=1}^{n}
\lvert0\cdots010\cdots0\rangle, ∣ W n ⟩ = n 1 k = 1 ∑ n ∣ 0 ⋯ 010 ⋯ 0 ⟩ ,
where the 1 1 1 is in slot k k k . A one-qubit reduction is
ρ k = n − 1 n ∣ 0 ⟩ ⟨ 0 ∣ + 1 n ∣ 1 ⟩ ⟨ 1 ∣ . \rho_k
=
\frac{n-1}{n}
\lvert0\rangle\langle0\rvert
+
\frac1n
\lvert1\rangle\langle1\rvert. ρ k = n n − 1 ∣ 0 ⟩ ⟨ 0 ∣ + n 1 ∣ 1 ⟩ ⟨ 1 ∣ .
Across a split with m m m qubits on one side,
∣ W n ⟩ = m n ∣ W m ⟩ A ∣ 0 ⟩ A ˉ ⊗ n − m + n − m n ∣ 0 ⟩ A ⊗ m ∣ W n − m ⟩ A ˉ . \lvert W_n\rangle
=
\sqrt{\frac{m}{n}}\,
\lvert W_m\rangle_A
\lvert0\rangle_{\bar A}^{\otimes n-m}
+
\sqrt{\frac{n-m}{n}}\,
\lvert0\rangle_A^{\otimes m}
\lvert W_{n-m}\rangle_{\bar A}. ∣ W n ⟩ = n m ∣ W m ⟩ A ∣ 0 ⟩ A ˉ ⊗ n − m + n n − m ∣ 0 ⟩ A ⊗ m ∣ W n − m ⟩ A ˉ .
The Schmidt probabilities are
m n , n − m n . \frac{m}{n},
\qquad
\frac{n-m}{n}. n m , n n − m .
Tracing out one qubit gives
ρ r e s t = n − 1 n ∣ W n − 1 ⟩ ⟨ W n − 1 ∣ + 1 n ∣ 0 ⟩ ⊗ n − 1 ⟨ 0 ∣ ⊗ n − 1 . \rho_{\mathrm{rest}}
=
\frac{n-1}{n}
\lvert W_{n-1}\rangle\langle W_{n-1}\rvert
+
\frac1n
\lvert0\rangle^{\otimes n-1}
\langle0\rvert^{\otimes n-1}. ρ rest = n n − 1 ∣ W n − 1 ⟩ ⟨ W n − 1 ∣ + n 1 ∣ 0 ⟩ ⊗ n − 1 ⟨ 0 ∣ ⊗ n − 1 .
Every pure bipartite state in finite-dimensional Hilbert spaces has a Schmidt decomposition
∣ Ψ ⟩ = ∑ k = 1 r λ k ∣ u k ⟩ A ∣ v k ⟩ B , \lvert\Psi\rangle
=
\sum_{k=1}^{r}
\lambda_k
\lvert u_k\rangle_A
\lvert v_k\rangle_B, ∣ Ψ ⟩ = k = 1 ∑ r λ k ∣ u k ⟩ A ∣ v k ⟩ B ,
with
λ k > 0 , ∑ k λ k 2 = 1. \lambda_k>0,
\qquad
\sum_k\lambda_k^2=1. λ k > 0 , k ∑ λ k 2 = 1.
The sets { ∣ u k ⟩ A } \{\lvert u_k\rangle_A\} {∣ u k ⟩ A } and { ∣ v k ⟩ B } \{\lvert v_k\rangle_B\} {∣ v k ⟩ B } are orthonormal, and r r r is the Schmidt rank.
Equivalently,
r = rank ρ A = rank ρ B = rank C , r
=
\operatorname{rank}\rho_A
=
\operatorname{rank}\rho_B
=
\operatorname{rank}C, r = rank ρ A = rank ρ B = rank C ,
where C C C is the coefficient matrix in a product basis.
The reduced states are
ρ A = ∑ k λ k 2 ∣ u k ⟩ ⟨ u k ∣ , ρ B = ∑ k λ k 2 ∣ v k ⟩ ⟨ v k ∣ . \rho_A
=
\sum_k
\lambda_k^2
\lvert u_k\rangle\langle u_k\rvert,
\qquad
\rho_B
=
\sum_k
\lambda_k^2
\lvert v_k\rangle\langle v_k\rvert. ρ A = k ∑ λ k 2 ∣ u k ⟩ ⟨ u k ∣ , ρ B = k ∑ λ k 2 ∣ v k ⟩ ⟨ v k ∣ .
A pure bipartite state is a product state exactly when
r = 1. r=1. r = 1.
The nonzero eigenvalues of ρ A \rho_A ρ A and ρ B \rho_B ρ B are the same:
spec > 0 ( ρ A ) = spec > 0 ( ρ B ) = { λ k 2 } . \operatorname{spec}_{>0}(\rho_A)
=
\operatorname{spec}_{>0}(\rho_B)
=
\{\lambda_k^2\}. spec > 0 ( ρ A ) = spec > 0 ( ρ B ) = { λ k 2 } .
The von Neumann entropy is
S ( ρ ) = − Tr ( ρ log ρ ) . S(\rho)
=
-\operatorname{Tr}(\rho\log\rho). S ( ρ ) = − Tr ( ρ log ρ ) .
If ρ \rho ρ has eigenvalues p α p_\alpha p α ,
S ( ρ ) = − ∑ α p α log p α , S(\rho)
=
-\sum_\alpha
p_\alpha\log p_\alpha, S ( ρ ) = − α ∑ p α log p α ,
with the convention 0 log 0 = 0 0\log0=0 0 log 0 = 0 .
For subsystem reductions,
S A = S ( ρ A ) , S B = S ( ρ B ) . S_A
=
S(\rho_A),
\qquad
S_B
=
S(\rho_B). S A = S ( ρ A ) , S B = S ( ρ B ) .
For a pure bipartite state,
S A = S ( ρ A ) , S B = S ( ρ B ) , S_A
=
S(\rho_A),
\qquad
S_B
=
S(\rho_B), S A = S ( ρ A ) , S B = S ( ρ B ) ,
and
S A = S B = − ∑ k λ k 2 log λ k 2 . S_A=S_B
=
-\sum_k
\lambda_k^2\log\lambda_k^2. S A = S B = − k ∑ λ k 2 log λ k 2 .
For a maximally entangled state of Schmidt rank d d d ,
S A = log d . S_A=\log d. S A = log d .
For mixed joint states, S A S_A S A is local mixedness rather than a standalone entanglement measure.
The Renyi entropy of order α > 0 \alpha>0 α > 0 , α ≠ 1 \alpha\ne1 α = 1 , is
S α ( ρ ) = 1 1 − α log Tr ( ρ α ) . S_\alpha(\rho)
=
\frac{1}{1-\alpha}
\log
\operatorname{Tr}(\rho^\alpha). S α ( ρ ) = 1 − α 1 log Tr ( ρ α ) .
The second Renyi entropy is
S 2 ( ρ ) = − log Tr ( ρ 2 ) . S_2(\rho)
=
-\log
\operatorname{Tr}(\rho^2). S 2 ( ρ ) = − log Tr ( ρ 2 ) .
The purity is
Tr ( ρ 2 ) . \operatorname{Tr}(\rho^2). Tr ( ρ 2 ) .
For a pure state, Tr ( ρ 2 ) = 1 \operatorname{Tr}(\rho^2)=1 Tr ( ρ 2 ) = 1 . For a mixed reduced state, Tr ( ρ A 2 ) < 1 \operatorname{Tr}(\rho_A^2)<1 Tr ( ρ A 2 ) < 1 .
For two distinguishable particles on a line,
H 12 = L 2 ( R ) ⊗ L 2 ( R ) ≅ L 2 ( R 2 ) . \mathcal H_{12}
=
L^2(\mathbb R)\otimes L^2(\mathbb R)
\cong
L^2(\mathbb R^2). H 12 = L 2 ( R ) ⊗ L 2 ( R ) ≅ L 2 ( R 2 ) .
A product wavefunction has the form
Ψ ( x 1 , x 2 ) = ψ ( x 1 ) ϕ ( x 2 ) . \Psi(x_1,x_2)
=
\psi(x_1)\phi(x_2). Ψ ( x 1 , x 2 ) = ψ ( x 1 ) ϕ ( x 2 ) .
The reduced density-kernel of particle 1 is
ρ 1 ( x , x ′ ) = ∫ − ∞ ∞ d y Ψ ( x , y ) Ψ ∗ ( x ′ , y ) . \rho_1(x,x')
=
\int_{-\infty}^{\infty}dy\,
\Psi(x,y)\Psi^*(x',y). ρ 1 ( x , x ′ ) = ∫ − ∞ ∞ d y Ψ ( x , y ) Ψ ∗ ( x ′ , y ) .
For a single oscillator mode,
q = a + a † 2 , p = a − a † i 2 , [ q , p ] = i . q
=
\frac{a+a^\dagger}{\sqrt2},
\qquad
p
=
\frac{a-a^\dagger}{i\sqrt2},
\qquad
[q,p]=i. q = 2 a + a † , p = i 2 a − a † , [ q , p ] = i .
For a unitary mode-basis change,
∣ χ α ⟩ = ∑ i U i α ∣ φ i ⟩ , b α † = ∑ i U i α a i † . \lvert\chi_\alpha\rangle
=
\sum_i U_{i\alpha}\lvert\varphi_i\rangle,
\qquad
b_\alpha^\dagger
=
\sum_i U_{i\alpha}a_i^\dagger. ∣ χ α ⟩ = i ∑ U i α ∣ φ i ⟩ , b α † = i ∑ U i α a i † .
A balanced single-excitation two-mode state is
∣ Ψ + ⟩ = 1 2 ( ∣ 1 A , 0 B ⟩ + ∣ 0 A , 1 B ⟩ ) , \lvert\Psi_+\rangle
=
\frac{1}{\sqrt2}
\left(
\lvert1_A,0_B\rangle
+
\lvert0_A,1_B\rangle
\right), ∣ Ψ + ⟩ = 2 1 ( ∣ 1 A , 0 B ⟩ + ∣ 0 A , 1 B ⟩ ) ,
with mode entropy S A = log 2 S_A=\log2 S A = log 2 .
The two-mode squeezed vacuum is
∣ T M S V ( r ) ⟩ = 1 cosh r ∑ n = 0 ∞ ( tanh r ) n ∣ n A , n B ⟩ . \lvert\mathrm{TMSV}(r)\rangle
=
\frac{1}{\cosh r}
\sum_{n=0}^{\infty}
(\tanh r)^n
\lvert n_A,n_B\rangle. ∣ TMSV ( r )⟩ = cosh r 1 n = 0 ∑ ∞ ( tanh r ) n ∣ n A , n B ⟩ .
For two identical particles, the exchange operator satisfies
P 12 2 = I . P_{12}^2=I. P 12 2 = I .
Bosonic states satisfy
P 12 ∣ Ψ ⟩ = ∣ Ψ ⟩ , P_{12}\lvert\Psi\rangle
=
\lvert\Psi\rangle, P 12 ∣ Ψ ⟩ = ∣ Ψ ⟩ ,
while fermionic states satisfy
P 12 ∣ Ψ ⟩ = − ∣ Ψ ⟩ . P_{12}\lvert\Psi\rangle
=
-\lvert\Psi\rangle. P 12 ∣ Ψ ⟩ = − ∣ Ψ ⟩ .
For two orthonormal one-particle states ∣ α ⟩ \lvert\alpha\rangle ∣ α ⟩ and ∣ β ⟩ \lvert\beta\rangle ∣ β ⟩ ,
∣ α , β ⟩ S = 1 2 ( ∣ α ⟩ 1 ∣ β ⟩ 2 + ∣ β ⟩ 1 ∣ α ⟩ 2 ) , \lvert\alpha,\beta\rangle_S
=
\frac{1}{\sqrt2}
\bigl(
\lvert\alpha\rangle_1\lvert\beta\rangle_2
+
\lvert\beta\rangle_1\lvert\alpha\rangle_2
\bigr), ∣ α , β ⟩ S = 2 1 ( ∣ α ⟩ 1 ∣ β ⟩ 2 + ∣ β ⟩ 1 ∣ α ⟩ 2 ) ,
and
∣ α , β ⟩ A = 1 2 ( ∣ α ⟩ 1 ∣ β ⟩ 2 − ∣ β ⟩ 1 ∣ α ⟩ 2 ) . \lvert\alpha,\beta\rangle_A
=
\frac{1}{\sqrt2}
\bigl(
\lvert\alpha\rangle_1\lvert\beta\rangle_2
-
\lvert\beta\rangle_1\lvert\alpha\rangle_2
\bigr). ∣ α , β ⟩ A = 2 1 ( ∣ α ⟩ 1 ∣ β ⟩ 2 − ∣ β ⟩ 1 ∣ α ⟩ 2 ) .
The N N N -particle symmetrizer and antisymmetrizer are
Π S ( N ) = 1 N ! ∑ π ∈ S N U ( π ) , \Pi_S^{(N)}
=
\frac{1}{N!}
\sum_{\pi\in S_N}
U(\pi), Π S ( N ) = N ! 1 π ∈ S N ∑ U ( π ) ,
and
Π A ( N ) = 1 N ! ∑ π ∈ S N sgn ( π ) U ( π ) . \Pi_A^{(N)}
=
\frac{1}{N!}
\sum_{\pi\in S_N}
\operatorname{sgn}(\pi)U(\pi). Π A ( N ) = N ! 1 π ∈ S N ∑ sgn ( π ) U ( π ) .
Pauli exclusion follows from antisymmetry:
∣ α , α ⟩ A = 0. \lvert\alpha,\alpha\rangle_A=0. ∣ α , α ⟩ A = 0.
If s = ⟨ α ∣ β ⟩ s=\langle\alpha\vert\beta\rangle s = ⟨ α ∣ β ⟩ is not zero, the normalized two-particle combinations are
∣ α , β ⟩ ± = ∣ α ⟩ 1 ∣ β ⟩ 2 ± ∣ β ⟩ 1 ∣ α ⟩ 2 2 ( 1 ± ∣ s ∣ 2 ) , \lvert\alpha,\beta\rangle_\pm
=
\frac{
\lvert\alpha\rangle_1\lvert\beta\rangle_2
\pm
\lvert\beta\rangle_1\lvert\alpha\rangle_2
}{
\sqrt{2(1\pm\lvert s\rvert^2)}
}, ∣ α , β ⟩ ± = 2 ( 1 ± ∣ s ∣ 2 ) ∣ α ⟩ 1 ∣ β ⟩ 2 ± ∣ β ⟩ 1 ∣ α ⟩ 2 ,
provided the denominator is nonzero.
For orthonormal spin-orbitals φ 1 , … , φ N \varphi_1,\ldots,\varphi_N φ 1 , … , φ N , the Slater determinant is
Ψ ( q 1 , … , q N ) = 1 N ! det [ φ j ( q i ) ] i , j = 1 N . \Psi(q_1,\ldots,q_N)
=
\frac{1}{\sqrt{N!}}
\det
\bigl[
\varphi_j(q_i)
\bigr]_{i,j=1}^{N}. Ψ ( q 1 , … , q N ) = N ! 1 det [ φ j ( q i ) ] i , j = 1 N .
For a factorized two-particle state
Ψ ( q 1 , q 2 ) = ψ ( x 1 , x 2 ) χ ( s 1 , s 2 ) , \Psi(q_1,q_2)
=
\psi(\mathbf x_1,\mathbf x_2)\chi(s_1,s_2), Ψ ( q 1 , q 2 ) = ψ ( x 1 , x 2 ) χ ( s 1 , s 2 ) ,
the total exchange parity is
η t o t a l = η s p a c e η s p i n . \eta_{\mathrm{total}}
=
\eta_{\mathrm{space}}\eta_{\mathrm{spin}}. η total = η space η spin .
For a one-particle Hilbert space h \mathcal h h ,
F B ( h ) = ⨁ N = 0 ∞ Sym N h , \mathcal F_B(\mathcal h)
=
\bigoplus_{N=0}^{\infty}
\operatorname{Sym}^N\mathcal h, F B ( h ) = N = 0 ⨁ ∞ Sym N h ,
and
F F ( h ) = ⨁ N = 0 ∞ ∧ N h . \mathcal F_F(\mathcal h)
=
\bigoplus_{N=0}^{\infty}
\wedge^N\mathcal h. F F ( h ) = N = 0 ⨁ ∞ ∧ N h .
The vacuum sector is
Sym 0 h ≅ ∧ 0 h ≅ C , \operatorname{Sym}^0\mathcal h
\cong
\wedge^0\mathcal h
\cong
\mathbb C, Sym 0 h ≅ ∧ 0 h ≅ C ,
with normalized vector ∣ 0 ⟩ \lvert0\rangle ∣ 0 ⟩ .
Bosonic occupations satisfy
n i = 0 , 1 , 2 , … . n_i=0,1,2,\ldots . n i = 0 , 1 , 2 , … .
Fermionic occupations satisfy
n i ∈ { 0 , 1 } . n_i\in\{0,1\}. n i ∈ { 0 , 1 } .
The total occupation is
N = ∑ i n i . N
=
\sum_i n_i. N = i ∑ n i .
Number states satisfy
N i ∣ n 1 , n 2 , … ⟩ = n i ∣ n 1 , n 2 , … ⟩ . N_i
\lvert n_1,n_2,\ldots\rangle
=
n_i
\lvert n_1,n_2,\ldots\rangle. N i ∣ n 1 , n 2 , … ⟩ = n i ∣ n 1 , n 2 , … ⟩ .
For a bosonic mode i i i ,
a i † ∣ … , n i , … ⟩ B = n i + 1 ∣ … , n i + 1 , … ⟩ B , a_i^\dagger
\lvert\ldots,n_i,\ldots\rangle_B
=
\sqrt{n_i+1}\,
\lvert\ldots,n_i+1,\ldots\rangle_B, a i † ∣ … , n i , … ⟩ B = n i + 1 ∣ … , n i + 1 , … ⟩ B ,
and
a i ∣ … , n i , … ⟩ B = n i ∣ … , n i − 1 , … ⟩ B . a_i
\lvert\ldots,n_i,\ldots\rangle_B
=
\sqrt{n_i}\,
\lvert\ldots,n_i-1,\ldots\rangle_B. a i ∣ … , n i , … ⟩ B = n i ∣ … , n i − 1 , … ⟩ B .
The canonical commutation relations are
[ a i , a j † ] = δ i j I , [ a i , a j ] = 0 , [ a i † , a j † ] = 0. [a_i,a_j^\dagger]
=
\delta_{ij}I,
\qquad
[a_i,a_j]
=0,
\qquad
[a_i^\dagger,a_j^\dagger]
=0. [ a i , a j † ] = δ ij I , [ a i , a j ] = 0 , [ a i † , a j † ] = 0.
The normalized bosonic occupation state is
∣ n 1 , n 2 , … ⟩ B = ∏ i ( a i † ) n i n i ! ∣ 0 ⟩ . \lvert n_1,n_2,\ldots\rangle_B
=
\prod_i
\frac{(a_i^\dagger)^{n_i}}{\sqrt{n_i!}}
\lvert0\rangle. ∣ n 1 , n 2 , … ⟩ B = i ∏ n i ! ( a i † ) n i ∣ 0 ⟩ .
For fermions, fix a mode ordering and define
S i = ∑ k < i n k . S_i
=
\sum_{k<i}n_k. S i = k < i ∑ n k .
Creation acts as
c i † ∣ n 1 , … , 0 i , … ⟩ F = ( − 1 ) S i ∣ n 1 , … , 1 i , … ⟩ F , c_i^\dagger
\lvert n_1,\ldots,0_i,\ldots\rangle_F
=
(-1)^{S_i}
\lvert n_1,\ldots,1_i,\ldots\rangle_F, c i † ∣ n 1 , … , 0 i , … ⟩ F = ( − 1 ) S i ∣ n 1 , … , 1 i , … ⟩ F ,
and annihilation acts as
c i ∣ n 1 , … , 1 i , … ⟩ F = ( − 1 ) S i ∣ n 1 , … , 0 i , … ⟩ F . c_i
\lvert n_1,\ldots,1_i,\ldots\rangle_F
=
(-1)^{S_i}
\lvert n_1,\ldots,0_i,\ldots\rangle_F. c i ∣ n 1 , … , 1 i , … ⟩ F = ( − 1 ) S i ∣ n 1 , … , 0 i , … ⟩ F .
If the attempted creation would double-occupy the mode, or the attempted annihilation acts on an empty mode, the result is zero.
The canonical anticommutation relations are
{ c i , c j † } = δ i j I , { c i , c j } = 0 , { c i † , c j † } = 0. \{c_i,c_j^\dagger\}
=
\delta_{ij}I,
\qquad
\{c_i,c_j\}
=0,
\qquad
\{c_i^\dagger,c_j^\dagger\}
=0. { c i , c j † } = δ ij I , { c i , c j } = 0 , { c i † , c j † } = 0.
Pauli exclusion in operator form is
( c i † ) 2 = 0. (c_i^\dagger)^2=0. ( c i † ) 2 = 0.
For bosons,
N i = a i † a i . N_i
=
a_i^\dagger a_i. N i = a i † a i .
For fermions,
N i = c i † c i . N_i
=
c_i^\dagger c_i. N i = c i † c i .
In both cases,
N i ∣ n 1 , n 2 , … ⟩ = n i ∣ n 1 , n 2 , … ⟩ . N_i
\lvert n_1,n_2,\ldots\rangle
=
n_i
\lvert n_1,n_2,\ldots\rangle. N i ∣ n 1 , n 2 , … ⟩ = n i ∣ n 1 , n 2 , … ⟩ .
The total number operator is
N = ∑ i N i . N
=
\sum_i N_i. N = i ∑ N i .
For bosons,
[ N i , a j † ] = δ i j a j † , [ N i , a j ] = − δ i j a j . [N_i,a_j^\dagger]
=
\delta_{ij}a_j^\dagger,
\qquad
[N_i,a_j]
=
-\delta_{ij}a_j. [ N i , a j † ] = δ ij a j † , [ N i , a j ] = − δ ij a j .
For fermions,
[ N i , c j † ] = δ i j c j † , [ N i , c j ] = − δ i j c j . [N_i,c_j^\dagger]
=
\delta_{ij}c_j^\dagger,
\qquad
[N_i,c_j]
=
-\delta_{ij}c_j. [ N i , c j † ] = δ ij c j † , [ N i , c j ] = − δ ij c j .
For a fermionic mode,
N i 2 = N i . N_i^2=N_i. N i 2 = N i .
Let d i , d i † d_i,d_i^\dagger d i , d i † denote bosonic or fermionic mode operators, with the appropriate algebra.
A one-body operator A A A with matrix elements
A i j = ⟨ φ i ∣ A ∣ φ j ⟩ A_{ij}
=
\langle\varphi_i\vert A\vert\varphi_j\rangle A ij = ⟨ φ i ∣ A ∣ φ j ⟩
is represented by
A ^ = ∑ i j A i j d i † d j . \widehat A
=
\sum_{ij}
A_{ij}
d_i^\dagger d_j. A = ij ∑ A ij d i † d j .
A two-body interaction with matrix elements V i j ; k l V_{ij;kl} V ij ; k l has the typical number-conserving form
V ^ = 1 2 ∑ i j k l V i j ; k l d i † d j † d l d k , \widehat V
=
\frac12
\sum_{ijkl}
V_{ij;kl}
d_i^\dagger d_j^\dagger d_l d_k, V = 2 1 ij k l ∑ V ij ; k l d i † d j † d l d k ,
with the ordering convention chosen consistently, especially for fermions.
A common many-particle Hamiltonian with one-body and two-body terms is
H = ∑ i j h i j d i † d j + 1 2 ∑ i j k l V i j ; k l d i † d j † d l d k . H
=
\sum_{ij}
h_{ij}d_i^\dagger d_j
+
\frac12
\sum_{ijkl}
V_{ij;kl}
d_i^\dagger d_j^\dagger d_l d_k. H = ij ∑ h ij d i † d j + 2 1 ij k l ∑ V ij ; k l d i † d j † d l d k .
A field operator expanded in a one-particle basis is
ψ ( x ) = ∑ i φ i ( x ) d i , ψ † ( x ) = ∑ i φ i ∗ ( x ) d i † . \psi(\mathbf x)
=
\sum_i
\varphi_i(\mathbf x)d_i,
\qquad
\psi^\dagger(\mathbf x)
=
\sum_i
\varphi_i^*(\mathbf x)d_i^\dagger. ψ ( x ) = i ∑ φ i ( x ) d i , ψ † ( x ) = i ∑ φ i ∗ ( x ) d i † .
For bosons,
[ ψ ( x ) , ψ † ( y ) ] = δ ( x − y ) . [\psi(\mathbf x),\psi^\dagger(\mathbf y)]
=
\delta(\mathbf x-\mathbf y). [ ψ ( x ) , ψ † ( y )] = δ ( x − y ) .
For fermions,
{ ψ ( x ) , ψ † ( y ) } = δ ( x − y ) . \{\psi(\mathbf x),\psi^\dagger(\mathbf y)\}
=
\delta(\mathbf x-\mathbf y). { ψ ( x ) , ψ † ( y )} = δ ( x − y ) .
For normal ordering , bosonic mode operators obey
: a i a j † : = a j † a i , a i a j † = : a i a j † : + δ i j . :a_i a_j^\dagger:
=
a_j^\dagger a_i,
\qquad
a_i a_j^\dagger
=
:a_i a_j^\dagger:
+\delta_{ij}. : a i a j † := a j † a i , a i a j † =: a i a j † : + δ ij .
For fermionic mode operators,
: c i c j † : = − c j † c i , c i c j † = : c i c j † : + δ i j . :c_i c_j^\dagger:
=
-c_j^\dagger c_i,
\qquad
c_i c_j^\dagger
=
:c_i c_j^\dagger:
+\delta_{ij}. : c i c j † := − c j † c i , c i c j † =: c i c j † : + δ ij .
With the empty-vacuum contraction convention used in Wick’s theorem preview ,
C ( a i , a j † ) = δ i j , C ( a i † , a j ) = 0 , C(a_i,a_j^\dagger)
=
\delta_{ij},
\qquad
C(a_i^\dagger,a_j)=0, C ( a i , a j † ) = δ ij , C ( a i † , a j ) = 0 ,
and similarly
C ( c i , c j † ) = δ i j , C ( c i † , c j ) = 0 , C(c_i,c_j^\dagger)
=
\delta_{ij},
\qquad
C(c_i^\dagger,c_j)=0, C ( c i , c j † ) = δ ij , C ( c i † , c j ) = 0 ,
with fermionic signs supplied by the operator swaps.
For two distinguishable particles on a line,
H 12 = L 2 ( R ) ⊗ L 2 ( R ) ≅ L 2 ( R 2 ) . \mathcal H_{12}
=
L^2(\mathbb R)\otimes L^2(\mathbb R)
\cong
L^2(\mathbb R^2). H 12 = L 2 ( R ) ⊗ L 2 ( R ) ≅ L 2 ( R 2 ) .
A position-space two-particle state is normalized by
∫ d x 1 d x 2 ∣ Ψ ( x 1 , x 2 ) ∣ 2 = 1. \int dx_1\,dx_2\,
\lvert\Psi(x_1,x_2)\rvert^2
=
1. ∫ d x 1 d x 2 ∣ Ψ ( x 1 , x 2 ) ∣ 2 = 1.
Product wavefunctions have
Ψ ( x 1 , x 2 ) = ψ ( x 1 ) ϕ ( x 2 ) . \Psi(x_1,x_2)
=
\psi(x_1)\phi(x_2). Ψ ( x 1 , x 2 ) = ψ ( x 1 ) ϕ ( x 2 ) .
For a pure two-particle state, the reduced density kernel of particle 1 is
ρ 1 ( x , x ′ ) = ∫ d y Ψ ( x , y ) Ψ ∗ ( x ′ , y ) . \rho_1(x,x')
=
\int dy\,
\Psi(x,y)\Psi^*(x',y). ρ 1 ( x , x ′ ) = ∫ d y Ψ ( x , y ) Ψ ∗ ( x ′ , y ) .
The diagonal is the marginal position density:
ρ 1 ( x , x ) = ∫ d y ∣ Ψ ( x , y ) ∣ 2 . \rho_1(x,x)
=
\int dy\,
\lvert\Psi(x,y)\rvert^2. ρ 1 ( x , x ) = ∫ d y ∣ Ψ ( x , y ) ∣ 2 .
For N N N bosonic modes, define the quadrature vector
R = ( q 1 , p 1 , … , q N , p N ) T , R
=
(q_1,p_1,\ldots,q_N,p_N)^T, R = ( q 1 , p 1 , … , q N , p N ) T ,
with
[ R j , R k ] = i Ω j k , Ω = ⨁ ℓ = 1 N ( 0 1 − 1 0 ) . [R_j,R_k]
=
i\Omega_{jk},
\qquad
\Omega
=
\bigoplus_{\ell=1}^{N}
\begin{pmatrix}
0 & 1\\
-1 & 0
\end{pmatrix}. [ R j , R k ] = i Ω j k , Ω = ℓ = 1 ⨁ N ( 0 − 1 1 0 ) .
The Gaussian covariance matrix is
V j k = 1 2 ⟨ Δ R j Δ R k + Δ R k Δ R j ⟩ , Δ R j = R j − ⟨ R j ⟩ . V_{jk}
=
\frac12
\left\langle
\Delta R_j\Delta R_k+\Delta R_k\Delta R_j
\right\rangle,
\qquad
\Delta R_j=R_j-\langle R_j\rangle. V j k = 2 1 ⟨ Δ R j Δ R k + Δ R k Δ R j ⟩ , Δ R j = R j − ⟨ R j ⟩ .
Physical covariance matrices satisfy
V + i 2 Ω ≥ 0. V+\frac{i}{2}\Omega
\ge
0. V + 2 i Ω ≥ 0.
An N N N -mode Gaussian Wigner function is
W ( R ) = 1 ( 2 π ) N det V exp [ − 1 2 ( R − d ) T V − 1 ( R − d ) ] . W(R)
=
\frac{1}{(2\pi)^N\sqrt{\det V}}
\exp\!\left[
-\frac12
(R-d)^T V^{-1}(R-d)
\right]. W ( R ) = ( 2 π ) N det V 1 exp [ − 2 1 ( R − d ) T V − 1 ( R − d ) ] .
Pure Gaussian states have symplectic eigenvalues ν j = 1 / 2 \nu_j=1/2 ν j = 1/2 and
det V = 2 − 2 N . \det V
=
2^{-2N}. det V = 2 − 2 N .
With the squeezing convention used in this volume,
S 1 ( r ) = exp [ r 2 ( a 2 − a † 2 ) ] , S_1(r)
=
\exp\!\left[
\frac r2
\left(
a^2-a^{\dagger 2}
\right)
\right], S 1 ( r ) = exp [ 2 r ( a 2 − a † 2 ) ] ,
and
S 1 † ( r ) q S 1 ( r ) = e − r q , S 1 † ( r ) p S 1 ( r ) = e r p . S_1^\dagger(r)qS_1(r)
=
e^{-r}q,
\qquad
S_1^\dagger(r)pS_1(r)
=
e^r p. S 1 † ( r ) q S 1 ( r ) = e − r q , S 1 † ( r ) p S 1 ( r ) = e r p .
The single-mode squeezed-vacuum covariance matrix is
V s q ( r ) = 1 2 ( e − 2 r 0 0 e 2 r ) . V_{\rm sq}(r)
=
\frac12
\begin{pmatrix}
e^{-2r} & 0\\
0 & e^{2r}
\end{pmatrix}. V sq ( r ) = 2 1 ( e − 2 r 0 0 e 2 r ) .
The two-mode squeeze operator is
S 2 ( r ) = exp [ r ( a A † a B † − a A a B ) ] , S_2(r)
=
\exp\!\left[
r
\left(
a_A^\dagger a_B^\dagger-a_Aa_B
\right)
\right], S 2 ( r ) = exp [ r ( a A † a B † − a A a B ) ] ,
and
∣ T M S V ( r ) ⟩ = 1 cosh r ∑ n = 0 ∞ ( tanh r ) n ∣ n A , n B ⟩ . \lvert\mathrm{TMSV}(r)\rangle
=
\frac{1}{\cosh r}
\sum_{n=0}^{\infty}
(\tanh r)^n
\lvert n_A,n_B\rangle. ∣ TMSV ( r )⟩ = cosh r 1 n = 0 ∑ ∞ ( tanh r ) n ∣ n A , n B ⟩ .
For two particles on a line,
[ X i , P j ] = i ℏ δ i j . [X_i,P_j]
=
i\hbar\,\delta_{ij}. [ X i , P j ] = i ℏ δ ij .
The EPR collective observables are
X 1 − X 2 , P 1 + P 2 , X_1-X_2,
\qquad
P_1+P_2, X 1 − X 2 , P 1 + P 2 ,
and they commute:
[ X 1 − X 2 , P 1 + P 2 ] = 0. [X_1-X_2,P_1+P_2]
=
0. [ X 1 − X 2 , P 1 + P 2 ] = 0.
A formal simultaneous generalized eigenfunction is
Ψ x 0 , P ( x 1 , x 2 ) = N exp [ i P 2 ℏ ( x 1 + x 2 ) ] δ ( x 1 − x 2 − x 0 ) . \Psi_{x_0,P}(x_1,x_2)
=
\mathcal N
\exp\!\left[
\frac{iP}{2\hbar}(x_1+x_2)
\right]
\delta(x_1-x_2-x_0). Ψ x 0 , P ( x 1 , x 2 ) = N exp [ 2ℏ i P ( x 1 + x 2 ) ] δ ( x 1 − x 2 − x 0 ) .
It obeys
( X 1 − X 2 ) Ψ x 0 , P = x 0 Ψ x 0 , P , ( P 1 + P 2 ) Ψ x 0 , P = P Ψ x 0 , P , (X_1-X_2)\Psi_{x_0,P}
=
x_0\Psi_{x_0,P},
\qquad
(P_1+P_2)\Psi_{x_0,P}
=
P\Psi_{x_0,P}, ( X 1 − X 2 ) Ψ x 0 , P = x 0 Ψ x 0 , P , ( P 1 + P 2 ) Ψ x 0 , P = P Ψ x 0 , P ,
but it is not normalizable in L 2 ( R 2 ) L^2(\mathbb R^2) L 2 ( R 2 ) .
A finite Gaussian approximation is
Ψ σ , Σ = 1 π σ Σ exp [ − ( x 1 − x 2 − x 0 ) 2 4 σ 2 − ( x 1 + x 2 ) 2 4 Σ 2 + i P 2 ℏ ( x 1 + x 2 ) ] . \Psi_{\sigma,\Sigma}
=
\frac{1}{\sqrt{\pi\sigma\Sigma}}
\exp\!\left[
-\frac{(x_1-x_2-x_0)^2}{4\sigma^2}
-\frac{(x_1+x_2)^2}{4\Sigma^2}
+\frac{iP}{2\hbar}(x_1+x_2)
\right]. Ψ σ , Σ = π σ Σ 1 exp [ − 4 σ 2 ( x 1 − x 2 − x 0 ) 2 − 4 Σ 2 ( x 1 + x 2 ) 2 + 2ℏ i P ( x 1 + x 2 ) ] .
With this convention,
Var ( X 1 − X 2 ) = σ 2 , Var ( P 1 + P 2 ) = ℏ 2 Σ 2 . \operatorname{Var}(X_1-X_2)
=
\sigma^2,
\qquad
\operatorname{Var}(P_1+P_2)
=
\frac{\hbar^2}{\Sigma^2}. Var ( X 1 − X 2 ) = σ 2 , Var ( P 1 + P 2 ) = Σ 2 ℏ 2 .
For dimensionless two-mode quadratures,
q j = a j + a j † 2 , p j = a j − a j † i 2 , q_j
=
\frac{a_j+a_j^\dagger}{\sqrt2},
\qquad
p_j
=
\frac{a_j-a_j^\dagger}{i\sqrt2}, q j = 2 a j + a j † , p j = i 2 a j − a j † ,
the two-mode squeezed vacuum satisfies, in a common phase convention,
Var ( q A − q B ) = e − 2 r , Var ( p A + p B ) = e − 2 r . \operatorname{Var}(q_A-q_B)
=
e^{-2r},
\qquad
\operatorname{Var}(p_A+p_B)
=
e^{-2r}. Var ( q A − q B ) = e − 2 r , Var ( p A + p B ) = e − 2 r .
The ideal EPR limit corresponds to singular limits such as σ → 0 \sigma\to0 σ → 0 , Σ → ∞ \Sigma\to\infty Σ → ∞ , or r → ∞ r\to\infty r → ∞ .
The standard Bell pair is
∣ Φ + ⟩ = 1 2 ( ∣ 00 ⟩ + ∣ 11 ⟩ ) . \lvert\Phi^+\rangle
=
\frac{1}{\sqrt2}
\left(
\lvert00\rangle+\lvert11\rangle
\right). ∣ Φ + ⟩ = 2 1 ( ∣ 00 ⟩ + ∣ 11 ⟩ ) .
Its one-qubit reduced states are
ρ A = ρ B = 1 2 I , \rho_A
=
\rho_B
=
\frac12 I, ρ A = ρ B = 2 1 I ,
and its entanglement entropy is
S ( ρ A ) = 1 bit . S(\rho_A)
=
1
\quad
\text{bit}. S ( ρ A ) = 1 bit .
This is one ebit of bipartite pure-state entanglement. In the quantum-information bridge , the standard resource shorthand is
1 ebit + 2 classical bits ⟶ 1 qubit transmission 1\ \text{ebit}
+
2\ \text{classical bits}
\longrightarrow
1\ \text{qubit transmission} 1 ebit + 2 classical bits ⟶ 1 qubit transmission
for teleportation, and
1 ebit + 1 qubit transmission ⟶ 2 classical bits 1\ \text{ebit}
+
1\ \text{qubit transmission}
\longrightarrow
2\ \text{classical bits} 1 ebit + 1 qubit transmission ⟶ 2 classical bits
for superdense coding. These are protocol resource statements, not local signaling mechanisms.
For the Bell state
∣ Φ + ⟩ = ∣ 00 ⟩ + ∣ 11 ⟩ 2 , \lvert\Phi^+\rangle
=
\frac{
\lvert00\rangle+\lvert11\rangle
}{\sqrt2}, ∣ Φ + ⟩ = 2 ∣ 00 ⟩ + ∣ 11 ⟩ ,
the one-qubit reduced states are
ρ A = ρ B = 1 2 I . \rho_A
=
\rho_B
=
\frac12 I. ρ A = ρ B = 2 1 I .
Local outcome probabilities are computed from the reduced state:
p ( a ) = Tr A ( ρ A E a ) . p(a)
=
\operatorname{Tr}_A(\rho_A E_a). p ( a ) = Tr A ( ρ A E a ) .
An ordinary trace-preserving operation on the remote subsystem cannot change the unconditioned local reduced state:
Tr B [ ( I A ⊗ E B ) ( ρ A B ) ] = ρ A . \operatorname{Tr}_B
\left[
(I_A\otimes\mathcal E_B)(\rho_{AB})
\right]
=
\rho_A. Tr B [ ( I A ⊗ E B ) ( ρ A B ) ] = ρ A .
Conditional states can change after a remote outcome is known:
ρ A ∣ b = σ A ∣ b p ( b ) . \rho_{A\vert b}
=
\frac{\sigma_{A\vert b}}{p(b)}. ρ A ∣ b = p ( b ) σ A ∣ b .
A standard CHSH expression has the local-hidden-variable bound
∣ E ( a , b ) + E ( a , b ′ ) + E ( a ′ , b ) − E ( a ′ , b ′ ) ∣ ≤ 2 , \left|
E(a,b)
+
E(a,b')
+
E(a',b)
-
E(a',b')
\right|
\le
2, ∣ E ( a , b ) + E ( a , b ′ ) + E ( a ′ , b ) − E ( a ′ , b ′ ) ∣ ≤ 2 ,
while suitable quantum states and measurements can reach 2 2 2\sqrt2 2 2 . See the foundations bridge for the distinction between entanglement, steering, Bell nonlocality, no-signaling, and measurement update.
A low-gain pair source can be written schematically as
∣ Ψ ⟩ ≈ ∣ 0 ⟩ + ϵ ∑ μ , ν f μ ν a μ † b ν † ∣ 0 ⟩ + O ( ϵ 2 ) . \lvert\Psi\rangle
\approx
\lvert0\rangle
+
\epsilon
\sum_{\mu,\nu}
f_{\mu\nu}\,
a_\mu^\dagger b_\nu^\dagger
\lvert0\rangle
+
O(\epsilon^2). ∣ Ψ ⟩ ≈ ∣ 0 ⟩ + ϵ μ , ν ∑ f μν a μ † b ν † ∣ 0 ⟩ + O ( ϵ 2 ) .
After selecting the one-pair component, signal-idler entanglement is present when the normalized amplitude F μ ν F_{\mu\nu} F μν is not factorizable:
∣ ψ 2 ⟩ = ∑ μ , ν F μ ν ∣ 1 μ ⟩ s ∣ 1 ν ⟩ i . \lvert\psi_2\rangle
=
\sum_{\mu,\nu}
F_{\mu\nu}
\lvert1_\mu\rangle_s
\lvert1_\nu\rangle_i. ∣ ψ 2 ⟩ = μ , ν ∑ F μν ∣ 1 μ ⟩ s ∣ 1 ν ⟩ i .
A polarization Bell-like state is
∣ Φ ϕ ⟩ p o l = ∣ H ⟩ s ∣ H ⟩ i + e i ϕ ∣ V ⟩ s ∣ V ⟩ i 2 . \lvert\Phi_\phi\rangle_{\mathrm{pol}}
=
\frac{
\lvert H\rangle_s\lvert H\rangle_i
+
e^{i\phi}
\lvert V\rangle_s\lvert V\rangle_i
}{\sqrt2}. ∣ Φ ϕ ⟩ pol = 2 ∣ H ⟩ s ∣ H ⟩ i + e i ϕ ∣ V ⟩ s ∣ V ⟩ i .
A single-photon path-entangled state is
∣ ψ ⟩ = ∣ 1 ⟩ A ∣ 0 ⟩ B + e i ϕ ∣ 0 ⟩ A ∣ 1 ⟩ B 2 . \lvert\psi\rangle
=
\frac{
\lvert1\rangle_A\lvert0\rangle_B
+
e^{i\phi}
\lvert0\rangle_A\lvert1\rangle_B
}{\sqrt2}. ∣ ψ ⟩ = 2 ∣ 1 ⟩ A ∣ 0 ⟩ B + e i ϕ ∣ 0 ⟩ A ∣ 1 ⟩ B .
Homodyne detection measures phase-dependent quadratures such as
q θ = a e − i θ + a † e i θ 2 . q_\theta
=
\frac{
a e^{-i\theta}
+
a^\dagger e^{i\theta}
}{\sqrt2}. q θ = 2 a e − i θ + a † e i θ .
For the two-mode squeezed vacuum convention used above,
Var ( q A − q B ) = e − 2 r , Var ( p A + p B ) = e − 2 r . \operatorname{Var}(q_A-q_B)
=
e^{-2r},
\qquad
\operatorname{Var}(p_A+p_B)
=
e^{-2r}. Var ( q A − q B ) = e − 2 r , Var ( p A + p B ) = e − 2 r .
See the quantum-optics bridge for mode labels, down-conversion cautions, polarization and path entanglement, and homodyne interpretation.
For a spatial, site, orbital, or mode region A A A ,
H ≅ H A ⊗ H A ˉ , ρ A = Tr A ˉ ρ . \mathcal H
\cong
\mathcal H_A\otimes\mathcal H_{\bar A},
\qquad
\rho_A
=
\operatorname{Tr}_{\bar A}\rho. H ≅ H A ⊗ H A ˉ , ρ A = Tr A ˉ ρ .
For a pure global state, the entanglement entropy across the A ∣ A ˉ A\vert\bar A A ∣ A ˉ cut is
S A = − Tr ( ρ A log ρ A ) . S_A
=
-\operatorname{Tr}(\rho_A\log\rho_A). S A = − Tr ( ρ A log ρ A ) .
In the many-body bridge , common scaling shorthands are
S A ∼ α ∣ ∂ A ∣ area law , S_A
\sim
\alpha\,\lvert\partial A\rvert
\quad
\text{area law}, S A ∼ α ∣ ∂ A ∣ area law ,
and
S A ∝ ∣ A ∣ volume law . S_A
\propto
\lvert A\rvert
\quad
\text{volume law}. S A ∝ ∣ A ∣ volume law .
The entanglement Hamiltonian is defined, up to additive constants, by
ρ A = e − H E Z E , Z E = Tr ( e − H E ) . \rho_A
=
\frac{e^{-H_E}}{Z_E},
\qquad
Z_E
=
\operatorname{Tr}(e^{-H_E}). ρ A = Z E e − H E , Z E = Tr ( e − H E ) .
For a matrix product state with bond dimension χ \chi χ across a cut,
S A ≤ log χ . S_A
\le
\log\chi. S A ≤ log χ .
For two-dimensional topologically ordered gapped phases, a common preview form is
S ( A ) = α ∣ ∂ A ∣ − γ + ⋯ , γ = log D S(A)
=
\alpha\,\lvert\partial A\rvert
-
\gamma
+
\cdots,
\qquad
\gamma
=
\log\mathcal D S ( A ) = α ∣ ∂ A ∣ − γ + ⋯ , γ = log D
in standard cases, where D \mathcal D D is the total quantum dimension. Extracting γ \gamma γ requires region combinations that cancel nonuniversal boundary terms.
For a regulated spatial region A A A in a field theory,
ρ A = Tr A ˉ ∣ 0 ⟩ ⟨ 0 ∣ , S A = − Tr ( ρ A log ρ A ) . \rho_A
=
\operatorname{Tr}_{\bar A}
\lvert0\rangle\langle0\rvert,
\qquad
S_A
=
-\operatorname{Tr}(\rho_A\log\rho_A). ρ A = Tr A ˉ ∣ 0 ⟩ ⟨ 0 ∣ , S A = − Tr ( ρ A log ρ A ) .
In d > 1 d>1 d > 1 spatial dimensions, a common leading cutoff behavior is
S A ∼ α Area ( ∂ A ) ϵ d − 1 + ⋯ . S_A
\sim
\alpha\,
\frac{\operatorname{Area}(\partial A)}{\epsilon^{d-1}}
+
\cdots. S A ∼ α ϵ d − 1 Area ( ∂ A ) + ⋯ .
For a one-dimensional conformal field theory interval of length ℓ \ell ℓ ,
S A = c 3 log ℓ ϵ + c o n s t . S_A
=
\frac{c}{3}
\log\frac{\ell}{\epsilon}
+
\mathrm{const}. S A = 3 c log ϵ ℓ + const .
The modular Hamiltonian shorthand is
K A = − log ρ A , ρ A = e − K A , K_A
=
-\log\rho_A,
\qquad
\rho_A
=
e^{-K_A}, K A = − log ρ A , ρ A = e − K A ,
up to additive constants if a separate normalization is used. Relative entropy can be written as
D ( ρ ∥ σ ) = Δ ⟨ K σ ⟩ − Δ S ≥ 0 , K σ = − log σ . D(\rho\Vert\sigma)
=
\Delta\langle K_\sigma\rangle
-
\Delta S
\ge
0,
\qquad
K_\sigma=-\log\sigma. D ( ρ ∥ σ ) = Δ ⟨ K σ ⟩ − Δ S ≥ 0 , K σ = − log σ .
In holographic theories, the simplest static classical relation is
S A = Area ( γ A ) 4 G N . S_A
=
\frac{\operatorname{Area}(\gamma_A)}{4G_N}. S A = 4 G N Area ( γ A ) .
These are preview shorthands. See Entanglement in QFT Preview for regulator, region, algebra, modular-Hamiltonian, and holography cautions.
In an orthonormal spin-orbital basis, an N N N -electron determinant is
∣ Φ ⟩ = a i 1 † a i 2 † ⋯ a i N † ∣ 0 ⟩ . \lvert\Phi\rangle
=
a_{i_1}^\dagger
a_{i_2}^\dagger
\cdots
a_{i_N}^\dagger
\lvert0\rangle. ∣ Φ ⟩ = a i 1 † a i 2 † ⋯ a i N † ∣ 0 ⟩ .
A correlated finite-basis state can be expanded as
∣ Ψ ⟩ = ∑ I C I ∣ Φ I ⟩ . \lvert\Psi\rangle
=
\sum_I C_I\lvert\Phi_I\rangle. ∣ Ψ ⟩ = I ∑ C I ∣ Φ I ⟩ .
Relative to a reference determinant ∣ Φ 0 ⟩ \lvert\Phi_0\rangle ∣ Φ 0 ⟩ , common excitation notation is
∣ Φ i a ⟩ = a a † a i ∣ Φ 0 ⟩ , ∣ Φ i j a b ⟩ = a a † a b † a j a i ∣ Φ 0 ⟩ . \lvert\Phi_i^a\rangle
=
a_a^\dagger a_i
\lvert\Phi_0\rangle,
\qquad
\lvert\Phi_{ij}^{ab}\rangle
=
a_a^\dagger a_b^\dagger
a_j a_i
\lvert\Phi_0\rangle. ∣ Φ i a ⟩ = a a † a i ∣ Φ 0 ⟩ , ∣ Φ ij ab ⟩ = a a † a b † a j a i ∣ Φ 0 ⟩ .
For a spatial orbital p p p with local basis
∣ 0 ⟩ , ∣ ↑ ⟩ , ∣ ↓ ⟩ , ∣ ↑ ↓ ⟩ , \lvert0\rangle,
\qquad
\lvert\uparrow\rangle,
\qquad
\lvert\downarrow\rangle,
\qquad
\lvert\uparrow\downarrow\rangle, ∣ 0 ⟩ , ∣ ↑ ⟩ , ∣ ↓ ⟩ , ∣ ↑↓ ⟩ ,
the one-orbital entropy is
s p = − Tr ( ρ p log ρ p ) , ρ p = Tr p ˉ ρ . s_p
=
-\operatorname{Tr}(\rho_p\log\rho_p),
\qquad
\rho_p
=
\operatorname{Tr}_{\bar p}\rho. s p = − Tr ( ρ p log ρ p ) , ρ p = Tr p ˉ ρ .
For orbitals p , q p,q p , q , the orbital mutual information is
I ( p : q ) = s p + s q − s p q , s p q = − Tr ( ρ p q log ρ p q ) . I(p:q)
=
s_p+s_q-s_{pq},
\qquad
s_{pq}
=
-\operatorname{Tr}(\rho_{pq}\log\rho_{pq}). I ( p : q ) = s p + s q − s pq , s pq = − Tr ( ρ pq log ρ pq ) .
These are basis-dependent mode-entanglement diagnostics. See the quantum-chemistry bridge for cautions about antisymmetry, electron correlation, orbital basis choices, and spin entanglement.
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.
L. E. Ballentine, Quantum Mechanics: A Modern Development , 2nd ed., World Scientific, 2014.
A. Einstein, B. Podolsky, and N. Rosen, “Can Quantum-Mechanical Description of Physical Reality Be Considered Complete?”, Physical Review 47 , 777-780, 1935.
M. Srednicki, “Entropy and Area,” Physical Review Letters 71 , 666-669, 1993, doi:10.1103/PhysRevLett.71.666.
M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information , Cambridge University Press, 2010.
C. H. Bennett and S. J. Wiesner, “Communication via one- and two-particle operators on Einstein-Podolsky-Rosen states”, Physical Review Letters 69 , 2881-2884, 1992.
C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres, and W. K. Wootters, “Teleporting an unknown quantum state via dual classical and Einstein-Podolsky-Rosen channels”, Physical Review Letters 70 , 1895-1899, 1993.
S. L. Braunstein and P. van Loock, “Quantum information with continuous variables”, Reviews of Modern Physics 77 , 513-577, 2005.
C. Weedbrook, S. Pirandola, R. Garcia-Patron, N. J. Cerf, T. C. Ralph, J. H. Shapiro, and S. Lloyd, “Gaussian quantum information”, Reviews of Modern Physics 84 , 621-669, 2012.
A. Serafini, Quantum Continuous Variables: A Primer of Theoretical Methods , CRC Press, 2017.
R. Simon, “Peres-Horodecki Separability Criterion for Continuous Variable Systems”, Physical Review Letters 84 , 2726-2729, 2000.
L.-M. Duan, G. Giedke, J. I. Cirac, and P. Zoller, “Inseparability Criterion for Continuous Variable Systems”, Physical Review Letters 84 , 2722-2725, 2000.
D. F. Walls and G. J. Milburn, Quantum Optics , 2nd ed., Springer, 2008.
A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems , McGraw-Hill, 1971.
A. Altland and B. Simons, Condensed Matter Field Theory , 2nd ed., Cambridge University Press, 2010.
L. Amico, R. Fazio, A. Osterloh, and V. Vedral, “Entanglement in many-body systems”, Reviews of Modern Physics 80 , 517-576, 2008, doi:10.1103/RevModPhys.80.517.
J. Eisert, M. Cramer, and M. B. Plenio, “Area laws for the entanglement entropy”, Reviews of Modern Physics 82 , 277-306, 2010, doi:10.1103/RevModPhys.82.277.
U. Schollwoeck, “The density-matrix renormalization group in the age of matrix product states”, Annals of Physics 326 , 96-192, 2011, doi:10.1016/j.aop.2010.09.012.
J. Rissler, R. M. Noack, and S. R. White, “Measuring orbital interaction using quantum information theory”, Chemical Physics 323 , 519-531, 2006, doi:10.1016/j.chemphys.2005.10.018.
G. K.-L. Chan and S. Sharma, “The density matrix renormalization group in quantum chemistry”, Annual Review of Physical Chemistry 62 , 465-481, 2011, doi:10.1146/annurev-physchem-032210-103338.
L. Mandel and E. Wolf, Optical Coherence and Quantum Optics , Cambridge University Press, 1995.
D. F. Walls and G. J. Milburn, Quantum Optics , 2nd ed., Springer, 2008.
P. G. Kwiat, K. Mattle, H. Weinfurter, A. Zeilinger, A. V. Sergienko, and Y. Shih, “New high-intensity source of polarization-entangled photon pairs”, Physical Review Letters 75 , 4337-4341, 1995, doi:10.1103/PhysRevLett.75.4337.
A. Einstein, B. Podolsky, and N. Rosen, “Can Quantum-Mechanical Description of Physical Reality Be Considered Complete?”, Physical Review 47 , 777-780, 1935, doi:10.1103/PhysRev.47.777.
J. F. Clauser, M. A. Horne, A. Shimony, and R. A. Holt, “Proposed Experiment to Test Local Hidden-Variable Theories”, Physical Review Letters 23 , 880-884, 1969, doi:10.1103/PhysRevLett.23.880.