Superoperators (or supermatrices) are linear operators acting on density matrices (or vectorized density matrices). They generalize unitary operators: while unitaries preserve trace and positivity only for some operations, superoperators are general linear maps.

The master equation is naturally expressed in superoperator form:

$$\frac{d\rho}{dt} = \mathcal{L}[\rho]$$

where $\mathcal{L}$ is the Liouvillian superoperator.

The Liouvillian is:

$$\mathcal{L}[\rho] = -\frac{i}{\hbar}[H, \rho] + \sum_k \left( L_k \rho L_k^\dagger - \frac{1}{2}\{L_k^\dagger L_k, \rho\}\right)$$

QuTiP constructs it via liouvillian():

from qutip import *
 
H = 0.5 * sigmaz()
c_ops = [0.1 * sigmam()]
L = liouvillian(H, c_ops)
 
# L is a Qobj with dims [[2, 2], [2, 2]]
# It acts on density matrices (vectorized)

Superoperators act on vectorized density matrices: $|\rho\rangle\rangle = \text{vec}(\rho)$ where elements are stacked column-wise. The map becomes:

$$|\dot{\rho}\rangle\rangle = L |\rho\rangle\rangle$$

This is a matrix equation, amenable to linear algebra.

Superoperators are useful for analyzing spectral properties of master equations, computing eigenmodes, and advanced open system theory.