skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A Large-scale Approach to Modeling Molecular Biosignatures: The Diatomics
Abstract This work presents the first steps to modeling synthetic rovibrational spectra for all molecules of astrophysical interest using a new approach implemented in the Prometheus code. The goal is to create a new comprehensive source of first-principles molecular spectra, thus bridging the gap for missing data to help drive future high-resolution studies. Our primary application domain is for molecules identified as signatures of life in planetary atmospheres (biosignatures), but our approach is general and can be applied to other systems. In this work we evaluate the accuracy of our method by studying four diatomic molecules, H 2 , O 2 , N 2 , and CO, all of which have well-known spectra. Prometheus uses the transition-optimised shifted Hermite (TOSH) theory to account for anharmonicity for the fundamental ν = 0 → ν = 1 band, along with thermal-profile modeling for the rotational transitions. To this end, we expand TOSH theory to enable the modeling of rotational constants. We show that our simple model achieves results that are a better approximation of the real spectra than those produced through an harmonic approach. We compare our results with high-resolution HITRAN and ExoMol spectral data. We find that modeling accuracy tends to diminish for rovibrational transition away from the band origin, thus highlighting the need for the theory to be further adapted.  more » « less
Award ID(s):
1927130
PAR ID:
10358255
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
925
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
57
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract In this work, we present a new approach to produce spectroscopic constants and model first-principles synthetic spectra for all molecules of astrophysical interest. We have generalized our previous diatomic molecule simulation framework, employing transition-optimized shifted Hermite (TOSH) theory, thereby enabling the modeling of polyatomic rotational constants for molecules with three or more atoms. These capabilities are now provided by our new code Epimetheus. As a first validation of our approach, we confront our predictions and assess their accuracy against the well-studied triatomic molecule ozone 666 (16O3), in addition to eight of its potential isotopomers: ozone 668 (16O16O18O), 686 (16O18O16O), 667 (16O16O17O), 676 (16O17O16O), 688 (16O18O18O), 868 (18O16O18O), 888 (18O3), and 777 (17O3). We then assess the accuracy of these rotational constants using the Epimetheus data in our code Pandora, and generate synthetic molecular spectra. The ozone spectra presented here are purely infrared and not Raman. Epimetheus builds upon the work from our previous code Prometheus, which used the TOSH theory to account for anharmonicity for the fundamentalν = 0 → ν = 1 band, going further to now account for triatomic molecules. This is combined with thermal profile modeling for the rotational transitions. We have found that this extended method performs well, typically approximating the spectroscopic constants with errors of less than 2%. Some issues do arise depending on the symmetry group of the ozone isotopomer. From these spectroscopic constants and using our own spectral modeling code, we show that we can provide the data to produce appreciable molecular spectra, which are good approximations until high-resolution studies can be done. 
    more » « less
  2. High-resolution direct absorption infrared spectra of metastable cis-formic acid (HCOOH) trapped in a cis-well resonance behind a 15 kcal/mol barrier are reported for the first time, with the energetically unstable conformer produced in a supersonic slit plasma expansion of trans-formic acid/H 2 mixtures. We present a detailed high-resolution rovibrational analysis for cis-formic acid species in the OH stretch ( ν 1 ) fundamental, providing first precision vibrational band origin, rotational constants, and term values, which in conjunction with ab initio calculations at the couple-cluster with single, double, and perturbative triple [CCSD(T)]/ANOn (n = 0, 1, 2) level support the experimental assignments and establish critical points on the potential energy surface for internal rotor trans-to-cis isomerization. Relative intensities for a- and b-type transitions observed in the spectra permit the transition dipole moment components to be determined in the body fixed frame and prove to be in good agreement with ab initio CCSD(T) theoretical estimates but in poor agreement with simple bond-dipole predictions. The observed signal dependence on H 2 in the discharge suggests the presence of a novel H atom radical chemical mechanism for strongly endothermic “up-hill” internal rotor isomerization between trans- and cis-formic acid conformers. 
    more » « less
  3. Two-dimensional infrared (2DIR) spectroscopy has become an established method for generating vibrational spectra in condensed phase samples composed of mixtures that yield heavily congested infrared and Raman spectra. These condensed phase 2DIR spectrometers can provide very high temporal resolution (<1 ps), but the spectral resolution is generally insufficient for resolving rotational peaks in gas phase spectra. Conventional (1D) rovibrational spectra of gas phase molecules are often plagued by severe spectral congestion, even when the sample is not a mixture. Spectral congestion can obscure the patterns in rovibrational spectra that are needed to assign peaks in the spectra. A method for generating high resolution 2DIR spectra of gas phase molecules has now been developed and tested using methane as the sample. The 2D rovibrational patterns that are recorded resemble an asterisk with a center position that provides the frequencies of both of the two coupled vibrational levels. The ability to generate easily recognizable 2D rovibrational patterns, regardless of temperature, should make the technique useful for a wide range of applications that are otherwise difficult or impossible when using conventional 1D rovibrational spectroscopy. 
    more » « less
  4. Highly reactive benzyl radicals are generated by electron dissociative attachment to benzyl chloride doped into a neon–hydrogen–helium discharge and immediately cooled to T rot = 15 K in a high density, supersonic slit expansion environment. The sub-Doppler spectra are fit to an asymmetric-top rotational Hamiltonian, thereby yielding spectroscopic constants for the ground ( v = 0) and first excited ( v = 1, ν 3 , ν 4 ) vibrational levels of the ground electronic state. The rotational constants obtained for the ground state are in good agreement with previous laser induced fluorescence measurements (LIF), with vibrational band origins ( ν 3 = 3073.2350 ± 0.0006 cm −1 , ν 4 = 3067.0576 ± 0.0006 cm −1 ) in agreement with anharmonically corrected density functional theory calculations. To assist in detection of benzyl radical in the interstellar medium, we have also significantly improved the precision of the ground state rotational constants through combined analysis of the ground state IR and LIF combination differences. Of dynamical interest, there is no evidence in the sub-Doppler spectra for tunneling splittings due to internal rotation of the CH 2 methylene subunit, which implies a significant rotational barrier consistent with partial double bond character in the CC bond. This is further confirmed with high level ab initio calculations at the CCSD(T)-f12b/ccpVdZ-f12 level, which predict a zero-point energy corrected barrier to internal rotation of Δ E tun ≈ 11.45 kcal mol −1 or 4005 cm −1 . In summary, the high-resolution infrared spectra are in excellent agreement with simple physical organic chemistry pictures of a strongly resonance-stabilized benzyl radical with a nearly rigid planar structure due to electron delocalization around the aromatic ring. 
    more » « less
  5. Quantum state-resolved spectroscopy was recently achieved for C60 molecules when cooled by buffer gas collisions and probed with a midinfrared frequency comb. This rovibrational quantum state resolution for the largest molecule on record is facilitated by the remarkable symmetry and rigidity of C60, which also present new opportunities and challenges to explore energy transfer between quantum states in this many-atom system. Here we combine state-specific optical pumping, buffer gas collisions, and ultrasensitive intracavity nonlinear spectroscopy to initiate and probe the rotation-vibration energy transfer and relaxation. This approach provides the first detailed characterization of C60 collisional energy transfer for a variety of collision partners, and determines the rotational and vibrational inelastic collision cross sections. These results compare well with our theoretical modeling of the collisions, and establish a route towards quantum state control of a new class of unprecedentedly large molecules. 
    more » « less