Icy Grains Research Articles

First detection of gas-phase ammonia in a planet-forming disk

Nitrogen chemistry in protoplanetary disks and the freeze-out on dust particles is key to understand the formation of nitrogen bearing species in early solar system analogs. So far, ammonia has not been detected beyond the snowline in protoplanetary disks. We aim to find gas-phase ammonia in a protoplanetary disk and characterize its abundance with respect to water vapor. Using HIFI on the Herschel Space Observatory we detect, for the first time, the ground-state rotational emission of ortho-NH$_3$ in a protoplanetary disk, around TW Hya. We use detailed models of the disk's physical structure and the chemistry of ammonia and water to infer the amounts of gas-phase molecules of these species. We explore two radial distributions ( confined to $<$60 au like the millimeter-sized grains) and two vertical distributions (near the midplane where water is expected to photodesorb off icy grains) to describe the (unknown) location of the molecules. These distributions capture the effects of radial drift and vertical settling of ice-covered grains. We use physical-chemical models to reproduce the fluxes with assuming that water and ammonia are co-spatial. We infer ammonia gas-phase masses of 0.7-11.0 $\times$10$^{21}$ g. For water, we infer gas-phase masses of 0.2-16.0 $\times$10$^{22}$ g. This corresponds to NH$_3$/H$_2$O abundance ratios of 7\%-84\%, assuming that water and ammonia are co-located. Only in the most compact and settled adopted configuration is the inferred NH$_3$/H$_2$O consistent with interstellar ices and solar system bodies of $\sim$ 5\%-10\%. Volatile release in the midplane may occur via collisions between icy bodies if the available surface for subsequent freeze-out is significantly reduced, e.g., through growth of small grains into pebbles or larger.

The radial dependence of pebble accretion rates: A source of diversity in planetary systems

Context. The classical "planetesimal" accretion scenario for the formation of planets has recently evolved with the idea that "pebbles", centimeter- to meter-sized icy grains migrating in protoplanetary disks, can control planetesimal and/or planetary growth. Aims. We investigate how pebble accretion depends on disk properties and affects the formation of planetary systems Methods. We construct analytical models of pebble accretion onto planetary embryos that consistently account for the mass and orbital evolution of the pebble flow and reflect disk structure. Results. We derive simple formulas for pebble accretion rates in the so-called "settling" regime for planetary embryos with more than 100 km in size. For relatively smaller embryos or in outer disk regions, the accretion mode is 3D, meaning that the thickness of the pebble flow must be taken into account, and resulting in an accretion rate that is independent of the embryo mass. For larger embryos or in inner regions, the accretion is in a 2D mode, i.e., the pebble disk may be considered to be infinitely thin. We show that the radial dependence of the pebble accretion rata is different (even the sign of the power-law exponent changes) for different disk conditions such as the disk heating source (viscous heating or stellar irradiation), the drag law (Stokes or Epstein, and weak or strong coupling), and in the 2D or 3D accretion modes. We also discuss the effect of the sublimation/destruction of icy pebbles inside the snow line. Conclusions. Pebble accretion easily produces a large diversity of planetary systems. In other words, to correctly infer the results of planet formation through pebble accretion, detailed prescriptions of disk evolution and pebble growth, sublimation/destruction and migration are required.

ORIGIN OF MOLECULAR OXYGEN IN COMET 67P/CHURYUMOV–GERASIMENKO

ABSTRACT Molecular oxygen has been detected in the coma of comet 67P/Churyumov–Gerasimenko with abundances in the 1%–10% range by the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis-Double Focusing Mass Spectrometer 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 (PSN), 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 PSN 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 comet 67P/Churyumov-Gerasimenko’s coma while the latter scenario requires that clathrates formed from ISM icy grains that crystallized when entering the PSN.

Isomerization and fragmentation of acetonitrile upon interaction with N(4S) atoms: the chemistry of nitrogen in dense molecular clouds

We experimentally show that the reaction between ground state nitrogen atoms N(4S) and acetonitrile CH3CN can lead to two distinct chemical pathways that are both thermally activated at very low temperatures. First is CH3CN isomerization which produces CH3NC and H2CCNH. Second is CH3CN decomposition which produces HNC and CH3CNH+CN− fragments, with the possible release of H2. Our results reveal that the mobility of N(4S)-atoms is stimulated in the 3–11 K temperature range, and that its subsequent encounter with one acetonitrile molecule is sufficient for the aforementioned reactions to occur without the need for additional energy to be supplied to the CH3CN + N(4S) system. These findings shed more light on the nitrogen chemistry that can possibly take place in dense molecular clouds, which until now was thought to only involve high-energy processes and therefore be unlikely to occur in such cold and dark interstellar regions. The reaction pathways we propose in this study have very important astrochemical implications, as it was shown recently that the atomic nitrogen might be more abundant, in many interstellar icy grain mantles, than previously thought. Also, these reaction pathways can now be considered within dense molecular clouds, and possibly affect the branching ratios for N-bearing molecules computed in astrochemical modelling.

Photochemistry and Astrochemistry: Photochemical Pathways to Interstellar Complex Organic Molecules.

The interstellar medium is characterized by a rich and diverse chemistry. Many of its complex organic molecules are proposed to form through radical chemistry in icy grain mantles. Radicals form readily when interstellar ices (composed of water and other volatiles) are exposed to UV photons and other sources of dissociative radiation, and if sufficiently mobile the radicals can react to form larger, more complex molecules. The resulting complex organic molecules (COMs) accompany star and planet formation and may eventually seed the origins of life on nascent planets. Experiments of increasing sophistication have demonstrated that known interstellar COMs as well as the prebiotically interesting amino acids can form through ice photochemistry. We review these experiments and discuss the qualitative and quantitative kinetic and mechanistic constraints they have provided. We finally compare the effects of UV radiation with those of three other potential sources of radical production and chemistry in interstellar ices: electrons, ions, and X-rays.

UNIFIED MICROSCOPIC–MACROSCOPIC MONTE CARLO SIMULATIONS OF COMPLEX ORGANIC MOLECULE CHEMISTRY IN COLD CORES

ABSTRACT The recent discovery of methyl formate and dimethyl ether in the gas phase of cold cores with temperatures as cold as 10 K challenges our previous astrochemical models concerning the formation of complex organic molecules (COMs). The strong correlation between the abundances and distributions of methyl formate and dimethyl ether further shows that current astrochemical models may be missing important chemical processes in cold astronomical sources. We investigate a scenario in which COMs and the methoxy radical can be formed on dust grains via a so-called chain reaction mechanism, in a similar manner to CO2. A unified gas-grain microscopic–macroscopic Monte Carlo approach with both normal and interstitial sites for icy grain mantles is used to perform the chemical simulations. Reactive desorption with varying degrees of efficiency is included to enhance the nonthermal desorption of species formed on cold dust grains. In addition, varying degrees of efficiency for the surface formation of methoxy are also included. The observed abundances of a variety of organic molecules in cold cores can be reproduced in our models. The strong correlation between the abundances of methyl formate and dimethyl ether in cold cores can also be explained. Nondiffusive chemical reactions on dust grain surfaces may play a key role in the formation of some COMs.

The (impossible?) formation of acetaldehyde on the grain surfaces: insights from quantum chemical calculations

Abstract Complex Organic Molecules (COMs) have been detected in the interstellar medium (ISM). However, it is not clear whether their synthesis occurs on the icy surfaces of interstellar grains or via a series of gas-phase reactions. As a test case of the COMs synthesis in the ISM, we present new quantum chemical calculations on the formation of acetaldehyde (CH3CHO) from the coupling of the HCO and CH3 radicals, both in gas phase and on water ice surfaces. The binding energies of HCO and CH3 on the amorphous water ice were also computed (2333 and 734 K, respectively). Results indicate that, in gas phase, the products could be either CH3CHO, CH4 + CO, or CH3OCH, depending on the relative orientation of the two radicals. However, on the amorphous water ice, only the CH4 + CO product is possible due to the geometrical constraints imposed by the water ice surface. Therefore, acetaldehyde cannot be synthesized by the CH3 + HCO coupling on the icy grains. We discuss the implications of these results and other cases, such as ethylene glycol and dimethyl ether, in which similar situations can occur, suggesting that formation of these molecules on the grain surfaces might be unlikely.

PENTACARBON DIOXIDE (C5O2) FORMATION AND ITS ROLE AS A TRACER OF SOLAR SYSTEM EVOLUTION

ABSTRACT Carbon monoxide is the second most abundant molecule on icy grains in the interstellar medium. These grains are under the influence of ionizing radiation, which induces the chemical reaction within the ice. Here we report the first observation of subliming pentacarbon dioxide (C5O2) after irradiation of pure carbon monoxide ice with energetic electrons. Our results show that pentacarbon dioxide is a stable reaction product in a carbon monoxide matrix that survives the sublimation in star-forming regions at sublimation temperatures of 175 K. Along with carbon suboxide (C3O2), this molecule can serve as a powerful tracer of the temperature history of formerly carbon monoxide rich ices in molecular clouds and star-forming regions.

A panoptic model for planetesimal formation and pebble delivery

The journey from dust particle to planetesimal involves physical processes acting on scales ranging from micrometers (the sticking and restructuring of aggregates) to hundreds of astronomical units (the size of the turbulent protoplanetary nebula). Considering these processes simultaneously is essential when studying planetesimal formation. We develop a novel, global, semi-analytical model for the evolution of the mass-dominating dust particles in a turbulent protoplanetary disk that takes into account the evolution of the dust surface density while preserving the essential characteristics of the porous coagulation process. This panoptic model is used to study the growth from submicron to planetesimal sizes in disks around Sun-like stars. For highly porous ices, unaffected by collisional fragmentation and erosion, rapid growth to planetesimal sizes is possible in a zone stretching out to ${\sim}10\mathrm{~AU}$ for massive disks. When porous coagulation is limited by erosive collisions, the formation of planetesimals through direct coagulation is not possible, but the creation of a large population of aggregates with Stokes numbers close to unity might trigger the streaming instability (SI). However, we find that reaching conditions necessary for SI is difficult and limited to dust-rich disks, (very) cold disks, or disks with weak turbulence. Behind the snow-line, porosity-driven aggregation of icy grains results in rapid (${\sim}10^{4}\mathrm{~yr}$) formation of planetesimals. If erosive collisions prevent this, SI might be triggered for specific disk conditions. The numerical approach introduced in this work is ideally suited for studying planetesimal formation and pebble delivery simultaneously and will help build a coherent picture of the start of the planet formation process.

DIRECT IMAGING OF THE WATER SNOW LINE AT THE TIME OF PLANET FORMATION USING TWO ALMA CONTINUUM BANDS

Molecular snow lines in protoplanetary disks have been studied theoretically for decades because of their importance in shaping planetary architectures and compositions. The water snow line lies in the planet formation region at < 10 AU, and so far its location has been estimated only indirectly from spatially-unresolved spectroscopy. This work presents a proof-of-concept method to directly image the water snow line in protoplanetary disks through its physical and chemical imprint in the local dust properties. We adopt a physical disk model that includes dust coagulation, fragmentation, drift, and a change in fragmentation velocities of a factor 10 between dry silicates and icy grains as found by laboratory work. We find that the presence of a water snow line leads to a sharp discontinuity in the radial profile of the dust emission spectral index {\alpha}_mm, due to replenishment of small grains through fragmentation. We use the ALMA simulator to demonstrate that this effect can be observed in protoplanetary disks using spatially-resolved ALMA images in two continuum bands. We explore the model dependence on the disk viscosity and find that the spectral index reveals the water snow line for a wide range of conditions, with opposite trends when the emission is optically thin rather than thick. If the disk viscosity is low ({\alpha}_visc < 10^-3) the snow line produces a ring-like structure with a minimum at {\alpha}_mm ~ 2 in the optically thick regime, possibly similar to what has been measured with ALMA in the innermost region of the HL Tau disk.

TRAPPIST photometry and imaging monitoring of comet C/2013 R1 (Lovejoy): Implications for the origin of daughter species

We report the results of the narrow band photometry and imaging monitoring of comet C/2013 R1 (Lovejoy) with the robotic telescope TRAPPIST (La Silla observatory). We gathered around 400 images over 8 months pre- and post-perihelion between September 12, 2013 and July 6, 2014. We followed the evolution of the OH, NH, CN, C3 , and C2 production rates computed with the Haser model as well as the evolution of the dust production. All five gas species display an asymmetry about perihelion, the rate of brightening being steeper than the rate of fading. The study of the coma morphology reveals gas and dust jets which indicate one or several active zone(s) on the nucleus. The dust, C2 , and C3 morphologies present some similarities while the CN morphology is different. OH and NH are enhanced in the tail direction. The study of the evolution of the comet activity shows that the OH, NH, and C2 production rates evolution with the heliocentric distance is correlated to the dust evolution. The CN and, to a lesser extent, the C3 do not display such a correlation with the dust. These evidences and the comparison with parent species production rates indicate that C2 and C3 on one side and OH and NH on the other side could be -at least partially- released from organic-rich grains and icy grains. On the contrary, all evidences point to HCN being the main parent of CN in this comet.

The CO2 abundance in Comets C/2012 K1 (PanSTARRS), C/2012 K5 (LINEAR), and 290P/Jäger as measured with Spitzer

Carbon dioxide is one of the most abundant ices present in comets and is therefore important for understanding cometary composition and activity. We present analysis of observations of CO2 and [OI] emission in three comets to measure the CO2 abundance and evaluate the possibility of employing observations of [OI] emission in comets as a proxy for CO2. We obtained NIR imaging sensitive to CO2 of comets C/2012 K1 (PanSTARRS), C/2012 K5 (LINEAR), and 290P/Jäger with the IRAC instrument on Spitzer. We acquired observations of [OI] emission in these comets with the ARCES echelle spectrometer mounted on the 3.5-m telescope at Apache Point Observatory and observations of OH with the Swift observatory (PanSTARRS) and with Keck HIRES (Jäger). The CO2/H2O ratios derived from the Spitzer images are 12.6±1.3% (PanSTARRS), 28.9±3.6% (LINEAR), and 31.3±4.2% (Jäger). These abundances are derived under the assumption that contamination from CO emission is negligible. The CO2 abundance for PanSTARRS is close to the average abundance measured in comets at similar heliocentric distance to date, while the abundances measured for LINEAR and Jäger are significantly larger than the average abundance. From the coma morphology observed in PanSTARRS and the assumed gas expansion velocity, we derive a rotation period for the nucleus of about 9.2h. Comparison of H2O production rates derived from ARCES and Swift data, as well as other observations, suggest the possibility of sublimation from icy grains in the inner coma. We evaluate the possibility that the [OI] emission can be employed as a proxy for CO2 by comparing CO2/H2O ratios inferred from the [OI] lines to those measured directly by Spitzer. We find that for PanSTARRS we can reproduce the observed CO2 abundance to an accuracy of ∼20%. For LINEAR and Jäger, we were only able to obtain upper limits on the CO2 abundance inferred from the [OI] lines. These upper limits are consistent with the CO2 abundances measured by Spitzer.

Low-temperature chemistry between water and hydroxyl radicals: H/D isotopic effects

Sets of systematic laboratory experiments are presented – combining Ultra High Vacuum cryogenic and plasma-line deposition techniques – that allow us to compare H/D isotopic effects in the reaction of H2O (D2O) ice with the hydroxyl radical OD (OH). The latter is known to play a key role as intermediate species in the solid-state formation of water on icy grains in space. The main finding of our work is that the reaction H2O + OD → OH + HDO occurs and that this may affect the HDO/H2O abundances in space. The opposite reaction D2O + OH → OD + HDO is much less effective, and also given the lower D2O abundances in space not expected to be of astronomical relevance. The experimental results are extended to the other four possible reactions between hydroxyl and water isotopes and are subsequently used as input for Kinetic Monte Carlo simulations. This way we interpret our findings in an astronomical context, qualitatively testing the influence of the reaction rates.

ION IRRADIATION OF H2-LADEN POROUS WATER-ICE FILMS: IMPLICATIONS FOR INTERSTELLAR ICES

To understand the effects of cosmic-ray (CR) impacts on interstellar icy grains immersed in H2 gas, we have irradiated porous water-ice films loaded with H2 with 100 keV H+. The ice films were exposed to H2 gas at different pressures following deposition and during irradiation. A net H2 loss is observed during irradiation due to competition between ion-induced sputtering and gas adsorption. The initial H2 loss cross-section, 4(1) × 10−14 cm2, was independent of film thickness, H2, and proton fluxes. In addition to sputtering, irradiation also closes nanopores, trapping H2 in the film with binding that exceeds physical absorption energies. As a result, 2%–7% H2 is retained in the ice following irradiation to high fluences. We find that the trapped H2 concentration increases with decreasing Φ, the ratio of ion to H2 fluxes, suggesting that as high as 8% solid H2 can be trapped in interstellar ice by CR or stellar wind impacts.

DOUBLE DCO+RINGS REVEAL CO ICE DESORPTION IN THE OUTER DISK AROUND IM LUP

In a protoplanetary disk, a combination of thermal and non-thermal desorption processes regulate where volatiles are liberated from icy grain mantles into the gas phase. Non-thermal desorption should result in volatile-enriched gas in disk-regions where complete freeze-out is otherwise expected. We present ALMA observations of the disk around the young star IM Lup in 1.4 mm continuum, C18O 2-1, H13CO+ 3-2 and DCO+ 3-2 emission at ~0".5 resolution. The images of these dust and gas tracers are clearly resolved. The DCO+ line exhibits a striking pair of concentric rings of emission that peak at radii of ~0".6 and 2" (~90 and 300 AU, respectively). Based on disk chemistry model comparison, the inner DCO+ ring is associated with the balance of CO freeze-out and thermal desorption due to a radial decrease in disk temperature. The outer DCO+ ring is explained by non-thermal desorption of CO ice in the low-column-density outer disk, repopulating the disk midplane with cold CO gas. The CO gas then reacts with abundant H2D+ to form the observed DCO+ outer ring. These observations demonstrate that spatially resolved DCO+ emission can be used to trace otherwise hidden cold gas reservoirs in the outmost disk regions, opening a new window onto their chemistry and kinematics.

Formation, destruction and chemical influences of water ice: A review of recent laboratory results

AbstractWater ice is the dominant constituent of icy grain mantles in the interstellar medium, and as such one of the most abundant species during all stages of star and planet formation. Its formation through atom addition reactions on grain surfaces, its destruction through different desorption channels, and its influence on the chemistry and desorption efficiencies of other species in icy grain mantles have all been the objects of intense study. This contribution reviews our current understanding of these processes, and the laboratory experiments that have been instrumental in establishing the existing paradigm.

Linking the Origin of Asteroids to Planetesimal Formation in the Solar Nebula

AbstractThe asteroids (more precisely: objects of the main asteroid belt) and Kuiper Belt objects (more precisely: objects of the cold classical Kuiper Belt) are leftovers of the building material for our earth and all other planets in our solar system from more than 4.5 billion years ago. At the time of their formation those were typically 100 km large objects. They were called planetesimals, built up from icy and dusty grains. In our current paradigm of planet formation it was turbulent flows and metastable flow patterns, like zonal flows and vortices, that concentrated mm to cm sized icy dust grains in sufficient numbers that a streaming instability followed by a gravitational collapse of these particle clump was triggered. The entire picture is sometimes referred to as gravoturbulent formation of planetesimals. What was missing until recently, was a physically motivated prediction on the typical sizes at which planetesimals should form via this process. Our numerical simulations in the past had only shown a correlation between numerical resolution and planetesimal size and thus no answer was possible (Johansen et al.2011). But with the lastest series of simulations on JUQUEEN (Stephan &amp; Doctor 2015), covering all the length scales down to the physical size of actual planetesimals, we were able to obtain values for the turbulent particle diffusion as a function of the particle load in the gas. Thus, we have all necessary data at hand to feed a 'back of the envelope' calculation that predicts the size of planetesimals as result of a competition between gravitational concentration and turbulent diffusion. Using the diffusion values obtained in the numerical simulations it predicts planetesimal sizes on the order of 100 km, which suprisingly coincides with the measured data from both asteroids (Bottke et al.2005) as well from Kuiper Belt objects (Nesvorny et al.2011).

Laboratory constraints on ice formation, restructuring and desorption

AbstractIces form on the surfaces of interstellar and circumstellar dust grains though freeze-out of molecules and atoms from the gas-phase followed by chemical reactions. The composition, chemistry, structure and desorption properties of these ices regulate two important aspects of planet formation: the locations of major condensation fronts in protoplanetary disks (i.e. snow lines) and the formation efficiencies of complex organic molecules in astrophysical environments. The latter regulates the availability of prebiotic material on nascent planets. With ALMA it is possible to directly observe both (CO) snowlines and complex organics in protoplanetary disks. The interpretation of these observations requires a detailed understanding of the fundamental ice processes that regulate the build-up, evolution and desorption of icy grain mantles. This proceeding reviews how experiments on thermal CO and N2 ice desorption, UV photodesorption of CO ice, and CO diffusion in H2O ice have been used to guide and interpret astrochemical observations of snowlines and complex molecules.

THE DISTRIBUTION OF DEUTERATED FORMALDEHYDE WITHIN ORION-KL

We report the first high angular resolution imaging (3.4\arcsec $\times$ 3.0\arcsec) of deuterated formaldehyde (HDCO) toward Orion--KL, carried out with the Submillimeter Array (SMA). We find that the spatial distribution of the formaldehyde emission systematically differs from that of methanol: while methanol is found towards the inner part of the region, HDCO is found in colder gas that wraps around the methanol emission on four sides. The HDCO/H$_2$CO ratios are determined to be 0.003--0.009 within the region, up to an order of magnitude higher than the D/H measured for methanol. These findings strengthen the previously suggested hypothesis that there are differences in the chemical pathways leading to HDCO (via deuterated gas phase chemistry) and deuterated methanol (through conversion of formaldehyde into methanol on the surface of icy grain mantles).

Shedding light on the formation of the pre-biotic molecule formamide with ASAI

Formamide (NH2CHO) has been proposed as a pre-biotic precursor with a key role in the emergence of life on Earth. While this molecule has been observed in space, most of its detections correspond to high-mass star-forming regions. Motivated by this lack of investigation in the low-mass regime, we searched for formamide, as well as isocyanic acid (HNCO), in 10 low- and intermediate-mass pre-stellar and protostellar objects. The present work is part of the IRAM Large Programme ASAI (Astrochemical Surveys At IRAM), which makes use of unbiased broadband spectral surveys at millimetre wavelengths. We detected HNCO in all the sources and NH2CHO in five of them. We derived their abundances and analysed them together with those reported in the literature for high-mass sources. For those sources with formamide detection, we found a tight and almost linear correlation between HNCO and NH2CHO abundances, with their ratio being roughly constant -between 3 and 10- across 6 orders of magnitude in luminosity. This suggests the two species are chemically related. The sources without formamide detection, which are also the coldest and devoid of hot corinos, fall well off the correlation, displaying a much larger amount of HNCO relative to NH2CHO. Our results suggest that, while HNCO can be formed in the gas phase during the cold stages of star formation, NH2CHO forms most efficiently on the mantles of dust grains at these temperatures, where it remains frozen until the temperature rises enough to sublimate the icy grain mantles. We propose hydrogenation of HNCO as a likely formation route leading to NH2CHO.