Photodissociation of interstellar ArH+

Aims. Following the recent detection of 36ArH+ in the Crab nebula spectrum, we have computed the photodissociation rate of ArH+ in order to constrain the physical processes at work in this environment. Methods. Photodissociation cross sections of ArH+ are computed in an ab initio approach including explicit account of spin-orbit coupling. Results. We report the photodissociation cross section of ArH+ as a function of wavelength. Photodissociation probabilities are derived for different impinging radiation fields.The photodissociation probability of for a very small unshielded cloud surrounded on all sides by the unshielded InterStellar Radiation Field (ISRF) model described by Draine (1978) is equal to 9.9e-12 s-1 and 1.9e-9 s-1 in the Crab nebula conditions. The dependence on the visual extinction is obtained by using the Meudon Photon Dominated Region (PDR) code and corresponding analytical fits are provided. Conclusions. These data will help to produce a realistic chemical network to interpret the observations. Photodissociation of ArH+ is found to be moderate and the presence of this molecular ion is mainly dependent on the molecular fractionComment: 11 pages, 6 Figures, Accepted in Astronomy Astrophysic

Asymptotic behavior of the quadratic variation of the sum of two Hermite processes of consecutive orders

Hermite processes are self--similar processes with stationary increments which appear as limits of normalized sums of random variables with long range dependence. The Hermite process of order $1$ is fractional Brownian motion and the Hermite process of order $2$ is the Rosenblatt process. We consider here the sum of two Hermite processes of order $q\geq 1$ and $q+1$ and of different Hurst parameters. We then study its quadratic variations at different scales. This is akin to a wavelet decomposition. We study both the cases where the Hermite processes are dependent and where they are independent. In the dependent case, we show that the quadratic variation, suitably normalized, converges either to a normal or to a Rosenblatt distribution, whatever the order of the original Hermite processes

Collisional excitation of singly deuterated ammonia NH2_2D by H2_2

The availability of collisional rate coefficients with H$_2$ is a pre-requisite for interpretation of observations of molecules whose energy levels are populated under non local thermodynamical equilibrium conditions. In the current study, we present collisional rate coefficients for the NH$_2$D / para--H$_2$($J_2 = 0,2$) collisional system, for energy levels up to $J_\tau = 7_7$ ($E_u$$\sim$735 K) and for gas temperatures in the range $T = 5-300$K. The cross sections are obtained using the essentially exact close--coupling (CC) formalism at low energy and at the highest energies, we used the coupled--states (CS) approximation. For the energy levels up to $J_\tau = 4_2$ ($E_u$$\sim$215 K), the cross sections obtained through the CS formalism are scaled according to a few CC reference points. These reference points are subsequently used to estimate the accuracy of the rate coefficients for higher levels, which is mainly limited by the use of the CS formalism. Considering the current potential energy surface, the rate coefficients are thus expected to be accurate to within 5\% for the levels below $J_\tau = 4_2$, while we estimate an accuracy of 30\% for higher levels

CH2D+, the Search for the Holy Grail

CH2D+, the singly deuterated counterpart of CH3+, offers an alternative way to mediate formation of deuterated species at temperatures of several tens of K, as compared to the release of deuterated species from grains. We report a longstanding observational search for this molecular ion, whose rotational spectroscopy is not yet completely secure. We summarize the main spectroscopic properties of this molecule and discuss the chemical network leading to the formation of CH2D+, with explicit account of the ortho/para forms of H2, H3+ and CH3+. Astrochemical models support the presence of this molecular ion in moderately warm environments at a marginal level.Comment: 25 pages, 6 Figures Accepted in Journal of Physical Chemistry A. "Oka Festschrift: Celebrating 45 years of Astrochemistry

The IRAM-30m line survey of the Horsehead PDR: III. High abundance of complex (iso-)nitrile molecules in UV-illuminated gas

Complex (iso-)nitrile molecules, such as CH3CN and HC3N, are relatively easily detected in our Galaxy and in other galaxies. We constrain their chemistry through observations of two positions in the Horsehead edge: the photo-dissociation region (PDR) and the dense, cold, and UV-shielded core just behind it. We systematically searched for lines of CH3CN, HC3N, C3N, and some of their isomers in our sensitive unbiased line survey at 3, 2, and 1mm. We derived column densities and abundances through Bayesian analysis using a large velocity gradient radiative transfer model. We report the first clear detection of CH3NC at millimeter wavelength. We detected 17 lines of CH3CN at the PDR and 6 at the dense core position, and we resolved its hyperfine structure for 3 lines. We detected 4 lines of HC3N, and C3N is clearly detected at the PDR position. We computed new electron collisional rate coefficients for CH3CN, and we found that including electron excitation reduces the derived column density by 40% at the PDR position. While CH3CN is 30 times more abundant in the PDR than in the dense core, HC3N has similar abundance at both positions. The isomeric ratio CH3NC/CH3CN is 0.15+-0.02. In the case of CH3CN, pure gas phase chemistry cannot reproduce the amount of CH3CN observed in the UV-illuminated gas. We propose that CH3CN gas phase abundance is enhanced when ice mantles of grains are destroyed through photo-desorption or thermal-evaporation in PDRs, and through sputtering in shocks. (abridged)Comment: Accepted for publication in Astronomy & Astrophysic

Collisional excitation of doubly and triply deuterated ammonia ND2_2H and ND3_3 by H2_2

The availability of collisional rate coefficients is a prerequisite for an accurate interpretation of astrophysical observations, since the observed media often harbour densities where molecules are populated under non--LTE conditions. In the current study, we present calculations of rate coefficients suitable to describe the various spin isomers of multiply deuterated ammonia, namely the ND$_2$H and ND$_3$ isotopologues. These calculations are based on the most accurate NH$_3$--H$_2$ potential energy surface available, which has been modified to describe the geometrical changes induced by the nuclear substitutions. The dynamical calculations are performed within the close--coupling formalism and are carried out in order to provide rate coefficients up to a temperature of $T$ = 50K. For the various isotopologues/symmetries, we provide rate coefficients for the energy levels below $\sim$ 100 cm$^{-1}$. Subsequently, these new rate coefficients are used in astrophysical models aimed at reproducing the NH$_2$D, ND$_2$H and ND$_3$ observations previously reported towards the prestellar cores B1b and 16293E. We thus update the estimates of the corresponding column densities and find a reasonable agreement with the previous models. In particular, the ortho--to--para ratios of NH$_2$D and NHD$_2$ are found to be consistent with the statistical ratios

Ionization fraction and the enhanced sulfur chemistry in Barnard 1

Barnard B1b has revealed as one of the most interesting globules from the chemical and dynamical point of view. It presents a rich molecular chemistry characterized by large abundances of deuterated and complex molecules. Furthermore, it hosts an extremely young Class 0 object and one candidate to First Hydrostatic Core (FHSC). Our aim was to determine the cosmic ray ionization rate and the depletion factors in this extremely young star forming region. We carried out a spectral survey towards Barnard 1b as part of the IRAM Large program ASAI using the IRAM 30-m telescope at Pico Veleta (Spain). This provided a very complete inventory of neutral and ionic C-, N- and S- bearing species with, up to our knowledge, the first secure detections of the deuterated ions DCS+ and DOCO+. We used a state-of-the-art pseudo-time-dependent gas-phase chemical model to determine the value of the cosmic ray ionization rate and the depletion factors. The observational data were well fitted with $\zeta_{H_2}$ between 3E-17 s$^{-1}$ and 1E-16 s$^{-1}$. Elemental depletions were estimated to be ~10 for C and O, ~1 for N and ~25 for S. Barnard B1b presents similar depletions of C and O than those measured in pre-stellar cores. The depletion of sulfur is higher than that of C and O but not as extreme as in cold cores. In fact, it is similar to the values found in some bipolar outflows, hot cores and photon-dominated regions. Several scenarios are discussed to account for these peculiar abundances. We propose that it is the consequence of the initial conditions (important outflows and enhanced UV fields in the surroundings) and a rapid collapse (~0.1 Myr) that permits to maintain most S- and N-bearing species in gas phase to great optical depths. The interaction of the compact outflow associated with B1b-S with the surrounding material could enhance the abundances of S-bearing molecules, as well.Comment: Paper accepted in Astronomy and Astrophysics; 28 pags, 21 figure

A model for atomic and molecular interstellar gas: The Meudon PDR code

We present the revised ``Meudon'' model of Photon Dominated Region (PDR code), presently available on the web under the Gnu Public Licence at: http://aristote.obspm.fr/MIS. General organisation of the code is described down to a level that should allow most observers to use it as an interpretation tool with minimal help from our part. Two grids of models, one for low excitation diffuse clouds and one for dense highly illuminated clouds, are discussed, and some new results on PDR modelisation highlighted.Comment: accepted in ApJ sup