In the practical implementations of this approximation, such as in the recently proposed spin-LSC, the nuclear DOFs were proposed to be updated with the Liouvillian in Eq. 94. This should be viewed as an independent approximation, in addition to the linearized approximation expressed in Eq. 93. Future investigations will be carried to develop new propagation schemes taking into account the two Liouvillian components for a trajectory based method.