Problem 1. Let
If we now define two projection operators/matrices:
Calculate the measurement probability of qubit 1 corresponding to
Solution.
xxxxxxxxxxbegin using SymPyendxxxxxxxxxxbegin halfsqrt = 1/sympy.sqrt(2) ⋅ = * ⊗(X,Y) = kron(X,Y)end;xxxxxxxxxxbegin b₀ = [1; 0] b₁ = [0; 1] ψ = halfsqrt ⋅ (b₀ ⊗ b₀ + b₁ ⊗ b₁) ψ′= ψ.adjoint()end;xxxxxxxxxxbegin P̂₁ = sympy.Matrix([ 1 0; 0 0]) P̂₂ = sympy.Matrix([ 0 0; 0 1]) I₂ = sympy.eye(2,2)end;Notice that
The measurment probability of qubit 1 having an outcome corresponding to
xxxxxxxxxx𝓞 = P̂₁ ⊗ I₂xxxxxxxxxx𝓞′ = 𝓞.adjoint();xxxxxxxxxxp₁ = ψ′ ⋅ ( 𝓞′ ⋅ 𝓞 ) ⋅ ψ;which is
.
The resulting state from the projective measurement is then:
xxxxxxxxxxϕ = 1/√p₁[1] ⋅ 𝓞 ⋅ ψ |> sympy.MatrixThis resulting state is no longer entangled (i.e. it collapsed in the parlance of measurment interpretation). The measurement probability on the projective measurment state for the outcome of qubit 2 corresponding to the operaor
x
begin ϕ′ = ϕ.adjoint() p₂ = ϕ′ ⋅ ( 𝓞′ ⋅ 𝓞 ) ⋅ ϕend;