Molecular clouds are the cold regions of the Milky Way where stars form. They are enriched by rather complex molecules. Many of these molecules are believed to be synthesized on the icy surfaces of the interstellar submicron-sized dust grains that permeate the Galaxy. At 10 K thermal desorption is ineffcient and, therefore, why these molecules are found in the cold gas has tantalized astronomers for years. The assumption of the current models, called chemical desorption, is that the molecule formation energy released by the chemical reaction at the grain surface is partially absorbed by the grain and the remaining one causes the ejection of the newly formed molecule into the gas. Here we report an accurate ab-initio molecular dynamics simulations aimed to study the fate of the energy released by the first reaction of the H addition chain on CO, CO + H $\rightarrow$ HCO, occurring on a crystalline ice surface model. We show that about 90% of the HCO formation energy is injected towards the ice in the first picosecond, leaving HCO with an energy content (10-15 kJ mol$^{-1}$) more than a factor two lower than its adsorption energy (30 kJ mol$^{-1}$). As a result, in agreement with laboratory experiments, we conclude that chemical desorption is ineffcient for this specific system, namely H + CO on crystalline ice. We suspect this behavior to be quite general when dealing with hydrogen bonds, which are responsible of both the cohesive energy of the ice mantle and the interaction with adsorbates, as the HCO radical, even though ad hoc simulations are needed to draw specific conclusions on other systems.

The SOLIS (Seeds Of Life In Space) IRAM/NOEMA Large Program aims at studying a set of crucial complex organic molecules in a sample of sources, with well-known physical structure, covering the various phases of Solar-type star formation. One representative object of the transition from the prestellar core to the protostar phases has been observed toward the Very Low Luminosity Object (VeLLO) called L1521F. This type of source is important to study to make the link between prestellar cores and Class 0 sources and also to constrain the chemical evolution during the process of star formation. Two frequency windows (81.6-82.6 GHz and 96.65-97.65 GHz) were used to observe the emission from several complex organics toward the L1521F VeLLO. Only 2 transitions of methanol (A+, E2) have been detected in the narrow window centered at 96.7 GHz (with an upper limit on E1) in a very compact emission blob (~7'' corresponding to ~1000au) toward the NE of the L1521F protostar. The CS 2-1 transition is also detected within the WideX bandwidth. Consistently, with what has been found in prestellar cores, the methanol emission appears ~1000au away from the dust peak. The location of the methanol blob coincides with one of the filaments previously reported in the literature. The Tex of the gas inferred from methanol is (10$\pm$2) K, while the H2 gas density (estimated from the detected CS 2-1 emission and previous CS 5-4 ALMA obs.) is a factor >25 higher than the density in the surrounding environment (n(H2) >10$^{7}$ cm$^{-3}$). From its compactness, low excitation temperature and high gas density, we suggest that the methanol emission detected with NOEMA is either a cold and dense shock-induced blob, recently formed ($\leq$ few hundred years) by infalling gas or a cold and dense fragment that may have just been formed as a result of the intense gas dynamics found within the L1521F VeLLO system.

Formation of interstellar complex organic molecules is currently thought to be dominated by the barrierless coupling between radicals on the interstellar icy grain surfaces. Previous standard DFT results on the reactivity between CH$_3$ and HCO on amorphous water surfaces, showed that formation of CH$_4$ + CO by H transfer from HCO to CH$_3$ assisted by water molecules of the ice was the dominant channel. However, the adopted description of the electronic structure of the biradical (i.e., CH$_3$/HCO) system was inadequate (without the broken-symmetry (BS) approach). In this work, we revisit the original results by means of BS-DFT both in gas phase and with one water molecule simulating the role of the ice. Results indicate that adoption of BS-DFT is mandatory to describe properly biradical systems. In the presence of the single water molecule, the water-assisted H transfer exhibits a high energy barrier. In contrast, CH$_3$CHO formation is found to be barrierless. However, direct H transfer from HCO to CH$_3$ to give CO and CH$_4$ presents a very low energy barrier, hence being a potential competitive channel to the radical coupling and indicating, moreover, that the physical insights ofthe original work remain valid.

Interstellar complex organic molecules (iCOMs) can be loosely defined as chemical compounds with at least six atoms in which at least one is carbon. The observations of iCOMs in star-forming regions have shown that they contain an important fraction of carbon in a molecular form, which can be used to synthesize more complex, even biotic molecules. Hence, iCOMs are major actors in the increasing molecular complexity in space and they might have played a role in the origin of terrestrial life. Understand-ing how iCOMs are formed is relevant for predicting the ultimate organic chemistry reached in the interstellar medium. One possibility is that they are synthesized on the interstellar grain icy surfaces, via recombination of previously formed radicals. The present work focuses on the reactivity of HCO with CH3/NH2 on the grain icy sur-faces, investigated by means of quantum chemical simulations. The goal is to carry outa systematic study using different computational approaches and models for the icy surfaces. Specifically, DFT computations have been bench-marked with CASPT2 and CCSD(T) methods, and the ice mantles have been mimicked with cluster models of 1, 2, 18 and 33 H2O molecules, in which different reaction sites have been considered. Our results indicate that the HCO + CH3/NH2 reactions, if they actually occur, have two major competitive channels: the formation of iCOMs CH3CHO/NH2CHO, or the formation of CO + CH4/NH3. These two channels are either barrierless or presentrelatively low ($\leq$ 10 kJ/mol equal to about 1200 K) energy barriers. Finally, we briefly discuss the astrophysical implications of these findings.

Interstellar formamide (NH2CHO) has recently attracted significant attention due to its potential role as a molecular building block in the formation of precursor biomolecules relevant for the origin of life. Its formation, whether on the surfaces of the interstellar grains or in the gas phase, is currently debated. The present article presents new theoretical quantum chemical computations on possible NH2CHO formation routes in water-rich amorphous ices, simulated by a 33-H2O-molecule cluster. We have considered three possible routes. The first one refers to a scenario used in several current astrochemical models, that is, the radical-radical association reaction between NH2 and HCO. Our calculations show that formamide can indeed be formed, but in competition with formation of NH3 and CO through a direct H transfer process. The final outcome of the NH2 + HCO reactivity depends on the relative orientation of the two radicals on the ice surface. We then analyzed two other possibilities, suggested here for the first time: reaction of either HCN or CN with water molecules of the ice mantle. The reaction with HCN has been found to be characterized by large energy barriers and, therefore, cannot occur under the interstellar ice conditions. On the contrary, the reaction with the CN radical can occur, possibly leading through multiple steps to the formation of NH2CHO. For this reaction, water molecules of the ice act as catalytic active sites since they help the H transfers involved in the process, thus reducing the energy barriers (compared to the gas-phase analogous reaction). Additionally, we apply a statistical model to estimate the reaction rate coefficient when considering the cluster of 33-H2O-molecules as an isolated moiety with respect to the surrounding environment, i.e., the rest of the ice.

Carbon monoxide (CO) is the second most abundant molecule in the gas-phase of the interstellar medium. In dense molecular clouds, it is also present in the solid-phase as a constituent of the mixed water-dominated ices covering dust grains. Its presence in the solid-phase is inferred from its infrared (IR) signals. In experimental observations of solid CO/water mixed samples, its IR frequency splits into two components, giving rise to a blue- and a redshifted band. However, in astronomical observations, the former has never been observed. Several attempts have been carried out to explain this peculiar behaviour, but the question still remains open. In this work, we resorted to pure quantum mechanical simulations in order to shed some light on this problem. We adopted different periodic models simulating the CO/H$_2$O ice system, such as single and multiple CO adsorption on water ice surfaces, CO entrapped into water cages and proper CO:H$_2$O mixed ices. We also simulated pure solid CO. The detailed analysis of our data revealed how the quadrupolar character of CO and the dispersive forces with water ice determine the energetic of the CO/H$_2$O ice interaction, as well as the CO spectroscopic behaviour. Our data suggest that the blueshifted peak can be assigned to CO interacting {\it via} the C atom with dangling H atoms of the water ice, while the redshifted one can actually be the result of CO involved in different reciprocal interactions with the water matrix. We also provide a possible explanation for the lack of the blueshifted peak in astronomical spectra. Our aim is not to provide a full account of the various interstellar ices, but rather to elucidate the sensitivity of the CO spectral features to different water ice environments.

Towards the pre-stellar core L1544, the methanol (CH$_3$OH) emission forms an asymmetric ring around the core centre, where CH$_3$OH is mostly in solid form, with a clear peak 4000~au to the north-east of the dust continuum peak. As part of the NOEMA Large Project SOLIS (Seeds of Life in Space), the CH$_3$OH peak has been spatially resolved to study its kinematics and physical structure and to investigate the cause behind the local enhancement. We find that methanol emission is distributed in a ridge parallel to the main axis of the dense core. The centroid velocity increases by about 0.2~km~s$^{-1}$ and the velocity dispersion increases from subsonic to transonic towards the central zone of the core, where the velocity field also shows complex structure. This could be indication of gentle accretion of material onto the core or interaction of two filaments, producing a slow shock. We measure the rotational temperature and show that methanol is in local thermodynamic equilibrium (LTE) only close to the dust peak, where it is significantly depleted. The CH$_3$OH column density, $N_{tot}({\rm CH_3OH})$, profile has been derived with non-LTE radiative transfer modelling and compared with chemical models of a static core. The measured $N_{tot}({\rm CH_3OH})$ profile is consistent with model predictions, but the total column densities are one order of magnitude lower than those predicted by models, suggesting that the efficiency of reactive desorption or atomic hydrogen tunnelling adopted in the model may be overestimated; or that an evolutionary model is needed to better reproduce methanol abundance.

Complex organic molecules have been observed for decades in the interstellar medium. Some of them might be considered as small bricks of the macromolecules at the base of terrestrial life. It is hence particularly important to understand organic chemistry in Solar-like star forming regions. In this article, we present a new observational project: SOLIS (Seeds Of Life In Space). This is a Large Project at the IRAM-NOEMA interferometer, and its scope is to image the emission of several crucial organic molecules in a sample of Solar-like star forming regions in different evolutionary stage and environments. Here, we report the first SOLIS results, obtained from analysing the spectra of different regions of the Class 0 source NGC1333-IRAS4A, the protocluster OMC-2 FIR4, and the shock site L1157-B1. The different regions were identified based on the images of formamide (NH2CHO) and cyanodiacetylene (HC5N) lines. We discuss the observed large diversity in the molecular and organic content, both on large (3000-10000 au) and relatively small (300-1000 au) scales. Finally, we derive upper limits to the methoxy fractional abundance in the three observed regions of the same order of magnitude of that measured in few cold prestellar objects, namely ~10^-12-10^-11 with respect to H2 molecules.

Context: Modern versions of the Miller-Urey experiment claim that formamide (NH$_2$CHO) could be the starting point for the formation of metabolic and genetic macromolecules. Intriguingly, formamide is indeed observed in regions forming Solar-type stars as well as in external galaxies. Aims: How NH$_2$CHO is formed has been a puzzle for decades: our goal is to contribute to the hotly debated question of whether formamide is mostly formed via gas-phase or grain surface chemistry. Methods: We used the NOEMA interferometer to image NH$_2$CHO towards the L1157-B1 blue-shifted shock, a well known interstellar laboratory, to study how the components of dust mantles and cores released into the gas phase triggers the formation of formamide. Results: We report the first spatially resolved image (size $\sim$ 9", $\sim$ 2300 AU) of formamide emission in a shocked region around a Sun-like protostar: the line profiles are blueshifted and have a FWHM $\simeq$ 5 km s$^{-1}$. A column density of $N_{\rm NH_2CHO}$ = 8 $\times$ 10$^{12}$ cm$^{-1}$, and an abundance (with respect to H-nuclei) of 4 $\times$ 10$^{-9}$ are derived. We show a spatial segregation of formamide with respect to other organic species. Our observations, coupled with a chemical modelling analysis, indicate that the formamide observed in L1157-B1 is formed by gas-phase chemical process, and not on grain surfaces as previously suggested. Conclusions: The SOLIS interferometric observations of formamide provide direct evidence that this potentially crucial brick of life is efficiently formed in the gas-phase around Sun-like protostars.

The interstellar delivery of carbon atoms locked into molecules might be one of the key ingredients for the emergence of life. Cyanopolyynes are carbon chains delimited at their two extremities by an atom of hydrogen and a cyano group, so that they might be excellent reservoirs of carbon. The simplest member, HC3N, is ubiquitous in the galactic interstellar medium and found also in external galaxies. Thus, understanding the growth of cyanopolyynes in regions forming stars similar to our Sun, and what affects it, is particularly relevant. In the framework of the IRAM/NOEMA Large Program SOLIS (Seeds Of Life In Space), we have obtained a map of two cyanopolyynes, HC3N and HC5N, in the protocluster OMC2-FIR4. Because our Sun is thought to be born in a rich cluster, OMC2-FIR4 is one of the closest and best known representatives of the environment in which the Sun may have been born. We find a HC3N/HC5N abundance ratio across the source in the range ~ 1 - 30, with the smallest values (< 10) in FIR5 and in the Eastern region of FIR4. The ratios < 10 can be reproduced by chemical models only if: (1) the cosmic-ray ionisation rate $\zeta$ is ~ $4 \times 10^{-14}$ s$^{-1}$; (2) the gaseous elemental ratio C/O is close to unity; (3) O and C are largely depleted. The large $\zeta$ is comparable to that measured in FIR4 by previous works and was interpreted as due to a flux of energetic (> 10 MeV) particles from embedded sources. We suggest that these sources could lie East of FIR4 and FIR5. A temperature gradient across FIR4, with T decreasing by about 10 K, could also explain the observed change in the HC3N/HC5N line ratio, without the need of a cosmic ray ionisation rate gradient. However, even in this case, a high constant cosmic-ray ionisation rate (of the order of $10^{-14}$ s$^{-1}$) is necessary to reproduce the observations.