“…In this method, the exponential time evolution operator is divided into several parts to reduce the numerical error. The wave packet at time t + ∆t is calculated as: ψ(R, r, θ, t + ∆t) = e −i V(R,r,θ)∆t/2ℏ e −i T(θ)∆t/2ℏ e −i T(R,r)∆t/ℏ e −i T(θ)∆t/2ℏ × e −i V(R,r,θ)∆t/2ℏ ψ(R, r, θ, t) + O(∆t) 3 , (60) where O(∆t) 3 is the error of third order, V(R, r, θ) is the potential energy of the triatomic system, T(θ) and T(R, r) are the angular and radial KE operators, respectively.…”