@@ -93,16 +93,16 @@ e = basis(3, 1) # excited state
9393g = basis(3, 2) # ground state
9494
9595# operators
96- ge = tensor(qeye(N), g * e') # |g><e|
97- ue = tensor(qeye(N), u * e') # |u><e|
98- uu = tensor(qeye(N), u * u') # |u><u|
99- gg = tensor(qeye(N), g * g') # |g><g|
96+ g_e = tensor(qeye(N), g * e') # |g><e|
97+ u_e = tensor(qeye(N), u * e') # |u><e|
98+ u_u = tensor(qeye(N), u * u') # |u><u|
99+ g_g = tensor(qeye(N), g * g') # |g><g|
100100a = tensor(destroy(N), qeye(3))
101101
102102# Hamiltonian
103103g = 5 # coupling strength
104- H0 = -g * (ge ' * a + a' * ge )
105- H1 = ue ' + ue
104+ H0 = -g * (g_e ' * a + a' * g_e )
105+ H1 = u_e ' + u_e
106106f1(p, t) = 9 * exp(-(t / 5)^2)
107107H_t = QobjEvo(
108108 (
@@ -116,8 +116,8 @@ H_t = QobjEvo(
116116γ0 = 6 # Atomic decay rate
117117c_ops = [
118118 sqrt(κ) * a,
119- sqrt(4 * γ0 / 9) * ge , # Use Rb branching ratio of 4/9 e -> g
120- sqrt(5 * γ0 / 9) * ue # 5/9 e -> u
119+ sqrt(4 * γ0 / 9) * g_e , # Use Rb branching ratio of 4/9 e -> g
120+ sqrt(5 * γ0 / 9) * u_e # 5/9 e -> u
121121]
122122
123123tlist = LinRange(0, 4, 200) # Define time vector
@@ -126,8 +126,8 @@ tlist = LinRange(0, 4, 200) # Define time vector
126126# Build observables
127127e_ops = [
128128 a' * a,
129- uu ,
130- gg
129+ u_u ,
130+ g_g
131131]
132132
133133# solve dynamics
0 commit comments