S2 has been observed for decades in comets, including comet 67P/Churyumov-Gerasimenko. Despite the fact that this molecule appears ubiquitous in these bodies, the nature of its source remains unknown. In this study, we assume that S2 is formed by irradiation (photolysis and/or radiolysis) of S-bearing molecules embedded in the icy grain precursors of comets, and that the cosmic ray flux simultaneously creates voids in ices within which the produced molecules can accumulate. We investigate the stability of S2 molecules in such cavities, assuming that the surrounding ice is made of H2S or H2O. We show that the stabilization energy of S2 molecules in such voids is close to that of the H2O ice binding energy, implying that they can only leave the icy matrix when this latter sublimates. Because S2 has a short lifetime in the vapor phase, we derive that its formation in grains via irradiation must occur only in low density environments such as the ISM or the upper layers of the protosolar nebula, where the local temperature is extremely low. In the first case, comets would have agglomerated from icy grains that remained pristine when entering the nebula. In the second case, comets would have agglomerated from icy grains condensed in the protosolar nebula and that would have been efficiently irradiated during their turbulent transport towards the upper layers of the disk. Both scenarios are found consistent with the presence of molecular oxygen in comets.

The nature of the icy material accreted by comets during their formation in the outer regions of the protosolar nebula is a major open question in planetary science. Some scenarios of comet formation predict that these bodies agglomerated from crystalline ices condensed in the protosolar nebula. Concurrently, alternative scenarios suggest that comets accreted amorphous ice originating from the interstellar cloud or from the very distant regions of the protosolar nebula. On the basis of existing laboratory and modeling data, we find that the N$_2$/CO and Ar/CO ratios measured in the coma of the Jupiter family comet 67P/Churyumov-Gerasimenko by the ROSINA instrument aboard the European Space Agency's Rosetta spacecraft match those predicted for gases trapped in clathrates. If these measurements are representative of the bulk N$_2$/CO and Ar/CO ratios in 67P/Churyumov-Gerasimenko, it implies that the ices accreted by the comet formed in the nebula and do not originate from the interstellar medium, supporting the idea that the building blocks of outer solar system bodies have been formed from clathrates and possibly from pure crystalline ices. Moreover, because 67P/Churyumov-Gerasimenko is impoverished in Ar and N$_2$, the volatile enrichments observed in Jupiter's atmosphere cannot be explained solely via the accretion of building blocks with similar compositions and require an additional delivery source. A potential source may be the accretion of gas from the nebula that has been progressively enriched in heavy elements due to photoevaporation.

Molecular oxygen has been detected in the coma of comet 67P/Churyumov-Gerasimenko with abundances in the 1-10% range by the ROSINA-DFMS instrument on board the Rosetta spacecraft. Here we find that the radiolysis of icy grains in low-density environments such as the presolar cloud may induce the production of large amounts of molecular oxygen. We also show that molecular oxygen can be efficiently trapped in clathrates formed in the protosolar nebula, and that its incorporation as crystalline ice is highly implausible because this would imply much larger abundances of Ar and N2 than those observed in the coma. Assuming that radiolysis has been the only O2 production mechanism at work, we conclude that the formation of comet 67P/Churyumov-Gerasimenko is possible in a dense and early protosolar nebula in the framework of two extreme scenarios: (1) agglomeration from pristine amorphous icy grains/particles formed in ISM and (2) agglomeration from clathrates that formed during the disk's cooling. The former scenario is found consistent with the strong correlation between O2 and H2O observed in 67P/C-G's coma while the latter scenario requires that clathrates formed from ISM icy grains that crystallized when entering the protosolar nebula.

On November 12, 2014, the ESA/Rosetta descent module Philae landed on the Abydos site of comet 67P/Churyumov-Gerasimenko. Aboard this module, the Ptolemy mass spectrometer measured a CO/CO2 ratio of 0.07 +/- 0.04 which differs substantially from the value obtained in the coma by the Rosetta/ROSINA instrument, suggesting a heterogeneity in the comet nucleus. To understand this difference, we investigated the physico-chemical properties of the Abydos subsurface leading to CO/CO2 ratios close to that observed by Ptolemy at the surface of this region. We used a comet nucleus model that takes into account different water ice phase changes (amorphous ice, crystalline ice and clathrates), as well as diffusion of molecules throughout the pores of the matrix. The input parameters of the model were optimized for the Abydos site and the ROSINA CO/CO2 measured ratio is assumed to correspond to the bulk value in the nucleus. We find that all considered structures of water ice are able to reproduce the Ptolemy observation with a time difference not exceeding ~50 days, i.e. lower than ~2% on 67P/Churyumov-Gerasimenko's orbital period. The suspected heterogeneity of 67P/Churyumov-Gerasimenko's nucleus is also found possible only if it is constituted of crystalline ices. If the icy phase is made of amorphous ice or clathrates, the difference between Ptolemy and ROSINA's measurements would rather originate from the spatial variations in illumination on the nucleus surface. An eventual new measurement of the CO/CO2 ratio at Abydos by Ptolemy could be decisive to distinguish between the three water ice structures.

Bimolecular diffusion coefficients are important parameters used by atmospheric models to calculate altitude profiles of minor constituents in an atmosphere. Unfortunately, laboratory measurements of these coefficients were never conducted at temperature conditions relevant to the atmosphere of Titan. Here we conduct a detailed uncertainty analysis of the bimolecular diffusion coefficient parameters as applied to Titan's upper atmosphere to provide a better understanding of the impact of uncertainty for this parameter on models. Because temperature and pressure conditions are much lower than the laboratory conditions in which bimolecular diffusion parameters were measured, we apply a Bayesian framework, a problem-agnostic framework, to determine parameter estimates and associated uncertainties. We solve the Bayesian calibration problem using the open-source QUESO library which also performs a propagation of uncertainties in the calibrated parameters to temperature and pressure conditions observed in Titan's upper atmosphere. Our results show that, after propagating uncertainty through the Massman model, the uncertainty in molecular diffusion is highly correlated to temperature and we observe no noticeable correlation with pressure. We propagate the calibrated molecular diffusion estimate and associated uncertainty to obtain an estimate with uncertainty due to bimolecular diffusion for the methane molar fraction as a function of altitude. Results show that the uncertainty in methane abundance due to molecular diffusion is in general small compared to eddy diffusion and the chemical kinetics description. However, methane abundance is most sensitive to uncertainty in molecular diffusion above 1200 km where the errors are nontrivial and could have important implications for scientific research based on diffusion models in this altitude range.