Search

We present the spectroscopic characterization of cyclopropenethione (c-C3H2S) in the laboratory and detect it in space using the Green Bank Telescope (GBT) Observations of TMC-1: Hunting Aromatic Molecules (GOTHAM) survey. The detection of this molecule - the missing link in understanding the C3H2S isomeric family in TMC-1 - completes the detection of all 3 low-energy isomers of C3H2S as both CH2CCS and HCCCHS have been previously detected in this source. The total column density of this molecule (N_T of 5.72+2.65/-1.61x10^10 cm^-2 at an excitation temperature of 4.7+1.3/-1.1 K) is smaller than both CH2CCS and HCCCHS and follows nicely the relative dipole principle (RDP), a kinetic rule-of-thumb for predicting isomer abundances which suggests that, all other chemistry among a family of isomers being the same, the member with the smallest dipole should be the most abundant. The RDP now holds for the astronomical abundance ratios of both the S-bearing and O-bearing counterparts observed in TMC-1; however, CH2CCO continues to elude detection in any astronomical source., Comment: 10 Pages, 7 Figures, 11 Tables

This work aims to spectroscopically characterize and provide for the first time direct experimental frequencies of the ground vibrational state and two excited states of the simplest alkynyl thiocyanate (HCCSCN) for astrophysical use. Both microwave (8-16~GHz) and millimeter wave regions (50-120~GHz) of the spectrum have been measured and analyzed in terms of Watson's semirigid rotor Hamiltonian. A total of 314 transitions were assigned to the ground state of HCCSCN and a first set of spectroscopic constants have been accurately determined. Spectral features of the molecule were then searched for in Sgr B2(N), NGC 6334I, G+0.693-0.027 and TMC-1 molecular clouds. Upper limits to the column density are provided., Comment: 9 Pages, 6 Figures, 3 Tables

Polycyclic aromatic hydrocarbons (PAHs) are among the most ubiquitous compounds in the universe, accounting for up to ~25% of all interstellar carbon. Since most unsubstituted PAHs do not possess permanent dipole moments, they are invisible to radio astronomy. Constraining their abundances relies on the detection of polar chemical proxies, such as aromatic nitriles. We report the detection of 2- and 4-cyanopyrene, isomers of the recently detected 1-cyanopyrene. We find that these isomers are present in an abundance ratio of ~2:1:2, which mirrors the number of equivalent sites available for CN addition. We conclude that there is evidence that the cyanopyrene isomers formed by direct CN addition to pyrene under kinetic control in hydrogen-rich gas at 10 K and discuss constraints on the H/CN ratio for PAHs in TMC-1., Comment: Submitted version to comply with licensing agreements

Polycyclic aromatic hydrocarbons (PAHs) are expected to be the most abundant class of organic molecules in space. Their interstellar lifecycle is not well understood, and progress is hampered by difficulties detecting individual PAH molecules. Here, we present the discovery of CN-functionalized pyrene, a 4-ring PAH, in the dense cloud TMC-1 using the 100-m Green Bank Telescope. We derive an abundance of 1-cyanopyrene of ~1.52 x $10^{12}$ cm$^{-2}$, and from this estimate that the un-substituted pyrene accounts for up to ~0.03-0.3% of the carbon budget in the dense interstellar medium which trace the birth sites of stars and planets. The presence of pyrene in this cold (~10 K) molecular cloud agrees with its recent measurement in asteroid Ryugu where isotopic clumping suggest a cold, interstellar origin. The direct link to the birth site of our solar system is strengthened when we consider the solid state pyrene content in the pre-stellar materials compared to comets, which represent the most pristine material in the solar system. We estimate that solid state pyrene can account for 1% of the carbon within comets carried by this one single organic molecule. The abundance indicates pyrene is an "island of stability" in interstellar PAH chemistry and suggests a potential cold molecular cloud origin of the carbon carried by PAHs that is supplied to forming planetary systems, including habitable worlds such as our own., Comment: Version of manuscript revised to comply with licensing requirements

We used new high spectral resolution observations of propynal (HCCCHO) towards TMC-1 and in the laboratory to update the spectral line catalog available for transitions of HCCCHO - specifically at frequencies lower than 30 GHz which were previously discrepant in a publicly available catalog. The observed astronomical frequencies provided high enough spectral resolution that, when combined with high-resolution (~2 kHz) measurements taken in the laboratory, a new, consistent fit to both the laboratory and astronomical data was achieved. Now with a nearly exact (<1 kHz) frequency match to the J=2-1 and 3-2 transitions in the astronomical data, using a Markov chain Monte Carlo (MCMC) analysis, a best fit to the total HCCCHO column density of 7.28+4.08/-1.94 x 10^12 cm^-2 was found with a surprisingly low excitation temperature of just over 3 K. This column density is around a factor of 5 times larger than reported in previous studies. Finally, this work highlights that care is needed when using publicly available spectral catalogs to characterize astronomical spectra. The availability of these catalogs is essential to the success of modern astronomical facilities and will only become more important as the next generation of facilities come online., Comment: 14 pages, 5 figures, 5 Tables, 1 Appendix

We demonstrate for the first time that Galactic cosmic rays with energies as high as 1e10 eV can trigger a cascade of low-energy (< 20 eV) secondary electrons that could be a significant contributor to the interstellar synthesis of prebiotic molecules whose delivery by comets, meteorites, and interplanetary dust particles may have kick-started life on Earth. We explore the relative importance of low-energy (< 20 eV) secondary electrons--agents of radiation chemistry--and low-energy (< 10 eV), non-ionizing photons--instigators of photochemistry. Our calculations indicate fluxes of 100 electrons/cm2/s for low-energy secondary electrons produced within interstellar ices due to incident attenuated Galactic cosmic-ray (CR) protons. Consequently, in certain star-forming regions where internal high-energy radiation sources produce ionization rates that are observed to be a thousand times greater than the typical interstellar Galactic ionization rate, the flux of low-energy secondary electrons should far exceed that of non-ionizing photons. Because reaction cross-sections can be several orders of magnitude larger for electrons than for photons, even in the absence of such enhancement our calculations indicate that secondary low-energy electrons are at least as significant as low-energy (< 10 eV) non-ionizing photons in the interstellar synthesis of prebiotic molecules. Most importantly, our results demonstrate the pressing need for explicitly incorporating low-energy electrons in current and future astrochemical simulations of cosmic ices. Such models are critically important for interpreting James Webb Space Telescope infrared measurements, which are currently being used to probe the origins of life by studying complex organic molecules found in ices near star-forming regions., Comment: 14 pages, 6 figures

Using data from the GOTHAM (GBT Observations of TMC-1: Hunting for Aromatic Molecules) survey, we report the first astronomical detection of the C10H- anion. The astronomical observations also provided the necessary data to refine the spectroscopic parameters of C10H-. From the velocity stacked data and the matched filter response, C10H- is detected at >9{\sigma} confidence level at a column density of 4.04e11 cm-2. A dedicated search for the C10H radical was also conducted towards TMC-1. In this case, the stacked molecular emission of C10H was detected at a ~3.2{\sigma} confidence interval at a column density of 2.02e11 cm-2. However, since the determined confidence level is currently <5{\sigma}, we consider the identification of C10H as tentative. The full GOTHAM dataset was also used to better characterize the physical parameters including column density, excitation temperature, linewidth, and source size for the C4H, C6H and C8H radicals and their respective anions, and the measured column densities were compared to the predictions from a gas/grain chemical formation model and from a machine learning analysis. Given the measured values, the C10H-/C10H column density ratio is ~2.0 - the highest value measured between an anion and neutral species to date. Such a high ratio is at odds with current theories for interstellar anion chemistry. For the radical species, both models can reproduce the measured abundances found from the survey; however, the machine learning analysis matches the detected anion abundances much better than the gas/grain chemical model, suggesting that the current understanding of the formation chemistry of molecular anions is still highly uncertain., Comment: 38 Pages, 24 Figures, 12 Tables, 8 Appendices

We present laboratory rotational spectroscopy of five isomers of cyanoindene (2-, 4-, 5-, 6-, and 7-cyanoindene) using a cavity Fourier-transform microwave spectrometer operating between 6-40 GHz. Based on these measurements, we report the detection of 2-cyanoindene (1H-indene-2-carbonitrile; 2-C$_9$H$_7$CN) in GOTHAM line survey observations of the dark molecular cloud TMC-1 using the Green Bank Telescope at centimeter wavelengths. Using a combination of Markov Chain Monte Carlo (MCMC), spectral stacking, and matched filtering techniques, we find evidence for the presence of this molecule at the 6.3$\sigma$ level. This provides the first direct observation of the ratio of a cyano-substituted polycyclic aromatic hydrocarbon (PAH) to its pure hydrocarbon counterpart, in this case indene, in the same source. We discuss the possible formation chemistry of this species, including why we have only detected one of the isomers in TMC-1. We then examine the overall hydrocarbon:CN-substituted ratio across this and other simpler species, as well as compare to those ratios predicted by astrochemical models. We conclude that while astrochemical models are not yet sufficiently accurate to reproduce absolute abundances of these species, they do a good job at predicting the ratios of hydrocarbon:CN-substituted species, further solidifying -CN tagged species as excellent proxies for their fully-symmetric counterparts., Comment: Accepted in ApJL

The origin of complex organic molecules (COMs) in young Class 0 protostars has been one of the major questions in astrochemistry and star formation. While COMs are thought to form on icy dust grains via gas-grain chemistry, observational constraints on their formation pathways have been limited to gas-phase detection. Sensitive mid-infrared spectroscopy with JWST enables unprecedented investigation of COM formation by measuring their ice absorption features. We present an overview of JWST/MIRI MRS spectroscopy and imaging of a young Class 0 protostar, IRAS 15398-3359, and identify several major solid-state absorption features in the 4.9-28 $\mu$m wavelength range. These can be attributed to common ice species, such as H$_2$O, CH$_3$OH, NH$_3$, and CH$_4$, and may have contributions from more complex organic species, such as C$_2$H$_5$OH and CH$_3$CHO. The MRS spectra show many weaker emission lines at 6-8 $\mu$m, which are due to warm CO gas and water vapor, possibly from a young embedded disk previously unseen. Finally, we detect emission lines from [Fe II], [Ne II], [S I], and H$_2$, tracing a bipolar jet and outflow cavities. MIRI imaging serendipitously covers the south-western (blue-shifted) outflow lobe of IRAS 15398-3359, showing four shell-like structures similar to the outflows traced by molecular emission at sub-mm wavelengths. This overview analysis highlights the vast potential of JWST/MIRI observations and previews scientific discoveries in the coming years., Comment: Accepted to ApJL. The methods of data reduction and spectral extraction are updated

Polycyclic Aromatic Hydrocarbons (PAHs) have long been invoked in the study of interstellar and protostellar sources, but the unambiguous identification of any individual PAH has proven elusive until very recently. As a result, the formation mechanisms for this important class of molecules remain poorly constrained. Here we report the first interstellar detection of a pure hydrocarbon PAH, indene (C$_9$H$_8$), as part of the GBT Observations of TMC-1: Hunting for Aromatic Molecules (GOTHAM) survey. This detection provides a new avenue for chemical inquiry, complementing the existing detections of CN-functionalized aromatic molecules. From fitting the GOTHAM observations, indene is found to be the most abundant organic ring detected in TMC-1 to date. And from astrochemical modeling with NAUTILUS, the observed abundance is greater than the model's prediction by several orders of magnitude suggesting that current formation pathways in astrochemical models are incomplete. The detection of indene in relatively high abundance implies related species such as cyanoindene, cyclopentadiene, toluene, and styrene may be detectable in dark clouds., Comment: 18 pages, 2 figures and 3 tables in the main text. 1 table and 1 figure in the Appendix. Published in The Astrophysical Journal Letters. Supplementary data available in the DataVerse entry provided in text

Ubiquitous unidentified infrared emission bands are seen in many astronomical sources. Although these bands are widely, if not unanimously, attributed to the collective emission from polycyclic aromatic hydrocarbons, no single species from this class has been detected in space. We present the discovery of two -CN functionalized polycyclic aromatic hydrocarbons, 1- and 2-cyanonaphthalene, in the interstellar medium aided by spectral matched filtering. Using radio observations with the Green Bank Telescope, we observe both bi-cyclic ring molecules in the molecular cloud TMC-1. We discuss potential in situ gas-phase formation pathways from smaller organic precursor molecules.

We report the discovery of two unsaturated organic species, trans-(E)-cyanovinylacetylene and vinylcyanoacetylene, using the second data release of the GOTHAM deep survey towards TMC-1 with the 100 m Green Bank Telescope. For both detections, we performed velocity stacking and matched filter analyses using Markov chain Monte Carlo simulations, and for trans-(E)-cyanovinylacetylene, three rotational lines were observed at low signal-to-noise (${\sim}$3$\sigma$). From this analysis, we derive column densities of $2\times10^{11}$ and $3\times10^{11}$ cm$^{-2}$ for vinylcyanoacetylene and trans-(E)-cyanovinylacetylene, respectively, and an upper limit of $<2\times10^{11}$ cm$^{-2}$ for trans-(Z)-cyanovinylacetylene. Comparisons with G3//B3LYP semi-empirical thermochemical calculations indicate abundances of the [H$_3$C$_5$N}] isomers are not consistent with their thermodynamic stability, and instead their abundances are mainly driven by dynamics. We provide discussion into how these species may be formed in TMC-1, with reference to related species like vinyl cyanide (CH$_2$=CHCN). As part of this discussion, we performed the same analysis for ethyl cyanide (CH$_3$CH$_2$CN), the hydrogenation product of CH$_2$=CHCN. This analysis provides evidence -- at 4.17$\sigma$ significance -- an upper limit to the column density of $<4\times10^{11}$ cm$^{-2}$; an order of magnitude lower than previous upper limits towards this source., Comment: Four figures, two tables. Five supplementary figures and four supplementary tables. Published in The Astrophysical Journal Letters

We propose a new model for treating solid-phase photoprocesses in interstellar ice analogues. In this approach, photoionization and photoexcitation are included in more detail, and the production of electronically-excited (suprathermal) species is explicitly considered. In addition, we have included non-thermal, non-diffusive chemistry to account for the low-temperature characteristic of cold cores. As an initial test of our method, we have simulated two previous experimental studies involving the UV irradiation of pure solid O$_2$. In contrast to previous solid-state astrochemical model calculations which have used gas-phase photoabsorption cross-sections, we have employed solid-state cross-sections in our calculations. This method allows the model to be tested using well-constrained experiments rather than poorly constrained gas-phase abundances in ISM regions. Our results indicate that inclusion of non-thermal reactions and suprathermal species allows for reproduction of low-temperature solid-phase photoprocessing that simulate interstellar ices within cold ($\sim$ 10 K) dense cores such as TMC-1., Comment: ApJ, accepted: 15 pages, 3 figures

Complex organic molecules (COMs) have been detected in a variety of interstellar sources. The abundances of these COMs in warming sources can be explained by syntheses linked to increasing temperatures and densities, allowing quasi-thermal chemical reactions to occur rapidly enough to produce observable amounts of COMs, both in the gas phase, and upon dust grain ice mantles. The COMs produced on grains then become gaseous as the temperature increases sufficiently to allow their thermal desorption. The recent observation of gaseous COMs in cold sources has not been fully explained by these gas-phase and dust grain production routes. Radiolysis chemistry is a possible non-thermal method of producing COMs in cold dark clouds. This new method greatly increases the modeled abundance of selected COMs upon the ice surface and within the ice mantle due to excitation and ionization events from cosmic ray bombardment. We examine the effect of radiolysis on three C$_{2}$H$_{4}$O$_{2}$ isomers -- methyl formate (HCOOCH$_3$), glycolaldehyde (HCOCH$_2$OH), and acetic acid (CH$_3$COOH) -- and a chemically similar molecule, dimethyl ether (CH$_3$OCH$_3$), in cold dark clouds. We then compare our modelled gaseous abundances with observed abundances in TMC-1, L1689B, and B1-b., Comment: 12 pages, 6 Figures, 5 tables, 33 numbered equations

Cosmic rays are ubiquitous in interstellar environments, and their bombardment of dust-grain ice mantles is a possible driver for the formation of complex, even prebiotic molecules. Yet, critical data that are essential for accurate modeling of this phenomenon, such as the average radii of cosmic-ray tracks in amorphous solid water (ASW) remain unconstrained. It is shown that cosmic ray tracks in ASW can be approximated as a cylindrical volume with an average radius that is mostly independent of the initial particle energy. Interactions between energetic ions and both a low-density amorphous (LDA) and high-density amorphous (HDA) ice target are simulated using the Geant4-DNA Monte Carlo toolkit, which allows for tracking secondary electrons down to subexcitation energies in the material. We find the peak track core radii, $r_\mathrm{cyl}$, for LDA and HDA ices to be 9.9 nm and 8.4 nm, respectively - somewhat less than double the value of 5 nm often assumed in astrochemical models., Comment: 19 pages, 7 figures, 1 table (accepted in ApJ)

Much like six-membered rings, five-membered rings are ubiquitous in organic chemistry, frequently serving as the building blocks for larger molecules, including many of biochemical importance. From a combination of laboratory rotational spectroscopy and a sensitive spectral line survey in the radio band toward the starless cloud core TMC-1, we report the astronomical detection of 1-cyano-1,3-cyclopentadiene, $c$-C$_5$H$_5$CN}, a highly polar, cyano derivative of cyclopentadiene, $c$-C$_5$H$_6$. The derived abundance of $c$-C$_5$H$_5$CN} is far greater than predicted from astrochemical models which well reproduce the abundance of many carbon chains. This finding implies either an important production mechanism or a large reservoir of aromatic material may need to be considered. The apparent absence of its closely-related isomer, 2-cyano-1,3-cyclopentadiene, may arise from its lower stability or be indicative of a more selective pathway for formation of the 1-cyano isomer, perhaps one starting from acyclic precursors. The absence of N-heterocycles such as pyrrole and pyridine is discussed in light of the astronomical finding., Comment: Accepted for publication in Nature Astronomy. 36 pages comprising five figures, two tables, five supplementary tables, and six supplementary figures

Benzonitrile ($c$-C$_6$H$_5$CN), a polar proxy for benzene ($c$-C$_6$H$_6$}), has the potential to serve as a highly convenient radio probe for aromatic chemistry, provided this ring can be found in other astronomical sources beyond the molecule-rich prestellar cloud TMC-1 where it was first reported by McGuire et al. in 2018. Here we present radio astronomical evidence of benzonitrile in four additional pre-stellar, and possibly protostellar, sources: Serpens 1A, Serpens 1B, Serpens 2, and MC27/L1521F. These detections establish benzonitrile is not unique to TMC-1; rather aromatic chemistry appears to be widespread throughout the earliest stages of star formation, likely persisting at least to the initial formation of a protostar. The abundance of benzonitrile far exceeds predictions from models which well reproduce the abundances of carbon chains, such as HC$_7$N, a cyanpolyyne with the same heavy atoms, indicating the chemistry responsible for planar carbon structures (as opposed to linear ones) in primordial sources is favorable but not well understood. The abundance of benzonitrile relative to carbon-chain molecules displays sizable variations between sources within the Taurus and Serpens clouds, implying the importance of physical conditions and initial elemental reservoirs of the clouds themselves., Comment: To Appear in Nature Astronomy, 53 Pages, 15 Figures, 12 Tables

As the inventory of interstellar molecules continues to grow, the gulf between small species, whose individual rotational lines can be observed with radio telescopes, and large ones, such as polycyclic aromatic hydrocarbons (PAHs) best studied in bulk via infrared and optical observations, is slowly being bridged. Understanding the connection between these two molecular reservoirs is critical to understanding the interstellar carbon cycle, but will require pushing the boundaries of how far we can probe molecular complexity while still retaining observational specificity. Toward this end, we present a method for detecting and characterizing new molecular species in single-dish observations toward sources with sparse line spectra. We have applied this method to data from the ongoing GOTHAM (GBT Observations of TMC-1: Hunting Aromatic Molecules) Green Bank Telescope (GBT) large program, discovering six new interstellar species. In this paper we highlight the detection of HC$_{11}$N, the largest cyanopolyyne in the interstellar medium., Comment: 42 pages, 36 figures, 9 tables, to appear in Nature Astronomy

We present an overview of the GOTHAM (GBT Observations of TMC-1: Hunting Aromatic Molecules) Large Program on the Green Bank Telescope. This and a related program were launched to explore the depth and breadth of aromatic chemistry in the interstellar medium at the earliest stages of star formation, following our earlier detection of benzonitrile ($c$-C$_6$H$_5$CN) in TMC-1. In this work, details of the observations, use of archival data, and data reduction strategies are provided. Using these observations, the interstellar detection of propargyl cyanide (HCCCH$_2$CN) is described, as well as the accompanying laboratory spectroscopy. We discuss these results, and the survey project as a whole, in the context of investigating a previously unexplored reservoir of complex, gas-phase molecules in pre-stellar sources. A series of companion papers describe other new astronomical detections and analyses., Comment: Accepted in the Astrophysical Journal Letters

We report an astronomical detection of HC$_4$NC for the first time in the interstellar medium with the Green Bank Telescope toward the TMC-1 molecular cloud with a minimum significance of $10.5 \sigma$. The total column density and excitation temperature of HC$_4$NC are determined to be $3.29^{+8.60}_{-1.20}\times 10^{11}$ cm$^{-2}$ and $6.7^{+0.3}_{-0.3}$ K, respectively, using the MCMC analysis. In addition to HC$_4$NC, HCCNC is distinctly detected whereas no clear detection of HC$_6$NC is made. We propose that the dissociative recombination of the protonated cyanopolyyne, HC$_5$NH$^+$, and the protonated isocyanopolyyne, HC$_4$NCH$^+$, are the main formation mechanisms for HC$_4$NC while its destruction is dominated by reactions with simple ions and atomic carbon. With the proposed chemical networks, the observed abundances of HC$_4$NC and HCCNC are reproduced satisfactorily., Comment: Accepted in the Astrophysical Journal Letters

We report a search for the diatomic hydrides SiH, PH, and FeH along the line of sight toward the chemically rich circumstellar envelopes of IRC+10216 and VY Canis Majoris. These molecules are thought to form in high temperature regions near the photospheres of these stars, and may then further react via gas-phase and dust-grain interactions leading to more complex species, but have yet to be constrained by observation. We used the GREAT spectrometer on SOFIA to search for rotational emission lines of these molecules in four spectral windows ranging from 600 GHz to 1500 GHz. Though none of the targeted species were detected in our search, we report their upper limit abundances in each source and discuss how they influence the current understanding of hydride chemistry in dense circumstellar media. We attribute the non-detections of these hydrides to their compact source sizes, high barriers of formation, and proclivity to react with other molecules in the winds., Comment: Accepted for publication in ApJ. 14 pages, 4 figures, 3 tables

In this work, we present the results of our investigation into the chemistry of Z- and E-cyanomethanimine (HNCHCN), both of which are possible precursors to the nucleobase adenine. Ab initio quantum chemical calculations for a number of reactions with atomic hydrogen were carried out. We find that the reaction H + Z/E-HNCHCN leading both to H-addition as well as H$_2$-abstraction proceed via similar short-range barriers with bimolecular rate coefficients on the order of $\sim10^{-17}$ cm$^{3}$ s$^{-1}$. These results were then incorporated into astrochemical models and used in simulations of the giant molecular cloud G+0.693. The calculated abundances obtained from these models were compared with previous observational data and found to be in good agreement, with a predicted [Z/E] ratio of $\sim3$ - somewhat smaller than the previously derived value of $6.1\pm2.4$. We find that the [Z/E] ratio in our simulations is due mostly to ion-molecule destruction rates driven by the different permanent dipoles of the two conformers. Based on these results, we propose a general rule-of-thumb for estimating the abundances of isomers in interstellar environments., Comment: The Astrophysical Journal, accepted: 17 pages, 9 figures, 7 tables

In this work, we reexamine sulfur chemistry occurring on and in the ice mantles of interstellar dust grains, and report the effects of two new modifications to standard astrochemical models; namely, (a) the incorporation of cosmic ray-driven radiation chemistry and (b) the assumption of fast, non-diffusive reactions for key radicals in the bulk. Results from our models of dense molecular clouds show that these changes can have a profound influence on the abundances of sulfur-bearing species in ice mantles, including a reduction in the abundance of solid-phase H$_2$S and HS, and a significant increase in the abundances of OCS, SO$_2$, as well as pure allotropes of sulfur, especially S$_8$. These pure-sulfur species - though nearly impossible to observe directly - have long been speculated to be potential sulfur reservoirs and our results represent possibly the most accurate estimates yet of their abundances in the dense ISM. Moreover, the results of these updated models are found to be in good agreement with available observational data. Finally, we examine the implications of our findings with regard to the as-yet-unknown sulfur reservoir thought to exist in dense interstellar environments., Comment: ApJ, accepted: 37 pages, 8 figures

A long standing problem in astrochemistry is the inability of many current models to account for missing sulfur content. Many relatively simple species that may be good candidates to sequester sulfur have not been measured experimentally at the high spectral resolution necessary to enable radioastronomical identification. On the basis of new laboratory data, we report searches for the rotational lines in the microwave, millimeter, and sub-millimeter regions of the sulfur-containing hydrocarbon HCCSH. This simple species would appear to be a promising candidate for detection in space owing to the large dipole moment along its $b$-inertial axis, and because the bimolecular reaction between two highly abundant astronomical fragments (CCH and SH radicals) may be rapid. An inspection of multiple line surveys from the centimeter to the far-infrared toward a range of sources from dark clouds to high-mass star-forming regions, however, resulted in non-detections. An analogous search for the lowest-energy isomer, H$_2$CCS, is presented for comparison, and also resulted in non-detections. Typical upper limits on the abundance of both species relative to hydrogen are $10^{-9}$-$10^{-10}$. We thus conclude that neither isomer is a major reservoir of interstellar sulfur in the range of environments studied. Both species may still be viable candidates for detection in other environments or at higher frequencies, providing laboratory frequencies are available., Comment: Accepted in the Astrophysical Journal

Shocks are a crucial probe for understanding the ongoing chemistry within ices on interstellar dust grains where many complex organic molecules (COMs) are believed to be formed. However, previous work has been limited to the initial liberation into the gas phase through non-thermal desorption processes such as sputtering. Here, we present results from the adapted three-phase gas-grain chemical network code NAUTILUS, with the inclusion of additional high-temperature reactions, non-thermal desorption, collisional dust heating, and shock-physics parameters. This enhanced model is capable of reproducing many of the molecular distributions and abundance ratios seen in our prior observations of the prototypical shocked-outflow L1157. In addition, we find that, among others, NH$_2$CHO, HCOOCH$_3$, and CH$_3$CHO have significant post-shock chemistry formation routes that differ from those of many other COMs observed in shocks. Finally, a number of selected species and phenomena are studied here with respect to their usefulness as shock tracers in various astrophysical sources., Comment: 20 pages, 10 figures, Accepted for publication in The Astrophysical Journal

To date, two isomers of H$_2$C$_3$O have been detected, namely, propynal (HCCCHO) and cylclopropenone (c-H$_2$C$_3$O). A third, propadienone (CH$_2$CCO), has thus far eluded observers despite the fact that it is the lowest in energy of the three. This previously noted result is in contradiction of the minimum energy principle, which posits that the abundances of isomers in interstellar environments can be predicted based on their relative stabilities - and suggests, rather, the importance of kinetic over thermodynamic effects in explaining the role of such species. Here, we report results of \textit{ab initio} quantum chemical calculations of the reaction between H and (a) HC$_3$O, (b) H$_2$C$_3$O (both propynal and propadienone), and (c) CH$_2$CHCO. We have found that, among all possible reactions between atomic hydrogen and either propadienone or propynal, only the destruction of propadienone is barrierless and exothermic. That this destruction pathway is indeed behind the non-detection of CH$_2$CCO is further suggested by our finding that the product of this process, the radical CH$_2$CHCO, can subsequently react barrierlessly with H to form propenal (CH$_2$CHCHO) which has, in fact, been detected in regions where the other two H$_2$C$_3$O isomers are observed. Thus, these results not only shed light on a previously unresolved astrochemical mystery, but also further highlight the importance of kinetics in understanding the abundances of interstellar molecules., Comment: ApJ, accepted: 14 pages, 2 figures

Many astrochemical models today explicitly consider the species that comprise the bulk of interstellar dust grain ice-mantles separately from those in the top few monolayers. Bombardment of these ices by ionizing radiation - whether in the form of cosmic rays, stellar winds, or radionuclide emission - represents an astrochemically viable means of driving a rich chemistry even in the bulk of the ice-mantle, now supported by a large body of work in laboratory astrophysics. In this study, using an existing rate equation-based astrochemical code modified to include a method of considering radiation chemistry recently developed by us, we attempted to simulate two such studies in which (a) pure O$_2$ ice at 5 K and, (b) pure H$_2$O ice at 16 K and 77 K, were bombarded by keV H$^+$ ions. Our aims are twofold: (1) to test the capability of our newly developed method to replicate the results of ice-irradiation experiments, and (2) to determine in such a well-constrained system how bulk chemistry is best handled using the same gas-grain codes that are used to model the interstellar medium (ISM). We find that our modified astrochemical model is able to reproduce both the abundance of O$_3$ in the 5 K pure O$_2$ ice, as well as both the abundance of H$_2$O$_2$ in the 16 K water ice and the previously noted decrease of hydrogen peroxide at higher temperatures. However, these results require the assumption that radicals and other reactive species produced via radiolysis react quickly and non-diffusively with neighbors in the ice., Comment: ApJ, accepted. 30 pages, 5 figures

One of the primary mechanisms for inferring the dynamical history of planets in our Solar System and in exoplanetary systems is through observation of elemental ratios (i.e. C/O). The ability to effectively use these observations relies critically on a robust understanding of the chemistry and evolutionary history of the observed abundances. Significant efforts have been devoted to this area from within astrochemistry circles, and these efforts should be supported going forward by the larger exoplanetary science community. In addition, the construction of a next-generation radio interferometer will be required to test many of these predictive models in situ, while simultaneously providing the resolution necessary to pinpoint the location of planets in formation., Comment: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)

We present the first results of a pilot program to conduct an ALMA Band 10 spectral line survey of the high-mass star-forming region NGC 6334I. The observations were taken in exceptional weather conditions (0.19 mm precipitable water) with typical system temperatures $T_{\rm{sys}}$ $<$950 K at $\sim$890 GHz. A bright, bipolar north-south outflow is seen in HDO and CS emission, driven by the embedded massive protostar MM1B. This has allowed, for the first time, a direct comparison of the thermal water in this outflow to the location of water maser emission from prior 22 GHz VLA observations. The maser locations are shown to correspond to the sites along the outflow cavity walls where high velocity gas impacts the surrounding material. We also compare our new observations to prior Herschel HIFI spectral line survey data of this field, detecting an order of magnitude more spectral lines (695 vs 65) in the ALMA data. We focus on the strong detections of the complex organic molecule glycolaldehyde (HC(O)CH$_2$OH) in the ALMA data that is not detected in the heavily beam-diluted HIFI spectra. Finally, we stress the need for dedicated THz laboratory spectroscopy to support and exploit future high-frequency molecular line observations with ALMA., Comment: Accepted in ApJ Letters

In this paper, we present preliminary results illustrating the effect of cosmic rays on solid-phase chemistry in models of both TMC-1 and several sources with physical conditions identical to TMC-1 except for hypothetically enhanced ionization rates. Using a recent theory for the addition of cosmic ray-induced reactions to astrochemical models, we calculated the radiochemical yields, called $G$ values, for the primary dust grain ice-mantle constituents. We show that the inclusion of this non-thermal chemistry can lead to the formation of complex organic molecules from simpler ice-mantle constituents, even under cold core conditions. In addition to enriching ice-mantles, we find that these new radiation-chemical processes can lead to increased gas-phase abundances as well, particularly for HOCO, NO$_2$, HC$_2$O, methyl formate (HCOOCH$_3$), and ethanol (CH$_3$CH$_2$OH). These model results imply that HOCO - and perhaps NO$_2$ - might be observable in TMC-1. Future detections of either of these two species in cold interstellar environments could provide strong support for the importance of cosmic ray-driven radiation chemistry. The increased gas-phase abundance of methyl formate can be compared with abundances achieved through other formation mechanisms such as pure gas-phase chemistry and three-body surface reactions., Comment: ApJ, accepted

Polycyclic aromatic hydrocarbons and polycyclic aromatic nitrogen heterocycles are thought to be widespread throughout the Universe, because these classes of molecules are probably responsible for the unidentified infrared bands, a set of emission features seen in numerous Galactic and extragalactic sources. Despite their expected ubiquity, astronomical identification of specific aromatic molecules has proven elusive. We present the discovery of benzonitrile ($c$-C$_6$H$_5$CN), one of the simplest nitrogen-bearing aromatic molecules, in the interstellar medium. We observed hyperfine-resolved transitions of benzonitrile in emission from the molecular cloud TMC-1. Simple aromatic molecules such as benzonitrile may be precursors for polycyclic aromatic hydrocarbon formation, providing a chemical link to the carriers of the unidentified infrared bands.

We report the detection of interstellar methoxymethanol (CH$_3$OCH$_2$OH) in ALMA Bands 6 and 7 toward the MM1 core in the high-mass star-forming region NGC 6334I at ~0.1" - 1" spatial resolution. A column density of 4(2) x $10^{18}$ cm$^{-2}$ at $T_{ex}$ = 200 K is derived toward MM1, ~34 times less abundant than methanol (CH$_3$OH), and significantly higher than predicted by astrochemical models. Probable formation and destruction pathways are discussed, primarily through the reaction of the CH$_3$OH photodissociation products, the methoxy (CH$_3$O) and hydroxymethyl (CH$_2$OH) radicals. Finally, we comment on the implications of these mechanisms on gas-phase vs grain-surface routes operative in the region, and the possibility of electron-induced dissociation of CH$_3$OH rather than photodissociation., Comment: Accepted in Astrophysical Journal Letters

We report the detection of the carbon-chain radical HC$_5$O for the first time in the interstellar medium toward the dark cloud TMC-1 using the 100 m Green Bank Telescope. We observe four hyperfine components of this radical in the $J = 17/2 \rightarrow 15/2$ rotational transition that originates from the $^2\Pi_{1/2}$ fine structure level of its ground state, and calculate an abundance of $n/n_{H_2}$ = $1.7\times 10^{-10}$, assuming an excitation temperature of $T_{ex} = 7$~K. No indication of HC$_3$O, HC$_4$O, HC$_6$O, is found in these or archival observations of the source, while we report tentative evidence for HC$_7$O. We compare calculated upper limits, and the abundance of HC$_5$O to predictions based on (1) the abundance trend of the analogous HC$_n$N family in TMC-1 and (2) a gas-grain chemical model. We find that the gas-grain chemical model well reproduces the observed abundance of HC$_5$O, as well as the upper limits of HC$_3$O, HC$_6$O, and HC$_7$O, but HC$_4$O is over produced. The prospects for astronomical detection of both shorter and longer HC$_n$O chains are discussed., Comment: Accepted in The Astrophysical Journal Letters

The collisions between high-energy ions and solids can result in significant physical and chemical changes to the material. These effects are potentially important for better understanding the chemistry of interstellar and planetary bodies, which are exposed to cosmic radiation and the solar wind, respectively, however, modeling such collisions on a detailed microscopic basis has thus far been largely unsuccessful. To that end, a new model, entitled CIRIS: the Chemistry of Ionizing Radiation in Solids, was created to calculate the physical and chemical effects of the irradiation of solid materials. With the new code, we simulate O2 ice irradiated with 100 keV protons. Our models are able to reproduce the measured ozone abundances of a previous experimental study, as well as independently predict the approximate thickness of the ice used in that work.

The chemistry of dense interstellar regions was analyzed using a time-dependent gas-grain astrochemical simulation and a new chemical network that incorporates deuterated chemistry taking into account nuclear spin-states for the hydrogen chemistry and its deuterated isotopologues. With this new network, the utility of the [HCO$^+$]/[DCO$^+$] abundance ratio as a probe of the cosmic ray ionization rate has been reexamined, with special attention paid to the effect of the initial value of the molecular hydrogen ortho-to-para ratio (OPR). After discussing the use of the probe for cold cores, we then compare our results with previous theoretical and observational results for a molecular cloud close to the supernova remnant W51C, which is thought to have an enhanced cosmic ray ionization rate $\zeta$ caused by the nearby $\gamma$-ray source. In addition, we attempt to use our approach to estimate the cosmic ray ionization rate for L1174, a dense core with an embedded star. Beyond the previously known sensitivity of [HCO$^+$]/[DCO$^+$] to $\zeta$, we demonstrate its additional dependence on the initial OPR and, secondarily, on the age of the source, its temperature, and its density. We conclude that the usefulness of the [HCO$^+$]/[DCO$^+$] abundance ratio to constrain the cosmic ray ionization rate in dense regions increases with source age and ionization rate as the ratio becomes far less sensitive to the initial value of the OPR., Comment: ApJ (accepted), 21 pages, 7 figures

Bell et al. (1997) reported the first detection of the cyanopolyyne HC$_{11}$N toward the cold dark cloud TMC-1; no subsequent detections have been reported toward any source. Additional observations of cyanopolyynes and other carbon-chain molecules toward TMC-1 have shown a log-linear trend between molecule size and column density, and in an effort to further explore the underlying chemical processes driving this trend, we have analyzed GBT observations of HC$_9$N and HC$_{11}$N toward TMC-1. Although we find an HC$_9$N column density consistent with previous values, HC$_{11}$N is not detected and we derive an upper limit column density significantly below that reported in Bell et al. Using a state-of-the-art chemical model, we have investigated possible explanations of non-linearity in the column density trend. Despite updating the chemical model to better account for ion-dipole interactions, we are not able to explain the non-detection of HC$_{11}$N, and we interpret this as evidence of previously unknown carbon-chain chemistry. We propose that cyclization reactions may be responsible for the depleted HC$_{11}$N abundance, and that products of these cyclization reactions should be investigated as candidate interstellar molecules., Comment: 10 pages, 7 figures, 3 tables, accepted for publication in MNRAS

L1157, a molecular dark cloud with an embedded Class 0 protostar possessing a bipolar outflow, is an excellent source for studying shock chemistry, including grain-surface chemistry prior to shocks, and post-shock, gas-phase processing. The L1157-B1 and B2 positions experienced shocks at an estimated ~2000 and 4000 years ago, respectively. Prior to these shock events, temperatures were too low for most complex organic molecules to undergo thermal desorption. Thus, the shocks should have liberated these molecules from the ice grain-surfaces en masse, evidenced by prior observations of SiO and multiple grain mantle species commonly associated with shocks. Grain species, such as OCS, CH3OH, and HNCO, all peak at different positions relative to species that are preferably formed in higher velocity shocks or repeatedly-shocked material, such as SiO and HCN. Here, we present high spatial resolution (~3") maps of CH3OH, HNCO, HCN, and HCO+ in the southern portion of the outflow containing B1 and B2, as observed with CARMA. The HNCO maps are the first interferometric observations of this species in L1157. The maps show distinct differences in the chemistry within the various shocked regions in L1157B. This is further supported through constraints of the molecular abundances using the non-LTE code RADEX (Van der Tak et al. 2007). We find the east/west chemical differentiation in C2 may be explained by the contrast of the shock's interaction with either cold, pristine material or warm, previously-shocked gas, as seen in enhanced HCN abundances. In addition, the enhancement of the HNCO abundance toward the the older shock, B2, suggests the importance of high-temperature O-chemistry in shocked regions., Comment: Accepted for publication in the Astrophysical Journal

Research over the past four decades has shown a rich variety of complex organic molecular content in some meteorites. This current study is an attempt to gain a better insight into the thermal conditions experienced by these molecules inside meteorites during atmospheric entry. In particular, we wish to understand possible chemical processes that can occur during entry and that might have had an effect on complex organic or prebiotic species that were delivered in this way to the early Earth. A simulation was written in Fortran to model heating by the shock generated during entry and the subsequent thermal diffusion inside the body of a meteorite. Experimental data was used for the thermal parameters of several types of meteorites, including iron-nickel and several classes of chondrites. A Sutton-Graves model of stagnation-point heating was used to calculate peak surface temperatures and an explicit difference formula was used to generate thermal diffusion profiles for both chondrites and iron-nickel type meteorites. Results from the simulation show pyrolytic temperature penetration to a depth of ca. 0.5 to 1 cm. Non-dissociative warming of meteorite interiors penetrates further to ca. 4 cm. These results support the findings that extraterrestrial delivery is a viable option for prebiotic molecular "seeding" of a planet., Comment: 22 pages, 4 figures, submitted to Astrobiology