The "+R decay" procedure (see Eq. A24 through Eq. A28) in Ehrenfest+R was designed to correct the decay of the diagonal density matrix elements and offdiagonal elements independently, with the goal that the combined Ehrenfest dynamics and +R decay dynamics will match the Lindblad decay of the density matrix elements. However, the procedure to adjust the diagonals of the density matrix by changing the magnitudes of the corresponding coefficients also changes the off-diagonal density matrix elements, causing unintended deviations in the dynamics. The off-diagonal relaxation decay rate in Ehrenfest+R, which is proposed to be (see Eq. 53 in Ref. 10)