标签: qutip

我一直在使用Lindblad方程对开放量子系统建模工作很长时间.汉密尔顿主义者如下:

然而,另外两个矩阵被添加到哈密顿量.其中一个的所有对角线项都等于-33.3333i,其他一切都为零.另一个是矩阵,第三个对角线项等于-0.033333i.

Lindblad方程是这样的:

其中L_i是矩阵(在列表中:[L1,L2,L3,L4,L5,L6,L7]).L_i的矩阵只是一个7x7矩阵,除了L_(ii)= 1之外全零.H是总哈密顿量, 是密度矩阵,和 是一个等于的常数 其中T是温度,k是玻尔兹曼常数,和 ,其中h是普朗克常数.(请注意,gamma是自然单位)

以下代码解决了Lindblad方程,因此计算密度矩阵.然后计算并绘制这个与时间的关系:

这被称为站点3人口. 被称为胸罩和 被称为ket.两者都是载体.在这种情况下,请参阅代码以了解其定义.

这是代码:

from qutip import Qobj, Options, mesolve
import numpy as np
import scipy
from math import *
import matplotlib.pyplot as plt

hamiltonian = np.array([
    [215, -104.1, 5.1, -4.3, 4.7, -15.1, -7.8],
    [-104.1, 220.0, 32.6, 7.1, 5.4, 8.3, 0.8],
    [5.1, 32.6, 0.0, -46.8, 1.0, -8.1, 5.1],
    [-4.3, 7.1, -46.8, 125.0, -70.7, -14.7, -61.5],
    [4.7, 5.4, 1.0, -70.7, 450.0, 89.7, -2.5],
    [-15.1, 8.3, …

我是 Python 新手，现在正在学习matplotlib使用动画功能。我在下面的qutip，因为我的学习教程。但是当我复制并粘贴qutip教程的示例代码时，它不起作用错误消息是 Axes3D object is not iterable.
所以，我想检查我创建的代码，但不知道问题是否是我的代码是另一回事。我不知道该怎么做，想知道为什么示例代码不起作用。

这是教程中的示例代码：

from qutip import *
from scipy import *
def qubit_integrate(w, theta, gamma1, gamma2, psi0, tlist):

    sx = sigmax(); sy = sigmay(); sz = sigmaz(); sm = sigmam()
    H = w * (cos(theta) * sz + sin(theta) * sx)
    c_op_list = []
    n_th = 0.5 # temperature
    rate = gamma1 * (n_th + 1)
    if rate > 0.0: c_op_list.append(sqrt(rate) * sm)
    rate = gamma1 …