python - Jaynes Cummigs 模特和 Qutip

Question

我在Qutip 手册中找到了用于可视化 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))