A “square-root” method for the density matrix and its applications to Lindblad operators

Attractive features of the present factorization method are complete positivity, the no higher than linear increase of the implementation effort with the number of states involved and the introduction of randomness only at the start of the process.