solving master equation with time dependent decay rates

[bandita1025]

'using QuantumToolbox
using CairoMakie

H = 0.5 * sigmax()
state0 = basis(2, 0)
tlist = LinRange(0, 10, 20)

gamma_t(p,t) = sin(p*t)
D = lindblad_dissipator(QuantumToolbox.sigmam())

L0 = QuantumToolbox.liouvillian(H)
L = QobjEvo([(L0, 1.0), (D, gamma_t)])

sol = mesolve(L, state0, tlist, e_ops=[sigmaz()])
sol.states'

i am trying to solve this lindbladian with a time dependent decay rate, howvevr i am not successful in running it, several errors, can you help me fix it

[ytdHuang]

@bandita1025

I add some comments below. For more details, please refer to our documentation https://qutip.org/QuantumToolbox.jl/v0.31.0/

H = 0.5 * sigmax()
state0 = basis(2, 0)
tlist = LinRange(0, 10, 20)

gamma_t(p,t) = sin(p[1]*t) # p should be a Vector or NamedTuple

D = lindblad_dissipator(sigmam())

L0 = QuantumToolbox.liouvillian(H)

# when you pass arguments into QobjEvo
# it must be a tuple, which contains pairs (inner tuples) of Qobj and coeficient_function
# if there is a time-independent term (your L0 here), you can just pass it without putting it in inner tuple (don't need to add 1.0)
L = QobjEvo((L0, (D, gamma_t)))

# pass the value of p into the solver for your gamma_t function using keyword argument params
coef = [0.123]
sol = mesolve(L, state0, tlist, e_ops=[sigmaz()], params = coef)
sol.states