Modelling of deuterium chemistry and its application to molecular clouds

We have developed new models of the chemistry of deuterium for investigating fractionation in interstellar molecular clouds. We have incorporated the latest information on reactions which affect deuteration, extended previous models to include S-D bonds for the first time and included the gasphase chemistry of some doubly-deuterated species. We present models for a wide range of physical parameters, including density, temperature, elemental abundances, and the freeze out of molecules on to dust grains. We discuss the detailed fractionation of particular species and show how fractionation can be used to probe the history of interstellar matter. The freeze out of molecules onto dust leads to significant enhancement in fractionation ratios and, in particular, to large fractionation in doubly-deuterated species.

SO and CS observations of molecular clouds - II. Analysis and modelling of the abundance ratios - probing O-2/CO with SO/CS?

We here analyse the observational SO and CS data presented in Nilsson ct al. (2000). The SO/CS integrated intensity ratio maps are presented for 19 molecular clouds, together with tables of relevant ratios at strategic positions, where we have also observed (SO)-S-34 and/or (CS)-S-34. The SO/CS abundance ratio as calculated from an LTE analysis is highly varying within and between the sources. Our isotopomer observations and Monte Carlo simulations verify that this is not an artifact due to optical depth problems. The variation of the maximum SO/CS abundance ratio between the clouds is 0.2-7. The largest variations within a cloud are found for the most nearby objects, possibly indicating resolution effects. We have also performed time dependent chemical simulations. We compare the simulations with our observed SO/CS abundance ratios and suggest a varying oxygen to carbon initial abundance, differing temporal evolution, density differences and X-ray sources associated with young stellar objects as possible explanations to the variations. In particular, the observed variation of the maximum SO/CS abundance ratio between the clouds can be explained by using initial O/C+ abundance ratios in the range 1.3-2.5. We finally derive a relationship between the SO/CS and O-2/CO abundance ratios, which may be used as a guide to find the most promising interstellar O-2 search targets.

Turbulent diffusion and mixing in interstellar dark clouds with adsorption onto grains

This paper reports the results of models of dark cloud chemistry incorporating a depth dependent density distribution with diffusive mixing and adsorption onto grains. The model is based on the approach taken by Xie et al. (1995), with the addition of grain accretion effects. Without diffusion, the central regions of the cloud freeze out in less than 10(7) years. Freeze-out time is dependent on density, so the diffuse outer region of the cloud remains abundant in gas for about an order of magnitude longer. We find that fairly small amounts of diffusive mixing can delay freeze-out at the centre of the model cloud for a time up to an order of magnitude greater than without diffusion, due to material diffusing inward from the edges of the cloud. The gas-phase lifetime of the cloud core can thus be increased by up to an order of magnitude or more by this process. We have run three different grain models with various diffusion coefficients to investigate the effects of changing the sticking parameters.

A 7000 yr perspective on volcanic ash clouds affecting northern Europe

The ash cloud resulting from the 2010 eruption of Eyjafjöll caused severe disruption to air travel across Europe but as a geological event, it is not unprecedented. Analysis of peat and lake sediments from northern Europe has revealed the presence of microscopic layers of Icelandic volcanic ash (tephra). These sedimentary records, together with historical records of Holocene ash falls, demonstrate that Icelandic volcanoes have generated substantial ash clouds that reached northern Europe many times. Here we present the first comprehensive compilation of sedimentary and historical records of ash-fall events in northern Europe, spanning the last 7000 years. Within this period ten tephra layers have been identified in the Faroe Islands, 14 in Great Britain, 11 in Germany, 38 in Scandinavia and 33 in Ireland. Seven ash fall events have been historically documented prior to the Eyjafjöll 2010 event. Ash fall events appear to be more frequent in the last 1500 years, but it is unclear whether this reflects a true increase in eruption frequency or dispersal, or is an artefact of the records themselves or the way they have been generated. In the last 1,000 years, volcanic ash clouds reached Northern Europe with a mean return interval of 53 ± 8 years (the range of return intervals is between 6 and 112 years). Modelling using the ash records for the last millennium indicates that for any 10 year period there is a 17% probability of tephra fallout event in Northern Europe. These values must be considered as conservative estimates due to the nature of tephra capture and preservation in the sedimentary record.

Desorption processes and the deuterium fractionation in molecular clouds

Several processes have been suggested as ways of returning accreted grain mantles to the gas, thus preventing the total removal of molecules from the gas phase in dark quiescent clouds. We attempt to distinguish between them by considering not only the calculated gas-phase abundances, but also the ratio of the abundances of deuterated species to non-deuterated species. We find that the D/H ratio in molecules is relatively model-independent, but that desorption due to the formation of H-2 on grains gives the best overall agreement with the observations.

Molecular abundances in the magellanic clouds III. LIRS 36, a star-forming region in the Small Magellanic Cloud

Detections of CO, CS, SO, C2H, HCO+, HCN, HNC, H2CO, and C3H2 are reported from LIRS 36, a star-forming region in the Small Magellanic Cloud. (CO)-O-18, NO, CH3OH, and most notably CN have not been detected, while the rare isotopes (CO)-C-13 and, tentatively, (CS)-S-34 ar,seen. This is so far the most extensive molecular multiline study of an interstellar medium with a heavy element depletion exceeding a factor of four.

Molecular abundances in the magellanic clouds .1. A multiline study of five cloud cores

Nine H II regions of the LMC were mapped in (CO)-C-13(1-0) and three in (CO)-C-12(1-0) to study the physical properties of the interstellar medium in the Magellanic Clouds. For N113 the molecular core is found to have a peak position which differs from that of the associated H II region by 20 ''. Toward this molecular core the (CO)-C-12 and (CO)-C-13 peak T-MB line temperatures of 7.3 K and 1.2 K are the highest so far found in the Magellanic Clouds. The molecular concentrations associated with N113, N44BC, N159HW, and N214DE in the LMC and LIRS 36 in the SMC were investigated in a variety of molecular species to study the chemical properties of the interstellar medium. I(HCO+)/I(HCN) and I(HCN)/I(HNC) intensity ratios as well as lower limits to the I((CO)-C-13)/I((CO)-O-18) ratio were derived for the rotational 1-0 transitions. Generally, HCO+ is stronger than HCN, and HCN is stronger than HNC. The high relative HCO+ intensities are consistent with a high ionization flux from supernovae remnants and young stars, possibly coupled with a large extent of the HCO+ emission region. The bulk of the HCN arises from relatively compact dense cloud cores. Warm or shocked gas enhances HCN relative to HNC. From chemical model calculations it is predicted that I(HCN)/I(HNC) close to one should be obtained with higher angular resolution (less than or similar to 30 '') toward the cloud cores. Comparing virial masses with those obtained from the integrated CO intensity provides an H-2 mass-to-CO luminosity conversion factor of 1.8 x 10(20) mol cm(-2) (K km s(-1))(-1) for N113 and 2.4 x 10(20) mol cm(-2) (K km s(-1))(-1) for N44BC. This is consistent with values derived for the Galactic disk.

Molecular abundances in the Magellanic Clouds .2. Deuterated species in the LMC

The first definite discoveries of extragalactic deuterium are reported. DCO+ has been detected in three and DCN has been measured in one star-forming region of the Large Magellanic Cloud (LMC). While the HCO+/DCO+ abundance ratios are found to be 19 +/- 3, 24 +/- 4, and 67 +/- 18 for N113, N44BC and N159HW, respectively, a HCN/DCN abundance ratio of 23 +/- 5 is obtained for N113. These results are consistent with a gas temperature of about 20 K and a D/H ratio of about 1.5 x 10(-5), consistent with that observed in the Galaxy. If the cloud temperature is closer to 30 K, then a D/H ratio is required to be up to an order of magnitude larger. Because this ratio provides a lower limit to the primordial D/H ratio, it indicates that the baryon mass density alone is unable to close the universe.

THE CHEMISTRY OF COMPLEX-MOLECULES IN INTERSTELLAR CLOUDS

This paper is concerned with the chemical evolution of large molecules in interstellar clouds. We consider the chemistry and ionisation balance of large polycyclic aromatic hydrocarbon (PAH) type molecules in diffuse clouds and show that certain PAH molecules can be doubly ionised by the interstellar ultraviolet radiation field. If recombination of the dications so produced with electrons is dissociative rather than radiative, then PAHs are rapidly destroyed. PAHs which can only be singly ionised have much smaller recombination energies and can be long lasting in these regions. This type of property may be very important in selecting the PAH species which can populate the general interstellar medium and account for certain of the diffuse bands observed in optical spectra. Destruction of PAH molecules via formation of dications may be responsible for the weakening of the diffuse bands observed in regions of high UV flux.

THE EFFECT OF RAPID NEUTRAL-NEUTRAL REACTIONS ON CHEMICAL-MODELS OF DENSE INTERSTELLAR CLOUDS

Recent experiments on rapid neutral-neutral reactions involving the radical CN at low temperature and the neutral C atom at room temperature suggest that atom-neutral and radical-neutral reactions may be generally more rapid at low temperature than hitherto thought. We have included a variety of rapid neutral-neutral reactions in our gas-phase chemical models of quiescent, dense interstellar clouds. We find the calculated abundances of many molecules to be greatly changed from previous values. In particular, the peak 'early-time' abundances of organic molecules are reduced.

OBSERVATIONS OF DEUTERATED CYANOACETYLENE IN DARK CLOUDS

We have observed DC3N and HC3N in a number of cold dust clouds in order to derive the degree of deuterium fractionation. We find that the ratio of DC3N to HC3N is large, at about 0.05 or more, and discuss the implications of this result for the synthesis of cyanoacetylene. The observations are most readily interpreted if the deuteration of HC3N is linked to that of cyclic C3H2, which is also observed to exhibit a large degree of deuterium fractionation. HC3N deuteration levels comparable with those we observed are found to he just compatible with the mechanism suggested by Howe & Millar, but with adjusted rate coefficients. Freeze-out on to grain surfaces is also considered, but produces widespread deuterium enhancement in many species. contrary to observed levels.

ALTERNATIVE ROUTES TO DEUTERATION IN DARK CLOUDS

We investigate new approaches to the deuteration of C3H2, HC3N and HC5N in dark clouds, following the suggestion that protonated HC3N might form different isomers, a linear structure (HC3NH+) being the most stable. We consider the effect of linear HC3NH+ and HC5NH+ on the formation of HC3N and HC5N, and find that deuteration ratios at approximately 10 K are reduced, in the case of HC3N to values significantly below those observed, such that a deuteration mechanism other than direct deuteron transfer is probably required for cold clouds.

ON THE ORIGIN OF NH IN DIFFUSE INTERSTELLAR CLOUDS

The observation by Meyer & Roch of NH in the interstellar clouds towards zeta Per and HD 27778 cannot be explained with conventional gas-phase chemistry models. A simple non-equilibrium model for the zeta Per cloud, which incorporates the grain-surface production of NH and OH or, alternatively, NH3 and H2O, is able to reproduce the abundances of all observed species (except CH+) quite accurately. Moreover, chemical models which include grain-surface reactions can reproduce the observed abundance not only of NH but also of CN, which is efficiently formed at low temperatures, initiated by the reaction of NH with C+. Pure gas-phase models and cloud interface models, in which NH and CH+ are formed in a warm and tenuous environment, fail to explain the observed high abundance of CN. Hence the observation of NH in zeta Per and HD 27778 provides evidence for the presence of grain-surface reactions leading to molecules other than H-2. It is predicted that NH2 and NH3 should have abundances not much below that of NH if NH3 instead of NH is formed on grains. With or without surface reactions, the column densities of H2O and C2H are expected to be about 10(13) cm-2, and these molecules may be detectable in the zeta Per cloud.