where H s (R, P ) is expressed in Eq. 86. The analytic expressions of the structure constants f ijk are provided in Eq. A5. Note that we choose to propagate the nuclear DOFs following the state-independent Liouvillian L 0 that leads to Eq. 94, not considering the statedependent Liouvillian L n i which corresponds to Eq. 97.