Photoelectron spectra of benzene: Can path dependent diabatic surfaces provide unique observables?

While carrying out Beyond Born-Oppenheimer theory based diabatization, the solutions of adiabatic-to-diabatic transformation equations depend on the paths of integration over two-dimensional cross-sections of multi-dimensional space of nuclear degrees of freedom. It is shown that such path-dependent solutions leading to diabatic potential energy surface matrices computed along any two different paths are related through an orthogonal matrix, and thereby, those surface matrices should provide unique observables. While exploring the numerical validity of the theoretical framework, we construct diabatic Hamiltonians for the five low-lying electronic states (X̃2E1g, B̃2E2g, and C̃2A2u) of benzene radical cation (C6H6+) along three different approaches of contour integration over two dimensional nuclear planes constituted by seven non-adiabatically active normal modes. Three different diabatic surface matrices are further employed to generate the photoelectron spectra of the benzene molecule (C6H6). It is interesting to note that the spectral peak positions and intensity patterns for all three cases are almost close to each other and also exhibit very good agreement with the experimental results.