Lab 2 - Quantum States
Useful for working examples and problems with photon quantum states. You may notice some similarity to the Jones Calculus ;-)
import numpy as np
from qutip import *These are the polarization states:
H = Qobj([[1],[0]])
V = Qobj([[0],[1]])
P45 = Qobj([[1/np.sqrt(2)],[1/np.sqrt(2)]])
M45 = Qobj([[1/np.sqrt(2)],[-1/np.sqrt(2)]])
R = Qobj([[1/np.sqrt(2)],[-1j/np.sqrt(2)]])
L = Qobj([[1/np.sqrt(2)],[1j/np.sqrt(2)]])VDevices:
HWP - Half-wave plate axis at to the horizontal
LP - Linear polarizer, axis at
QWP - Quarter-wave plate, axis at
Note, these are functions so you need to call them with a specific value of theta.
def HWP(theta):
return Qobj([[np.cos(2*theta),np.sin(2*theta)],[np.sin(2*theta),-np.cos(2*theta)]]).tidyup()def LP(theta):
return Qobj([[np.cos(theta)**2,np.cos(theta)*np.sin(theta)],[np.sin(theta)*np.cos(theta),np.sin(theta)**2]]).tidyup()def QWP(theta):
return Qobj([[np.cos(theta)**2 + 1j*np.sin(theta)**2,
(1-1j)*np.sin(theta)*np.cos(theta)],
[(1-1j)*np.sin(theta)*np.cos(theta),
np.sin(theta)**2 + 1j*np.cos(theta)**2]]).tidyup()QWP(np.pi/4)Example 1) Check that the state is normalized¶
H.dag()*HTo show more information on an object, use the question mark after the function or object:
np.sin?Example 2) Converting from ket to bra:¶
psi = Qobj([[1+1j],[2-1j]])
psipsi.dag()psi.dag().dag()the .dag() python method computes the “daggar” or the complex transpose.