First, we generated a number of initial conditions (nuclear coordinates and momenta) near the ground-state minimum of the cis isomer using the Wigner sampling at T=0 K. Based on these initial conditions, we performed the on-the-fly Born-Oppenheimer molecular dynamics (BOMD) simulations at the |g⟩ and |e⟩ surfaces with 100 and 200 trajectories, respectively. The nuclear motions are propagated using the velocity-Verlet algorithm with a time step of dt = 0.5 fs. A total of 30 nuclear steps were run for the dynamics at the |g⟩ state and 700 nuclear steps were run for dynamics at the |e⟩ state. To enhance the sampling of the distribution of the geometries in the dynamics process, we further generate another 200 initial conditions using the Wigner sampling at T = 1000 K and performed the BOMD simulations at the |e⟩ state. A total of 100 nuclear propagation steps are run for each trajectory. We performed the same initial sampling and dynamics for the trans isomer. In addition, we created 200 conformations by rotating the central NN double bond from the cis isomer to the trans isomer and performed short-time (30 nuclear steps) BOMD at the |e⟩ state initialized from these conformations.