Skip to Content
Math401Math 401, Summer 2025: Freiwald research project notesMath 401, Topic 4: The quantum version of probabilistic concepts

Math401 Topic 4: The quantum version of probabilistic concepts

In mathematics, on often speaks of non-commutative instead of quantum constructions.

Note, in this section, we will see a lot of mixed used terms used in physics and mathematics. I will use italic to denote the terminology used in physics. It is safe to ignore them if you just care about the mathematics.

Section 1: Generalities about classical versus quantum systems

In classical physics, we assume that a system’s properties have well-defined values regardless of how we choose to measure them.

Basic terminology

Set of states

The preparation of a system builds a convex set of states as our initial condition for the system.

For a collection of NN system. Let procedure N1=λP1N_1=\lambda P_1 be a preparation procedure for state P1P_1, and N2=(1λ)P2N_2=(1-\lambda) P_2 be a preparation procedure for state P2P_2. The state of the collection is N=λN1+(1λ)N2N=\lambda N_1+(1-\lambda) N_2.

Set of effects

The set of effects is the set of all possible outcomes of a measurement. Ω={ω1,ω2,,ωn}\Omega=\{\omega_1, \omega_2, \ldots, \omega_n\}. Where each ωi\omega_i is an associated effect, or some query problems regarding the system. (For example, is outcome ωi\omega_i observed?)

Registration of outcomes

A pair of state and effect determines a probability Ei(P)=p(ωiP)E_i(P)=p(\omega_i|P). By the law of large numbers, this probability shall converge to N(ωi)/NN(\omega_i)/N as NN increases.

Quantum states, observables (random variables), and effects can be represented mathematically by linear operators on a Hilbert space.

Section 2: Examples of physical experiment in language of mathematics

Sten-Gernach experiment

State preparation: Silver tams are emitted from a thermal source and collimated to form a beam.

Measurement: Silver atoms interact with the field produced by the magnet and impinges on the class plate.

Registration: The impression left on the glass pace by the condensed silver atoms.

Section 3: Finite probability spaces in the language of Hilbert space and operators

Superposition is a linear combination of two or more states.

A quantum coin can be represented mathematically by linear combination of 0|0\rangle and 1|1\rangle.α0+β1HC2\alpha|0\rangle+\beta|1\rangle\in\mathscr{H}\cong\mathbb{C}^2.

For the rest of the material, we shall take the H\mathscr{H} to be vector space over C\mathbb{C}.

Definitions in classical probability under generalized probability theory

Definition of states (classical probability)

A state in classical probability is a probability distribution on the set of all possible outcomes. We can list as (p1,p2,,pn)(p_1,p_2,\cdots,p_n).

To each event AΩA\in \Omega, we associate the operator on H\mathscr{H} of multiplication by the indicator function PAMIA:fIAfP_A\coloneqq M_{\mathbb{I}_A}:f\mapsto \mathbb{I}_A f that projects onto the subspace of H\mathscr{H} corresponding to the event AA.

PA=k=1nakkkP_A=\sum_{k=1}^n a_k|k\rangle\langle k|

where ak{0,1}a_k\in\{0,1\}, and ak=1a_k=1 if and only if kAk\in A. Note that PA=PAP_A^*=P_A and PA2=PAP_A^2=P_A.

Definition of density operator (classical probability)

Let (p1,p2,,pn)(p_1,p_2,\cdots,p_n) be a probability distribution on XX, where pk0p_k\geq 0 and k=1npk=1\sum_{k=1}^n p_k=1. The density operator ρ\rho is defined by

ρ=k=1npkkk\rho=\sum_{k=1}^n p_k|k\rangle\langle k|

The probability of event AA relative to the probability distribution (p1,p2,,pn)(p_1,p_2,\cdots,p_n) becomes the trace of the product of ρ\rho and PAP_A.

Probρ(A)Tr(ρPA)\operatorname{Prob}_\rho(A)\coloneqq\text{Tr}(\rho P_A)

Definition of random variables (classical probability)

A random variable is a function f:XRf:X\to\mathbb{R}. It can also be written in operator form:

F=k=1nf(k)P{k}F=\sum_{k=1}^n f(k)P_{\{k\}}

The expectation of ff relative to the probability distribution (p1,p2,,pn)(p_1,p_2,\cdots,p_n) is given by

Eρ(f)=k=1npkf(k)=Tr(ρF)\mathbb{E}_\rho(f)=\sum_{k=1}^n p_k f(k)=\operatorname{Tr}(\rho F)

Note, by our definition of the operator F,PA,ρF,P_A,\rho (all diagonal operators) commute among themselves, which does not hold in general, in non-commutative (quantum) theory.

Section 4: Why we need generalized probability theory to study quantum systems

Story of light polarization and violation of Bell’s inequality.

Classical picture of light polarization and Bell’s inequality

An interesting story will be presented here.

Polarization of light

The light which comes through a polarizer is polarized in a certain direction. If we fixed the first filter and rotate the second filter, we will observe the intensity of the light will change.

The light intensity decreased with α\alpha (the angle between the two filters). The light should vanished when α=π/2\alpha=\pi/2.

Filter figure

By experimental measurement, the intensity of the light passing the first filter is half the beam intensity (Assume the original beam is completely unpolarized).

Then I1=I0/2I_1=I_0/2, and

I2=I0cos2αI_2=I_0\cos^2\alpha

Claim: there exist a smallest package of monochromatic light, which is a photon.

We can model the behavior of each individual photon passing through the filter with direction α\alpha with random variable PαP_\alpha. The Pα(ω)=1P_\alpha(\omega)=1 if the photon passes through the filter, and Pα(ω)=0P_\alpha(\omega)=0 if the photon does not pass through the filter.

Then, the probability of the photon passing through the two filters with direction α\alpha and β\beta is given by

E(PαPβ)=Prob(Pα=1 and Pβ=1)=12cos2(αβ)\mathbb{E}(P_\alpha P_\beta)=\operatorname{Prob}(P_\alpha=1 \text{ and } P_\beta=1)=\frac{1}{2}\cos^2(\alpha-\beta)

However, for system of 3 polarizing filters F1,F2,F3F_1,F_2,F_3, having direction α1,α2,α3\alpha_1,\alpha_2,\alpha_3. If we put them on the optical bench in pairs, Then we will have three random variables P1,P2,P3P_1,P_2,P_3.

Bell’s 3 variable inequality

Prob(P1=1,P3=0)Prob(P1=1,P2=0)+Prob(P2=1,P3=0)\operatorname{Prob}(P_1=1,P_3=0)\leq \operatorname{Prob}(P_1=1,P_2=0)+\operatorname{Prob}(P_2=1,P_3=0)

Proof

By the law of total probability, (The event that the photon passes through the first filter but not the third filter is the union of the event that the photon did not pass through the second filter and the event that the photon passed the second filter and did not pass through the third filter) we have

Prob(P1=1,P3=0)=Prob(P1=1,P2=0,P3=0)+Prob(P1=1,P2=1,P3=0)=Prob(P1=1,P2=0)Prob(P3=0)+Prob(P2=1,P3=0)Prob(P1=1)Prob(P1=1,P2=0)+Prob(P2=1,P3=0)\begin{aligned} \operatorname{Prob}(P_1=1,P_3=0)&=\operatorname{Prob}(P_1=1,P_2=0,P_3=0)+\operatorname{Prob}(P_1=1,P_2=1,P_3=0)\\ &=\operatorname{Prob}(P_1=1,P_2=0)\operatorname{Prob}(P_3=0)+\operatorname{Prob}(P_2=1,P_3=0)\operatorname{Prob}(P_1=1)\\ &\leq\operatorname{Prob}(P_1=1,P_2=0)+\operatorname{Prob}(P_2=1,P_3=0) \end{aligned}

However, according to our experimental measurement, for any pair of polarizers Fi,FjF_i,F_j, by the complement rule, we have

Prob(Pi=1,Pj=0)=Prob(Pi=1)Prob(Pi=1,Pj=1)=1212cos2(αiαj)=12sin2(αiαj)\begin{aligned} \operatorname{Prob}(P_i=1,P_j=0)&=\operatorname{Prob}(P_i=1)-\operatorname{Prob}(P_i=1,P_j=1)\\ &=\frac{1}{2}-\frac{1}{2}\cos^2(\alpha_i-\alpha_j)\\ &=\frac{1}{2}\sin^2(\alpha_i-\alpha_j) \end{aligned}

This leads to a contradiction if we apply the inequality to the experimental data.

12sin2(α1α3)12sin2(α1α2)+12sin2(α2α3)\frac{1}{2}\sin^2(\alpha_1-\alpha_3)\leq\frac{1}{2}\sin^2(\alpha_1-\alpha_2)+\frac{1}{2}\sin^2(\alpha_2-\alpha_3)

If α1=0,α2=π6,α3=π3\alpha_1=0,\alpha_2=\frac{\pi}{6},\alpha_3=\frac{\pi}{3}, then

12sin2(π3)12sin2(π6)+12sin2(π6π3)3818+183814\begin{aligned} \frac{1}{2}\sin^2(-\frac{\pi}{3})&\leq\frac{1}{2}\sin^2(-\frac{\pi}{6})+\frac{1}{2}\sin^2(\frac{\pi}{6}-\frac{\pi}{3})\\ \frac{3}{8}&\leq\frac{1}{8}+\frac{1}{8}\\ \frac{3}{8}&\leq\frac{1}{4} \end{aligned}

This is a contradiction, so Bell’s inequality is violated.

Other revised experiments (eg. Aspect’s experiment, Calcium entangled photon experiment) are also conducted and the inequality is still violated.

The true model of light polarization

The full description of the light polarization is given belows:

State of polarization of a photon: ψ=α0+β1\psi=\alpha|0\rangle+\beta|1\rangle, where 0|0\rangle and 1|1\rangle are the two orthogonal polarization states in C2\mathbb{C}^2.

Polarization filter (generalized 0,1 valued random variable): orthogonal projection PαP_\alpha on C2\mathbb{C}^2 corresponding to the direction α\alpha. (operator satisfies Pα=Pα=Pα2P_\alpha^*=P_\alpha=P_\alpha^2.)

The matrix representation of PαP_\alpha is given by

Pα=(cos2(α)cos(α)sin(α)cos(α)sin(α)sin2(α))P_\alpha=\begin{pmatrix} \cos^2(\alpha) & \cos(\alpha)\sin(\alpha)\\ \cos(\alpha)\sin(\alpha) & \sin^2(\alpha) \end{pmatrix}

Probability of a photon passing through the filter PαP_\alpha is given by Pαψ,ψ\langle P_\alpha\psi,\psi\rangle, this is cos2(α)\cos^2(\alpha) if we set ψ=0\psi=|0\rangle.

Since the probability of a photon passing through the three filters is not commutative, it is impossible to discuss Prob(P1=1,P3=0)\operatorname{Prob}(P_1=1,P_3=0) in the classical setting.

Section 5: The non-commutative (quantum) probability theory

Let H\mathscr{H} be a Hilbert space. H\mathscr{H} consists of complex-valued functions on a finite set Ω={1,2,,n}\Omega=\{1,2,\cdots,n\}. and that the functions (e1,e2,,en)(e_1,e_2,\cdots,e_n) form an orthonormal basis of H\mathscr{H}. We use Dirac notation k|k\rangle to denote the basis vector eke_k.

In classical settings, multiplication operators is now be the full space of bounded linear operators on H\mathscr{H}. (Denoted by B(H)\mathscr{B}(\mathscr{H}))

Let A,BFA,B\in\mathscr{F} be the set of all events in the classical probability settings. XX denotes the set of all possible outcomes.

A orthogonal projection on a Hilbert space is a projection operator satisfying P=PP^*=P and P2=PP^2=P. We denote the set of all orthogonal projections on H\mathscr{H} by P\mathscr{P}.

This can be found in linear algebra. Orthogonal projection 

Let P,QPP,Q\in\mathscr{P} be the event in non-commutative (quantum) probability space. R()R(\cdot) is the range of the operator. PP^\perp is the orthogonal complement of PP.

ClassicalClassical interpretationNon-commutative (Quantum)Non-commutative (Quantum) interpretation
ABA\subset BEvent AA is a subset of event BBPQP\leq QR(P)R(Q)R(P)\subseteq R(Q) Range of event PP is a subset of range of event QQ
ABA\cap BBoth event AA and BB happenedPQP\land Qprojection to R(P)R(Q)R(P)\cap R(Q) Range of event PP and event QQ happened
ABA\cup BAny of the event AA or BB happenedPQP\lor Qprojection to R(P)R(Q)R(P)\cup R(Q) Range of event PP or event QQ happened
XAX\subset A or AcA^cEvent AA did not happenPP^\perpprojectionR(P)R(P)^\perp Range of event PP is the orthogonal complement of range of event PP

In such setting, some rules of classical probability theory are not valid in quantum probability theory.

In classical probability theory, A(BC)=(AB)(AC)A\cap(B\cup C)=(A\cap B)\cup(A\cap C).

In quantum probability theory, P(QR)(PQ)(PR)P\land(Q\lor R)\neq(P\land Q)\lor(P\land R) in general.

Definitions of non-commutative (quantum) probability theory under generalized probability theory

Definition of states (non-commutative (quantum) probability theory)

A state on (B(H),P)(\mathscr{B}(\mathscr{H}),\mathscr{P}) is a map μ:P[0,1]\mu:\mathscr{P}\to[0,1] such that:

  1. μ(O)=0\mu(O)=0, where OO is the zero projection.
  2. If P1,P2,,PnP_1,P_2,\cdots,P_n are pairwise disjoint orthogonal projections, then μ(P1P2Pn)=i=1nμ(Pi)\mu(P_1\lor P_2\lor\cdots\lor P_n)=\sum_{i=1}^n\mu(P_i).

Where projections are disjoint if PiPj=PjPi=OP_iP_j=P_jP_i=O.

Definition of density operator (non-commutative (quantum) probability theory)

A density operator ρ\rho on the finite-dimensional Hilbert space H\mathscr{H} is:

  1. self-adjoint (A=AA^*=A, that is Ax,y=x,Ay\langle Ax,y\rangle=\langle x,Ay\rangle for all x,yHx,y\in\mathscr{H})
  2. positive semi-definite (all eigenvalues are non-negative)
  3. Tr(ρ)=1\operatorname{Tr}(\rho)=1.

If (ψ1,ψ2,,ψn)(|\psi_1\rangle,|\psi_2\rangle,\cdots,|\psi_n\rangle) is an orthonormal basis of H\mathscr{H} consisting of eigenvectors of ρ\rho, for the eigenvalue p1,p2,,pnp_1,p_2,\cdots,p_n, then pj0p_j\geq 0 and j=1npj=1\sum_{j=1}^n p_j=1.

We can write ρ\rho as

ρ=j=1npjψjψj\rho=\sum_{j=1}^n p_j|\psi_j\rangle\langle\psi_j|

(under basis ψj|\psi_j\rangle, it is a diagonal matrix with pjp_j on the diagonal)

Every basis of H\mathscr{H} can be decomposed to these forms.

Theorem: Born’s rule

Let ρ\rho be a density operator on H\mathscr{H}. then

μ(P)Tr(ρP)=j=1npjψjPψj\mu(P)\coloneqq\operatorname{Tr}(\rho P)=\sum_{j=1}^n p_j\langle\psi_j|P|\psi_j\rangle

Defines a probability measure on the space P\mathscr{P}.

[Proof ignored here]

Theorem: Gleason’s theorem (very important)

Let H\mathscr{H} be a Hilbert space over C\mathbb{C} or R\mathbb{R} of dimension n3n\geq 3. Let μ\mu be a state on the space P\mathscr{P} of projections on H\mathscr{H}. Then there exists a unique density operator ρ\rho such that

μ(P)=Tr(ρP)\mu(P)=\operatorname{Tr}(\rho P)

for all PPP\in\mathscr{P}. P\mathscr{P} is the space of all orthogonal projections on H\mathscr{H}.

[Proof ignored here]

Definition of random variable or Observables (non-commutative (quantum) probability theory)

It is the physical measurement of a system that we are interested in. (kinetic energy, position, momentum, etc.)

B(R)\mathscr{B}(\mathbb{R}) is the set of all Borel sets on R\mathbb{R}.

An random variable on the Hilbert space H\mathscr{H} is a projection valued map P:B(R)PP:\mathscr{B}(\mathbb{R})\to\mathscr{P}.

With the following properties:

  1. P()=OP(\emptyset)=O (the zero projection)
  2. P(R)=IP(\mathbb{R})=I (the identity projection)
  3. For any sequence A1,A2,,AnB(R)A_1,A_2,\cdots,A_n\in \mathscr{B}(\mathbb{R}). the following holds:
    (a) P(i=1nAi)=i=1nP(Ai)P(\bigcup_{i=1}^n A_i)=\bigvee_{i=1}^n P(A_i)
    (b) P(i=1nAi)=i=1nP(Ai)P(\bigcap_{i=1}^n A_i)=\bigwedge_{i=1}^n P(A_i)
    (c) P(Ac)=IP(A)P(A^c)=I-P(A) (d) If AjA_j are mutually disjoint (that is P(Ai)P(Aj)=P(Aj)P(Ai)=OP(A_i)P(A_j)=P(A_j)P(A_i)=O for iji\neq j), then P(j=1nAj)=j=1nP(Aj)P(\bigcup_{j=1}^n A_j)=\sum_{j=1}^n P(A_j)

Definition of probability of a random variable

For a system prepared in state ρ\rho, the probability of the random variable by the projection-valued measure PP is in the Borel set AA is Tr(ρP(A))\operatorname{Tr}(\rho P(A)).

Expectation of an random variable and projective measurement

Notice that if we set λ\lambda is observed with probability pλ=Tr(ρPλ)p_\lambda=\operatorname{Tr}(\rho P_\lambda). ρλsp(T)PλρPλ\rho'\coloneqq\sum_{\lambda\in sp(T)}P_\lambda \rho P_\lambda is a density operator.

Definition of expectation of operators

Let TT be a self-adjoint operator on H\mathscr{H}. The expectation of TT relative to the density operator ρ\rho is given by

Eρ(T)=Tr(ρT)=λsp(T)λTr(ρP(λ))\mathbb{E}_\rho(T)=\operatorname{Tr}(\rho T)=\sum_{\lambda\in sp(T)}\lambda \operatorname{Tr}(\rho P(\lambda))

if we set T=λsp(T)λPλT=\sum_{\lambda\in sp(T)}\lambda P_\lambda, then Eρ(T)=λsp(T)λTr(ρP(λ))\mathbb{E}_\rho(T)=\sum_{\lambda\in sp(T)}\lambda \operatorname{Tr}(\rho P(\lambda)).

The uncertainty principle

Let A,BA,B be two self-adjoint operators on H\mathscr{H}. Then we define the following two self-adjoint operators:

i[A,B]i(ABBA)i[A,B]\coloneqq i(AB-BA) ABAB+BA2A\circ B\coloneqq \frac{AB+BA}{2}

Note that ABA\circ B satisfies Jordan’s identity.

(AB)(AA)=A(B(AA))(A\circ B)\circ (A\circ A)=A\circ (B\circ (A\circ A))

Definition of variance

Given a state ρ\rho, the variance of AA is given by

Varρ(A)Eρ(A2)Eρ(A)2=Tr(ρA2)Tr(ρA)2\operatorname{Var}_\rho(A)\coloneqq\mathbb{E}_\rho(A^2)-\mathbb{E}_\rho(A)^2=\operatorname{Tr}(\rho A^2)-\operatorname{Tr}(\rho A)^2

Definition of covariance

Given a state ρ\rho, the covariance of AA and BB is given by the Jordan product of AA and BB.

Covρ(A,B)Eρ(AB)Eρ(A)Eρ(B)=Tr(ρAB)Tr(ρA)Tr(ρB)\operatorname{Cov}_\rho(A,B)\coloneqq\mathbb{E}_\rho(A\circ B)-\mathbb{E}_\rho(A)\mathbb{E}_\rho(B)=\operatorname{Tr}(\rho A\circ B)-\operatorname{Tr}(\rho A)\operatorname{Tr}(\rho B)

Robertson-Schrödinger uncertainty relation in finite dimensional Hilbert space

Let ρ\rho be a state on H\mathscr{H}, A,BA,B be two self-adjoint operators on H\mathscr{H}. Then

Varρ(A)Varρ(B)Covρ(A,B)2+14Eρ([A,B])2\operatorname{Var}_\rho(A)\operatorname{Var}_\rho(B)\geq\operatorname{Cov}_\rho(A,B)^2+\frac{1}{4}|\mathbb{E}_\rho([A,B])|^2

If ρ\rho is a pure state (ρ=ψψ\rho=|\psi\rangle\langle\psi| for some unit vector ψH|\psi\rangle\in\mathscr{H}), and the equality holds, then AA and BB are collinear (i.e. A=cBA=c B for some constant cRc\in\mathbb{R}).

When AA and BB commute, the classical inequality is recovered (Cauchy-Schwarz inequality).

Varρ(A)Varρ(B)Covρ(A,B)2\operatorname{Var}_\rho(A)\operatorname{Var}_\rho(B)\geq\operatorname{Cov}_\rho(A,B)^2

[Proof ignored here]

The uncertainty relation for unbounded symmetric operators

Definition of symmetric operator

An operator AA is symmetric if for all ϕ,ψH\phi,\psi\in\mathscr{H}, we have

Aϕ,ψ=ϕ,Aψ\langle A\phi,\psi\rangle=\langle\phi,A\psi\rangle

An example of symmetric operator is T(ψ)=idψdxT(\psi)=i\hbar\frac{d\psi}{dx}. If we let H=D(T)\mathscr{H}=\mathscr{D}(T), \hbar is the Planck constant.

D(T)\mathscr{D}(T) be the space of all square integrable, differentiable, and it’s derivative is also square integrable functions on R\mathbb{R}.

Definition of joint domain of two operators

Let (A,D(A)),(B,D(B))(A,\mathscr{D}(A)),(B,\mathscr{D}(B)) be two symmetric operators on their corresponding domains. The domain of ABAB is defined as

D(AB){ψD(B):BψD(A)}\mathscr{D}(AB)\coloneqq\{\psi\in\mathscr{D}(B):B\psi\in\mathscr{D}(A)\}

Since (AB)ψ=A(Bψ)(AB)\psi=A(B\psi), the variance of an operator AA relative to a pure state ρ=ψψ\rho=|\psi\rangle\langle\psi| is given by

Varρ(A)=Tr(ρA2)Tr(ρA)2=ψ,A2ψψ,Aψ2\operatorname{Var}_\rho(A)=\operatorname{Tr}(\rho A^2)-\operatorname{Tr}(\rho A)^2=\langle\psi,A^2\psi\rangle-\langle\psi,A\psi\rangle^2

If AA is symmetric, then Varρ(A)=Aψ,Aψψ,Aψ2\operatorname{Var}_\rho(A)=\langle A\psi,A\psi\rangle-\langle \psi, A\psi\rangle^2.

Robertson-Schrödinger uncertainty relation for unbounded symmetric operators

Let (A,D(A)),(B,D(B))(A,\mathscr{D}(A)),(B,\mathscr{D}(B)) be two symmetric operators on their corresponding domains. Then

Varρ(A)Varρ(B)Covρ(A,B)2+14Eρ([A,B])2\operatorname{Var}_\rho(A)\operatorname{Var}_\rho(B)\geq\operatorname{Cov}_\rho(A,B)^2+\frac{1}{4}|\mathbb{E}_\rho([A,B])|^2

If ρ\rho is a pure state (ρ=ψψ\rho=|\psi\rangle\langle\psi| for some unit vector ψH|\psi\rangle\in\mathscr{H}), and the equality holds, then AψA\psi and BψB\psi are collinear (i.e. Aψ=cBψA\psi=c B\psi for some constant cRc\in\mathbb{R}).

[Proof ignored here]

Summary of analog of classical probability theory and non-commutative (quantum) probability theory

Classical probabilityNon-commutative (Quantum) probability
Sample space Ω\Omega, cardinality Ω=n\vert\Omega\vert=n, example: Ω={0,1}\Omega=\{0,1\}Complex Hilbert space H\mathscr{H}, dimension dimH=n\dim\mathscr{H}=n, example: H=C2\mathscr{H}=\mathbb{C}^2
Common algebra of C\mathbb{C} valued functionsAlgebra of bounded operators B(H)\mathscr{B}(\mathscr{H})
ffˉf\mapsto \bar{f} complex conjugationPPP\mapsto P^* adjoint
Events: indicator functions of setsProjections: space of orthogonal projections PB(H)\mathscr{P}\subseteq\mathscr{B}(\mathscr{H})
functions ff such that f2=f=ff^2=f=\overline{f}orthogonal projections PP such that P=P=P2P^*=P=P^2
R\mathbb{R}-valued functions f=ff=\overline{f}self-adjoint operators A=AA=A^*
If1({λ})\mathbb{I}_{f^{-1}(\{\lambda\})} is the indicator function of the set f1({λ})f^{-1}(\{\lambda\})P(λ)P(\lambda) is the orthogonal projection to eigenspace
f=λRange(f)λIf1({λ})f=\sum_{\lambda\in \operatorname{Range}(f)}\lambda \mathbb{I}_{f^{-1}(\{\lambda\})}A=λsp(A)λP(λ)A=\sum_{\lambda\in \operatorname{sp}(A)}\lambda P(\lambda)
Probability measure μ\mu on Ω\OmegaDensity operator ρ\rho on H\mathscr{H}
Delta measure δω\delta_\omegaPure state ρ=ψψ\rho=\vert\psi\rangle\langle\psi\vert
μ\mu is non-negative measure and i=1nμ({i})=1\sum_{i=1}^n\mu(\{i\})=1ρ\rho is positive semi-definite and Tr(ρ)=1\operatorname{Tr}(\rho)=1
Expected value of random variable ff is Eμ(f)=i=1nf(i)μ({i})\mathbb{E}_{\mu}(f)=\sum_{i=1}^n f(i)\mu(\{i\})Expected value of operator AA is Eρ(A)=Tr(ρA)\mathbb{E}_\rho(A)=\operatorname{Tr}(\rho A)
Variance of random variable ff is Varμ(f)=i=1n(f(i)Eμ(f))2μ({i})\operatorname{Var}_\mu(f)=\sum_{i=1}^n (f(i)-\mathbb{E}_\mu(f))^2\mu(\{i\})Variance of operator AA is Varρ(A)=Tr(ρA2)Tr(ρA)2\operatorname{Var}_\rho(A)=\operatorname{Tr}(\rho A^2)-\operatorname{Tr}(\rho A)^2
Covariance of random variables ff and gg is Covμ(f,g)=i=1n(f(i)Eμ(f))(g(i)Eμ(g))μ({i})\operatorname{Cov}_\mu(f,g)=\sum_{i=1}^n (f(i)-\mathbb{E}_\mu(f))(g(i)-\mathbb{E}_\mu(g))\mu(\{i\})Covariance of operators AA and BB is Covρ(A,B)=Tr(ρAB)Tr(ρA)Tr(ρB)\operatorname{Cov}_\rho(A,B)=\operatorname{Tr}(\rho A\circ B)-\operatorname{Tr}(\rho A)\operatorname{Tr}(\rho B)
Composite system is given by Cartesian product of the sample spaces Ω1×Ω2\Omega_1\times\Omega_2Composite system is given by tensor product of the Hilbert spaces H1H2\mathscr{H}_1\otimes\mathscr{H}_2
Product measure μ1×μ2\mu_1\times\mu_2 on Ω1×Ω2\Omega_1\times\Omega_2Tensor product of space ρ1ρ2\rho_1\otimes\rho_2 on H1H2\mathscr{H}_1\otimes\mathscr{H}_2
Marginal distribution πv\pi_*vPartial trace Tr2(ρ)\operatorname{Tr}_2(\rho)

States of two dimensional system and the complex projective space (Bloch sphere)

Let v=eiθuv=e^{i\theta}u, then the space of pure states (ρ=uu\rho=|u\rangle\langle u|) is the complex projective space CP1\mathbb{C}P^1.

α=xi+iyi,β=x2+iy2\alpha=x_i+iy_i,\beta=x_2+iy_2 must satisfy α2+β2=1|\alpha|^2+|\beta|^2=1, that is x12+x22+y12+y22=1x_1^2+x_2^2+y_1^2+y_2^2=1.

The set of unit vectors in C2\mathbb{C}^2 is the unit sphere in R3\mathbb{R}^3.

So the space of pure states is the unit circle in R2\mathbb{R}^2.

Mapping between the space of pure states and the complex projective space

Any two dimensional pure state can be written as eiθue^{i\theta}u, where uu is a unit vector in R2\mathbb{R}^2. There exists a bijective map P:S2P1M2(C)P:S^2\to\mathscr{P}_1\subseteq M_2(\mathbb{C}) such that P(u)=uuP(u)=|u\rangle\langle u|.

P(x)=12(I+aσ)=12(1001)+ax2(0110)+ay2(0ii0)+az2(1001)P(\vec{x})=\frac{1}{2}(I+\vec{a}\cdot\vec{\sigma})=\frac{1}{2}\begin{pmatrix} 1&0\\ 0&1 \end{pmatrix}+\frac{a_x}{2}\begin{pmatrix} 0&1\\ 1&0 \end{pmatrix}+\frac{a_y}{2}\begin{pmatrix} 0&-i\\ i&0 \end{pmatrix}+\frac{a_z}{2}\begin{pmatrix} 1&0\\ 0&-1 \end{pmatrix}
Last updated on