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.

MARVEL analysis of high-resolution rovibrational spectra of 13 C 16 O 2 .

A set of empirical rovibrational energy levels, obtained through the MARVEL (measured active rotational-vibrational energy levels) procedure, is presented for the 13 C 16 O 2 isotopologue of carbon dioxide. This procedure begins with the collection and analysis of experimental rovibrational transitions from the literature, allowing for a comprehensive review of the literature on the high-resolution spectroscopy of 13 C 16 O 2 , which is also presented. A total of 60 sources out of more than 750 checked provided 14,101 uniquely measured and assigned rovibrational transitions in the wavenumber range of 579-13,735 cm - 1 . This is followed by a weighted least-squares refinement yielding the energy levels of the states involved in the measured transitions. Altogether 6318 empirical rovibrational energies have been determined for 13 C 16 O 2 . Finally, estimates have been given for the uncertainties of the empirical energies, based on the experimental uncertainties of the transitions. The detailed analysis of the lines and the spectroscopic network built from them, as well as the uncertainty estimates, all serve to pinpoint possible errors in the experimental data, such as typos, misassignment of quantum numbers, and misidentifications. Errors found in the literature data were corrected before including them in the final MARVEL dataset and analysis.

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.

Infrared Spectra of Small Radicals for Exoplanetary Spectroscopy: OH, NH, CN and CH: The State of Current Knowledge.

In this study, we present a current state-of-the-art review of middle-to-near IR emission spectra of four simple astrophysically relevant molecular radicals-OH, NH, CN and CH. The spectra of these radicals were measured by means of time-resolved Fourier transform infrared spectroscopy in the 700-7500 cm-1 spectral range and with 0.07-0.02 cm-1 spectral resolution. The radicals were generated in a glow discharge of gaseous mixtures in a specially designed discharge cell. The spectra of short-lived radicals published here are of great importance, especially for the detailed knowledge and study of the composition of exoplanetary atmospheres in selected new planets. Today, with the help of the James Webb telescope and upcoming studies with the help of Plato and Ariel satellites, when the investigated spectral area is extended into the infrared spectral range, it means that detailed knowledge of the infrared spectra of not only stable molecules but also the spectra of short-lived radicals or ions, is indispensable. This paper follows a simple structure. Each radical is described in a separate chapter, starting with historical and actual theoretical background, continued by our experimental results and concluded by spectral line lists with assigned notation.

A variational model for the hyperfine resolved spectrum of VO in its ground electronic state.

A variational model for the infra-red spectrum of vanadium monoxide (VO) is presented, which aims to accurately predict the hyperfine structure within the VO XΣ-4 electronic ground state. To give the correct electron spin splitting of the XΣ-4 state, electron spin dipolar interaction within the ground state and the spin-orbit coupling between XΣ-4 and two excited states, AΠ4 and 1Σ+2, are calculated ab initio alongside hyperfine interaction terms. Four hyperfine coupling terms are explicitly considered: Fermi-contact interaction, electron spin-nuclear spin dipolar interaction, nuclear spin-rotation interaction, and nuclear electric quadrupole interaction. These terms are included as part of a full variational solution of the nuclear-motion Schrödinger equation performed using program Duo, which is used to generate both hyperfine-resolved energy levels and spectra. To improve the accuracy of the model, ab initio curves are subject to small shifts. The energy levels generated by this model show good agreement with the recently derived empirical term values. This and other comparisons validate both our model and the recently developed hyperfine modules in Duo.

New physical insights: Formamide discharge decomposition and the role of fragments in the formation of large biomolecules.

In this work we present a time-resolved FTIR spectroscopic study on kinetics of atomic and molecular species, specifically CO, CN radical, N2, HCN and CO2 generated in a glow discharge of formamide-nitrogen-water mixture in a helium buffer gas. Radicals such as NH, CH and OH have been proven to be fundamental stones of subsequent chemical reactions having a crucial role in a prebiotic synthesis of large organic molecules. This work contains three main goals. Firstly, we present our time-resolved spectra of formamide decomposition products and discuss the mechanism of collisional excitations between specific species. Secondly, according to our time resolution, we demonstrate and explain the band shape of CO's first overtone and the energy transfer between excited nitrogen and CO, present in our spectra. Lastly, we present theoretical results for the non-LTE modelling of the spectra using bi-temperature approach and a 1D harmonic Franck-Condon approach for the multi-molecule spectra of the formamide decomposition process in the 1800-5600 cm-1 spectral range.

VIB5 database with accurate ab initio quantum chemical molecular potential energy surfaces.

High-level ab initio quantum chemical (QC) molecular potential energy surfaces (PESs) are crucial for accurately simulating molecular rotation-vibration spectra. Machine learning (ML) can help alleviate the cost of constructing such PESs, but requires access to the original ab initio PES data, namely potential energies computed on high-density grids of nuclear geometries. In this work, we present a new structured PES database called VIB5, which contains high-quality ab initio data on 5 small polyatomic molecules of astrophysical significance (CH3Cl, CH4, SiH4, CH3F, and NaOH). The VIB5 database is based on previously used PESs, which, however, are either publicly unavailable or lacking key information to make them suitable for ML applications. The VIB5 database provides tens of thousands of grid points for each molecule with theoretical best estimates of potential energies along with their constituent energy correction terms and a data-extraction script. In addition, new complementary QC calculations of energies and energy gradients have been performed to provide a consistent database, which, e.g., can be used for gradient-based ML methods.

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.

Rovibronic spectroscopy of PN from first principles.

We report an ab initio study on the rovibronic spectroscopy of the closed-shell diatomic molecule phosphorous mononitride, PN. The study considers the nine lowest electronic states, X 1Σ+, A 1Π, C 1Σ-, D 1Δ, E 1Σ-, a 3Σ+, b 3Π, d 3Δ and e 3Σ- using high level electronic structure theory and accurate nuclear motion calculations. The ab initio data cover 9 potential energy, 14 spin-orbit coupling, 7 electronic angular momentum coupling, 9 electric dipole moment and 8 transition dipole moment curves. The Duo nuclear motion program is used to solve the coupled nuclear motion Schrödinger equations for these nine electronic states and to simulate rovibronic absorption spectra of 31P14N for different temperatures, which are compared to available spectroscopic studies. Lifetimes for all states are calculated and compared to previous results from the literature. The calculated lifetime of the A1Π state shows good agreement with an experimental value from the literature, which is an important quality indicator for the ab initio A-X transition dipole moment.

Electric quadrupole transitions in carbon dioxide.

Recent advances in high sensitivity spectroscopy have made it possible, in combination with accurate theoretical predictions, to observe, for the first time, very weak electric quadrupole transitions in a polar polyatomic molecule of water. Here, we present accurate theoretical predictions of the complete quadrupole rovibrational spectrum of a non-polar molecule CO2, important in atmospheric and astrophysical applications. Our predictions are validated by recent cavity enhanced absorption spectroscopy measurements and are used to assign few weak features in the recent ExoMars Atmospheric Chemistry Suite mid-infrared spectroscopic observations of the Martian atmosphere. Predicted quadrupole transitions appear in some of the mid-infrared CO2 and water vapor transparency regions, making them important for detection and characterization of the minor absorbers in water- and CO2-rich environments, such as those present in the atmospheres of Earth, Venus, and Mars.

Theoretical rovibronic spectroscopy of the calcium monohydroxide radical (CaOH).

The rovibronic (rotation-vibration-electronic) spectrum of the calcium monohydroxide radical (CaOH) is of interest to studies of exoplanet atmospheres and ultracold molecules. Here, we theoretically investigate the Ã2Π-X̃2Σ+ band system of CaOH using high-level ab initio theory and variational nuclear motion calculations. New potential energy surfaces (PESs) are constructed for the X̃2Σ+ and Ã2Π electronic states along with Ã-X̃ transition dipole moment surfaces (DMSs). For the ground X̃2Σ+ state, a published high-level ab initio PES is empirically refined to all available experimental rovibrational energy levels up to J = 15.5, reproducing the observed term values with a root-mean-square error of 0.06 cm-1. Large-scale multireference configuration interaction calculations using quintuple-zeta quality basis sets are employed to generate the Ã2Π state PESs and Ã-X̃ DMSs. Variational calculations consider both Renner-Teller and spin-orbit coupling effects, which are essential for a correct description of the spectrum of CaOH. Computed rovibronic energy levels of the Ã2Π state, line list calculations up to J = 125.5, and an analysis of Renner-Teller splittings in the ν2 bending mode of CaOH are discussed.

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.

Modelling the non-local thermodynamic equilibrium spectra of silylene (SiH2).

This paper sets out a robust methodology for modelling spectra of polyatomic molecules produced in reactive or dissociative environments, with vibrational populations outside local thermal equilibrium (LTE). The methodology is based on accurate, extensive ro-vibrational line lists containing transitions with high vibrational excitations and relies on the detailed ro-vibrational assignments. The developed methodology is applied to model non-LTE IR and visible spectra of silylene (SiH2) produced in a decomposition of disilane (Si2H6), a reaction of technological importance. Two approaches for non-LTE vibrational populations of the product SiH2 are introduced: a simplistic 1D approach based on the Harmonic approximation and a full 3D model incorporating accurate vibrational wavefunctions of SiH2 computed variationally with the TROVE (Theoretical ROVibrational Energy) program. We show how their non-LTE spectral signatures can be used to trace different reaction channels of molecular dissociations.

Five carbon- and nitrogen-bearing species in a hot giant planet's atmosphere.

The atmospheres of gaseous giant exoplanets orbiting close to their parent stars (hot Jupiters) have been probed for nearly two decades1,2. They allow us to investigate the chemical and physical properties of planetary atmospheres under extreme irradiation conditions3. Previous observations of hot Jupiters as they transit in front of their host stars have revealed the frequent presence of water vapour4 and carbon monoxide5 in their atmospheres; this has been studied in terms of scaled solar composition6 under the usual assumption of chemical equilibrium. Both molecules as well as hydrogen cyanide were found in the atmosphere of HD 209458b5,7,8, a well studied hot Jupiter (with equilibrium temperature around 1,500 kelvin), whereas ammonia was tentatively detected there9 and subsequently refuted10. Here we report observations of HD 209458b that indicate the presence of water (H2O), carbon monoxide (CO), hydrogen cyanide (HCN), methane (CH4), ammonia (NH3) and acetylene (C2H2), with statistical significance of 5.3 to 9.9 standard deviations per molecule. Atmospheric models in radiative and chemical equilibrium that account for the detected species indicate a carbon-rich chemistry with a carbon-to-oxygen ratio close to or greater than 1, higher than the solar value (0.55). According to existing models relating the atmospheric chemistry to planet formation and migration scenarios3,11,12, this would suggest that HD 209458b formed far from its present location and subsequently migrated inwards11,13. Other hot Jupiters may also show a richer chemistry than has been previously found, which would bring into question the frequently made assumption that they have solar-like and oxygen-rich compositions.

A spectroscopic model for the low-lying electronic states of NO.

The rovibronic structure of A2Σ+, B2Π, and C2Π states of nitric oxide (NO) is studied with the aim of producing comprehensive line lists for its near ultraviolet spectrum. Empirical energy levels for the three electronic states are determined using a combination of the empirical measured active rotation-vibration energy level (MARVEL) procedure and ab initio calculations, and the available experimental data are critically evaluated. Ab initio methods that deal simultaneously with the Rydberg-like A2Σ+ and C2Π and the valence B2Π state are tested. Methods of modeling the sharp avoided crossing between the B2Π and C2Π states are tested. A rovibronic Hamiltonian matrix is constructed using the variational nuclear motion program Duo whose eigenvalues are fitted to the MARVEL. The matrix also includes coupling terms obtained from the refinement of the ab initio potential energy and spin-orbit coupling curves. Calculated and observed energy levels agree well with each other, validating the applicability of our method and providing a useful model for this open shell system.

Treating linear molecules in calculations of rotation-vibration spectra.

In this article, a numerical implementation of the exact kinetic energy operator (KEO) for triatomic molecules (symmetric of XY2-type and asymmetric of YXZ-type) is presented. The implementation is based on the valence coordinates with the bisecting (XY2-type molecules) and bond-vector (YXZ) embeddings and includes the treatment of the singularity at linear geometry. The KEO is represented in a sum-of-product form. The singularity caused by the undetermined angle at the linear configuration is resolved with the help of the associated Legendre and Laguerre polynomials used as parameterized bending basis functions in the finite basis set representation. The exact KEO implementation is combined with the variational solver theoretical rovibrational energies, equipped with a general automatic symmetry-adaptation procedure and efficient basis step contraction schemes, providing a powerful computational solver of triatomic molecules for accurate computations of highly excited ro-vibrational spectra. The advantages of different basis set choices are discussed. Examples of specific applications for computing hot spectra of linear molecules are given.

Detection of electric-quadrupole transitions in water vapour near 5.4 and 2.5 µm.

Nowadays, the spectroscopic databases used for the modeling of Earth and planetary atmospheres provide only electric-dipole transitions for polyatomic molecules (H2O, CO2, N2O, CH4, O3…). Very recently, electric-quadrupole transitions have been detected in the high sensitivity cavity ring down spectrum (CRDS) of water vapour near 1.3 µm [A. Campargue et al., Phys. Rev. Res., 2020, 2, 023091, DOI: 10.1103/PhysRevResearch.2.023091]. This discovery paved the way to systematic searches of quadrupole transitions in water vapor and other polyatomic molecules. In the present work, on the basis of high accuracy ab initio predictions, H216O quadrupole lines are detected for the first time in the 5.4 µm and 2.5 µm regions where they are predicted to have their largest intensities (up to 10-26 cm per molecule). A total of twelve quadrupole lines are identified in two high sensitivity Fourier transform spectra recorded with a 1064 m path length. Ten lines in the 4030-4150 cm-1 region are assigned to the ν3 band while the lines near 1820 and 1926 cm-1 belong to the ν2 band. The derived line intensities which are largely above the dipole intensity cut-off of the standard spectroscopic databases, agree nicely with the theoretical predictions. We thus conclude that the calculated line list of quadrupole transitions, validated by the present measurements, should be incorporated in the spectroscopic databases.

Spectroscopy of YO from first principles.

We report an ab initio study on the spectroscopy of the open-shell diatomic molecule yttrium oxide, YO. The study considers the six lowest doublet states, X2Σ+, A'2Δ, A2Π, B2Σ+, C2Π, D2Σ+, and a few higher-lying quartet states using high levels of electronic structure theory and accurate nuclear motion calculations. The coupled cluster singles, doubles, and perturbative triples, CCSD(T), and multireference configuration interaction (MRCI) methods are employed in conjunction with a relativistic pseudopotential on the yttrium atom and a series of correlation-consistent basis sets ranging in size from triple-ζ to quintuple-ζ quality. Core-valence correlation effects are taken into account and complete basis set limit extrapolation is performed for CCSD(T). Spin-orbit coupling is included through the use of both MRCI state-interaction with spin-orbit (SI-SO) approach and four-component relativistic equation-of-motion CCSD calculations. Using the ab initio data for bond lengths ranging from 1.0 to 2.5 Å, we compute 6 potential energy, 12 spin-orbit, 8 electronic angular momentum, 6 electric dipole moment and 12 transition dipole moment (4 parallel and 8 perpendicular) curves which provide a complete description of the spectroscopy of the system of six lowest doublet states. The Duo nuclear motion program is used to solve the coupled nuclear motion Schrödinger equation for these six electronic states. The spectra of 89Y16O simulated for different temperatures are compared with several available high resolution experimental studies; good agreement is found once minor adjustments are made to the electronic excitation energies.

Variationally Computed IR Line List for the Methyl Radical CH3.

We present the first variational calculation of a hot-temperature ab initio line list for the CH3 radical. It is based on a high-level ab initio potential energy surface and dipole moment surface of CH3 in the ground electronic state. The ro-vibrational energy levels and Einstein A coefficients were calculated using the general-molecule variational approach implemented in the computer program TROVE. Vibrational energies and vibrational intensities are found to be in very good agreement with the available experimental data. The line list comprises 9 127 123 ro-vibrational states ( J ≤ 40) and 2 058 655 166 transitions, covering the wavenumber range up to 10 000 cm-1 and should be suitable for temperatures up to T = 1500 K.

Theoretical rotation-vibration spectroscopy of cis- and trans-diphosphene (P2H2) and the deuterated species P2HD.

Growing astronomical interest in phosphorous (P) chemistry is stimulating the search for new interstellar P-bearing molecules, a task requiring detailed knowledge of the microwave and infrared molecular spectrum. In this work, we present comprehensive rotation-vibration line lists of the cis- and trans-isomers of diphosphene (P2H2). The line lists have been generated using robust, first-principles methodologies based on newly computed, high-level ab initio potential energy and dipole moment surfaces. Transitions are considered between states with energies up to 8000 cm-1 and total angular momentum J ≤ 25. These are the first-ever line lists to be reported for P2H2, and they should significantly facilitate future spectroscopic characterization of this system. The deuterated species trans-P2HD and the effect of its dynamic dipole moment on the rovibrational spectrum are also discussed.