ABSTRACT:We explore the workability of a parallelized algorithm of time-dependent discrete variable representation (TDDVR) methodology formulated by involving "classical" trajectories on each DOF of a multi-mode multi-state Hamiltonian to reproduce the population dynamics, photoabsorption spectra and nuclear dynamics of the benzene radical cation. To perform such dynamics, we have used a realistic model Hamiltonian consists of five lowest electronic states (X 2 E 1g , B 2 E 2g , C 2 A 2u , D 2 E 1u , and E 2 B 2u ) which are interconnected through several conical intersections with nine vibrational modes. The calculated nuclear dynamics and photoabsorption spectra with the advent of our parallelized TDDVR approach show excellent agreement with the results obtained by multiconfiguration time-dependent Hartree method and experimental findings, respectively. The major focus of this article is to demonstrate how the "classical" trajectories for the different modes and the "classical" energy functional for those modes on each surface can enlight the time-dependent feature of nuclear density and its' nodal structure.