Python qutip用法(舉例介紹)
Python qutip用法(舉例介紹)
一、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 = spin_q_function(result.states[60], theta, phi)
Q4, THETA4, PHI4 = spin_q_function(result.states[90], theta, phi)
#在四個子圖中分別畫出四個狀態下的husimi q函式
fig = plt.figure(dpi=150,constrained_layout=1)
ax1 = fig.add_subplot(221,projection='3d')
ax2 = fig.add_subplot(222,projection='3d')
ax3 = fig.add_subplot(223,projection='3d')
ax4 = fig.add_subplot(224,projection='3d')
plot_spin_distribution_3d(Q1, THETA1, PHI1,fig=fig,ax=ax1)
plot_spin_distribution_3d(Q2, THETA2, PHI2,fig=fig,ax=ax2)
plot_spin_distribution_3d(Q3, THETA3, PHI3,fig=fig,ax=ax3)
plot_spin_distribution_3d(Q4, THETA4, PHI4,fig=fig,ax=ax4)
for ax in [ax1,ax2,ax3,ax4]:
ax.view_init(0.5*np.pi, 0)
ax.axis('off')#不顯示座標軸
fig.show()
執行結果:
二、原子與光場相互作用
from qutip import *
import numpy as np
import matplotlib.pyplot as plt
alpha=1#相干光的引數alpha
n=2#原子數
j = n/2
psi0 = tensor(coherent(10,alpha),spin_coherent(j, 0, 0))#設定系統的初態
a=destroy(10)#光場的湮滅算符
a_plus=a.dag()#光場的產生算符
Jp=destroy(n+1).dag()#原子的升算符
J_=destroy(n+1)#原子的降算符
Jx=(Jp+J_)/2#原子的Jx算符
Jy=(Jp-J_)/(2j)#原子的Jy算符,這裡的j是虛數單位
Jz=(Jp*J_-J_*Jp)/2#原子的Jz算符
H=tensor(a,Jp)+tensor(a_plus,J_)#系統的哈密頓量
tlist=np.linspace(0,10,1000)#時間列表
result=mesolve(H,psi0,tlist)#態隨時間的演化
fig=plt.figure()
ax1 = fig.add_subplot(221)
ax2 = fig.add_subplot(222)
ax3 = fig.add_subplot(223)
ax4 = fig.add_subplot(224)
ax1.plot(tlist,expect(tensor(qeye(10),Jx),result.states))#Jx的平均值隨時間變化圖
ax2.plot(tlist,expect(tensor(qeye(10),Jy),result.states))#Jy的平均值隨時間變化圖
ax3.plot(tlist,expect(tensor(qeye(10),Jz),result.states))#Jz的平均值隨時間變化圖
ax4.plot(tlist,expect(tensor(qeye(10),Jx**2+Jy**2+Jz*2),result.states))#J平方的平均值隨時間變化圖
fig.subplots_adjust(top=None,bottom=None,left=None,right=None,wspace=0.4,hspace=0.4)#設定子圖間距
fig.show()
執行結果:
相關文章
- python BeautifulSoup用法介紹Python
- $.getJSON()用法程式碼例項介紹JSON
- MongoDB Sharding Balancer介紹和設定方法舉例MongoDB
- CONSTRAINT的用法舉例AI
- Travis CI 配置檔案 .travis.yml 的語法介紹和一些用法舉例
- cpp的json的用法舉例JSON
- 【mysql】SUBSTRING_INDEX 用法舉例MySqlIndex
- css url()用法介紹CSS
- getElementsByClassName()方法用法介紹
- css vm用法介紹CSS
- getCurrentPosition用法介紹
- js通過例項介紹關於閉包的用法JS
- mysql_config_editor用法舉例MySql
- MySQL query_cache_type的DEMAND引數介紹和使用舉例MySql
- 簡單介紹Python中異常處理用法Python
- MySQL 5.7 NOT EXISTS用法介紹MySql
- jQuery css()方法用法介紹jQueryCSS
- javascript中加號(+)用法介紹JavaScript
- jQuery(html,[ownerDocument])用法介紹jQueryHTML
- replaceChild()函式用法介紹函式
- Object.isSealed()用法介紹Object
- require.js用法介紹UIJS
- 網路測試DOS命令詳細介紹及使用方法舉例
- jquery九大選擇器的用法舉例jQuery
- linux LVM 磁碟管理 基本用法舉例LinuxLVM
- C++ 列舉型別介紹C++型別
- css em單位用法介紹CSS
- jQuery filter() 用法簡單介紹jQueryFilter
- css transition屬性用法介紹CSS
- PostgreSQL資料rotate用法介紹SQL
- javascript的this用法簡單介紹JavaScript
- js WebSocket用法簡單介紹JSWeb
- javascript arguments用法簡單介紹JavaScript
- onerror事件用法簡單介紹Error事件
- Android.mk 用法介紹Android
- spam和saint的用法介紹AI
- 舉例說明在組合語言中,"[]"的用法組合語言
- 番茄鐘的團隊用法——以會議舉例