Ab initio simulation of spin-vibronic spectra of methoxy radical

Despite the fact that experimental and theoretical work on the spectrum of methoxy has stretched from the microwave to the ultraviolet and proceeded for nearly 50 years, parts of the spectrum have remained a challenge to simulate theoretically and make reliable line-by-line assignments. The spectral complexity arises because the radical has a non-zero electron spin and significant vibronic coupling between the two electronic components of the ground state due to the presence of a conical intersection. This work describes a completely ab initio effort to understand and assign the spin-vibronic levels of the X̃2E state from 0 to above 3000 cm−1, a region that includes the fundamental transitions of the C–H symmetric and asymmetric stretches that have not previously been identified uniquely. A potential energy surface for methoxy was calculated at the equation-of-motion (EOM)-coupled cluster singles, doubles, and triples (CCSDT)/atomic natural orbital (ANO1) level of theory. Subsequently, this potential energy surface was fit to a quartic power series expansion of all nine vibrational normal coordinates (as determined at the minimum of the conical intersection) by the use of a machine-learning-based algorithm. After the addition of spin–orbit coupling, the spin-vibronic problem was solved using both the Krylov–Schur and Lanczos algorithms with the SOCJT3 software to converge eigenvalues up to 3500 cm−1 and their eigenvectors. The latter were used, in conjunction with the calculated dipole moment and its derivatives (calculated using finite differences at the EOM-CCSDT/ANO1 level), to determine spectral intensities for the spin-vibronic spectra. The calculated transition frequencies and intensities were used to simulate and assign the observed transitions of the spin-vibronic spectra of the radical. The credibility of the assignments and their significance is discussed in detail.