Код: Выделить всё
from qutip.nonmarkov.heom import HEOMSolver
# Define the bath and coupling
gamma = 0.05 # Damping rate
temperature = 0.1 # in units where kB = 1
bath = BosonicBath(gamma, temperature)
# Setup the solver
solver = HEOMSolver(H, [bath])
# Solve the dynamics
tlist = np.linspace(0, 10, 100)
result = solver.run(psi0, tlist)
Подробнее здесь: https://stackoverflow.com/questions/784 ... -nonmarkov