Python qutip用法(举例介绍)

最新推荐文章于 2024-08-16 09:51:56 发布

置顶

nejssd

最新推荐文章于 2024-08-16 09:51:56 发布

阅读量3k

点赞数 9

分类专栏：科学计算量子力学 qutip 文章标签： python 量子力学

本文链接：https://blog.csdn.net/nejssd/article/details/108896035

版权

Python qutip用法(举例介绍)

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

from qutip import spin_coherent,jmat,spin_q_function,plot_spin_distribution_3d,mesolve
import numpy as np
import matplotlib.pyplot as plt

n=30#原子数
j = n//2
psi0 = spin_coherent(j, np.pi/2, 0)#设置初态,自旋相干态
Jz=jmat(j,"z")
H=Jz**2#系统的哈密顿量

tlist=np.linspace(0,1,100)#时间列表
result=mesolve(H,psi0,tlist)#态随时间的演化

theta=np.linspace(0, np.pi, 100)
phi=np.linspace(0, 2*np.pi, 100)
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 = spin_q_function(result.states[60], theta, phi)
Q4, THETA4