Numerical convergence of the Sinc discrete variable representation for solving molecular vibrational states with a conical intersection in adiabatic representation

Within the Born-Oppenheimer (BO) approximation, nuclear motions of a molecule are often envisioned to occur on an adiabatic potential energy surface (PES). However, this single PES picture should be reconsidered if a conical intersection (CI) is present, although the energy is well below the CI. The presence of the CI results in two additional terms in the nuclear Hamiltonian in the adiabatic presentation, i.e., the diagonal BO correction (DBOC) and the geometric phase (GP), which are divergent at the CI. At the same time, there are cusps in the adiabatic PESs. Thus usually it is regarded that there is numerical difficulty in a quantum dynamics calculation for treating CI in the adiabatic representation. A popular numerical method in nuclear quantum dynamics calculations is the Sinc discrete variable representation (DVR) method. We examine the numerical accuracy of the Sinc DVR method for solving the Schrodinger equation of a two dimensional model of two electronic states with a CI in both the adiabatic and diabatic representation. The results suggest that the Sinc DVR method is capable of giving reliable results in the adiabatic representation with usual density of the grid points, without special treatment of the divergence of the DBOC and the GP. The numerical uncertainty is not worse than that after the introduction of an arbitrary vector potential for accounting the GP, whose accurate form usually is not easy to obtain.