python示例中讲解_Python中qutip用法示例详解

前言

QuTip是用于模拟开放量子系统动力学的开源库。QuTip库依赖于的Numpy、Scipy和Cython的数值包。此外,matplotlib提供了图形输出。http://qutip.org/。

python安装比较容易，需要选择一个版本，python2或python3，稍微麻烦的是Scipy。

一、N原子系综自旋概率分布

from qutip import *

import numpy as np

import matplotlib.pyplot as plt

n=2#原子数

j = n//2

psi0 = spin_coherent(j, np.pi/3, 0)#设置系统的初态为自旋相干态

Jp=destroy(2*j+1).dag()#升算符

J_=destroy(2*j+1)#降算符

Jz=(Jp*J_-J_*Jp)/2#Jz

H=Jz**2#系统的哈密顿量

tlist=np.linspace(0,3,100)#时间列表

result=mesolve(H,psi0,tlist)#态随时间的演化

theta=np.linspace(0, np.pi, 50)

phi=np.linspace(0, 2*np.pi, 50)

#分别计算四个状态下的 husimi q函数

Q1, THETA1, PHI1 = spin_q_function(result.states[0], theta, phi)

Q2, THETA2, PHI2 = spin_q_function(result.states[30], theta, phi)

Q3, THETA3, PHI3