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))