Spectroscopic investigation of [Al,N,C,O] refractory molecules

Spectroscopy and Photochemistry of [Al, N, C, O, H]: Connectivity to Aluminium-Bearing Species in the Universe

Aluminium is one of the most abundant metals in the universe and impacts the evolution of various astrophysical environments. Currently detected Al-bearing molecules represent only a small fraction of the aluminium budget, suggesting that aluminium may reside in other species. AlO and AlOH molecules are abundant in the oxygen-rich supergiant stars such as VY Canis Majoris, a stellar molecular factory with 60+ molecules including the prebiotic NC-bearing species. Additional Al-bearing molecules with N, C, O, and H may form in O-rich environments with radiation-accelerated chemistry. Here, we present spectroscopic identification of novel aluminium-bearing molecules composed of [Al, N, C, O, H] and [Al, N, C, O] from the reactions of Al atoms and HNCO in solid argon matrix, which are potential Al-bearing molecules in space. Photoinduced transformations among six [Al, N, C, O, H] isomers and three [Al, N, C, O] isomers, along with their dissociation reactions forming the known interstellar species, have been disclosed. These results provide new insight into the chemical network of astronomically detected Al-bearing species in space.

Spectroscopic Study of [Mg, H, N, C, O] Species: Implications for the Astrochemical Magnesium Chemistry

Magnesium is an abundant metal element in space, and magnesium chemistry has vital importance in the evolution of interstellar medium (ISM) and circumstellar regions, such as the asymptotic giant branch star IRC+10216 where a variety of Mg compounds bearing H, C, N, and O have been detected and proposed as the important components in the gas-phase molecular clouds and solid-state dust grains. Herein, we report the formation and infrared spectroscopic characterization of the Mg-bearing molecules HMg, [Mg, N, C], [Mg, H, N, C], [Mg, N, C, O], and [Mg, H, N, C, O] from the reactions of Mg/Mg+ and the prebiotic isocyanic acid (HNCO) in the solid neon matrix. Based on their thermal diffusion and photochemical behavior, a complex reactivity landscape involving association, decomposition, and isomerization reactions of these Mg-bearing molecules is developed, which can not only help understand the chemical processes of the magnesium (iso)cyanides in astrochemistry but also provide implications on the presence of magnesium (iso)cyanates in the ISM and the chemical model for the dust grain surface reactions. It also provides a new paradigm of the key intermediate nature of the cationic complexes in the formation of neutral interstellar species.

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.

F12+EOM Quartic Force Fields for Rovibrational Predictions of Electronically Excited States

Quartic force fields (QFFs) constructed using a sum of ground-state CCSD(T)-F12b energies with EOM-CCSD excitation energies are proposed for computation of spectroscopic properties of electronically excited states. This is dubbed the F12+EOM approach and is shown to provide similar accuracy to previous methodologies at lower computational cost. Using explicitly correlated F12 approaches instead of canonical CCSD(T), as in the corresponding (T)+EOM approach, allows for 70-fold improvement in computational time. The mean percent difference between the two methods for anharmonic vibrational frequencies is only 0.10%. A similar approach is also developed herein which accounts for core correlation and scalar relativistic effects, named F12cCR+EOM. The F12+EOM and F12cCR+EOM approaches both match to within 2.5% mean absolute error of experimental fundamental frequencies. These new methods should help in clarifying astronomical spectra by assigning features to vibronic and vibrational transitions of small astromolecules when such data are not available experimentally.

Spectroscopic characterization of [H, Cl, S, O] molecular system: Potential candidate for detection in Venus atmosphere

Accurate spectroscopic parameters have been obtained for the identification of the [H, Cl, S, O] molecular system in the Venus atmosphere using computational methods. These calculations employed both standard and explicitly correlated coupled cluster techniques. All isomers possess C1 symmetry, with HOSCl being the most stable isomer. Only HOSCl and trigonal-HSOCl isomers are thermodynamically stable relative to the first dissociation limit HCl + SO. Fundamental modes of the lowest three isomers exhibit many anharmonic resonances, resulting in complex spectra. All isomers are found to be stable in the visible region as the calculation of vertical energy transition indicates. No electronic states were found to strongly absorb in the near UV-vis region.

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.

Spectral Signatures of Hydrogen Thioperoxide (HOSH) and Hydrogen Persulfide (HSSH): Possible Molecular Sulfur Sinks in the Dense ISM

For decades, sulfur has remained underdetected in molecular form within the dense interstellar medium (ISM), and somewhere a molecular sulfur sink exists where it may be hiding. With the discovery of hydrogen peroxide (HOOH) in the ISM in 2011, a natural starting point may be found in sulfur-bearing analogs that are chemically similar to HOOH: hydrogen thioperoxide (HOSH) and hydrogen persulfide (HSSH). The present theoretical study couples the accuracy in the anharmonic fundamental vibrational frequencies from the explicitly correlated coupled cluster theory with the accurate rotational constants provided by canonical high-level coupled cluster theory to produce rovibrational spectra for use in the potential observation of HOSH and HSSH. The ν6 mode for HSSH at 886.1 cm−1 is within 0.2 cm−1 of the gas-phase experiment, and the B0 rotational constant for HSSH of 6979.5 MHz is within 9.0 MHz of the experimental benchmarks, implying that the unknown spectral features (such as the first overtones and combination bands) provided herein are similarly accurate. Notably, a previous experimentally-attributed 2ν1 mode, at 7041.8 cm−1, has been reassigned to the ν1+ν5 combination band based on the present work’s ν1+ν5 value at 7034.3 cm−1. The most intense vibrational transitions for each molecule are the torsions, with HOSH having a more intense transition of 72 km/mol compared to HSSH’s intensity of 14 km/mol. Furthermore, HOSH has a larger net dipole moment of 1.60 D compared to HSSH’s 1.15 D. While HOSH may be the more likely candidate of the two for possible astronomical observation via vibrational spectroscopy due to the notable difference in their intensities, both HSSH and HOSH have large enough net dipole moments to be detectable by rotational spectroscopy to discover the role these molecules may have as possible molecular sulfur sinks in the dense ISM.

Theoretical investigation of the structures, stabilities, and vibrational and rotational spectroscopic parameters of linear HOMgNC and HMgNCO molecules by density functional theory and coupled-cluster method

Linear HOMgNC and HMgNCO molecules: two appropriate candidates for interstellar observation and experimental preparation.

(T)+EOM Quartic Force Fields for Theoretical Vibrational Spectroscopy of Electronically Excited States

(T)+EOM quartic force fields (QFFs) are proposed for ab initio rovibrational properties of electronically excited states of small molecules. The (T)+EOM method is a simple treatment of the potential surface of the excited state using a composite energy from the CCSD(T) energy for the ground-state configuration and the EOM-CCSD excitation energy for the target state. The method is benchmarked with two open-shell species, HOO and HNF, and two closed-shell species, HNO and HCF. A (T)+EOM QFF with a complete basis set extrapolation (C) and corrections for core correlation (cC) and scalar relativity (R), dubbed (T)+EOM/CcCR, achieves a mean absolute error (MAE) as low as 1.6 cm^-1 for the à ²A' state of HOO versus an established benchmark QFF with CCSD(T)-F12b/cc-pVTZ-F12 (F12-TZ) for this variationally accessible electronically excited state. The MAE for anharmonic frequencies for (T)+EOM/CcCR versus F12-TZ for HNF is 7.5 cm^-1. The closed-shell species are compared directly with the experiment, where a simpler (T)+EOM QFF using the aug-cc-pVTZ basis set compares more favorably than the more costly (T)+EOM/CcCR, suggesting a possible influence of decreasing accuracy with basis set size. Scans along internal coordinates are also provided which show reasonable modeling of the potential surface by (T)+EOM compared to benchmark QFFs computed for variationally accessible electronic states. The agreement between (T)+EOM/CcCR with F12-TZ and CcCR benchmarks is also shown to be quite accurate for rotational constants and geometries, with an MAE of 0.008 MHz for the rotational constants of (T)+EOM/CcCR versus CcCR for à ²A' HOO and agreement within 0.003 Å for bond lengths.

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.

Anharmonic Frequencies and Spectroscopic Constants of OAlOH and AlOH: Strong Bonding but Unhindered Motion

The astrophysical buildup of premineral nanocrystals from atoms to the smallest network-covalent solids will require observations of various small molecules containing the most common elements in minerals including aluminum and oxygen. The present work utilizes high-level quantum chemical quartic force field (QFF) approaches to produce anharmonic vibrational frequencies and spectroscopic constants for such species. The computed B_eff for the astrochemically known AlOH molecule at 15780.5 MHz is a mere 40 MHz above the experimental value implying that the B_eff for OAlOH at 5580.9 MHz is similarly accurate. The additional 7.31 D dipole moment in OAlOH implies that this molecule is a viable target for interstellar observation. Unlike the other anharmonic vibrational frequencies reported in this work, the Al-O-H bending frequencies in both AlOH and OAlOH are poorly described in the present QFF results. However, this failing actually highlights the fact that these bends are exceptionally floppy yet with counterintuitive exceedingly strong bonding. The Al-O bond energies are 128.2 and 107.2 kcal/mol, respective of AlOH and OAlOH, while the barriers to linearity are meager 16.6 and 380.7 cm^-1 (0.1 and 1.1 kcal/mol).