The Ortho-to-Para Ratio of Ammonia in the L1157 Outflow

Theoretical Distribution of the Ammonia Binding Energy at Interstellar Icy Grains: A New Computational Framework

The binding energies (BE) of molecules on the interstellar grains are crucial in the chemical evolution of the interstellar medium (ISM). Both temperature-programmed desorption (TPD) laboratory experiments and quantum chemistry computations have often provided, so far, only single values of the BE for each molecule. This is a severe limitation, as the ices enveloping the grain mantles are structurally amorphous, giving rise to a manifold of possible adsorption sites, each with different BEs. However, the amorphous ice nature prevents the knowledge of structural details, hindering the development of a common accepted atomistic icy model. In this work, we propose a computational framework that closely mimics the formation of the interstellar grain mantle through a water by water accretion. On that grain, an unbiased random (but well reproducible) positioning of the studied molecule is then carried out. Here we present the test case of NH₃, a ubiquitous species in the molecular ISM. We provide the BE distribution computed by a hierarchy approach, using the semiempirical xTB-GFN2 as a low-level method to describe the whole icy cluster in combination with the B97D3 DFT functional as a high-level method on the local zone of the NH₃ interaction. The final ZPE-corrected BE is computed at the ONIOM(DLPNO-CCSD(T)//B97D3:xTB-GFN2) level, ensuring the best cost/accuracy ratio. The main peak of the predicted NH₃ BE distribution is in agreement with experimental TPD and computed data in the literature. A second broad peak at very low BE values is also present, which has never been detected before. It may provide the solution to a longstanding puzzle about the presence of gaseous NH₃ also observed in cold ISM objects.

Herschel/HIFI observations of the circumstellar ammonia lines in IRC+10216

CONTEXT

A discrepancy exists between the abundance of ammonia (NH3) derived previously for the circumstellar envelope (CSE) of IRC+10216 from far-IR submillimeter rotational lines and that inferred from radio inversion or mid-infrared (MIR) absorption transitions.

AIMS

To address the discrepancy described above, new high-resolution far-infrared (FIR) observations of both ortho- and para-NH3 transitions toward IRC+10216 were obtained with Herschel, with the goal of determining the ammonia abundance and constraining the distribution of NH3 in the envelope of IRC+10216.

METHODS

We used the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel to observe all rotational transitions up to the J = 3 level (three ortho- and six para-NH3 lines). We conducted non-LTE multilevel radiative transfer modelling, including the effects of near-infrared (NIR) radiative pumping through vibrational transitions. The computed emission line profiles are compared with the new HIFI data, the radio inversion transitions, and the MIR absorption lines in the ν2 band taken from the literature.

RESULTS

We found that NIR pumping is of key importance for understanding the excitation of rotational levels of NH3. The derived NH3 abundances relative to molecular hydrogen were (2.8 ± 0.5) × 10-8 for ortho-NH3 and [Formula: see text] for para-NH3, consistent with an ortho/para ratio of 1. These values are in a rough agreement with abundances derived from the inversion transitions, as well as with the total abundance of NH3 inferred from the MIR absorption lines. To explain the observed rotational transitions, ammonia must be formed near to the central star at a radius close to the end of the wind acceleration region, but no larger than about 20 stellar radii (1σ confidence level).