Analysis of measured high-resolution doublet rovibronic spectra and related line lists of CH and OH

Predissociation dynamics of the hydroxyl radical (OH) based on a five-state spectroscopic model

Multi-reference configuration interaction potential energy curves (PECs) and spin-orbit couplings for the X 2Π, A 2Σ+, 1 2Σ-, 1 4Σ-, and 1 4Π states of OH are computed and refined against empirical energy levels and transitions to produce a spectroscopic model. Predissociation lifetimes are determined by discretizing continuum states in the variational method nuclear motion calculation by restricting the calculation to a finite range of internuclear separations. Varying this range gives a series of avoided crossings between quasi-bound states associated with the A 2Σ+ and continuum states, from which predissociation lifetimes are extracted. 424 quasi-bound A 2Σ+ state rovibronic energy levels are analyzed, and 374 predissociation lifetimes are produced, offering good coverage of the predissociation region. Agreement with measured lifetimes is satisfactory, and a majority of computed results were within experimental uncertainty. A previously unreported A 2Σ+ state predissociation channel that goes via X 2Π is identified in the calculations. A Python package, binSLT, produced to calculate predissociation lifetimes, associated line broadening parameters, and lifetime uncertainties is made available. The PECs and other curves from this work will be used to produce a rovibronic ExoMol line list and temperature-dependent photodissociation cross sections for the hydroxyl radical.

Numerical Equivalence of Diabatic and Adiabatic Representations in Diatomic Molecules

The (time-independent) Schrödinger equation for atomistic systems is solved by using the adiabatic potential energy curves (PECs) and the associated adiabatic approximation. In cases where interactions between electronic states become important, the associated nonadiabatic effects are taken into account via derivative couplings (DDRs), also known as nonadiabatic couplings (NACs). For diatomic molecules, the corresponding PECs in the adiabatic representation are characterized by avoided crossings. The alternative to the adiabatic approach is the diabatic representation obtained via a unitary transformation of the adiabatic states by minimizing the DDRs. For diatomics, the diabatic representation has zero DDR and nondiagonal diabatic couplings ensue. The two representations are fully equivalent and so should be the rovibronic energies and wave functions, which result from the solution of the corresponding Schrödinger equations. We demonstrate (for the first time) the numerical equivalence between the adiabatic and diabatic rovibronic calculations of diatomic molecules using the ab initio curves of yttrium oxide (YO) and carbon monohydride (CH) as examples of two-state systems, where YO is characterized by a strong NAC, while CH has a strong diabatic coupling. Rovibronic energies and wave functions are computed using a new diabatic module implemented in the variational rovibronic code Duo. We show that it is important to include both the diagonal Born-Oppenheimer correction and nondiagonal DDRs. We also show that the convergence of the vibronic energy calculations can strongly depend on the representation of nuclear motion used and that no one representation is best in all cases.