In order to derive the closed formula of the EOMs in Eq. 102 with variables {Θ n , ϕ n }, we express the Hamiltonian in terms of these variables. Using Eq. 75 and Eq. 77, the mapping Hamiltonian in Eq. 58 (with replacing R to R upon nuclear Wigner transform) can be expressed as