Here, E represents the relevant electronic state energy in the neutral ground state geometry, and Q j is the j-th normal mode coordinate. Thus, g j , g + j , g - j in Eq.( ) are determined by calculating the energy gradients of the neutral singlet excited state, the cationic state, and the anionic state, respectively.