and likewise for the backward mapping variables. The Lindblad decay procedure described above is combined with the regular propagation of the mapping variables by PLDM. The overall propagation of the density matrix estimator during a time step (similar to the previous work as described in Ref. 47) can thus be summarized as