We used the EMoCA survey data to search for isocyanides in Sgr B2(N2) and their corresponding cyanide analogs. We then used the coupled three-phase chemical kinetics code MAGICKAL to simulate their chemistry. Several new species, and over 100 new reactions have been added to the network. In addition, a new single-stage simultaneous collapse/warm-up model has been implemented, thus eliminating the need for the previous two-stage models. A variable, visual extinction-dependent $\zeta$ was also incorporated into the model and tested. We report the tentative detection of CH$_3$NC and HCCNC in Sgr B2(N2), which represents the first detection of both species in a hot core of Sgr B2. Our updated chemical models can reproduce most observed NC:CN ratios reasonably well depending on the physical parameters chosen. The model that performs best has an extinction-dependent cosmic-ray ionization rate that varies from ~2 $\times$ 10$^{-15}$ s$^{-1}$ at the edge of the cloud to ~1 $\times$ 10$^{-16}$ s$^{-1}$ in the center. Models with higher extinction-dependent $\zeta$ than this model generally do not agree as well, nor do models with a constant $\zeta$ greater than the canonical value of 1.3 $\times$ 10$^{-17}$ s$^{-1}$ throughout the source. Radiative transfer models are run using results of the best-fit chemical model. Column densities produced by the radiative transfer models are significantly lower than those determined observationally. Inaccuracy in the observationally determined density and temperature profiles is a possible explanation. Excitation temperatures are well reproduced for the true ``hot core'' molecules, but are more variable for other molecules such as HC$_3$N, for which fewer lines exist in ALMA Band 3.

Context. Propyne (CH$_3$CCH) has been detected in a variety of environments, from Galactic star-forming regions to extragalactic sources. Such molecules are excellent tracers of the physical conditions in star-forming regions. Aims. This study explores the emission of CH$_3$CCH in the low-mass protostellar binary, IRAS 16293$-$2422, examining the spatial scales traced by this molecule, as well as its formation and destruction pathways. Methods. ALMA observations from the Protostellar Interferometric Line Survey (PILS) are used to determine the abundances and excitation temperatures of CH$_3$CCH towards both protostars, exploring spatial scales from 70 to 2400 au. The three-phase chemical kinetics model MAGICKAL is also used, to explore the chemical reactions of this molecule. Results. CH$_3$CCH is detected towards both IRAS 16293A and IRAS 16293B and is found to trace the hot corino component around each source in the PILS dataset. Eighteen transitions above 3$\sigma$ are detected, enabling robust excitation temperatures and column densities to be determined in each source. In IRAS 16293A, an excitation temperature of 90 K and a column density of 7.8$\times$10$^{15}$ cm$^{-2}$ best fits the spectra. In IRAS 16293B, an excitation temperature of 100 K and 6.8$\times$10$^{15}$ cm$^{-2}$ best fits the spectra. The chemical modelling finds that in order to reproduce the observed abundances, both gas-phase and grain-surface reactions are needed. Conclusions. CH$_3$CCH is a molecule whose brightness and abundance in many different regions can be utilised to provide a benchmark of molecular variation with the physical properties of star-forming regions. It is essential when making such comparisons, that the abundances are determined with a good understanding of the spatial scale of the emitting region, to ensure that accurate abundances are derived.

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.

The relative column densities of the structural isomers methyl formate, glycolaldehyde, and acetic acid are derived for a dozen positions towards the massive star-forming regions MM1 and MM2 in the NGC 6334I complex, which are separated by $\sim$4000 AU. Relative column densities of these molecules are also gathered from the literature for 13 other star-forming regions. In this combined dataset, a clear bi-modal distribution is observed in the relative column densities of glycolaldehyde and methyl formate. No such distribution is evident with acetic acid. The two trends are comprised of star-forming regions with a variety of masses, suggesting that there must be some other common parameter that is heavily impacting the formation of glycolaldehyde. This is indicative of some demonstrable differentiation in these cores; studying the abundances of these isomers may provide a clue as to the integral chemical processes ongoing in a variety of protostellar environments.

As the number of complex organic molecules (COMs) detected in the interstellar medium increases, it becomes important to place meaningful constraints on the formation pathways of these species. The molecular cloud SgrB2(N) is host to several hot molecular cores in the early stage of star formation, where a great variety of COMs are detected in the gas phase. Because of its exposure to the extreme conditions of the the Galactic center region, SgrB2(N) is one of the best targets to study the impact of environmental conditions on the production of COMs. Our main goal is to characterize the physico-chemical evolution of SgrB2(N)'s sources in order to explain their chemical differences and constrain their environmental conditions. The chemical composition of SgrB2(N)'s hot cores, N2, N3, N4, and N5 is derived by modeling their 3mm emission spectra extracted from the EMoCA imaging spectral line survey performed with ALMA. We derive the density distribution in the envelope of the sources based on the masses computed from the ALMA dust continuum emission maps. We use the radiative transfer code RADMC-3D to compute temperature profiles based on the COM rotational temperatures derived from population diagrams. We use published results of 3D RMHD simulations of high-mass star formation to estimate the time evolution of the sources properties. We employ the chemical code MAGICKAL to compute time-dependent chemical abundances in the sources and investigate how physical properties and environmental conditions influence the production of COMs. We find that chemical models with a cosmic-ray ionization rate of 7e-16s-1 best reproduce the abundances with respect to methanol of ten COMs observed toward SgrB2(N2-N5). We also show that COMs still form efficiently on dust grains with minimum dust temperatures in the prestellar phase as high as 15K, but that minimum temperatures higher than 25K are excluded.

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.

Methyl isocyanide (CH$_3$NC) is the isocyanide with the largest number of atoms confirmed in the interstellar medium (ISM), but it is not an abundant molecule, having only been detected towards a handful of objects. Conversely, its isomer, methyl cyanide (CH$_3$CN), is one of the most abundant complex organic molecules detected in the ISM, with detections in a variety of low- and high-mass sources. We use ALMA observations from the Protostellar Interferometric Line Survey (PILS) to search for methyl isocyanide and compare its abundance with that of its isomer methyl cyanide. We use a new line catalogue from the Cologne Database for Molecular Spectroscopy (CDMS) to identify methyl isocyanide lines. We also model the chemistry with an updated version of the three-phase chemical kinetics model {\em MAGICKAL}, presenting the first chemical modelling of methyl isocyanide to date. We detect methyl isocyanide for the first time in a solar-type protostar, IRAS 16293$-$2422 B, and present upper limits for its companion protostar, IRAS 16293$-$2422 A. Methyl isocyanide is found to be at least 20 times more abundant in source B compared to source A, with a CH$_3$CN/CH$_3$NC abundance ratio of 200 in IRAS 16293--2422 B and >5517 in IRAS 16293$-$2422 A. We also present the results of a chemical model of methyl isocyanide chemistry in both sources, and discuss the implications on methyl isocyanide formation mechanisms and the relative evolutionary stages of both sources.

Searches for the prebiotically-relevant cyanamide (NH$_2$CN) towards solar-type protostars have not been reported in the literature. We here present the first detection of this species in the warm gas surrounding two solar-type protostars, using data from the Atacama Large Millimeter/Submillimeter Array Protostellar Interferometric Line Survey (PILS) of IRAS 16293-2422 B and observations from the IRAM Plateau de Bure Interferometer of NGC1333 IRAS2A. We furthermore detect the deuterated and $^{13}$C isotopologues of NH$_2$CN towards IRAS 16293-2422 B. This is the first detection of NHDCN in the interstellar medium. Based on a local thermodynamic equilibrium analysis, we find that the deuteration of cyanamide ($\sim$ 1.7%) is similar to that of formamide (NH$_2$CHO), which may suggest that these two molecules share NH$_2$ as a common precursor. The NH$_2$CN/NH$_2$CHO abundance ratio is about 0.2 for IRAS 16293-2422 B and 0.02 for IRAS2A, which is comparable to the range of values found for Sgr B2. We explored the possible formation of NH$_2$CN on grains through the NH$_2$ + CN reaction using the chemical model MAGICKAL. Grain-surface chemistry appears capable of reproducing the gas-phase abundance of NH$_2$CN with the correct choice of physical parameters.

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.

Back-diffusion is the phenomenon by which random walkers revisit binding sites on a lattice. This phenomenon must occur on interstellar dust particles, slowing down dust-grain reactions, but it is not accounted for by standard rate-equation models. Microscopic kinetic Monte Carlo models have been used to investigate the effect of back-diffusion on reaction rates on interstellar dust grains. Grain morphology, size, and grain-surface coverage were varied and the effects of these variations on the magnitude of the back-diffusion effect were studied for the simple H+H reaction system. This back-diffusion effect is seen to reduce reaction rates by a maximum factor of ~5 for the canonical grain of 10$^6$ binding sites.The resulting data were fit to logarithmic functions that can be used to reproduce the effects of back-diffusion in rate-equation models.

Rate-equation models are a widely-used and inexpensive tool for the simulation of interstellar chemistry under a range of physical conditions. However, their application to grain-surface chemical systems necessitates a number of simplifying assumptions, due to the requirement to treat only the total population of each species, using averaged rates, rather than treating each surface particle as an independent entity. While the outputs from rate-equation models are strictly limited to such population information, the inputs -- in the form of the averaged rates that control the time-evolution of chemical populations -- can be guided by the results from more exact simulation methods. Here, we examine the effects of back-diffusion, wherein particles diffusing on a surface revisit binding sites on the lattice, slowing the total reaction rate. While this effect has been studied for two-particle systems, its influence at greater surface coverage of reactants has not been explored. Results from two Monte Carlo kinetics models (one a 2-D periodic lattice, the other the surface of a three dimensionally-realized grain) were used to develop a means to incorporate the grain-surface back-diffusion effect into rate-equation methods. The effects of grain size, grain morphology, and surface coverage on the magnitude of the back-diffusion effect were studied for the simple H+H reaction system. The results were fit with expressions that can be easily incorporated into astrochemical rate-equation models to reproduce accurately the effects of back-diffusion on grain-surface reaction rates. Back-diffusion reduces reaction rates by a maximum factor of around 5 for the canonical grain of $\sim$10$^6$ surface sites, but this falls to unity at close to full surface coverage.