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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

A deep search for the potential glycine precursor hydroxylamine (NH$_2$OH) using the Caltech Submillimeter Observatory (CSO) at $\lambda = 1.3$ mm and the Combined Array for Research in Millimeter-wave Astronomy (CARMA) at $\lambda = 3$ mm is presented toward the molecular outflow L1157, targeting the B1 and B2 shocked regions. We report non-detections of NH$_2$OH in both sources. We a perform non-LTE analysis of CH$_3$OH observed in our CSO spectra to derive kinetic temperatures and densities in the shocked regions. Using these parameters, we derive upper limit column densities of NH$_2$OH of $\leq1.4 \times 10^{13}$~cm$^{-2}$ and $\leq1.5 \times 10^{13}$~cm$^{-2}$ toward the B1 and B2 shocks, respectively, and upper limit relative abundances of $N_{NH_2OH}/N_{H_2} \leq1.4 \times 10^{-8}$ and $\leq1.5 \times 10^{-8}$, respectively.

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.

Recently, Lattelais et al. (2009) have interpreted aggregated observations of molecular isomers to suggest that there exists a "minimum energy principle'', such that molecular formation will favor more stable molecular isomers for thermodynamic reasons. To test the predictive power of this principle, we have fully characterized the spectra of the three isomers of C$_{3}$H$_{2}$O toward the well known molecular region Sgr B2(N). Evidence for the detection of the isomers cyclopropenone (c-C$_{3}$H$_{2}$O) and propynal (HCCCHO) is presented, along with evidence for the non-detection of the lowest zero-point energy isomer, propadienone (CH$_2$CCO). We interpret these observations as evidence that chemical formation pathways, which may be under kinetic control, have a more pronounced effect on final isomer abundances than thermodynamic effects such as the minimum energy principle.