Performance of density functionals for calculating barrier heights of chemical reactions relevant to astrophysics

The performance of 39 different LDA, GGA, meta-GGA, and hybrid density functionals has been evaluated, for calculating forward and reverse barrier heights of 10 gas-phase reactions involving hydrogen. The reactions are all relevant to astrochemistry. Special focus is put on the applicability of DFT for calculating the rates of corresponding surface hydrogenation reactions that are relevant to the chemistry of ice-coated interstellar grains. General trends in the performance of the density functionals for reactions involving H atoms, H-2, and OH are discussed. The OH+CO reaction is shown to be a very problematic case for DFT. The best overall performance is found for the hybrid density functionals, such as MPW1K, B97-1, B97-2, and B1B95. For several reactions, the HCTH GGA functionals and the VS98 and OLAP3 meta-GGA functionals also give results that are almost as good as those of the hybrid functionals.

Molecular line emission from a protoplanetary disk irradiated externally by a nearby massive star

Star formation often occurs within or nearby stellar clusters. Irradiation by nearby massive stars can photoevaporate protoplanetary disks around young stars (so-called proplyds) which raises questions regarding the ability of planet formation to take place in these environments. We investigate the two-dimensional physical and chemical structure of a protoplanetary disk surrounding a low-mass (T Tauri) star which is irradiated by a nearby massive O-type star to determine the survivability and observability of molecules in proplyds. Compared with an isolated star-disk system, the gas temperature ranges from a factor of a few (in the disk midplane) to around two orders of magnitude (in the disk surface) higher in the irradiated disk. Although the UV flux in the outer disk, in particular, is several orders of magnitude higher, the surface density of the disk is sufficient for effective shielding of the disk midplane so that the disk remains predominantly molecular in nature. We also find that non-volatile molecules, such as HCN and H2O, are able to freeze out onto dust grains in the disk midplane so that the formation of icy planetesimals, e.g., comets, may also be possible in proplyds. We have calculated the molecular line emission from the disk assuming LTE and determined that multiple transitions of atomic carbon, CO (and isotopologues, 13CO and C18O), HCO+, CN, and HCN may be observable with ALMA, allowing characterization of the gas column density, temperature, and optical depth in proplyds at the distance of Orion (˜400 pc).

Photodetachment as destruction mechanism for CN- and C3N- anions in circumstellar envelopes

Absolute photodetachment cross sections of two anions of astrophysical importance CN- and C3N- were measured to be (1.18 +- (0.03)_stat (0.17)_sys) * 10^-17 cm^2 and (1.43 +- (0.14)_stat (0.37)_sys) * 10^-17 cm^2 respectively at the ultraviolet wavelength of 266 nm (4.66 eV). These relatively large values of the cross sections imply that photodetachment can play a major role in the destruction mechanisms of these anions particularly in photon-dominated regions. We have therefore carried out model calculations using the newly measured cross sections to investigate the abundance of these molecular anions in the cirumstellar envelope of the carbon-rich star IRC+10216. The model predicts the relative importance of the various mechanisms of formation and destruction of these species in different regions of the envelope. UV photodetachment was found to be the major destruction mechanism for both CN- and C3N- anions in those regions of the envelope, where they occur in peak abundance. It was also found that photodetachment plays a crucial role in the degradation of these anions throughout the circumstellar envelope.

Complex organic molecules in protoplanetary disks

Context. Protoplanetary disks are vital objects in star and planet formation, possessing all the material, gas and dust, which may form a planetary system orbiting the new star. Small, simple molecules have traditionally been detected in protoplanetary disks; however, in the ALMA era, we expect the molecular inventory of protoplanetary disks to significantly increase.

Aims. We investigate the synthesis of complex organic molecules (COMs) in protoplanetary disks to put constraints on the achievable chemical complexity and to predict species and transitions which may be observable with ALMA.

Methods. We have coupled a 2D steady-state physical model of a protoplanetary disk around a typical T Tauri star with a large gas-grain chemical network including COMs. We compare the resulting column densities with those derived from observations and perform ray-tracing calculations to predict line spectra. We compare the synthesised line intensities with current observations and determine those COMs which may be observable in nearby objects. We also compare the predicted grain-surface abundances with those derived from cometary comae observations.

Results. We find COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances similar to 10(-6)-10(-4) that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, similar to 10(-12)-10(-7). Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards T Tauri star-disk systems. There is poor agreement with HC3(N) lines observed towards LkCa 15 and GO Tau and we discuss possible explanations for these discrepancies. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging, even with ALMA "Full Science" capabilities. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.

Modelling the warm H(2) infrared emission of the Helix nebula cometary knots

Molecular hydrogen emission is commonly observed in planetary nebulae. Images taken in infrared H(2) emission lines show that at least part of the molecular emission is produced inside the ionized region. In the best studied case, the Helix nebula, the H(2) emission is produced inside cometary knots (CKs), comet-shaped structures believed to be clumps of dense neutral gas embedded within the ionized gas. Most of the H(2) emission of the CKs seems to be produced in a thin layer between the ionized diffuse gas and the neutral material of the knot, in a mini-photodissociation region (mini-PDR). However, PDR models published so far cannot fully explain all the characteristics of the H(2) emission of the CKs. In this work, we use the photoionization code AANGABA to study the H(2) emission of the CKs, particularly that produced in the interface H(+)/H(0) of the knot, where a significant fraction of the H(2) 1-0 S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such a region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H(2) infrared lines. We find that the temperature derived from H(2) observations, even of a single knot, will depend very strongly on the observed transitions, with much higher temperatures derived from excited levels. We also proposed that the separation between the H alpha and [N II] peak emission observed in the images of CKs may be an effect of the distance of the knot from the star, since for knots farther from the central star the [N II] line is produced closer to the border of the CK than H alpha.

Formation of SO, SO+, and S-2 by radiative association

Rate coefficients for radiative association of SO, SO+, and S-2 are estimated. For temperatures ranging from 300 to 14,000 K, the direct radiative association rate coefficients are found to vary with temperature from 1.73 x 10(-19) to 7.29 x 10(-19) cm(3) s(-1) and from 1.49 x 10(-21) to 3.70 x 10(-19) cm(3) s(-1) for S-2 and SO, respectively. The rate coefficients for formation through the inverse predissociation for S-2 are found to vary from 3.59 x 10(-18) to 1.44 x 10(-20) cm(3) s(-1). For SO+, the direct rate coefficient varies rapidly with temperature from 3.62 x 10(-27) cm(3) s(-1) at 2000 K to 2.34 x 10(-20) cm(3) s(-1) at 14,000 K. The direct radiative association rate coefficients increase with the increase in temperature, but the inverse predissociation rate coefficients decrease with the increase in temperature.

Síntese dos isótopos do monóxido de carbono no meio interestelar

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Molecular chemistry and the missing mass problem in planetary nebulae

Context. Detections of molecular lines, mainly from H-2 and CO, reveal molecular material in planetary nebulae. Observations of a variety of molecules suggest that the molecular composition in these objects differs from that found in interstellar clouds or in circumstellar envelopes. The success of the models, which are mostly devoted to explain molecular densities in specific planetary nebulae, is still partial however. Aims. The present study aims at identifying the influence of stellar and nebular properties on the molecular composition of planetary nebulae by means of chemical models. A comparison of theoretical results with those derived from the observations may provide clues to the conditions that favor the presence of a particular molecule. Methods. A self-consistent photoionization numerical code was adapted to simulate cold molecular regions beyond the ionized zone. The code was used to obtain a grid of models and the resulting column densities are compared with those inferred from observations. Results. Our models show that the inclusion of an incident flux of X-rays is required to explain the molecular composition derived for planetary nebulae. We also obtain a more accurate relation for the N(CO)/N(H-2) ratio in these objects. Molecular masses obtained by previous works in the literature were then recalculated, showing that these masses can be underestimated by up to three orders of magnitude. We conclude that the problem of the missing mass in planetary nebulae can be solved by a more accurate calculation of the molecular mass.

Effect of Nontronite Smectite Clay on the Chemical Evolution of Several Organic Molecules under Simulated Martian Surface Ultraviolet Radiation Conditions

Most of the phyllosilicates detected at the surface of Mars today are probably remnants of ancient environments that sustained long-term bodies of liquid water at the surface or subsurface and were possibly favorable for the emergence of life. Consequently, phyllosilicates have become the main mineral target in the search for organics on Mars. But are phyllosilicates efficient at preserving organic molecules under current environmental conditions at the surface of Mars? We monitored the qualitative and quantitative evolutions of glycine, urea, and adenine in interaction with the Fe3+-smectite clay nontronite, one of the most abundant phyllosilicates present at the surface of Mars, under simulated martian surface ultraviolet light (190-400 nm), mean temperature (218 +/- 2 K), and pressure (6 +/- 1 mbar) in a laboratory simulation setup. We tested organic-rich samples that were representative of the evaporation of a small, warm pond of liquid water containing a high concentration of organics. For each molecule, we observed how the nontronite influences its quantum efficiency of photodecomposition and the nature of its solid evolution products. The results reveal a pronounced photoprotective effect of nontronite on the evolution of glycine and adenine; their efficiencies of photodecomposition were reduced by a factor of 5 when mixed at a concentration of 2.6x10(-2) mol of molecules per gram of nontronite. Moreover, when the amount of nontronite in the sample of glycine was increased by a factor of 2, the gain of photoprotection was multiplied by a factor of 5. This indicates that the photoprotection provided by the nontronite is not a purely mechanical shielding effect but is also due to stabilizing interactions. No new evolution product was firmly identified, but the results obtained with urea suggest a particular reactivity in the presence of nontronite, leading to an increase of its dissociation rate. Key Words: Martian surface-Organic chemistry-Photochemistry-Astrochemistry-Nontronite-Phyllosilicates. Astrobiology 15, 221-237.

A laboratory study of water ice erosion by low-energy ions

Water ice covers the surface of various objects in the outer Solar system.Within the heliopause, surface ice is constantly bombarded and sputtered by energetic particles from the solar wind and magnetospheres. We report a laboratory investigation of the sputtering yield of water ice when irradiated at 10 K by 4 keV singly (13C+, N+, O+, Ar+) and doubly charged ions (13C2+, N2+, O2+). The experimental values for the sputtering yields are in good agreement with the prediction of a theoretical model. There is no significant difference in the yield for singly and doubly charged ions. Using these yields, we estimate the rate of water ice erosion in the outer Solar system objects due to solar wind sputtering. Temperature-programmed desorption of the ice after irradiation with 13C+ and 13C2+ demonstrated the formation of 13CO and 13CO2, with 13CO being the dominant formed species.

Investigating the cosmic-ray ionization rate in the galactic interstellar medium through observations of H3+

Observations of H3+ in the Galactic diffuse interstellar medium (ISM) have led to various surprising results, including the conclusion that the cosmic-ray ionization rate (zeta_2) is about 1 order of magnitude larger than previously thought. The present survey expands the sample of diffuse cloud sight lines with H3+ observations to 50, with detections in 21 of those. Ionization rates inferred from these detections are in the range (1.7+-1.0)x10^-16 s^-1 < zeta_2 < (10.6+-6.8)x10^-16 s^-1 with a mean value of zeta_2=(3.3+-0.4)x10^-16 s^-1. Upper limits (3sigma) derived from non-detections of H3+ are as low as zeta_2 < 0.4x10^-16 s^-1. These low upper-limits, in combination with the wide range of inferred cosmic-ray ionization rates, indicate variations in zeta_2 between different diffuse cloud sight lines. Calculations of the cosmic-ray ionization rate from theoretical cosmic-ray spectra require a large flux of low-energy (MeV) particles to reproduce values inferred from observations. Given the relatively short range of low-energy cosmic rays --- those most efficient at ionization --- the proximity of a cloud to a site of particle acceleration may set its ionization rate. Variations in zeta_2 are thus likely due to variations in the cosmic-ray spectrum at low energies resulting from the effects of particle propagation. To test this theory, H3+ was observed in sight lines passing through diffuse molecular clouds known to be interacting with the supernova remnant IC 443, a probable site of particle acceleration. Where H3+ is detected, ionization rates of zeta_2=(20+-10)x10^-16 s^-1 are inferred, higher than for any other diffuse cloud. These results support both the concept that supernova remnants act as particle accelerators, and the hypothesis that propagation effects are responsible for causing spatial variations in the cosmic-ray spectrum and ionization rate. Future observations of H3+ near other supernova remnants and in sight lines where complementary ionization tracers (OH+, H2O+, H3O+) have been observed will further our understanding of the subject.

On the temperature dependence of the rate coefficient of formation of C2+ from C + CH+

[EN] We carry out quasi-classical trajectory caculations for theC + CH+ → C2+ + H reaction on an ad hoc computed high-level ab initio potential energy surface. Thermal rate coefficients at the temperatures of relevance in cold interstellar clouds are derived and compared with the assumed, temperature-independent estimates publicly available in kinetic databases KIDA and UDfA. For a temperature of 10 K the database value overestimates by a factor of two the one obtained by us (thus improperly enhancing the destruction route of CH+ in astrochemical kinetic models) which is seen to double in the temperature range 5–300 K with a sharp increase in the first 50 K. The computed values are fitted via the popular Arrhenius–Kooij formula and best-fitting parameters α = 1:32 X 10-9 cm3s-1, β = 0:10 and γ = 2:19 K to be included in the online mentioned databases are provided. Further investigation shows that the temperature dependence of the thermal rate coefficient better conforms to the recently proposed so-called ‘deformed Arrhenius’ law by Aquilanti and Mundim.