Anharmonic Frequencies of (MO)2 and Related Hydrides for M = Mg, Al, Si, P, S, Ca, and Ti and Heuristics for Predicting Anharmonic Corrections of Inorganic Oxides

The fundamental vibrational frequencies and spectroscopic constants of the C2O2H2 isomers: molecules known in simulated interstellar ice analogues

While trans-glyoxal may not be easily observable in astronomical sources through either IR or radioastronomy due to its C2h symmetry, its cis conformer along with the cyc-H2COCO epoxide isomer should be ready targets for astrochemical detection. The present quantum chemical study shows that not only are both molecular isomers strongly polar, they also have notable IR features and low isomerisation energies of 4.1 kcal mol-1 and 10.7 kcal mol-1, respectively. These three isomers along with two other C2O2H2 isomers have had their full set of fundamental vibrational frequencies and spectroscopic constants characterised herein. These isomers have previously been shown to occur in simulated astrophysical ices making them worthy targets of astronomical search. Furthermore, the hybrid quartic force field (QFF) approach utilized herein to produce the needed spectral data has a mean absolute percent error compared to the experimentally-available, gas phase fundamental vibrational frequencies of 0.6% and rotational constants to better than 0.1%. The hybrid QFF is defined from explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples level [CCSD(T)-F12b] including core electron correlation and a canonical CCSD(T) relativity correction for the harmonic (quadratic) terms in the QFF and simple CCSD(T)-F12b/cc-pVDZ energies for the cubic and quartic terms, the so-called "F12-TcCR+DZ QFF." This method is producing spectroscopically-accurate predictions for both fundamental vibrational frequencies and principal spectroscopic constants. Hence, the values computed in this work should be notably accurate and, hence, exceptionally useful to the spectroscopy and astrochemistry communities.

Picking up Good Vibrations through Quartic Force Fields and Vibrational Perturbation Theory

Quartic force fields (QFFs) define sparse potential energy surfaces (compared to semiglobal surfaces) that are the cheapest and easiest means of computing anharmonic vibrational frequencies, especially when utilized with second-order vibrational perturbation theory (VPT2). However, flat and shallow potential surfaces are exceedingly difficult for QFFs to treat through a combination of numerical noise in the often numerically computed derivatives and in competing energy factors in the composite energies often utilized to provide high-level spectroscopic predictions. While some of these issues can be alleviated with analytic derivatives, hybrid QFFs, and intelligent choices in coordinate systems, the best practice is for predicting good molecular vibrations via QFFs is to understand what they cannot do, and this manuscript documents such cases where QFFs may fail.

Performance of EOM-CCSD(T)(a)*-Based Quartic Force Fields in Computing Fundamental, Anharmonic Vibrational Frequencies of Molecular Electronically Excited States with Application to the Ã1A″ State of :CCH2 (Vinylidene)

Highly accurate anharmonic vibrational frequencies of electronically excited states are not as easily computed as their ground electronic state counterparts, but recently developed approximate triple excited state methods may be changing that. One emerging excited state method is equation of motion coupled cluster theory at the singles and doubles level with perturbative triples computed via the (a)* formalism, EOMEE-CCSD(T)(a)*. One of the most employed means for the ready computation of vibrational anharmonic frequencies for ground electronic states is second-order vibrational perturbation theory (VPT2), a theory based on quartic force fields (QFFs),fourth-order Taylor series expansions of the potential portion of the internuclear Watson Hamiltonian. The combination of these two is herein benchmarked for its performance for use as a means of computing rovibrational spectra of electronically excited states. Specifically, the EOMEE-CCSD(T)(a)* approach employing a complete basis set extrapolation along with core electron inclusion and relativity (the so-called "CcCR" approach) defining the QFF produces anharmonic fundamental vibrational frequencies within 2.83%, on the average, of reported gas-phase experimentally assigned values for the test set including the A ~ 1 A ″ states of HCF, HCCl, HSiF, HNO, and HPO. However, some states have exceptional accuracy in the fundamentals, most notably for ν2 of A ~ 1 A ″ HCCl in which the CcCR QFF value is within 1.8 cm-1 at 927.9 cm-1 (or 0.2%) of the experiment. Additionally, this approach produces rotational constants to, on the absolute average, within 0.41% of available experimental data, showcasing notable accuracy in the computation of rovibronic spectral data. Furthermore, utilizing a hybrid approach composed of harmonic CcCR force constants along with a set of simple EOMEE-CCSD(T)(a)*/aug-cc-pVQZ QFF cubic and quartic force constants is faster than using pure CcCR and better represents those modes that suffer from numerical instability in the anharmonic portion of the QFF, implying that this so-called "CcCR + QZ" QFF approach may be the best for future applications. Finally, complete, rovibrational spectral data are provided for A ~ 1 A 2 :CCH2, a molecule of potential astrochemical interest, in order to aid in its potential future experimental rovibronic characterization.

Quantum Chemistry and Astrochemistry: A Match Made in the Heavens

Quantum chemistry can uniquely answer astrochemical questions that no other technique can provide. Computations can be parallelized, automated, and left to run continuously providing exceptional molecular throughput that cannot be done through experimentation. Additionally, the granularity of the individual computations that are required of potential energy surfaces, reaction mechanism pathways, or other quantum chemically derived observables produces a unique mosaic that make up the larger whole. These pieces can be dissected for their individual contributions or evaluated in an ad hoc fashion for each of their roles in generating the larger whole. No other scientific approach is capable of reporting such fine-grained insights. Quantum chemistry also works from a bottom-up approach in providing properties directly from the desired molecule instead of a top-down perspective as required of experiment where molecules have to be linked to observed phenomena. Furthermore, modern quantum chemistry is well within the range of "chemical accuracy" and is approaching "spectroscopic accuracy." As such, the seemingly difficult questions asked by astrochemistry that would not be asked initially for any other application require quantum chemical reference data. While the results of quantum chemical computations are needed to interpret astrochemical observation, modeling, or laboratory experimentation, such hard questions, regardless of the original need to answer them, produce unique solutions. While questions in astrochemistry often require novel developments in and implementations of quantum chemistry as outlined herein, the applications of these solutions will stretch beyond astrochemistry and may yet impact fields much closer to Earth.

Spectroscopy and Photochemistry of Aluminum-Bearing Species in the Universe

ConspectusMetal-bearing molecules impact the chemical and physical environment of many astronomical sources such as the circumstellar envelopes of large asymptotic giant branch and red supergiant stars, the interstellar medium, and planetary atmospheres (e.g., ablation of ∼20 tons per day into the Earth's upper atmosphere). In recent decades, the number of successfully detected metal-containing molecules has increased via rotational spectroscopic observations, which are driven by theoretical and experimental investigations. Following formation, the ultimate fate of each species (stabilization, dissociation, etc.) is determined by its electronic structure and electronic spectroscopic properties as it encounters the pervasive radiation fields in the vacuum of space. Studying these properties can evince the possibility of detection and predict the impact each molecule has on its surrounding environment. Aluminum, one of the most abundant elements and metals, is distributed throughout the universe as a constituent of gas-phase molecules (e.g., AlO, AlOH, AlCl, etc.) as well as condensed onto solid dust grains such as Al2O3. Free gas-phase aluminum-bearing molecules are synthesized by nonthermal equilibrium processes such as shocks and pulsations near the stellar photosphere or via the reaction of molecules on the surface of dust grains. Recent investigations in our research group utilizing quantum chemical methods, such as coupled cluster (CCSD(T) and CCSD(T)-F12) and multireference configuration interaction (MRCI) with large basis sets, have explored a wide breadth of spectroscopy and photochemistry of small (triatomic and tetratomic) aluminum-bearing molecules, including Al-H, Al-C, Al-N, Al-O, Al-Si, Al-P, and Al-S bonds, among others. The ground-state spectroscopy (rotational and vibrational) of various aluminum-bearing molecules is discussed in the context of experimental and observational detection potentials. These detection potentials depend on various factors, such as the magnitude of the permanent dipole moment (PDM) and the population of states yielding transition frequencies in detectable ranges. Many aluminum-bearing molecules possess large PDMs and may be prime candidates for astronomical and laboratory detection. Within this discussion, interesting aspects of the ground-state molecular orbital configuration of OAlNO are shown to lead to an uncommon triplet ground state. Additionally, the electronic absorption spectrum of the quasi-isoenergetic ground-state isomers of AlOSO is discussed as a sensitive method for detecting this species and differentiating between the two isomers. Finally, photochemical mechanisms key to the production of AlO and AlOH in low-density regions and the destruction of AlCO and AlOC are also discussed in order to understand the radiation-induced formation and destruction of these molecules.

Spectroscopic Constants and Anharmonic Vibrational Frequencies of C(O)OC, c-C2O2 and Their Silicon-Containing Analogues

Comets are likely to contain various carbon oxide molecules potentially including C(O)OC and c-C2O2 on their surfaces and comae, as well as their silicon-substituted analogues possibly playing a role in the formation of interstellar dust grains. In this work, high-level quantum chemical data are provided to support such potential future astrophysical detection through the generation of predicted rovibrational data. Laboratory-based chemistry would also benefit from such aforementioned computational benchmarking considering these molecules' historic computational and experimental elusiveness. Coupled-cluster singles, doubles, and perturbative triples, the F12b formalism, and the cc-pCVTZ-F12 basis set garner the rapid, yet highly trusted F12-TcCR level of theory leveraged presently. This current work points to all four molecules' strong IR activity, coupled with large intensities, thus suggesting the potential for JWST detection. Although Si(O)OSi possesses a permanent dipole moment significantly larger than those of the other molecules of present interest, the significant abundance of the potential precursor carbon monoxide suggests that the dicarbon dioxide molecules may yet be observable in the microwave region of the electromagnetic spectrum. Thus, this present work details the likely existence and detectability of these four cyclic molecules, providing updated implications compared to previous work performed both experimentally and computationally.

pbqff: Push-Button Quartic Force Fields

pbqff is an open-source program for fully automating the production of quartic force fields (QFFs) and their corresponding anharmonic spectroscopic data. Rather than being a monolithic piece of code, it consists of several key modules including a generic interface to quantum chemistry codes and, notably, queuing systems; a molecular point group symmetry library; an internal-to-Cartesian coordinate conversion module; a module for the ordinary least-squares fitting of potential energy surfaces; and an improved second-order rotational and vibrational perturbation theory package for asymmetric and symmetric tops that handles type-1 and -2 Fermi resonances, Fermi resonance polyads, and Coriolis resonances. All of these pieces are written in Rust, a modern, safe, and performant programming language that has much to offer for scientific programming. This work introduces pbqff and its surrounding ecosystem, in addition to reporting new anharmonic vibrational data for c-(C)C₃H₂ and describing how the components of pbqff can be leveraged in other projects.

Complete, Theoretical Rovibronic Spectral Characterization of the Carbon Monoxide, Water, and Formaldehyde Cations

New high-level ab initio quartic force field (QFF) methods are explored which provide spectroscopic data for the electronically excited states of the carbon monoxide, water, and formaldehyde cations, sentinel species for expanded, recent cometary spectral analysis. QFFs based on equation-of-motion ionization potential (EOM-IP) with a complete basis set extrapolation and core correlation corrections provide assignment for the fundamental vibrational frequencies of the A˜2B1 and B˜2A1 states of the formaldehyde cation; only three of these frequencies have experimental assignment available. Rotational constants corresponding to these vibrational excitations are also provided for the first time for all electronically excited states of both of these molecules. EOM-IP-CCSDT/CcC computations support tentative re-assignment of the ν1 and ν3 frequencies of the B˜2B2 state of the water cation to approximately 2409.3 cm-1 and 1785.7 cm-1, respectively, due to significant disagreement between experimental assignment and all levels of theory computed herein, as well as work by previous authors. The EOM-IP-CCSDT/CcC QFF achieves agreement to within 12 cm-1 for the fundamental vibrational frequencies of the electronic ground state of the water cation compared to experimental values and to the high-level theoretical benchmarks for variationally-accessible states. Less costly EOM-IP based approaches are also explored using approximate triples coupled cluster methods, as well as electronically excited state QFFs based on EOM-CC3 and the previous (T)+EOM approach. The novel data, including vibrationally corrected rotational constants for all states studied herein, provided by these computations should be useful in clarifying comet evolution or other remote sensing applications in addition to fundamental spectroscopy.

Theoretical Rovibrational Spectroscopy of Magnesium Tricarbide-Multireference Character Thwarts a Full Analysis of All Isomers

Magnesium tricarbide isomers are studied herein with coupled cluster theory and multireference configuration interaction to support their possible detection in astrochemical environments such as the circumstellar envelope surrounding the star IRC +10216 or in terrestrial laboratories. Magnesium-bearing species may abound in the interstellar medium (ISM), but only eight (MgNC, MgCN, HMgNC, MgC₂H, MgC₃N, MgC₄H, MgC₅N, and MgC₆H) have been directly identified thus far. Several possible isomers for the related MgC₃ system are explored in their singlet and triplet spin multiplicities. Overall, this work offers quantum chemical insight of rovibrational spectroscopic data for MgC₃ using quartic force fields (QFFs) based on the CCSD(T) and CCSD(T)-F12 levels of theory at the complete basis set (CBS) limit. Additional corrections with small basis set CCSDT(Q) and scalar relativistic effects are also included in the analysis. Salient multireference character is found in the singlet diamond electronic state, which makes a definitive assignment of the ground state challenging. Nevertheless, coupled cluster-based composite energies and multireference configuration interaction both predict that the ¹A₁ diamond isomer is 1.6-2.2 kcal mol^-1 lower in energy than the ³A₁ diamond isomer. Furthermore, highly accurate binding energies of various isomers MgC₃ are provided for comparison to photodetachment experiments. Dipole moments along with harmonic infrared intensities will guide efforts for astronomical and spectroscopic characterization.

Polycyclic Aliphatic Hydrocarbons: Is Tetrahedrane Present in UIR Spectra?

The smallest Platonic hydrocarbon, tetrahedrane, has been subject to frequent theoretical and experimental study for 50 years, but its infrared spectrum and synthetic pathway remain a mystery. The recent partial...

Anharmonic Vibrational Frequencies of Water Borane and Associated Molecules

Water borane (BH3OH2) and borinic acid (BH2OH) have been proposed as intermediates along the pathway of hydrogen generation from simple reactants: water and borane. However, the vibrational spectra for neither water borane nor borinic acid has been investigaged experimentally due to the difficulty of isolating them in the gas phase, making accurate quantum chemical predictions for such properties the most viable means of their determination. This work presents theoretical predictions of the full rotational and fundamental vibrational spectra of these two potentially application-rich molecules using quartic force fields at the CCSD(T)-F12b/cc-pCVTZ-F12 level with additional corrections included for the effects of scalar relativity. This computational scheme is further benchmarked against the available gas-phase experimental data for the related borane and HBO molecules. The differences are found to be within 3 cm−1 for the fundamental vibrational frequencies and as close as 15 MHz in the B0 and C0 principal rotational constants. Both BH2OH and BH3OH2 have multiple vibrational modes with intensities greater than 100 km mol−1, namely ν2 and ν4 in BH2OH, and ν1, ν3, ν4, ν9, and ν13 in BH3OH2. Finally, BH3OH2 has a large dipole moment of 4.24 D, which should enable it to be observable by rotational spectroscopy, as well.

Fundamental Vibrational Frequencies and Spectroscopic Constants of Substituted Cyclopropenylidene (c-C3HX, X = F, Cl, CN)

The recent detection of ethynyl-functionalized cyclopropenylidene (c-C₃HC₂H) has initiated the search for other functional forms of cyclopropenylidene (c-C₃H₂) in space. There is existing gas-phase rotational spectroscopic data for cyano-cyclopropenylidene (c-C₃HCN), but the present work provides the first anharmonic vibrational spectral data for that molecule, as well as the first full set of both rotational and vibrational spectroscopic data for fluoro- and chloro-cyclopropenylidenes (c-C₃HF and c-C₃HCl). All three molecules have fundamental vibrational frequencies with substantial infrared intensities. Namely, c-C₃HCN has a moderately intense fundamental frequency at 1244.4 cm^-1, while c-C₃HF has two large intensity modes at 1765.4 and 1125.3 cm^-1 and c-C₃HCl again has two large intensity modes at 1692.0 and 1062.5 cm^-1. All of these frequencies are well within the spectral range covered by the high-resolution EXES instrument on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA). Further, all three molecules have dipole moments of around 3.0 D in line with c-C₃H₂, enabling them to be observed by pure rotational spectroscopy, as well. Thus, the rovibrational spectral data presented herein should assist with future laboratory studies of functionalized cyclopropenylidenes and may lead to their interstellar or circumstellar detection.

Anharmonic fundamental vibrational frequencies and spectroscopic constants of the potential HSO2 radical astromolecule

The recent report that HSO₂ is likely kinetically favored over the HOSO thermodynamic product in hydrogen addition to sulfur dioxide in simulated Venusian atmospheric conditions has led to the need for reference rotational, vibrational, and rovibrational spectral data for this molecule. While matrix-isolation spectroscopy has been able to produce vibrational frequencies for some of the vibrational modes, the full infrared to microwave spectrum of 1 ²A' HSO₂ is yet to be generated. High-level quantum chemical computations show in this work that the >2.5 D dipole moment of this radical makes it a notable target for possible radioastronomical observation. Additionally, the high intensity antisymmetric S-O stretch is computed here to be 1298.3 cm^-1, a 13.9 cm^-1 blueshift up from H₂ matrix analysis. In any case, the full set of rotational and spectroscopic constants and anharmonic fundamental vibrational frequencies is provided in this work in order to help characterize HSO₂ and probe its kinetic favorability.

Anharmonic vibrational frequencies of ammonia borane (BH3NH3)

The fundamental vibrational frequency of the B-N stretch in BH₃NH₃ has eluded gas-phase experimental observation for decades. This work offers a theoretical anharmonic prediction of this mode to be 644 cm^-1, using a Cartesian quartic force field at the CCSD(T)-F12/cc-pVTZ-F12 level of theory. The other fundamental frequencies reported herein have a mean absolute error of only 5 cm^-1 from the seven available gas-phase experimental frequencies, making the anharmonic vibrational frequencies and rotational constants the most accurate computational data available for BH₃NH₃ to date. The inclusion of Fermi, Coriolis, and Darling-Dennison resonances is a major source of this accuracy, with the non-resonance-corrected frequencies having a mean absolute error of 10 cm^-1. In particular, the inclusion of the 2ν₆ = ν₅ type 1 Fermi resonance increases the B-N stretching frequency by 14 cm^-1 compared to previous work. Ammonia borane also represents one of the largest molecules ever studied by quartic force fields, making this work an important step in extending the breadth of application for these theoretical rovibrational techniques.

Rovibrational Quantum Chemical Treatment of Inorganic and Organometallic Astrochemicals

ConspectusOur two groups have both independently and collaboratively been pushing quantum-chemical techniques to produce highly accurate predictions of anharmonic vibrational frequencies and spectroscopic constants for molecules containing atoms outside of the typical upper p block. Methodologies employ composite approaches, relying on various levels of coupled cluster theory-most often at the singles, doubles, and perturbative triples level-and quartic force field constructions of the potential portion of the intramolecular Watson Hamiltonian. Such methods are known to perform well for organic species, and we have extended this to molecules containing atoms outside of this realm.One notable atom that has received much attention in this application is magnesium. Mg is the second-most-abundant element in the Earth's mantle, and while molecules containing this element are among the confirmed astrochemicals, its further atomic abundance in the galaxy implies that many more molecules (both purely inorganic and organometallic) containing element 12 exist in astrophysical regions in chemical sizes between those of atoms and dust-sized nanocrystals. Our approach discussed herein is producing quality benchmarks and predicting novel data for magnesium-bearing molecules.The story is similar for Al and Si, which are also notably abundant in both rocky bodies and the universe at large. While Na, Sc, and Cu may not be as abundant as Mg, Al, and Si, molecules containing Na and transition metals have also previously been reported to be detected beyond the Earth. Consequently, the need to produce spectral reference data for molecules containing such atoms is growing. While several experimental groups (including, notably, the groups in Arizona, Boston, and France/Spain) have clearly led the way in detection of inorganic/organometallic molecules in space, computational support and even rational design can provide novel avenues for the detection of molecules containing atoms not typically studied in most laboratories. The application of quantum chemistry to other elements beyond carbon and its cronies at the top right of the periodic table promises a better understanding of the observable universe. It will also provide novel and fundamental chemical insights pushing the "central science" into new molecular territory.