Jaynes Cummigs模型和Qutip

Question

我在Qutip manual中找到了一个简单的程序，用于可视化JCM的周期性行为。

我的代码是

from qutip import *
from pylab import *
import matplotlib.pyplot as plt    
import numpy as np

N=10
a=tensor(destroy(N),qeye(2))
sz=tensor(destroy(N),sigmaz())
s=tensor(qeye(N),destroy(2))
wc=wq=2*np.pi
g=np.pi
psi0=tensor(basis(N,0),basis(2,0))
H = wc*a.dag()*a-0.5*wq*sz+0.5*g*(a*s.dag()+a.dag()*s)
tlist=np.linspace(0,3,100)
out=mesolve(H,psi0,tlist,[],[a.dag()*a])
plot(tlist,out.expect[0])
plt.show()

在图片上我应该看到鼻窦线，但是在此程序中，我只有一条窄线。怎么了？

Answer 1

答案在物理学上： 对于激发原子：

psi0=tensor(basis(N,0),basis(2,1))