使用 qutip.mesolve() 计算具有时间因变量的时间演化哈密顿量

如何解决使用 qutip.mesolve() 计算具有时间因变量的时间演化哈密顿量

我试图找到一些初始状态的时间演化与时间相关的 hailtonian，其中 z 变量（时间演化的梯度）也是时间相关的。通常，当 z 变量为常数且哈密顿量定义如下时，可以使用 qutip.mesolve() 函数：

def hamiltonian_t(t,args):
          """ evaluate the hamiltonian at time t. """
          H0 = args[0]
          Hp = args[1]
          z  = args[2]
          return (1-t*z)*H0 + t*z*Hp 

这是一个从 H0 演化到 Hp 的哈密顿量。然后我将 qt.mesolve(hamiltonian_t,psi_init,tlist,[],H_args) 与 args = (H0,Hp,z) 一起使用。

但是对于与时间相关的 z 变量的当前问题，我尝试将哈密顿量定义如下：

def hamiltonian_t(tz,args):
          """ evaluate the hamiltonian at time t. """
          H0 = args[0]
          Hp = args[1]
          t  = tz[0]
          z  = tz[1]
          return (1-t*z)*H0 + t*z*Hp 

这里 args = (H0,Hp) 表示一些 H0 和 Hp。然后我尝试使用 qt.mesolve(hamiltonian_t,tz,H_args) 来计算 tz=[(tlist[0],z[0]),...,(tlist[-1],z[-1])] 的时间演变。我在 TypeError: 'float' object is not subscriptable 行收到错误 t = tz[0]。我不明白为什么会出现此错误，因为 tz 的输入是元组列表而不是浮点数。然后 mesolve 函数应该获取每个列表元素并计算哈密顿量，至少我是这么认为的。

当我对某些 H0、Hp 和 z 使用 args = (H0,z) 并且我使用 qt.mesolve(hamiltonian_t,H_args) 时，此函数通常用于查找时间演化，但我的问题涉及 az 变量，该变量对于 tlist 中的每个 t 都不同.有人可以帮我正确定义这个函数，以便 qutip.mesolve() 可以使用 tlist 中每个时间对应的 z 并计算时间演化吗？