Low temperature scattering with the R-matrix method: argon-argon scattering

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.

A Method for the Variational Calculation of Hyperfine-Resolved Rovibronic Spectra of Diatomic Molecules

An algorithm for the calculation of hyperfine structure and spectra of diatomic molecules based on the variational nuclear motion is presented. The hyperfine coupling terms considered are Fermi-contact, nuclear spin-electron spin dipole-dipole, nuclear spin-orbit, nuclear spin-rotation, and nuclear electric quadrupole interactions. Initial hyperfine-unresolved wave functions are obtained for a given set of potential energy curves and associated couplings by a variation solution of the nuclear-motion Schrödinger equation. Fully hyperfine-resolved parity-conserved rovibronic Hamiltonian matrices for a given final angular momentum, F, are constructed and then diagonalized to give hyperfine-resolved energies and wave functions. Electric transition dipole moment curves can then be used to generate a hyperfine-resolved line list by applying rigorous selection rules. The algorithm is implemented in Duo, which is a general program for calculating spectra of diatomic molecules. This approach is tested for NO and MgH, and the results are compared to experiment and shown to be consistent with those given by the well-used effective Hamiltonian code PGOPHER.

A method for calculating temperature-dependent photodissociation cross sections and rates

The destruction of molecules by photodissociation plays a major role in many radiation-rich environments, including the evolution of the atmospheres of exoplanets, which often exist close to UV-rich stars. Most current photodissociation calculations and databases assume T = 0 K, which is inadequate for hot exoplanets and stars. A method is developed for computing photodissociation spectra of diatomic molecules as a function of temperature exploiting bound state variational nuclear motion program Duo and post-processing program ExoCross. Discrete transition intensities are spread out to represent a continuous photodissociation spectrum either by Gaussian smoothing or by averaging calculations over a range of different grid sizes. Our approach is tested on four different chemical species (HCl, HF, NaCl and BeH+), showing its ability to reproduce photodissociation cross sections and rates computed with other approaches and experiment. The temperature dependence of photodissociation cross sections and rates is studies showing strong temperature variation of the photodissociation cross sections.

Quantum and quasiclassical dynamical simulations for the Ar2H+ on a new global analytical potential energy surface

A new analytical potential energy surface (PES) has been constructed for the Ar₂H⁺ system from a dataset consisting of a large number of ab initio energies computed using the coupled-cluster singles, doubles and perturbative triples method and aug-cc-pVQZ basis set. The long-range interaction is added to the diatomic potentials using a standard long range expansion form to better describe the asymptotic regions. The vibrational states for the most stable structures of the Ar₂H⁺ system have been calculated, and few low lying states are assigned to quantum numbers. Reactive scattering studies have been performed for the Ar + Ar'H⁺ → Ar' + ArH⁺ proton exchange reaction on the newly generated PES. Reaction probability, cross sections, and rate constants are calculated for the Ar + Ar'H⁺(v = 0, j = 0) collisions within 0.01 eV-0.6 eV of relative translational energy using exact quantum dynamical simulations as well as quasiclassical trajectory (QCT) calculations. The effect of vibrational excitation of the reactants is also explored for the reaction. State averaged rate constants are calculated for the proton exchange reaction at different temperatures using the QCT method. The mechanistic pathways for the reaction are understood by analyzing the quasiclassical trajectories.

General mathematical formulation of scattering processes in atom-diatomic collisions in the RmatReact methodology

Accurately modelling cold and ultracold reactive collisions occurring over deep potential wells, such as [Formula: see text], requires the development of new theoretical and computational methodologies. One potentially useful framework is the R-matrix method adopted widely for electron-molecule collisions which has more recently been applied to non-reactive heavy-particle collisions such as Ar-Ar. The existing treatment of non-reactive elastic and inelastic scattering needs to be substantially extended to enable modelling of reactive collisions: this is the subject of this paper. Herein, we develop the general mathematical formulation for non-reactive elastic and inelastic scattering, photoassociation, photodissociation, charge exchange and reactive scattering using the R-matrix method. Of particular note is that the inner region, of central importance to calculable R-matrix methodologies, must be finite in all scattering coordinates rather than a single scattering coordinate as for non-reactive scattering. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H₃⁺, H₅⁺ and beyond'.