A variational calculation of vibrational levels of vinyl radical

2500 vibronic eigenstates of the NO3 radical

The nitrate radical NO3 plays an important role in atmospheric chemistry, yet many aspects of its coupled and anharmonic vibronic structure remain elusive. Here, using an accurate, coupled full-dimensional diabatic potential that includes five electronic states, we revisit the vibronic spectrum associated with the electronic state. Using recently developed tensor network state methods, we are able to compute more than 2500 vibronic states, thereby increasing the number of computed full-dimensional states by a factor of 50, compared to previous work. While we obtain good agreement with experiment for most of the assigned vibronic levels, for several others, we observe striking disagreement. Further, for the antisymmetric bending motion we find remarkably large symmetry-induced level splittings that are larger than the zero-order reference. We discuss non-negligible nonadiabatic effects and show that the Born-Oppenheimer approximation leads to significant errors in the spectrum.

Tunneling splittings using modified WKB method in Cartesian coordinates: The test case of vinyl radical

Modified WKB theory for calculating tunneling splittings in symmetric multi-well systems in full dimensionality is re-derived using Cartesian coordinates. It is explicitly shown that the theory rests on the wavefunction that is exact for harmonic potentials. The theory was applied to calculate tunneling splittings in vinyl radical and some of its deuterated isotopologues in their vibrational ground states and the low-lying vibrationally excited states and compared to exact variational results. The exact results are reproduced within a factor of 2 in most states. Remarkably, all large enhancements of tunneling splittings relative to the ground state, up to three orders in magnitude in some excited mode combinations, are well reproduced. It is also shown that in the asymmetrically deuterated vinyl radical, the theory correctly predicts the states that are localized in a single well and the delocalized tunneling states. Modified WKB theory on the minimum action path is computationally inexpensive and can also be applied without modification to much larger systems in full dimensionality; the results of this test case serve to give insight into the expected accuracy of the method.

High-dimensional neural network potentials for accurate vibrational frequencies: the formic acid dimer benchmark

In recent years, machine learning potentials (MLP) for atomistic simulations have attracted a lot of attention in chemistry and materials science. Many new approaches have been developed with the primary aim to transfer the accuracy of electronic structure calculations to large condensed systems containing thousands of atoms. In spite of these advances, the reliability of modern MLPs in reproducing the subtle details of the multi-dimensional potential-energy surface is still difficult to assess for such systems. On the other hand, moderately sized systems enabling the application of tools for thorough and systematic quality-control are nowadays rarely investigated. In this work we use benchmark-quality harmonic and anharmonic vibrational frequencies as a sensitive probe for the validation of high-dimensional neural network potentials. For the case of the formic acid dimer, a frequently studied model system for which stringent spectroscopic data became recently available, we show that high-quality frequencies can be obtained from state-of-the-art calculations in excellent agreement with coupled cluster theory and experimental data.

An ab initio study of the photodissociation of the vinyl radical

Photodissociation of the vinyl radical through pathways CH₂CH → CH₂C + H, CH₂CH → CHCH + H, and CH₂CH → CH₂ + CH is investigated by means of high-level ab initio calculations. Potential-energy curves (PECs) along the corresponding dissociating bond distance associated with the ground and several excited electronic states involved in the above fragmentation pathways, as well as the nonadiabatic couplings connecting the different states, are obtained. The findings of several experiments on vinyl photodissociation performed at different excitation wavelengths are analyzed and explained qualitatively in the light of the present PECs. A two-dimensional representation (consisting of radial and angular coordinates to represent one of the H atoms of the CH₂ group) is also used to calculate the electronic states. The surfaces obtained reflect a rich variety of conical intersections, exit barriers, and nonadiabatic couplings leading to predissociation in different regions of energy and of the two coordinates, suggesting a complex photodissociation dynamics of the CH₂CH → CHCH + H pathway, with rather different fragmentation mechanisms involved. The two-dimensional results also provide interesting information on the mechanism of in-plane hydrogen migration from the CH₂ group to the CH one through a high-lying transition state.

Using nondirect product Wigner D basis functions and the symmetry-adapted Lanczos algorithm to compute the ro-vibrational spectrum of CH4-H2O

By doing calculations on the methane-water van der Waals complex, we demonstrate that highly converged energy levels and wavefunctions can be obtained using Wigner D basis functions and the Symmetry-Adapted Lanczos (SAL) method. The Wigner D basis is a nondirect product basis and, therefore, efficient when the kinetic energy operator has accessible singularities. The SAL method makes it possible to exploit symmetry to label energy levels and reduce the cost of the calculation, without explicitly using symmetry-adapted basis functions. Line strengths are computed, and new bands are identified. In particular, we find unusually strong transitions between states associated with the isomers of the global minimum and the secondary minimum.

Spectroscopy and Scattering Studies Using Interpolated Ab Initio Potentials

The Born-Oppenheimer potential energy surface (PES) has come a long way since its introduction in the 1920s, both conceptually and in predictive power for practical applications. Nevertheless, nearly 100 years later-despite astonishing advances in computational power-the state-of-the-art first-principles prediction of observables related to spectroscopy and scattering dynamics is surprisingly limited. For example, the water dimer, (H₂O)₂, with only six nuclei and 20 electrons, still presents a formidable challenge for full-dimensional variational calculations of bound states and is considered out of reach for rigorous scattering calculations. The extremely poor scaling of the most rigorous quantum methods is fundamental; however, recent progress in development of approximate methodologies has opened the door to fairly routine high-quality predictions, unthinkable 20 years ago. In this review, in relation to the workflow of spectroscopy and/or scattering studies, we summarize progress and challenges in the component areas of electronic structure calculations, PES fitting, and quantum dynamical calculations.