As I understand from the official QuTiP guidlines, it is only capable of solving "cross-damping off" Master Equation in form:
$\dot{\rho(t)} = -\dfrac{i}{\hbar}[H(t),\rho(t)] + \sum\limits_n\dfrac{1}{2}\left[ C_n\rho(t)C_n^{\dagger} - \rho(t)C_n^{\dagger}C_n - C_n^{\dagger}C_n\rho(t) \right].$
Where $C_n$ are collapse operators. However when we need to consider cross-damping, the ME takes the form like:
$\dot{\rho(t)} = -\dfrac{i}{\hbar}[H(t),\rho(t)] + \sum\limits_{i,j}\dfrac{1}{2}\left[ C_i\rho(t)C_j^{\dagger} - \rho(t)C_i^{\dagger}C_j - C_i^{\dagger}C_j\rho(t) \right].$
Is QuTiP capable of solving this? And how one can manage to code this if it is?