Aminoacetonitrile, a molecule detected in the interstellar medium only towards the star-forming region Sagittarius B2 (Sgr B2) thus far, is considered an important prebiotic species. To date, observations were limited to ground state emission lines, whereas transitions from within vibrationally excited states remained undetected. We wanted to accurately determine the energies of the low-lying vibrational states of aminoacetonitrile, which are expected to be populated in Sgr B2(N1), the main hot core of Sgr B2(N). This step is fundamental in order to properly evaluate the vibration-rotation partition function of aminoacetonitrile as well as the line strengths of the rotational transitions of its vibrationally excited states. This is necessary to derive accurate column densities and secure the identification of these transitions in astronomical spectra. The far-infrared ro-vibrational spectrum of aminoacetonitrile has been recorded in absorption against a synchrotron source of continuum emission. Three bands, corresponding to the lowest vibrational modes of aminoacetonitrile, were observed in the frequency region below 500 cm$^{-1}$. The combined analysis of ro-vibrational and pure rotational data allowed us to prepare new spectral line catalogs for all the states under investigation. We used the imaging spectral line survey ReMoCA performed with ALMA to search for vibrationally excited aminoacetonitrile toward Sgr B2(N1). On the basis of these spectroscopic predictions, we report the interstellar detection of aminoacetonitrile in its $v_{11}=1$ and $v_{18}=1$ vibrational states toward Sgr B2(N1) in addition to emission in its vibrational ground state. The intensities of the identified $v_{11}=1$ and $v_{18}=1$ lines are consistent with the detected $v=0$ lines under LTE at a temperature of 200 K for an aminoacetonitrile column density of $1.1 \times 10^{17}$ cm$^{-2}$.

Glycolamide is a glycine isomer and also one of the simplest derivatives of acetamide (e.g., one hydrogen atom is replaced with a hydroxyl group), which is a known interstellar molecule. Using a battery of state of the art rotational spectroscopic techniques in the frequency and time domain, around 1500 transitions have been newly assigned. Based on the reliable frequency predictions, we report a radioastronomical search for glycolamide in the well known high-mass star forming region Sgr B2(N) using the ALMA imaging spectral line survey ReMoCA. We also searched for glycolamide toward Sgr B2(N) with the Effelsberg radio telescope. We report the nondetection of glycolamide toward this source with an abundance at least six and five times lower than that of acetamide and glycolaldehyde, respectively. Our astrochemical model suggests that glycolamide may be present in this source at a level just below the upper limit, which was derived from the EMoCA survey. We could also not detect the molecule in the region's extended molecular envelope, which was probed with the Effelsberg telescope. We find an upper limit to its column density that is similar to the column densities obtained earlier for acetamide and glycolaldehyde with the Green Bank Telescope.

A prevailing theory for the interstellar production of complex organic molecules (COMs) involves formation on warm dust-grain surfaces, via the diffusion and reaction of radicals produced through grain-surface photodissociation of stable molecules. However, some gas-phase O-bearing COMs, notably acetaldehyde(CH$_3$CHO), methyl formate(CH$_3$OCHO), and dimethyl ether(CH$_3$OCH$_3$), are now observed at very low temperatures, challenging the warm scenario. Here, we introduce a selection of new non-diffusive mechanisms into an astrochemical model, to account for the failure of the standard diffusive picture and to provide a more generalized scenario of COM formation on interstellar grains. New generic rate formulations are provided for cases where: (i) radicals are formed by reactions occurring close to another reactant, producing an immediate follow-on reaction; (ii) radicals are formed in an excited state, allowing them to overcome activation barriers to react with nearby stable molecules; (iii) radicals are formed through photo-dissociation close to a reaction partner, followed by immediate reaction. Each process occurs without the diffusion of large radicals. The new mechanisms significantly enhance cold COM abundances, successfully reproducing key observational results for prestellar core L1544. H-abstraction from grain-surface COMs, followed by recombination, plays a crucial role in amplifying chemical desorption into the gas phase. The UV-induced chemistry produces significant COM abundances in the bulk ices, which are retained on the grains and may persist to later stages. O$_2$ is also formed strongly in the mantle though photolysis, suggesting cometary O$_2$ could indeed be interstellar.

Star-forming regions show a rich and varied chemistry, including the presence of complex organic molecules - both in the cold gas distributed on large scales, and in the hot regions close to young stars where protoplanetary disks arise. Recent advances in observational techniques have opened new possibilities for studying this chemistry. In particular, the Atacama Large Millimeter/submillimeter Array (ALMA) has made it possible to study astrochemistry down to Solar System size scales, while also revealing molecules of increasing variety and complexity. In this review, we discuss recent observations of the chemistry of star-forming environments, with a particular focus on complex organic molecules, taking context from the laboratory experiments and chemical models that they have stimulated. The key takeaway points are: The physical evolution of individual sources plays a crucial role in their inferred chemical signatures, and remains an important area for observations and models to elucidate. Comparisons of the abundances measured toward different star-forming environments (high-mass versus low-mass, Galactic center versus Galactic disk) reveal a remarkable similarity, an indication that the underlying chemistry is relatively independent of variations in their physical conditions. Studies of molecular isotopologs in star-forming regions provide a link with measurements in our own Solar System, and thus may shed light on the chemical similarities and differences expected in other planetary systems.

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.

Methyl mercaptan (CH3SH) is a known interstellar molecule with abundances high enough that the detection of some of its minor isotopologues is promising. The present study aims at providing accurate spectroscopic parameters for the (13)CH3SH isotopologue to facilitate its identification in the interstellar medium at millimetre and submillimetre wavelengths. Through careful analysis of recent CH3SH spectra from 49-510 GHz and 1.1-1.5 THz recorded at natural isotopic composition, extensive assignments were possible not only for the ground torsional state of (13)CH3SH, but also in the first and second excited states. The torsion-rotation spectrum displays complex structure due to the large-amplitude internal rotation of the (13)CH3 group, similar to the main and other minor isotopic species of methyl mercaptan. The assigned transition frequencies have been fitted to within experimental error with a 52-parameter model employing the RAM36 programme. With predictions based on this fit, (13)CH3SH was searched for in spectra from the Atacama Large Millimetre/submillimetre Array (ALMA) towards the Galactic centre source Sgr B2(N2). Several transitions were expected to be observable, but all of them turned out to be severely blended with emission from other species, which prevents us from identifying (13)CH3SH in this source.

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.

The first computational model of solid-phase chemistry in cometary nuclear ices is presented. An astrochemical kinetics model, MAGICKAL, is adapted to trace the chemical evolution in multiple layers of cometary ice, over a representative period of 5 Gyr. Physical conditions are chosen appropriate for "cold storage" of the cometary nucleus in the outer Solar System, prior to any active phase. The chemistry is simulated at a selection of static temperatures in the range 5 - 60 K, while the ice is exposed to the interstellar radiation field, inducing a photochemistry in the outer ice layers that produces significant formation of complex organic molecules. A treatment for the chemistry resulting from cosmic-ray bombardment of the ices is also introduced into the model, along with a new formulation for low-temperature photochemistry. Production of simple and complex molecules to depth on the order of 10~m or more is achieved, with local fractional abundances comparable to observed values in many cases. The production of substantial amounts of O$_2$ (and H$_2$O$_2$) is found, suggesting that long-term processing by high-energy cosmic rays of cometary ices in situ, over a period on the order of 1 Gyr, may be sufficient to explain the large observed abundances of O$_2$, if the overall loss of material from the comet is limited to a depth on the order of 10 m. Entry into the inner solar system could produce a further enhancement in the molecular content of the nuclear ices that may be quantifiable using this modeling approach.

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.

We want to study the rotational spectra of CH$_3^{13}$C(O)CH$_3$ and $^{13}$CH$_3$C(O)CH$_3$ and search for them in Sagittarius B2(N2). We investigated the laboratory rotational spectrum of isotopically enriched CH$_3^{13}$C(O)CH$_3$ between 40 GHz and 910 GHz and of acetone between 36 GHz and 910 GHz in order to study $^{13}$CH$_3$C(O)CH$_3$ in natural isotopic composition. In addition, we searched for emission lines produced by these species in a molecular line survey of Sagittarius B2(N) carried out with ALMA. Discrepancies between predictions of the main isotopic species and the ALMA spectrum prompted us to revisit the rotational spectrum of this isotopolog. We assigned 9711 new transitions of CH$_3^{13}$C(O)CH$_3$ and 63 new transitions of $^{13}$CH$_3$C(O)CH$_3$ in the laboratory spectra. More than 1000 additional lines were assigned for the main isotopic species. We modeled the ground state data of all three isotopologs satisfactorily with the ERHAM program. We find that models of the torsionally excited states $v _{12} = 1$ and $v _{17} = 1$ of CH$_3$C(O)CH$_3$ improve only marginally. No transition of CH$_3^{13}$C(O)CH$_3$ is clearly detected toward the hot molecular core Sgr B2(N2). However, we report a tentative detection of $^{13}$CH$_3$C(O)CH$_3$ with a $^{12}$C/$^{13}$C isotopic ratio of 27 that is consistent with the ratio previously measured for alcohols in this source. Several dozens of transitions of both torsional states of the main isotopolog are detected as well. Our predictions of CH$_3^{13}$C(O)CH$_3$ and CH$_3$C(O)CH$_3$ are reliable into the terahertz region. The spectrum of $^{13}$CH$_3$C(O)CH$_3$ should be revisited in the laboratory with an enriched sample. Transitions pertaining to the torsionally excited states $v _{12} = 1$ and $v _{17} = 1$ of CH$_3$C(O)CH$_3$ could be identified unambiguously in Sagittarius B2(N2).

Urea, NH2C(O)NH2, is a molecule of great importance in organic chemistry and biology. Two searches for urea in the interstellar medium were reported in the past, but neither were conclusive. We want to take advantage of the increased sensitivity and angular resolution provided by ALMA to search for urea toward the hot cores embedded in the high-mass star forming region Sgr B2(N). We use the new spectral line survey ReMoCA performed toward Sgr B2(N) with ALMA in its observing cycle 4. The spectra are analyzed under the LTE approximation. We construct a full synthetic spectrum that includes all the molecules identified so far. We use new spectroscopic predictions for urea in its vibrational ground state and first vibrationally excited state to search for this complex organic molecule in the ReMoCA data set. We employ the gas-grain chemical kinetics model MAGICKAL to interpret the astronomical observations. We report the secure detection of urea toward the hot core Sgr B2(N1) at a position called N1S slightly offset from the continuum peak, which avoids obscuration by the dust.. We derive a column density of 2.7x10^16 cm-2 for urea, two orders of magnitude lower than formamide, and one order of magnitude below methyl isocyanate, acetamide, and N-methylformamide. The latter molecule is reliably identified toward N1S with 60 clearly detected lines, confirming an earlier claim of its tentative interstellar detection. We report the first interstellar detections of NH2CH18O and 15NH2CHO. We also report the nondetection of urea toward the secondary hot core Sgr B2(N2) with an abundance relative to the other four species at least one order of magnitude lower than toward the main hot core. Our chemical model roughly reproduces the relative abundances of formamide, methyl isocyanate, acetamide, and N-methylformamide, but it overproduces urea by at least one order of magnitude.

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.

Methyl mercaptan (CH$_3$SH) is an important sulfur-bearing species in the interstellar medium, terrestrial environment, and potentially in planetary atmospheres. The aim of the present study is to provide accurate spectroscopic parameters for the most abundant minor isotopolog CH$_3$$^{34}$SH to support radio astronomical observations at millimeter and submillimeter wavelengths. The rotational spectrum of CH$_3$$^{34}$SH, which is complicated by the large-amplitude internal rotation of the CH$_3$ group versus the $^{34}$SH frame, was investigated in the 49$-$510 GHz and 1.1$-$1.5 THz frequency ranges in natural isotopic abundance. The analysis of the spectrum was performed up to the second excited torsional state, and the obtained data were modeled with the RAM36 program. A fit within experimental accuracy was obtained with a RAM Hamiltonian model that uses 72 parameters. Predictions based on this fit are used to search for CH$_3$$^{34}$SH with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the hot molecular core Sgr B2(N2), but blends with emission lines of other species prevent its firm identification in this source.

The need to characterize ices coating dust grains in dense interstellar clouds arises from the importance of ice morphology in facilitating the diffusion and storage of radicals and reaction products in ices, a well-known place for the formation of complex molecules. Yet, there is considerable uncertainty about the structure of ISM ices, their ability to store volatiles and under what conditions. We measured the infrared absorption spectra of CO on the pore surface of porous amorphous solid water (ASW), and quantified the effective pore surface area of ASW. Additionally, we present results obtained from a Monte Carlo model of ASW in which the morphology of the ice is directly visualized and quantified. We found that 200 ML of ASW annealed to 20 K has a total pore surface area that is equivalent to 46 ML. This surface area decreases linearly with temperature to about 120 K. We also found that (1) dangling OH bonds only exist on the surface of pores; (2) almost all of the pores in the ASW are connected to the vacuum--ice interface, and are accessible for adsorption of volatiles from the gas phase; there are few closed cavities inside ASW at least up to a thickness of 200 ML; (3) the total pore surface area is proportional to the total 3-coordinated water molecules in the ASW in the temperature range 60--120 K. We also discuss the implications on the structure of ASW and surface reactions in the ice mantle in dense clouds.

We have extended the pure rotational investigation of the two isomers syn and anti vinyl mercaptan to the millimeter domain using a frequency-multiplication spectrometer. The species were produced by a radiofrequency discharge in 1,2-ethanedithiol. Additional transitions have been re-measured in the centimeter band using Fourier-transform microwave spectroscopy to better determine rest frequencies of transitions with low-$J$ and low-$K_a$ values. Experimental investigations were supported by quantum chemical calculations on the energetics of both the [C$_2$,H$_4$,S] and [C$_2$,H$_4$,O] isomeric families. Interstellar searches for both syn and anti vinyl mercaptan as well as vinyl alcohol were performed in the EMoCA (Exploring Molecular Complexity with ALMA) spectral line survey carried out toward Sagittarius (Sgr) B2(N2) with the Atacama Large Millimeter/submillimeter Array (ALMA). Highly accurate experimental frequencies (to better than 100 kHz accuracy) for both syn and anti isomers of vinyl mercaptan have been measured up to 250 GHz; these deviate considerably from predictions based on extrapolation of previous microwave measurements. Reliable frequency predictions of the astronomically most interesting millimeter-wave lines for these two species can now be derived from the best-fit spectroscopic constants. From the energetic investigations, the four lowest singlet isomers of the [C$_2$,H$_4$,S] family are calculated to be nearly isoenergetic, which makes this family a fairly unique test bed for assessing possible reaction pathways. Upper limits for the column density of syn and anti vinyl mercaptan are derived toward the extremely molecule-rich star-forming region Sgr B2(N2) enabling comparison with selected complex organic molecules.

The diffuse and translucent molecular clouds traced in absorption along the line of sight to strong background sources have so far been investigated mainly in the spectral domain because of limited angular resolution or small sizes of the background sources. We aim to resolve and investigate the spatial structure of molecular clouds traced by several molecules detected in absorption along the line of sight to SgrB2(N). We have used spectral line data from the EMoCA survey performed with ALMA, taking advantage of the high sensitivity and angular resolution. We identify, on the basis of the spectral analysis of c-C3H2 across the field of view, 15 main velocity components along the line of sight to SgrB2(N) and several components in the envelope of SgrB2. The c-C3H2 column densities reveal two categories of clouds. Clouds in Category I (3 kpc arm, 4 kpc arm, and some GC clouds) have smaller c-C3H2 column densities, smaller linewidths, and smaller widths of their column density PDFs than clouds in Category II (Scutum arm, Sgr arm, and other GC clouds). To investigate the spatial structure we derive opacity maps for the following molecules: c-C3H2, H13CO+, 13CO, HNC, HN13C, HC15N, CS, C34S, 13CS, SiO, SO, and CH3OH. These maps reveal that most molecules trace relatively homogeneous structures that are more extended than the field of view defined by the background continuum emission (about 15", that is 0.08-0.6pc depending on the distance). SO and SiO show more complex structures with smaller clumps of size ~5-8". Our analysis suggests that the driving of the turbulence is mainly solenoidal in the investigated clouds. On the basis of HCO+, we conclude that most line-of-sight clouds towards SgrB2 are translucent, including all clouds where complex organic molecules were recently detected. We also conclude that CCH and CH are good probes of H2 in both diffuse and translucent clouds.

In dense starless and protostellar cores, the relative abundance of deuterated species to their non-deuterated counterparts can become orders of magnitude greater than in the local interstellar medium. This enhancement proceeds through multiple pathways in the gas phase and on dust grains, where the chemistry is strongly dependent on the physical conditions. In this Chapter, we discuss how sensitive, high resolution observations with the ngVLA of emission from deuterated molecules will trace both the dense gas structure and kinematics on the compact physical scales required to track the gravitational collapse of star-forming cores and the subsequent formation of young protostars and circumstellar accretion regions. Simultaneously, such observations will play a critical role in tracing the chemical history throughout the various phases of star and planet formation. Many low-J transitions of key deuterated species, along with their undeuterated counterparts, lie within the 60-110 GHz frequency window, the lower end of which is largely unavailable with current facilities and instrumentation. The combination of sensitivity and angular resolution provided only by the ngVLA will enable unparalleled detailed studies of the physics and chemistry of the earliest stages of star formation.

Extraterrestrial amino acids, the chemical building blocks of the biopolymers that comprise life as we know it on Earth are present in meteoritic samples. More recently, glycine (NH$_2$CH$_2$COOH), the simplest amino acid, was detected by the Rosetta mission in comet 67P. Despite these exciting discoveries, our understanding of the chemical and physical pathways to the formation of (pre)biotic molecules is woefully incomplete. This is largely because our knowledge of chemical inventories during the different stages of star and planet formation is incomplete. It is therefore imperative to solidify our accounting of the chemical inventories, especially of critical yet low-abundance species, in key regions and to use this knowledge to inform, expand, and constrain chemical models of these reactions. This is followed naturally by a requirement to understand the spatial distribution and temporal evolution of this inventory. Here, we briefly outline a handful of particularly-impactful use cases in which the ngVLA will drive the field forward.

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.

Complex organic molecules are readily detected in the inner regions of the gaseous envelopes of forming protostars. In particular, molecules that contain nitrogen are interesting due to the role nitrogen plays in the development of life and the compact scales such molecules have been found to trace around forming protostars. The goal of this work is to determine the inventory of one family of nitrogen-bearing organic molecules, complex nitriles (molecules with a $-$CN functional group) towards two hot corino sources in the low-mass protostellar binary IRAS 16293$-$2422. This work explores the abundance differences between the two sources, the isotopic ratios, and the spatial extent derived from molecules containing the nitrile functional group. Using data from the Protostellar Interferometric Line Survey (PILS) obtained with ALMA we determine abundances and excitation temperatures for the detected nitriles. We also present a new method for determining the spatial structure of sources with high line density and large velocity gradients $-$ Velocity-corrected INtegrated emission (VINE) maps. We detect methyl cyanide (CH$_3$CN) as well as 5 of its isotopologues, including the detection of CHD$_2$CN which is the first detection in the ISM. We also detect ethyl cyanide (C$_2$H$_5$CN), vinyl cyanide (C$_2$H$_3$CN), and cyanoacetylene (HC$_3$N). We find that abundances are similar between IRAS 16293A and IRAS 16293B on small scales except for vinyl cyanide which is only detected towards the latter source. This suggests an important difference between the sources either in their evolutionary stage or warm-up timescales. We also detect a spatially double-peaked emission for the first time in molecular emission in the A source, suggesting that this source is showing structure related to a rotating toroid of material.

Massive young stellar objects in the Magellanic Clouds show infrared absorption features corresponding to significant abundances of CO, CO$_2$ and H$_2$O ice along the line of sight, with the relative abundances of these ices differing between the Magellanic Clouds and the Milky Way. CO ice is not detected towards sources in the Small Magellanic Cloud, and upper limits put its relative abundance well below sources in the Large Magellanic Cloud and the Milky Way. We use our gas-grain chemical code MAGICKAL, with multiple grain sizes and grain temperatures, and further expand it with a treatment for increased interstellar radiation field intensity to model the elevated dust temperatures observed in the MCs. We also adjust the elemental abundances used in the chemical models, guided by observations of HII regions in these metal-poor satellite galaxies. With a grid of models, we are able to reproduce the relative ice fractions observed in MC massive young stellar objects (MYSOs), indicating that metal depletion and elevated grain temperature are important drivers of the MYSO envelope ice composition. Magellanic Cloud elemental abundances have a sub-galactic C/O ratio, increasing H$_2$O ice abundances relative to the other ices; elevated grain temperatures favor CO$_2$ production over H$_2$O and CO. The observed shortfall in CO in the Small Magellanic Cloud can be explained by a combination of reduced carbon abundance and increased grain temperatures. The models indicate that a large variation in radiation field strength is required to match the range of observed LMC abundances. CH$_3$OH abundance is found to be enhanced in low-metallicity models, providing seed material for complex organic molecule formation in the Magellanic Clouds.

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.

HC$_3$N is an ubiquitous molecule in interstellar environments, from external galaxies, to Galactic interstellar clouds, star forming regions, and planetary atmospheres. Observations of its rotational and vibrational transitions provide important information on the physical and chemical structure of the above environments. We present the most complete global analysis of the spectroscopic data of HC$_3$N. We have recorded the high-resolution infrared spectrum from 450 to 1350 cm$^{-1}$, a region dominated by the intense $\nu_5$ and $\nu_6$ fundamental bands, located at 660 and 500 cm$^{-1}$, respectively, and their associated hot bands. Pure rotational transitions in the ground and vibrationally excited states have been recorded in the millimetre and sub-millimetre regions in order to extend the frequency range so far considered in previous investigations. All the transitions from the literature and from this work involving energy levels lower than 1000 cm$^{-1}$ have been fitted together to an effective Hamiltonian. Because of the presence of various anharmonic resonances, the Hamiltonian includes a number of interaction constants, in addition to the conventional rotational and vibrational l-type resonance terms. The data set contains about 3400 ro-vibrational lines of 13 bands and some 1500 pure rotational lines belonging to 12 vibrational states. More than 120 spectroscopic constants have been determined directly from the fit, without any assumption deduced from theoretical calculations or comparisons with similar molecules. An extensive list of highly accurate rest frequencies has been produced to assist astronomical searches and data interpretation. These improved data, have enabled a refined analysis of the ALMA observations towards Sgr B2(N2).

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.

Up to now, mostly relatively simple molecules have been detected in interstellar diffuse molecular clouds in our galaxy, but more complex species have been reported in the diffuse/translucent medium of a z = 0.89 spiral galaxy. We aim at searching for complex organic molecules (COMs) in diffuse molecular clouds along the line of sight to Sgr B2(N), taking advantage of the high sensitivity and angular resolution of the Atacama Large Millimeter/submillimeter Array (ALMA). We use data acquired as part of the EMoCA survey performed with ALMA. To analyse the absorption features of the molecules detected towards the ultracompact HII region K4 in Sgr B2(N), we calculate synthetic spectra for these molecules and fit their column densities, line widths, centroid velocities, and excitation temperatures. We report the detection of CH3OH, CH3CN, CH3CHO, HC3N, and NH2CHO in Galactic center (GC) diffuse clouds and CH3OH and CH3CN in a diffuse cloud in the Scutum arm. The chemical composition of one of the diffuse GC clouds is found to be similar to the one of the diffuse/translucent medium of the z=0.89 spiral galaxy. The chemical processes leading to chemical complexity in the diffuse molecular ISM appear to have remained similar since z=0.89. As proposed in previous studies, the presence of COMs in diffuse molecular clouds may result from a cyclical interstellar process of cloud contraction and expansion between diffuse and dense states.

Complex organic molecules (COMs) have been observed towards several low-mass young stellar objects (LYSOs). Small and heterogeneous samples have so far precluded conclusions on typical COM abundances, as well as the origin(s) of abundance variations between sources. We present observations towards 16 deeply embedded (Class 0/I) low-mass protostars using the IRAM 30m telescope. We detect CH$_2$CO, CH$_3$CHO, CH$_3$OCH$_3$, CH$_3$OCHO, CH$_3$CN, HNCO, and HC$_3$N towards 67%, 37%, 13%, 13%, 44%, 81%, and 75% of sources respectively. Median column densities derived using survival analysis range between 6.0x10$^{10}$ cm$^{-2}$ (CH$_3$CN) and 2.4x10$^{12}$ cm$^{-2}$ (CH$_3$OCH$_3$) and median abundances range between 0.48% (CH$_3$CN) and 16% (HNCO) with respect to CH$_3$OH. Column densities for each molecule vary by about one order of magnitude across the sample. Abundances with respect to CH$_3$OH are more narrowly distributed, especially for oxygen-bearing species. We compare observed median abundances with a chemical model for low-mass protostars and find fair agreement, although some modeling work remains to bring abundances higher with respect to CH$_3$OH. Median abundances with respect to CH$_3$OH in LYSOs are also found to be generally comparable to observed abundances in hot cores, hot corinos, and massive young stellar objects. Compared with comets, our sample is comparable for all molecules except HC$_3$N and CH$_2$CO, which likely become depleted at later evolutionary stages.

Methyl isocyanate (CH$_{3}$NCO) belongs to a select group of interstellar molecules considered to be relevant precursors in the formation of larger organic compounds, including those with peptide bonds. The molecule has only been detected in a couple of high-mass protostars and potentially on comets. A formation route on icy grains has been postulated for this molecule but experimental evidence is lacking. Here we extend the range of environments where methyl isocyanate is found, and unambiguously identify CH$_{3}$NCO through the detection of 43 unblended transitions in the ALMA Protostellar Interferometric Line Survey (PILS) of the low mass solar-type protostellar binary IRAS 16293-2422. The molecule is detected toward both components of the binary with a ratio HNCO/CH$_3$NCO $\sim$4--12. The isomers CH$_{3}$CNO and CH$_3$OCN are not identified, resulting in upper abundance ratios of CH$_{3}$NCO/CH$_{3}$CNO > 100 and CH$_{3}$NCO/CH$_3$OCN > 10. The resulting abundance ratios compare well with those found for related N-containing species toward high-mass protostars. To constrain its formation, a set of cryogenic UHV experiments is performed. VUV irradiation of CH$_{4}$:HNCO mixtures at 20 K strongly indicate that methyl isocyanate can be formed in the solid-state through CH$_{3}$ and (H)NCO recombinations. Combined with gas-grain models that include this reaction, the solid-state route is found to be a plausible scenario to explain the methyl isocyanate abundances found in IRAS 16293-2422.

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.

The SgrB2 molecular cloud contains several sites forming high-mass stars. SgrB2(N) is one of its main centers of activity. It hosts several compact and UCHII regions, as well as two known hot molecular cores (SgrB2(N1) and SgrB2(N2)), where complex organic molecules are detected. Our goal is to use the high sensitivity of ALMA to characterize the hot core population in SgrB2(N) and shed a new light on the star formation process. We use a complete 3 mm spectral line survey conducted with ALMA to search for faint hot cores in SgrB2(N). We report the discovery of three new hot cores that we call SgrB2(N3), SgrB2(N4), and SgrB2(N5). The three sources are associated with class II methanol masers, well known tracers of high-mass star formation, and SgrB2(N5) also with a UCHII region. The chemical composition of the sources and the column densities are derived by modelling the whole spectra under the assumption of LTE. The H2 column densities are computed from ALMA and SMA continuum emission maps. The H2 column densities of these new hot cores are found to be 16 up to 36 times lower than the one of the main hot core Sgr B2(N1). Their spectra have spectral line densities of 11 up to 31 emission lines per GHz, assigned to 22-25 molecules. We derive rotational temperatures around 140-180 K for the three new hot cores and mean source sizes of 0.4 for SgrB2(N3) and 1.0 for SgrB2(N4) and SgrB2(N5). SgrB2(N3) and SgrB2(N5) show high velocity wing emission in typical outflow tracers, with a bipolar morphology in their integrated intensity maps suggesting the presence of an outflow, like in SgrB2(N1). The associations of the hot cores with class II methanol masers, outflows, and/or UCHII regions tentatively suggest the following age sequence: SgrB2(N4), SgrB2(N3), SgrB2(N5), SgrB2(N1). The status of SgrB2(N2) is unclear. It may contain two distinct sources, a UCHII region and a very young hot core.

Using mm-wavelength data from ALMA, the EMoCA spectral line survey revealed the presence of both the straight-chain and branched forms of propyl cyanide (C$_3$H$_7$CN) toward the Galactic Center star-forming source Sgr B2(N2). This was the first interstellar detection of a branched aliphatic molecule. Through computational methods, we seek to explain the observed $i$:$n$ ratio for propyl cyanide, and to predict the abundances of the four different forms of the homologous nitrile, butyl cyanide (C$_4$H$_9$CN). We also investigate whether other molecules will show a similar degree of branching, by modeling alkanes up to pentane (C$_5$H$_{12}$). We use a coupled three-phase chemical kinetics model to simulate the chemistry of Sgr B2(N2), using an updated chemical network that includes grain-surface/ice-mantle formation routes for branched nitriles and alkanes. We use the EMoCA survey data to search for the straight-chain form of butyl cyanide toward Sgr B2(N2). The observed $i$:$n$ ratio for propyl cyanide is reproduced by the models. Butyl cyanide is predicted to show similar abundances to propyl cyanide, and to exhibit strong branching, with the $sec$ form clearly dominant over all others. The addition of CN to acetylene and ethene is found to be important to the production of vinyl, ethyl, propyl, and butyl cyanide. We report a non-detection of $n$-C$_4$H$_9$CN toward Sgr B2(N2), with an abundance at least 1.7 times lower than that of $n$-C$_3$H$_7$CN. This value is within the range predicted by the chemical models. The models indicate that the degree of branching rises with increasing molecular size. The efficiency of CN addition to unsaturated hydrocarbons boosts the abundances of nitriles in the model, and enhances the ratio of straight-to-branched molecule production. The predicted abundance of $s$-C$_4$H$_9$CN makes it a good candidate for future detection toward Sgr B2(N2).

N-methylformamide, CH3NHCHO, may be an important molecule for interstellar pre-biotic chemistry because it contains a peptide bond. The rotational spectrum of the most stable trans conformer of CH3NHCHO is complicated by strong torsion-rotation interaction due to the low barrier of the methyl torsion. We use two absorption spectrometers in Kharkiv and Lille to measure the rotational spectra over 45--630 GHz. The analysis is carried out using the Rho-axis method and the RAM36 code. We search for N-methylformamide toward the hot molecular core Sgr B2(N2) using a spectral line survey carried out with ALMA. The astronomical results are put into a broader astrochemical context with the help of a gas-grain chemical kinetics model. The laboratory data set for the trans conformer of CH3NHCHO consists of 9469 line frequencies with J <= 62, including the first assignment of the rotational spectra of the first and second excited torsional states. All these lines are fitted within experimental accuracy. We report the tentative detection of CH3NHCHO towards Sgr B2(N2). We find CH3NHCHO to be more than one order of magnitude less abundant than NH2CHO, a factor of two less abundant than CH3NCO, but only slightly less abundant than CH3CONH2. The chemical models indicate that the efficient formation of HNCO via NH + CO on grains is a necessary step in the achievement of the observed gas-phase abundance of CH3NCO. Production of CH3NHCHO may plausibly occur on grains either through the direct addition of functional-group radicals or through the hydrogenation of CH3NCO. Provided the detection of CH3NHCHO is confirmed, the only slight underabundance of this molecule compared to its more stable structural isomer acetamide and the sensitivity of the model abundances to the chemical kinetics parameters suggest that the formation of these two molecules is controlled by kinetics rather than thermal equilibrium.

One of the open questions in astrochemistry is how complex organic and prebiotic molecules are formed. Aims. Our aim is to start the process of compiling an inventory of oxygen-bearing complex organic molecules toward the solar-type Class 0 protostellar binary IRAS16293-2422 from an unbiased spectral survey with ALMA (PILS). Here we focus on the new detections of ethylene oxide (c-C$_2$H$_4$O), acetone (CH$_3$COCH$_3$), and propanal (C$_2$H$_5$CHO). Methods. With ALMA, we surveyed the spectral range from 329 to 363 GHz at 0.5$"$ (60 AU diameter) resolution. Using a simple model for the molecular emission in LTE, the excitation temperatures and column densities of each species were constrained. Results. We successfully detect propanal (44 lines), ethylene oxide (20 lines) and acetone (186 lines) toward one component of the protostellar binary, IRAS16293B. The high resolution maps demonstrate that the emission for all investigated species originates from the compact central region close to the protostar. This, along with a derived common excitation temperature of $\sim$ 125 K, is consistent with a coexistence of these molecules in the same gas. Conclusions. The observations mark the first detections of acetone, propanal and ethylene oxide toward a low-mass protostar. The relative abundance ratios of the two sets of isomers (CH$_3$COCH$_3$/C$_2$H$_5$CHO $\sim$ 8 and CH$_3$CHO/c-C$_2$H$_4$O $\sim$ 12) are comparable to previous observations toward high-mass protostars. The majority of observed abundance ratios from these results as well as those measured toward high-mass protostars are up to an order of magnitude above the predictions from chemical models. This may reflect either missing reactions or uncertain rates in the chemical networks. The physical conditions, such as temperatures or densities, used in the models, may not be applicable to solar-type protostars either.

We performed a spectral line survey called Exploring Molecular Complexity with ALMA (EMoCA) toward Sagittarius~B2(N) between 84.1 and 114.4 GHz with the Atacama Large Millimeter/submillimeter Array (ALMA) in its Cycles 0 and 1. We determined line intensities of n-propyl cyanide in the ground vibrational states of its gauche and anti conformers toward the hot molecular core Sagittarius B2(N2) which suggest that we should also be able to detect transitions pertaining to excited vibrational states. We wanted to determine spectroscopic parameters of low-lying vibrational states of both conformers of n-propyl cyanide to search for them in our ALMA data. We recorded laboratory rotational spectra of n-propyl cyanide in two spectral windows between 36 and 127 GHz. We searched for emission lines produced by these states in the ALMA spectrum of Sagittarius B2(N2). We modeled their emission and the emission of the ground vibrational states assuming local thermodynamic equilibrium (LTE). We have made extensive assignments of a- and b-type transitions of the four lowest vibrational states of the gauche conformer which reach J and Ka quantum numbers of 65 and 20, respectively. We assigned mostly a-type transitions for the anti conformer with J and Ka quantum numbers up to 48 and 24, respectively. Rotational and Fermi perturbations between two anti states allowed us to determine their energy difference. The resulting spectroscopic parameters enabled us to identify transitions of all four vibrational states of each conformer in our ALMA data. The emission features of all states, including the ground vibrational state, are well-reproduced with the same LTE modeling parameters, which gives us confidence in the reliability of the identifications, even for the states with only one clearly detected line.

The inner regions of the envelopes surrounding young protostars are characterised by a complex chemistry, with prebiotic molecules present on the scales where protoplanetary disks eventually may form. This paper introduces a systematic survey, "Protostellar Interferometric Line Survey (PILS)" of the Class 0 protostellar binary IRAS 16293-2422 using the Atacama Large Millimeter/submillimeter Array (ALMA). The survey covers the full frequency range from 329 to 363 GHz (0.8 mm) with additional targeted observations at 3.0 and 1.3 mm. More than 10,000 features are detected toward one component in the protostellar binary. Glycolaldehyde, its isomers, methyl formate and acetic acid, and its reduced alcohol, ethylene glycol, are clearly detected. For ethylene glycol both lowest state conformers, aGg' and gGg', are detected, the latter for the first time in the ISM. The abundance of glycolaldehyde is comparable to or slightly larger than that of ethylene glycol. In comparison to the Galactic Center, these two species are over-abundant relative to methanol, possibly an indication of formation at low temperatures in CO-rich ices. Both 13C and deuterated isotopologues of glycolaldehyde are detected, also for the first time ever in the ISM. For the deuterated species, a D/H ratio of approximately 5% is found with no differences between the deuteration in the different functional groups of glycolaldehyde. Measurements of the 13C-species lead to a 12C:13C ratio of approximately 30, lower than the typical ISM value. This low ratio may reflect an enhancement of 13CO in the ice due to either ion-molecule reactions in the gas before freeze-out or differences in the temperatures where 12CO and 13CO ices sublimate. The results reinforce the importance of low-temperature grain surface chemistry for the formation of prebiotic molecules seen here in the gas after sublimation of the entire ice mantle.

Formamide (NH$_2$CHO) has previously been detected in several star-forming regions and is thought to be a precursor for different prebiotic molecules. Its formation mechanism is still debated, however. Observations of formamide, related species and their isopotologues may provide useful clues to the chemical pathways leading to their formation. The Protostellar Interferometric Line Survey (PILS) represents an unbiased high angular resolution and sensitivity spectral survey of the low-mass protostellar binary IRAS 16293-2422 with the Atacama Large Millimeter/submillimeter Array (ALMA). We detect for the first time the three singly deuterated forms of NH$_2$CHO (NH$_2$CDO, cis- and trans-NHDCHO) as well as DNCO towards the component B of this binary source. The images reveal that the different isotopologues are all present in the same region. Based on observations of the $^{13}$C isotopologues of formamide and a standard $^{12}$C/$^{13}$C ratio, the deuterium fractionation is found to be similar for the three different forms with a value of about 2%. The DNCO/HNCO ratio is also comparable to the D/H ratio of formamide ($\sim$1%). These results are in agreement with the hypothesis that NH$_2$CHO and HNCO are chemically related through grain surface formation.

We have performed a spectral line survey called EMoCA toward Sagittarius B2(N) between 84 and 114 GHz with ALMA. Line intensities of the main isotopic species of ethyl cyanide and its singly 13C-substituted isotopomers observed toward the hot molecular core Sgr B2(N2) suggest that the doubly 13C-substituted isotopomers should be detectable also. We want to determine the spectroscopic parameters of all three doubly 13C-substituted isotopologues of ethyl cyanide to search for them in our ALMA data. We investigated the laboratory rotational spectra of the three species between 150 and 990 GHz. We searched for emission lines produced by these species in the ALMA spectrum of Sgr B2(N2). We modeled their emission as well as the emission of the 12C and singly 13C-substituted isotopologues assuming local thermodynamic equilibrium. We identified more than 5000 rotational transitions, pertaining to more than 3500 different transition frequencies, in the laboratory for each of the three isotopomers. The quantum numbers reach J ~ 115 and K_a ~ 35, resulting in accurate spectroscopic parameters and accurate rest frequency calculations beyond 1000 GHz for strong to moderately weak transitions of either isotopomer. All three species are unambiguously detected in our ALMA data. The 12C/13C column density ratio of the isotopomers with one 13C atom to the ones with two 13C atoms is about 25. Ethyl cyanide is the second molecule after methyl cyanide for which isotopologues containing two 13C atoms have been securely detected in the interstellar medium. The model of our ethyl cyanide data suggests that we should be able to detect vibrational satellites of the main species up to at least v_19 = 1 at 1130 K and up to v_13 + v_21 = 2 at 600 K for the isotopologues with one 13C atom in our present ALMA data. Such satellites may be too weak to be identified unambiguously for isotopologues with two 13C atoms.

Computational models of interstellar gas-grain chemistry have historically adopted a single dust-grain size of 0.1 micron, assumed to be representative of the size distribution present in the interstellar medium. Here, we investigate the effects of a broad grain-size distribution on the chemistry on dust-grain surfaces and the subsequent build-up of molecular ices on the grains, using a three-phase gas-grain chemical model of a quiescent dark cloud. We include an explicit treatment of the grain temperatures, governed both by the visual extinction of the cloud and the size of each individual grain-size population. We find that the temperature difference plays a significant role in determining the total bulk ice composition across the grain-size distribution, while the effects of geometrical differences between size populations appear marginal. We also consider collapse from a diffuse to a dark cloud, allowing dust temperatures to fall. Under the initial diffuse conditions, small grains are too warm to promote grain-mantle build-up, with most ices forming on the mid-sized grains. As collapse proceeds, the more abundant, smallest grains cool and become the dominant ice carriers; the large population of small grains means that this ice is distributed across many grains, with perhaps no more than 40 monolayers of ice each (versus several hundred assuming a single grain size). This effect may be important for the subsequent processing and desorption of the ice during the hot-core phase of star-formation, exposing a significant proportion of the ice to the gas phase, increasing the importance of ice-surface chemistry and surface-gas interactions.

Over the past five decades, radio astronomy has shown that molecular complexity is a natural outcome of interstellar chemistry, in particular in star forming regions. However, the pathways that lead to the formation of complex molecules are not completely understood and the depth of chemical complexity has not been entirely revealed. In addition, the sulfur chemistry in the dense interstellar medium is not well understood. We want to know the relative abundances of alkanethiols and alkanols in the Galactic Center source Sagittarius B2(N2), the northern hot molecular core in Sgr B2(N), whose relatively small line widths are favorable for studying the molecular complexity in space. We investigated spectroscopic parameter sets that were able to reproduce published laboratory rotational spectra of ethanethiol and studied effects that modify intensities in the predicted rotational spectrum of ethanol. We used the Atacama Large Millimeter Array (ALMA) in its Cycles~0 and 1 for a spectral line survey of Sagittarius B2(N) between 84 and 114.4 GHz. These data were analyzed by assuming local thermodynamic equilibrium (LTE) for each molecule. Our observations are supplemented by astrochemical modeling; a new network is used for the first time that includes reaction pathways for alkanethiols. The column density ratios involving methanol, ethanol, and methanethiol in Sgr B2(N2) are similar to values reported for Orion KL, but those involving ethanethiol are significantly different and suggest that the detection of ethanethiol reported toward Orion KL is uncertain. Our chemical model presently does not permit the prediction of sufficiently accurate column densities of alkanethiols or their ratios among alkanethiols and alkanols. Therefore, additional observational results are required to establish the level of C2H5SH in the dense and warm interstellar medium with certainty.

Deuteration is a powerful tracer of the history of the cold prestellar phase in star forming regions. Apart from methanol, little is known about deuterium fractionation of complex organic molecules in the interstellar medium, especially in high mass star forming regions. We use a spectral line survey performed with ALMA to search for deuterated complex organic molecules toward the hot molecular core Sgr B2(N2). Population diagrams and integrated intensity maps are constructed to fit rotational temperatures and emission sizes for each molecule. Column densities are derived by modelling the full spectrum under the LTE assumption. The results are compared to predictions of two astrochemical models that treat the deuteration process. We report the detection of CH2DCN toward Sgr B2(N2) with a deuteration level of 0.4%, and tentative detections of CH2DOH, CH2DCH2CN, the chiral molecule CH3CHDCN, and DC3N with levels in the range 0.05%-0.12%. A stringent deuteration upper limit is obtained for CH3OD (<0.07%). Upper limits in the range 0.5-1.8% are derived for the three deuterated isotopologues of vinyl cyanide, the four deuterated species of ethanol, and CH2DOCHO. Ethyl cyanide is less deuterated than methyl cyanide by at least a factor five. Except for methyl cyanide, the measured deuteration levels lie at least a factor four below the predictions of current astrochemical models. The deuteration levels in Sgr B2(N2) are also lower than in Orion KL by a factor of a few up to a factor ten. The discrepancy between the deuteration levels of Sgr B2(N2) and the predictions of chemical models, and the difference between Sgr B2(N2) and Orion KL may both be due to the higher kinetic temperatures that characterize the Galactic Center region compared to nearby clouds. Alternatively, they may result from a lower overall abundance of deuterium itself in the Galactic Center region by up to a factor ten.

The quest to detect prebiotic molecules in space, notably amino acids, requires an understanding of the chemistry involving nitrogen atoms. Hydroxylamine (NH$_2$OH) is considered a precursor to the amino acid glycine. Although not yet detected, NH$_2$OH is considered a likely target of detection with ALMA. We report on an experimental investigation of the formation of hydroxylamine on an amorphous silicate surface via the oxidation of ammonia. The experimental data are then fed into a simulation of the formation of NH$_2$OH in dense cloud conditions. On ices at 14 K and with a modest activation energy barrier, NH$_2$OH is found to be formed with an abundance that never falls below a factor 10 with respect to NH$_3$. Suggestions of conditions for future observations are provided.

Massive young stellar objects (MYSOs) with hot cores are classic sources of complex organic molecules. The origins of these molecules in such sources, as well as the small- and large-scale differentiation between nitrogen- and oxygen-bearing complex species, are poorly understood. We aim to use complex molecule abundances toward a chemically less explored class of MYSOs with weak hot organic emission lines to constrain the impact of hot molecular cores and initial ice conditions on the chemical composition toward MYSOs. We use the IRAM 30m and the Submillimeter Array to search for complex organic molecules over 8-16 GHz in the 1~mm atmospheric window toward three MYSOs with known ice abundances, but without luminous molecular hot cores. Complex molecules are detected toward all three sources at comparable abundances with respect to CH$_3$OH to classical hot core sources. The relative importance of CH$_3$CHO, CH$_3$CCH, CH$_3$OCH$_3$, CH$_3$CN, and HNCO differ between the organic-poor MYSOs and hot cores, however. Furthermore, the N-bearing molecules are generally concentrated toward the source centers, while most O- and C-bearing molecules are present both in the center and in the colder envelope. Gas-phase HNCO/CH$_3$OH ratios are tentatively correlated with the ratios of NH$_3$ ice over CH$_3$OH ice in the same lines of sight, which is consistent with new gas-grain model predictions. Hot cores are not required to form complex organic molecules, and source temperature and initial ice composition both seem to affect complex organic distributions toward MYSOs. To quantify the relative impact of temperature and initial conditions requires, however, a larger spatially resolved survey of MYSOs with ice detections.

The Chamaeleon clouds are excellent targets for low-mass star formation studies. Cha I and II are actively forming stars while Cha III shows no sign of ongoing star formation. We aim to determine the driving factors that have led to the very different levels of star formation activity in Cha I and III and examine the dynamical state and possible evolution of the starless cores within them. Observations were performed in various molecular transitions with APEX and Mopra. Five cores are gravitationally bound in Cha I and one in Cha III. The infall signature is seen toward 8-17 cores in Cha I and 2-5 cores in Cha III, which leads to a range of 13-28% of the cores in Cha I and 10-25% of the cores in Cha III that are contracting and may become prestellar. Future dynamical interactions between the cores will not be dynamically significant in either Cha I or III, but the subregion Cha I North may experience collisions between cores within ~0.7 Myr. Turbulence dissipation in the cores of both clouds is seen in the high-density tracers N2H+ 1-0 and HC3N 10-9. Evidence of depletion in the Cha I core interiors is seen in the abundance distributions of C17O, C18O, and C34S. Both contraction and static chemical models indicate that the HC3N to N2H+ abundance ratio is a good evolutionary indicator in the prestellar phase for both gravitationally bound and unbound cores. In the framework of these models, we find that the cores in Cha III and the southern part of Cha I are in a similar evolutionary stage and are less chemically evolved than the central region of Cha I. The measured HC3N/N2H+ abundance ratio and the evidence for contraction motions seen towards the Cha III starless cores suggest that Cha III is younger than Cha I Centre and that some of its cores may form stars in the future. The cores in Cha I South may on the other hand be transient structures. (abridged)

The largest non-cyclic molecules detected in the interstellar medium (ISM) are organic with a straight-chain carbon backbone. We report an interstellar detection of a branched alkyl molecule, iso-propyl cyanide (i-C3H7CN), with an abundance 0.4 times that of its straight-chain structural isomer. This detection suggests that branched carbon-chain molecules may be generally abundant in the ISM. Our astrochemical model indicates that both isomers are produced within or upon dust grain ice mantles through the addition of molecular radicals, albeit via differing reaction pathways. The production of iso-propyl cyanide appears to require the addition of a functional group to a non-terminal carbon in the chain. Its detection therefore bodes well for the presence in the ISM of amino acids, for which such side-chain structure is a key characteristic.

The first off-lattice Monte Carlo kinetics model of interstellar dust-grain surface chemistry is presented. The positions of all surface particles are determined explicitly, according to the local potential minima resulting from the pair-wise interactions of contiguous atoms and molecules, rather than by a pre-defined lattice structure. The model is capable of simulating chemical kinetics on any arbitrary dust-grain morphology, as determined by the user-defined positions of each individual dust-grain atom. A simple method is devised for the determination of the most likely diffusion pathways and their associated energy barriers for surface species. The model is applied to a small, idealized dust grain, adopting various gas densities and using a small chemical network. Hydrogen and oxygen atoms accrete onto the grain, to produce H2O, H2, O2 and H2O2. The off-lattice method allows the ice structure to evolve freely; ice mantle porosity is found to be dependent on the gas density, which controls the accretion rate. A gas density of 2 x 10^{4} cm^{-3}, appropriate to dark interstellar clouds, is found to produce a fairly smooth and non-porous ice mantle. At all densities, H2 molecules formed on the grains collect within the crevices that divide nodules of ice, and within micropores (whose extreme inward curvature produces strong local potential minima). The larger pores produced in the high-density models are not typically filled with H2. Direct deposition of water molecules onto the grain indicates that amorphous ices formed in this way may be significantly more porous than interstellar ices that are formed by surface chemistry.

Based on astrochemical theory, the complex molecular composition around high-mass YSOs should evolve from the outer envelope in toward the central hot region as a sequence of temperature dependent chemical pathways are activated in ices and in the gas-phase. The resulting complex molecules have been broadly classified into three generations dependent on the temperature (<25, >25, and >100 K) required for formation. We combine IRAM 30m and Submillimeter Array observations to explore the spatial distribution of organic molecules around the high-mass young stellar object NGC 7538 IRS9, whose weak complex molecule emission previously escaped detection, quantifying the emission and abundance profiles of key organic molecules as a function of distance from the central protostar. We find that emission from N-bearing organics and saturated O-bearing organics present large increases in emission around 8000 AU and R<3000 AU, while O-bearing molecules and hydrocarbons do not. The increase in flux from some complex molecules in the envelope, around 8000 AU or 25 K, is consistent with recent model predictions of an onset of complex ice chemistry at 20-30 K. The emission increase for some molecules at R<3000 AU suggests the presence of a weak hot core, where thermal ice evaporation and hot gas-phase chemistry drives the chemistry. Complex organics thus form at all radii and temperatures around this protostar, but the composition changes dramatically as the temperature increases, which is used to constrain the chemical generation(s) to which different classes of molecule belong.

A new chemical model is presented that simulates fully-coupled gas-phase, grain-surface and bulk-ice chemistry in hot cores. Glycine (NH2CH2COOH), the simplest amino acid, and related molecules such as glycinal, propionic acid and propanal, are included in the chemical network. Glycine is found to form in moderate abundance within and upon dust-grain ices via three radical-addition mechanisms, with no single mechanism strongly dominant. Glycine production in the ice occurs over temperatures ~40-120 K. Peak gas-phase glycine fractional abundances lie in the range 8 x 10^{-11} - 8 x 10^{-9}, occuring at ~200 K, the evaporation temperature of glycine. A gas-phase mechanism for glycine production is tested and found insignificant, even under optimal conditions. A new spectroscopic radiative-transfer model is used, allowing the translation and comparison of the chemical-model results with observations of specific sources. Comparison with the nearby hot-core source NGC 6334 IRS1 shows excellent agreement with integrated line intensities of observed species, including methyl formate. The results for glycine are consistent with the current lack of a detection of this molecule toward other sources; the high evaporation temperature of glycine renders the emission region extremely compact. Glycine detection with ALMA is predicted to be highly plausible, for bright, nearby sources with narrow emission lines. Photodissociation of water and subsequent hydrogen-abstraction from organic molecules by OH, and NH2, are crucial to the build-up of complex organic species in the ice. The inclusion of alternative branches within the network of radical-addition reactions appears important to the abundances of hot-core molecules; less favorable branching ratios may remedy the anomalously high abundance of glycolaldehyde predicted by this and previous models.

We investigate the molecular evolution and D/H abundance ratios that develop as star formation proceeds from a dense-cloud core to a protostellar core, by solving a gas-grain reaction network applied to a 1-D radiative hydrodynamic model with infalling fluid parcels. Spatial distributions of gas and ice-mantle species are calculated at the first-core stage, and at times after the birth of a protostar. Gas-phase methanol and methane are more abundant than CO at radii $r\lesssim 100$ AU in the first-core stage, but gradually decrease with time, while abundances of larger organic species increase. The warm-up phase, when complex organic molecules are efficiently formed, is longer-lived for those fluid parcels in-falling at later stages. The formation of unsaturated carbon chains (warm carbon-chain chemistry) is also more effective in later stages; C$^+$, which reacts with CH$_4$ to form carbon chains, increases in abundance as the envelope density decreases. The large organic molecules and carbon chains are strongly deuterated, mainly due to high D/H ratios in the parent molecules, determined in the cold phase. We also extend our model to simulate simply the chemistry in circumstellar disks, by suspending the 1-D infall of a fluid parcel at constant disk radii. The species CH$_3$OCH$_3$ and HCOOCH$_3$ increase in abundance in $10^4-10^5$ yr at the fixed warm temperature; both also have high D/H ratios.

High-angular resolution observations of dense molecular cores show that these cores can be clumpier at smaller scales, and that some of these clumps can also be unbound or transient. The use of chemical models of the evolution of the molecular gas provides a way to probe the physical properties of the clouds. We study the properties of the clump and inter-clump medium in the starless CS core in LDN 673 by carrying out a molecular line survey with the IRAM 30-m telescope toward two clumps and two inter-clump positions. We also observed the 1.2-mm continuum with the MAMBO-II bolometer at IRAM. The dust continuum map shows four condensations, three of them centrally peaked, coinciding with previously identified sub-millimetre sources. We confirm that the denser clump of the region, $n\sim3.6 \times10^5$\cmt, is also the more chemically evolved, and it could still undergo further fragmentation. The inter-clump medium positions are denser than previously expected, likely $n\sim1\times10^3$--1$\times10^4$\cmt\ due to contamination, and are chemically young, similar to the gas in the lower density clump position. We argue that the density contrast between these positions and their general young chemical age would support the existence of transient clumps in the lower density material of the core. We were also able to find reasonable fits of the observationally derived chemical abundances to models of the chemistry of transient clumps.

We investigate the formation and evolution of interstellar dust-grain ices under dark-cloud conditions, with a particular emphasis on CO2. We use a three-phase model (gas/surface/mantle) to simulate the coupled gas--grain chemistry, allowing the distinction of the chemically-active surface from the ice layers preserved in the mantle beneath. The model includes a treatment of the competition between barrier-mediated surface reactions and thermal-hopping processes. The results show excellent agreement with the observed behavior of CO2, CO and water ice in the interstellar medium. The reaction of the OH radical with CO is found to be efficient enough to account for CO2 ice production in dark clouds. At low visual extinctions, with dust temperatures ~12 K, CO2 is formed by direct diffusion and reaction of CO with OH; we associate the resultant CO2-rich ice with the observational polar CO2 signature. CH4 ice is well correlated with this component. At higher extinctions, with lower dust temperatures, CO is relatively immobile and thus abundant; however, the reaction of H and O atop a CO molecule allows OH and CO to meet rapidly enough to produce a CO:CO2 ratio in the range ~2--4, which we associate with apolar signatures. We suggest that the observational apolar CO2/CO ice signatures in dark clouds result from a strongly segregated CO:H2O ice, in which CO2 resides almost exclusively within the CO component. Observed visual-extinction thresholds for CO2, CO and H2O are well reproduced by depth-dependent models. Methanol formation is found to be strongly sensitive to dynamical timescales and dust temperatures.

(Abridged) Gas-phase complex organic molecules are commonly detected in the warm inner regions of protostellar envelopes. Recent models show that photochemistry in ices followed by desorption may explain the observed abundances. This study aims to experimentally quantify the broad-band UV-induced production rates of complex organics in CH3OH-rich ices at 20-70 K under ultra-high vacuum conditions. The reaction products are mainly identified by RAIRS and TPD experiments. Complex organics are readily formed in all experiments, both during irradiation and during a slow warm-up of the ices to 200 K after the UV lamp is turned off. The relative abundances of photoproducts depend on the UV fluence, the ice temperature, and whether pure CH3OH ice or CH3OH:CH4/CO ice mixtures are used. C2H6, CH3CHO, CH3CH2OH, CH3OCH3, HCOOCH3, HOCH2CHO and (CH2OH)2 are all detected in at least one experiment. The derived product-formation yields and their dependences on different experimental parameters, such as the initial ice composition, are used to estimate the CH3OH photodissociation branching ratios in ice and the relative diffusion barriers of the formed radicals. The experiments show that ice photochemistry in CH3OH ices is efficient enough to explain the observed abundances of complex organics around protostars and that ratios of complex molecules can be used to constrain their formation pathway.

Sensitive and high angular resolution ($\sim$ 0.4\arcsec) SO$_2$[22$_{2,20}$ $\to$ 22$_{1,21}$] and SiO[5$\to$4] line and 1.3 and 7 mm continuum observations made with the Submillimeter Array (SMA) and the Very Large Array (VLA) towards the young massive cluster W51 IRS2 are presented. We report the presence of a large (of about 3000 AU) and massive (40 M$_\odot$) dusty circumstellar disk and a hot gas molecular ring around a high-mass protostar or a compact small stellar system associated with W51 North. The simultaneous observations of the silicon monoxide molecule, an outflow gas tracer, further revealed a massive (200 M$_\odot$) and collimated ($\sim14^\circ$) outflow nearly perpendicular to the dusty and molecular structures suggesting thus the presence of a single very massive protostar with a bolometric luminosity of more than 10$^5$ L$_\odot$. A molecular hybrid LTE model of a Keplerian and infalling ring with an inner cavity and a central stellar mass of more than 60 M$_\odot$ agrees well with the SO$_2$[22$_{2,20}$ $\to$ 22$_{1,21}$] line observations. Finally, these results suggest that mechanisms, such as mergers of low- and intermediate- mass stars, might be not necessary for forming very massive stars.

In recent years, organic molecules of increasing complexity have been found toward the prolific Galactic center source Sagittarius B2. We wish to explore the degree of complexity that the interstellar chemistry can reach in star-forming regions. We carried out a complete line survey of the hot cores Sgr B2(N) and (M) with the IRAM 30 m telescope in the 3 mm range. We analyzed this spectral survey in the LTE approximation. We modeled the emission of all known molecules simultaneously, which allows us to search for less abundant, more complex molecules. We compared the derived column densities with the predictions of a coupled gas-phase and grain-surface chemical code. We report the first detection in space of ethyl formate (C2H5OCHO) and n-propyl cyanide (C3H7CN) toward Sgr B2(N). The abundances of ethyl formate and n-propyl cyanide relative to H2 are estimated to be 3.6e-9 and 1.0e-9, respectively. Our chemical modeling suggests that the sequential, piecewise construction of ethyl and n-propyl cyanide from their constituent functional groups on the grain surfaces is their most likely formation route. Ethyl formate is primarily formed on the grains by adding CH3 to functional-group radicals derived from methyl formate, although ethanol may also be a precursor. The detection in Sgr B2(N) of the next stage of complexity in two classes of complex molecule, esters and alkyl cyanides, suggests that greater complexity in other classes of molecule may be present in the interstellar medium. {Abridged}

Understanding grain-surface processes is crucial to interpreting the chemistry of the ISM. However, accurate surface chemistry models are computationally expensive and are difficult to integrate with gas-phase simulations. A new modified-rate method for solving grain-surface chemical systems is presented. Its purpose is accurately to model highly complex systems that can otherwise only be treated using the sometimes inadequate rate-equation approach. In contrast to previous rate-modification techniques, the functional form of the surface production rates was modified, and not simply the rate coefficient. This form is appropriate to the extreme "small-grain" limit, and can be verified using an analytical master-equation approach. Various further modifications were made to this basic form, to account for competition between processes, to improve estimates of surface occupation probabilities, and to allow a switch-over to the normal rate equations where these are applicable. The new method was tested against systems solved previously using exact techniques. Even the simplest method is quite accurate, and a great improvement over rate equations. Further modifications allow the master-equation results to be reproduced exactly for the methanol-producing system, within computational accuracy. Small discrepancies arise when non-zero activation energies are assumed for the methanol system, which result from complex reaction-competition processes that cannot be resolved easily without using exact methods. Inaccuracies in computed abundances are never greater than a few tens of percent, and typically of the order of one percent, in the most complex systems tested. Implementation of the method in simple networks, including hydrogen-only systems, is trivial, whilst the results are highly accurate.

Sakai et al. have observed long-chain unsaturated hydrocarbons and cyanopolyynes in the low-mass star-forming region L1527, and have attributed this result to a gas-phase ion-molecule chemistry, termed ``Warm Carbon Chain Chemistry'', which occurs during and after the evaporation of methane from warming grains. The source L1527 is an envelope surrounding a Class 0/I protostar with regions that possess a slightly elevated temperature of ~30 K. The molecules detected by Sakai et al. are typically associated only with dark molecular clouds, and not with the more evolved hot corino phase. In order to determine if L1527 is chemically distinct from a dark cloud, we compute models including various degrees of heating. The results indicate that the composition of L1527 is somewhat more likely to be due to ``Warm Carbon Chain Chemistry'' than to be a remnant of a colder phase. If so, the molecular products provide a signature of a previously uncharacterized early phase of low mass star formation, which can be characterized as a ``lukewarm'' corino. We also include predictions for other molecular species that might be observed toward candidate lukewarm corino sources. Although our calculations show that unsaturated hydrocarbons and cyanopolyynes can be produced in the gas phase as the grains warm up to 30 K, they also show that such species do not disappear rapidly from the gas as the temperature reaches 200 K, implying that such species might be detected in hot corinos and hot cores.

Gas-phase processes were long thought to be the key formation mechanisms for complex organic molecules in star-forming regions. However, recent experimental and theoretical evidence has cast doubt on the efficiency of such processes. Grain-surface chemistry is frequently invoked as a solution, but until now there have been no quantitative models taking into account both the high degree of chemical complexity and the evolving physical conditions of star-forming regions. Here, we introduce a new gas-grain chemical network, wherein a wide array of complex species may be formed by reactions involving radicals. The radicals we consider (H, OH, CO, HCO, CH3, CH3O, CH2OH, NH and NH2) are produced primarily by cosmic ray-induced photodissociation of the granular ices formed during the colder, earlier stages of evolution. The gradual warm-up of the hot core is crucial to the formation of complex molecules, allowing the more strongly-bound radicals to become mobile on grain surfaces. This type of chemistry is capable of reproducing the high degree of complexity seen in Sgr B2(N), and can explain the observed abundances and temperatures of a variety of previously detected complex organic molecules, including structural isomers. Many other complex species are predicted by this model, and several of these species may be detectable in hot cores. Differences in the chemistry of high- and low-mass star-formation are also addressed; greater chemical complexity is expected where evolution timescales are longer.

We present high angular resolution ($\sim$ 3$''$) and sensitive 1.3 mm continuum, cyanogen (CN) and vinyl cyanide (C$_2$H$_3$CN) line observations made with the Submillimeter Array (SMA) toward one of most highly obscured objects of the W51 IRS2 region, W51 North. We find that the CN line exhibits a pronounced inverse P-Cygni profile indicating that the molecular gas is infalling inwards this object with a mass accretion rate between 4 and 7 $\times$ 10$^{-2}$ M$_\odot$ yr$^{-1}$. The C$_2$H$_3$CN traces an east-west rotating molecular envelope that surrounds either a single obscured (proto)star with a kinematic mass of 40 M$_{\odot}$ or a small central cluster of B-type stars and that is associated with a compact high velocity bipolar outflow traced by H$_2$O masers and SiO molecular emission. We thus confirm that the W51 North region is part of the growing list of young massive star forming regions that have been associated with infalling motions and with large mass accretion rates ($\sim$ 10$^{-2}$ -- 10$^{-4}$), strengthening the evidence for massive stars forming with very high accretion rates sufficient to quench the formation of an UCHII region.

We investigate molecular evolution in a star-forming core that is initially a hydrostatic starless core and collapses to form a low-mass protostar. The results of a one-dimensional radiation-hydrodynamics calculation are adopted as a physical model of the core. We first derive radii at which CO and large organic species sublimate. CO sublimation in the central region starts shortly before the formation of the first hydrostatic core. When the protostar is born, the CO sublimation radius extends to 100 AU, and the region inside $\lesssim 10$ AU is hotter than 100 K, at which some large organic species evaporate. We calculate the temporal variation of physical parameters in infalling shells, in which the molecular evolution is solved using an updated gas-grain chemical model to derive the spatial distribution of molecules in a protostellar core. The shells pass through the warm region of $10 -100$ K in several $\times$ $10^4$ yr, and fall into the central star $\sim 100$ yr after they enter the region where $T \gtrsim 100$ K. We find that large organic species are formed mainly via grain-surface reactions at temperatures of $20 -40$ K and then desorbed into the gas-phase at their sublimation temperatures. Carbon-chain species can be formed by a combination of gas-phase reactions and grain-surface reactions following the sublimation of CH$_4$. Our model also predicts that CO$_2$ is more abundant in isolated cores, while gas-phase large organic species are more abundant in cores embedded in ambient clouds.

Aims: The gas-phase abundance of methanol in dark quiescent cores in the interstellar medium cannot be explained by gas-phase chemistry. In fact, the only possible synthesis of this species appears to be production on the surfaces of dust grains followed by desorption into the gas. Yet, evaporation is inefficient for heavy molecules such as methanol at the typical temperature of 10 K. It is necessary then to consider non-thermal mechanisms for desorption. But, if such mechanisms are considered for the production of methanol, they must be considered for all surface species. Methods: Our gas-grain network of reactions has been altered by the inclusion of a non-thermal desorption mechanism in which the exothermicity of surface addition reactions is utilized to break the bond between the product species and the surface. Our estimated rate for this process derives from a simple version of classical unimolecular rate theory with a variable parameter only loosely contrained by theoretical work. Results: Our results show that the chemistry of dark clouds is altered slightly at times up to 10^6 yr, mainly by the enhancement in the gas-phase abundances of hydrogen-rich species such as methanol that are formed on grain surfaces. At later times, however, there is a rather strong change. Instead of the continuing accretion of most gas-phase species onto dust particles, a steady-state is reached for both gas-phase and grain-surface species, with significant abundances for the former. Nevertheless, most of the carbon is contained in an undetermined assortment of heavy surface hydrocarbons. Conclusions: The desorption mechanism discussed here will be better constrained by observational data on pre-stellar cores, where a significant accretion of species such as CO has already occurred.

We investigate the chemical and observational implications of repetitive transient dense core formation in molecular clouds. We allow a transient density fluctuation to form and disperse over a period of 1 Myr, tracing its chemical evolution. We then allow the same gas immediately to undergo further such formation and dispersion cycles. The chemistry of the dense gas in subsequent cycles is similar to that of the first, and a limit cycle is reached quickly (2 - 3 cycles). Enhancement of hydrocarbon abundances during a specific period of evolution is the strongest indicator of previous dynamical history. The molecular content of the diffuse background gas in the molecular cloud is expected to be strongly enhanced by the core formation and dispersion process. Such enhancement may remain for as long as 0.5 Myr. The frequency of repetitive core formation should strongly determine the level of background molecular enhancement. We also convolve the emission from a synthesised dark cloud, comprised of ensembles of transient dense cores. We find that the dynamical history of the gas, and therefore the chemical state of the diffuse inter-core medium, may be determined if a sufficient sample of cores is present in an ensemble. Molecular ratios of key hydrocarbons with SO and SO2 are crucial to this distinction. Only surveys with great enough angular resolution to resolve individual cores, or very small groupings, are expected to show evidence of repetitive dynamical processing. The existence of non-equilibrium chemistry in the diffuse background may have implications for the initial conditions used in chemical models. Observed variations in the chemistries of diffuse and translucent regions may be explained by lines of sight which intersect a number of molecular cloud cores in various stages of evolution.

Aims: The production of saturated organic molecules in hot cores and corinos is not well understood. The standard approach is to assume that, as temperatures heat up during star formation, methanol and other species evaporate from grain surfaces and undergo a warm gas-phase chemistry at 100 K or greater to produce species such as methyl formate, dimethyl ether, and others. But a series of laboratory results shows that protonated ions, typical precursors to final products in ion-molecule schemes, tend to fragment upon dissociative recombination with electrons rather than just ejecting a hydrogen atom. Moreover, the specific proposed reaction to produce protonated methyl formate is now known not to occur at all. Methods: We utilize a gas-grain chemical network to probe the chemistry of the relatively ignored stage of hot core evolution during which the protostar switches on and the temperature of the surrounding gas and dust rises from 10 K to over 100 K. During this stage, surface chemistry involving heavy radicals becomes more important as surface hydrogen atoms tend to evaporate rather than react. Results: Our results show that complex species such as methyl formate, formic acid, and dimethyl ether can be produced in large abundance during the protostellar switch-on phase, but that both grain-surface and gas-phase processes help to produce most species. The longer the timescale for protostellar switch-on, the more important the surface processes.

We construct models of molecular clouds that are considered as ensembles of transient cores. Each core is assumed to develop in the background gas of the cloud, grow to high density and decay into the background. The chemistry in each core responds to the dynamical state of the gas and to the gas-dust interaction. Ices are deposited on the dust grains in the core's dense phase, and this material is returned to the gas as the core expands to low density. The cores of the ensemble number typically one thousand and are placed randomly in position within the cloud, and are assigned a random evolutionary phase. The models are used to generate molecular line contour maps of a typical dark cloud. These maps are found to represent extremely well the characteristic features of observed maps of the dark cloud L673, which has been observed at both low and high resolutions. The computed maps are found to exhibit the general morphology of the observed maps, and to generate similar sizes of emitting regions, molecular column densities, and the separations between peaks of emissions of various molecular species. The models give insight into the nature of molecular clouds and the dynamical processes occurring within them, and significantly constrain dynamical and chemical processes in the interstellar medium.