$P$-representation

The idea here is to convert the operator equation into a partial differential equation (PDE) for the $P$-representation of the reduced density operator $\rho$.