来自qutip的mesolve的密度矩阵元素与时间的函数

Question

我将qutip的mesolve用于系统的开放系统动力学。我需要将所有密度矩阵元素保存为时间的函数。 mesolve提供了仅返回任何给定运算符的期望值的选项。我想知道是否存在一种使密度矩阵的每个元素随时间变化而不必为每个元素编写运算符的方法。如果您有任何建议，请帮助。

Answer 1

如果未指定期望运算符，则将保存每个时间步长的整个密度矩阵。诸如out = mesolve(H, rho0, tlist)之类的东西，即未指定e_ops，将导致tlist中每个时间点的状态存储在out.states中。