When computing the population dynamics in a representation that is not the adiabatic states of V = Ĥpl , there is no unique way to calculate them in the TSH approach. In this work, we follow the estimator proposed by Subotnik et al, which shows more accurate results in our previous work for polariton dynamics. Below, we briefly introduce this estimator.