Problem 1. Let
If we now define two projection operators/matrices:
Calculate the measurement probability of qubit 1 corresponding to
Solution.
xxxxxxxxxx
begin
using SymPy
end
xxxxxxxxxx
begin
halfsqrt = 1/sympy.sqrt(2)
⋅ = *
⊗(X,Y) = kron(X,Y)
end;
xxxxxxxxxx
begin
b₀ = [1;
0]
b₁ = [0;
1]
ψ = halfsqrt ⋅ (b₀ ⊗ b₀ + b₁ ⊗ b₁)
ψ′= ψ.adjoint()
end;
xxxxxxxxxx
begin
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();
xxxxxxxxxx
p₁ = ψ′ ⋅ ( 𝓞′ ⋅ 𝓞 ) ⋅ ψ;
which is
.
The resulting state from the projective measurement is then:
xxxxxxxxxx
ϕ = 1/√p₁[1] ⋅ 𝓞 ⋅ ψ |> sympy.Matrix
This 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;