sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Wiesenfeld, L.

Normalized to: Wiesenfeld, L.

29 article(s) in total. 203 co-authors, from 1 to 19 common article(s). Median position in authors list is 6,0.

[1]  oai:arXiv.org:1706.00732  [pdf] - 1790546
Collisional excitation of NH3 by atomic and molecular hydrogen
Comments:
Submitted: 2017-06-01, last modified: 2018-11-28
We report extensive theoretical calculations on the rotation-inversion excitation of interstellar ammonia (NH3) due to collisions with atomic and molecular hydrogen (both para- and ortho-H2). Close-coupling calculations are performed for total energies in the range 1-2000 cm-1 and rotational cross sections are obtained for all transitions among the lowest 17 and 34 rotation-inversion levels of ortho- and para-NH3, respectively. Rate coefficients are deduced for kinetic temperatures up to 200 K. Propensity rules for the three colliding partners are discussed and we also compare the new results to previous calculations for the spherically symmetrical He and para-H2 projectiles. Significant differences are found between the different sets of calculations. Finally, we test the impact of the new rate coefficients on the calibration of the ammonia thermometer. We find that the calibration curve is only weakly sensitive to the colliding partner and we confirm that the ammonia thermometer is robust.
[2]  oai:arXiv.org:1811.00766  [pdf] - 1783997
Modelling the abundance structure of isocyanic acid (HNCO) toward the low-mass solar type protostar IRAS 16293-2422
Comments:
Submitted: 2018-11-02
Isocyanic acid (HNCO), the most stable of the simplest molecules containing the four main elements essential for organic chemistry, has been observed in several astrophysical environments such as molecular clouds, star-forming regions, external galaxies and comets. In this work, we model HNCO spectral line profiles toward the low-mass solar type protostar IRAS 16293$-2$422 observed with the ALMA interferometer, the IRAM, JCMT and APEX single-dish radio telescopes, and the HIFI instrument on board the Herschel Space Observatory. In star-forming environments, the HNCO emission is not always in Local Thermodynamical Equilibrium (LTE). A non-LTE radiative transfer approach is necessary to properly interpret the line profiles, and accurate collisional rate coefficients are needed. Here, we used the RADEX package with a completely new set of collisional quenching rates between HNCO and both ortho-H$_2$ and para-H$_2$ obtained from quantum chemical calculations yielding a novel potential energy surface in the rigid rotor approximation. We find that the lines profiles toward IRAS 16293$-$2422 are very well reproduced if we assume that the HNCO emission arises from a compact, dense and hot physical component associated with the hot corino, a warm component associated with the internal part of the protostellar envelope, and a cold and more extended component associated with the outer envelope. The derived HNCO abundances from our model agree well with those computed with the Nautilus chemical code.
[3]  oai:arXiv.org:1810.00266  [pdf] - 1791087
The Dense Gas Fraction in Galactic Center Clouds
Comments: 20 pages, 8 figures, 3 tables. ApJ accepted
Submitted: 2018-09-29
We present an analysis of gas densities in the central R=300 parsecs of the Milky Way, focusing on three clouds: GCM-0.02-0.07 (the 50 km/s cloud), GCM-0.13-0.08 (the 20 km/s cloud), and GCM0.25+0.01 (the "Brick"). Densities are determined using observations of the J=(3-2), (4-3), (5-4), (10-9), (18-17), (19-18), (21-20), and (24-23) transitions of the molecule HC3N. We find evidence of at least two excitation regimes for HC3N and constrain the low-excitation component to have a density less than 10^4 cm^-3 and the high-excitation component to have a density between 10^5 and 10^6 cm^-3. This is much less than densities of 10^7 cm^-3 that are found in Sgr B2, the most actively star-forming cloud in the Galactic center. This is consistent with the requirement of a higher density threshold for star formation in the Galactic center than is typical in the Galactic disk. We are also able to constrain the column density of each component in order to determine the mass fraction of 'dense' (n>10^5 cm^-3) gas for these clouds. We find that this is ~15% for all three clouds. Applying the results of our models to ratios of the (10-9) and (3-2) line across the entire central R=300 pc, we find that the fraction of gas with n>10^4 cm^-3 increases inward of a radius of ~140 pc, consistent with the predictions of recent models for the gas dynamics in this region. Our observations show that HC3N is an excellent molecule for probing the density structure of clouds in the Galactic center.
[4]  oai:arXiv.org:1710.10437  [pdf] - 1602745
Seeds Of Life In Space (SOLIS): The organic composition diversity at 300--1000 au scale in Solar-type star forming regions
Comments: Accepted for publication on The Astrophysical Journal
Submitted: 2017-10-28
Complex organic molecules have been observed for decades in the interstellar medium. Some of them might be considered as small bricks of the macromolecules at the base of terrestrial life. It is hence particularly important to understand organic chemistry in Solar-like star forming regions. In this article, we present a new observational project: SOLIS (Seeds Of Life In Space). This is a Large Project at the IRAM-NOEMA interferometer, and its scope is to image the emission of several crucial organic molecules in a sample of Solar-like star forming regions in different evolutionary stage and environments. Here, we report the first SOLIS results, obtained from analysing the spectra of different regions of the Class 0 source NGC1333-IRAS4A, the protocluster OMC-2 FIR4, and the shock site L1157-B1. The different regions were identified based on the images of formamide (NH2CHO) and cyanodiacetylene (HC5N) lines. We discuss the observed large diversity in the molecular and organic content, both on large (3000-10000 au) and relatively small (300-1000 au) scales. Finally, we derive upper limits to the methoxy fractional abundance in the three observed regions of the same order of magnitude of that measured in few cold prestellar objects, namely ~10^-12-10^-11 with respect to H2 molecules.
[5]  oai:arXiv.org:1708.04663  [pdf] - 1587137
Seeds of Life in Space (SOLIS) III. Formamide in protostellar shocks: evidence for gas-phase formation
Comments: 7 pages, 1 table, 5 figures, A&A Letters, in press
Submitted: 2017-08-15
Context: Modern versions of the Miller-Urey experiment claim that formamide (NH$_2$CHO) could be the starting point for the formation of metabolic and genetic macromolecules. Intriguingly, formamide is indeed observed in regions forming Solar-type stars as well as in external galaxies. Aims: How NH$_2$CHO is formed has been a puzzle for decades: our goal is to contribute to the hotly debated question of whether formamide is mostly formed via gas-phase or grain surface chemistry. Methods: We used the NOEMA interferometer to image NH$_2$CHO towards the L1157-B1 blue-shifted shock, a well known interstellar laboratory, to study how the components of dust mantles and cores released into the gas phase triggers the formation of formamide. Results: We report the first spatially resolved image (size $\sim$ 9", $\sim$ 2300 AU) of formamide emission in a shocked region around a Sun-like protostar: the line profiles are blueshifted and have a FWHM $\simeq$ 5 km s$^{-1}$. A column density of $N_{\rm NH_2CHO}$ = 8 $\times$ 10$^{12}$ cm$^{-1}$, and an abundance (with respect to H-nuclei) of 4 $\times$ 10$^{-9}$ are derived. We show a spatial segregation of formamide with respect to other organic species. Our observations, coupled with a chemical modelling analysis, indicate that the formamide observed in L1157-B1 is formed by gas-phase chemical process, and not on grain surfaces as previously suggested. Conclusions: The SOLIS interferometric observations of formamide provide direct evidence that this potentially crucial brick of life is efficiently formed in the gas-phase around Sun-like protostars.
[6]  oai:arXiv.org:1707.01384  [pdf] - 1585542
SOLIS II. Carbon-chain growth in the Solar-type protocluster OMC2-FIR4
Comments: 12 pages, 8 figures, accepted for publication in A&A
Submitted: 2017-07-05
The interstellar delivery of carbon atoms locked into molecules might be one of the key ingredients for the emergence of life. Cyanopolyynes are carbon chains delimited at their two extremities by an atom of hydrogen and a cyano group, so that they might be excellent reservoirs of carbon. The simplest member, HC3N, is ubiquitous in the galactic interstellar medium and found also in external galaxies. Thus, understanding the growth of cyanopolyynes in regions forming stars similar to our Sun, and what affects it, is particularly relevant. In the framework of the IRAM/NOEMA Large Program SOLIS (Seeds Of Life In Space), we have obtained a map of two cyanopolyynes, HC3N and HC5N, in the protocluster OMC2-FIR4. Because our Sun is thought to be born in a rich cluster, OMC2-FIR4 is one of the closest and best known representatives of the environment in which the Sun may have been born. We find a HC3N/HC5N abundance ratio across the source in the range ~ 1 - 30, with the smallest values (< 10) in FIR5 and in the Eastern region of FIR4. The ratios < 10 can be reproduced by chemical models only if: (1) the cosmic-ray ionisation rate $\zeta$ is ~ $4 \times 10^{-14}$ s$^{-1}$; (2) the gaseous elemental ratio C/O is close to unity; (3) O and C are largely depleted. The large $\zeta$ is comparable to that measured in FIR4 by previous works and was interpreted as due to a flux of energetic (> 10 MeV) particles from embedded sources. We suggest that these sources could lie East of FIR4 and FIR5. A temperature gradient across FIR4, with T decreasing by about 10 K, could also explain the observed change in the HC3N/HC5N line ratio, without the need of a cosmic ray ionisation rate gradient. However, even in this case, a high constant cosmic-ray ionisation rate (of the order of $10^{-14}$ s$^{-1}$) is necessary to reproduce the observations.
[7]  oai:arXiv.org:1604.05602  [pdf] - 1557673
Deuteration of ammonia in the starless core Ophiuchus/H-MM1
Comments: to appear in Astronomy & Astrophysics
Submitted: 2016-04-19, last modified: 2016-12-22
Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.4, NHD2/NH2D=0.2, and ND3/NHD2=0.06 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but predicts uniformly too low ortho/para-NH2D, and too large ortho/para-NHD2 ratios. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected. The present results suggest that while the deuteration of ammonia changes with physical conditions and time, the nuclear spin ratios of ammonia isotopologues do not probe the evolutionary stage of a cloud.
[8]  oai:arXiv.org:1610.00880  [pdf] - 1531856
History of the solar-type protostar IRAS16293-2422 as told by the cyanopolyynes
Comments: Accepted for publication in A&A
Submitted: 2016-10-04
Cyanopolyynes are chains of carbon atoms with an atom of hydrogen and a CN group on either side. They are detected almost everywhere in the ISM, as well as in comets. In the past, they have been used to constrain the age of some molecular clouds, since their abundance is predicted to be a strong function of time. We present an extensive study of the cyanopolyynes distribution in the solar-type protostar IRAS16293-2422 based on TIMASSS IRAM-30m observations. The goals are (i) to obtain a census of the cyanopolyynes in this source and of their isotopologues; (ii) to derive how their abundance varies across the protostar envelope; and (iii) to obtain constraints on the history of IRAS16293-2422. We detect several lines from HC3N and HC5N, and report the first detection of DC3N, in a solar-type protostar. We found that the HC3N abundance is roughly constant (~1.3x10^(-11)) in the outer cold envelope of IRAS16293-2422, and it increases by about a factor 100 in the inner region where Tdust>80K. The HC5N has an abundance similar to HC3N in the outer envelope and about a factor of ten lower in the inner region. The HC3N abundance derived in the inner region, and where the increase occurs, also provide strong constraints on the time taken for the dust to warm up to 80K, which has to be shorter than ~10^3-10^4yr. Finally, the cyanoacetylene deuteration is about 50\% in the outer envelope and <5$\% in the warm inner region. The relatively low deuteration in the warm region suggests that we are witnessing a fossil of the HC3N abundantly formed in the tenuous phase of the pre-collapse and then frozen into the grain mantles at a later phase. The accurate analysis of the cyanopolyynes in IRAS16293-2422 unveils an important part of its past story. It tells us that IRAS16293-2422 underwent a relatively fast (<10^5yr) collapse and a very fast (<10^3-10^4yr) warming up of the cold material to 80K.
[9]  oai:arXiv.org:1610.00438  [pdf] - 1490039
Theory of Gas Phase Scattering and Reactivity for Astrochemistry
Comments: 24 pages, 7 figures. White paper reporting a workshop held in Garching, Nov 23-Dec 4, 2015
Submitted: 2016-10-03
Because of the very peculiar conditions of chemistry in many astrophysical gases (low densities, mostly low temperatures, kinetics-dominated chemical evolution), great efforts have been devoted to study molecular signatures and chemical evolution. While experiments are being performed in many laboratories, it appears that the efforts directed towards theoretical works are not as strong. This report deals with the present status of chemical physics/physical chemistry theory, for the qualitative and quantitative understanding of kinetics of molecular scattering, being it reactive or inelastic. By gathering several types of expertise, from applied mathematics to physical chemistry, dialog is made possible, as a step towards new and more adapted theoretical frameworks, capable of meeting the theoretical, methodological and numerical challenges of kinetics-dominated gas phase chemistry in astrophysical environments. A state of the art panorama is presented, alongside present-day strengths and shortcomings. However, coverage is not complete, being limited in this report to actual attendance of the workshop. Some paths towards relevant progress are proposed.
[10]  oai:arXiv.org:1605.03786  [pdf] - 1422348
Collisional excitation of HC3N by para- and ortho-H2
Comments: 8 pages, 2 tables, 4 figures, accepted for publication in MNRAS
Submitted: 2016-05-12
New calculations for rotational excitation of cyanoacetylene by collisions with hydrogen molecules are performed to include the lowest 38 rotational levels of HC3N and kinetic temperatures to 300 K. Calculations are based on the interaction potential of Wernli et al. A&A, 464, 1147 (2007) whose accuracy is checked against spectroscopic measurements of the HC3N-H2 complex. The quantum coupled-channel approach is employed and complemented by quasi-classical trajectory calculations. Rate coefficients for ortho-H2 are provided for the first time. Hyperfine resolved rate coefficients are also deduced. Collisional propensity rules are discussed and comparisons between quantum and classical rate coefficients are presented. This collisional data should prove useful in interpreting HC3N observations in the cold and warm ISM, as well as in protoplanetary disks.
[11]  oai:arXiv.org:1601.02386  [pdf] - 1359192
Collisional excitation of doubly and triply deuterated ammonia ND$_2$H and ND$_3$ by H$_2$
Comments:
Submitted: 2016-01-11
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.
[12]  oai:arXiv.org:1503.04112  [pdf] - 1129041
Chemical composition of the circumstellar disk around AB Aurigae
Comments:
Submitted: 2015-03-13, last modified: 2015-05-11
Aims. Our goal is to determine the molecular composition of the circumstellar disk around AB Aurigae (hereafter, AB Aur). AB Aur is a prototypical Herbig Ae star and the understanding of its disk chemistry is of paramount importance to understand the chemical evolution of the gas in warm disks. Methods. We used the IRAM 30-m telescope to perform a sensitive search for molecular lines in AB Aur as part of the IRAM Large program ASAI (A Chemical Survey of Sun-like Star-forming Regions). These data were complemented with interferometric observations of the HCO+ 1-0 and C17O 1-0 lines using the IRAM Plateau de Bure Interferometer (PdBI). Single-dish and interferometric data were used to constrain chemical models. Results. Throughout the survey, several lines of CO and its isotopologues, HCO+, H2CO, HCN, CN and CS, were detected. In addition, we detected the SO 54-33 and 56-45 lines, confirming the previous tentative detection. Comparing to other T Tauri's and Herbig Ae disks, AB Aur presents low HCN 3-2/HCO+ 3-2 and CN 2-1/HCN 3-2 line intensity ratios, similar to other transition disks. AB Aur is the only protoplanetary disk detected in SO thus far. Conclusions. We modeled the line profiles using a chemical model and a radiative transfer 3D code. Our model assumes a flared disk in hydrostatic equilibrium. The best agreement with observations was obtained for a disk with a mass of 0.01 Msun , Rin=110 AU, Rout=550 AU, a surface density radial index of 1.5 and an inclination of 27 deg. The intensities and line profiles were reproduced within a factor of 2 for most lines. This agreement is reasonable taking into account the simplicity of our model that neglects any structure within the disk. However, the HCN 3-2 and CN 2-1 line intensities were predicted more intense by a factor of >10. We discuss several scenarios to explain this discrepancy.
[13]  oai:arXiv.org:1408.5757  [pdf] - 1216493
Collisional excitation of singly deuterated ammonia NH$_2$D by H$_2$
Comments:
Submitted: 2014-08-25, last modified: 2014-08-26
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.
[14]  oai:arXiv.org:1401.1136  [pdf] - 1202568
Weak maser emission of methyl formate toward Sagittarius B2(N) in the Green Bank Telescope PRIMOS Survey
Comments: 33 pages, 9 figures, accepted in The Astrophysical Journal (january 4 2014)
Submitted: 2014-01-06
A non-LTE radiative transfer treatment of cis-methyl formate (HCOOCH3) rotational lines is presented for the first time using a set of theoretical collisional rate coefficients. These coefficients have been computed in the temperature range 5-30 K by combining coupled-channel scattering calculations with a high accuracy potential energy surface for HCOOCH3-He. The results are compared to observations toward the Sagittarius B2(N) molecular cloud using the publicly available PRIMOS survey from the Green Bank Telescope. A total of 49 low-lying transitions of methyl formate, with upper levels below 25 K, are identified. These lines are found to probe a presumably cold (~30 K), moderately dense (~1e4 cm-3) and extended region surrounding Sgr B2(N). The derived column density of ~4e14 cm-2 is only a factor of ~10 larger than the column density of the trans conformer in the same source. Provided that the two conformers have the same spatial distribution, this result suggests that strongly non-equilibrium processes must be involved in their synthesis. Finally, our calculations show that all detected emission lines with a frequency below 30 GHz are (collisionally pumped) weak masers amplifying the continuum of Sgr B2(N). This result demonstrates the importance and generality of non-LTE effects in the rotational spectra of complex organic molecules at centimetre wavelengths.
[15]  oai:arXiv.org:1312.0282  [pdf] - 1944169
CII in the Interstellar Medium: Excitation by H2 Revisited
Comments:
Submitted: 2013-12-01
C$^+$ is a critical constituent of many regions of the interstellar medium, as it can be a major reservoir of carbon and, under a wide range of conditions, the dominant gas coolant. Emission from its 158$\mu$m fine structure line is used to trace the structure of photon dominated regions in the Milky Way and is often employed as a measure of the star formation rate in external galaxies. Under most conditions, the emission from the single [CII] line is proportional to the collisional excitation rate coefficient. We here used improved calculations of the deexcitation rate of [CII] by collisions with H$_2$ to calculate more accurate expressions for interstellar C$^+$ fine structure emission, its critical density, and its cooling rate. The collision rates in the new quantum calculation are $\sim$ 25% larger than those previously available, and narrow the difference between rates for excitation by atomic and molecular hydrogen. This results in [CII] excitation being quasi-independent of the molecular fraction and thus dependent only on the total hydrogen particle density. A convenient expression for the cooling rate at temperatures between 20 K and 400 K, assuming an LTE H$_2$ ortho to para ration is $\Lambda ({\rm LTE~OPR}) = \left(11.5 + 4.0\,e^{-100\,\mathrm K/T^{\rm kin}}\right)\;e^{-91.25\,\mathrm K/T^{\rm kin}}\,n ({\rm C}^{+})\,n({\rm H}_2)\times 10^{-24}\;{\rm ergs}~{\rm cm}^{-3}~{\rm s}^{-1}$. The present work should allow more accurate and convenient analysis of the [\CII] line emission and its cooling.
[16]  oai:arXiv.org:1311.2840  [pdf] - 1180659
The CHESS survey of the L1157-B1 bow-shock: high and low excitation water vapor
Comments: Accepted for publication in A&A. 12 pages, 9 figures, 4 tables
Submitted: 2013-11-12
Molecular outflows powered by young protostars strongly affect the kinematics and chemistry of the natal molecular cloud through strong shocks resulting in substantial modifications of the abundance of several species. As part of the "Chemical Herschel Surveys of Star forming regions" guaranteed time key program, we aim at investigating the physical and chemical conditions of H20 in the brightest shock region B1 of the L1157 molecular outflow. We observed several ortho- and para-H2O transitions using HIFI and PACS instruments on board Herschel, providing a detailed picture of the kinematics and spatial distribution of the gas. We performed a LVG analysis to derive the physical conditions of H2O shocked material, and ultimately obtain its abundance. We detected 13 H2O lines probing a wide range of excitation conditions. PACS maps reveal that H2O traces weak and extended emission associated with the outflow identified also with HIFI in the o-H2O line at 556.9 GHz, and a compact (~10") bright, higher-excitation region. The LVG analysis of H2O lines in the bow-shock show the presence of two gas components with different excitation conditions: a warm (Tkin~200-300 K) and dense (n(H2)~(1-3)x10^6 cm-3) component with an assumed extent of 10" and a compact (~2"-5") and hot, tenuous (Tkin~900-1400 K, n(H2)~10^3-10^4 cm-3) gas component, which is needed to account for the line fluxes of high Eu transitions. The fractional abundance of the warm and hot H2O gas components is estimated to be (0.7-2)x10^{-6} and (1-3)x10^{-4}, respectively. Finally, we identified an additional component in absorption in the HIFI spectra of H2O lines connecting with the ground state level, probably arising from the photodesorption of icy mantles of a water-enriched layer at the edges of the cloud.
[17]  oai:arXiv.org:1304.4804  [pdf] - 1166038
Rotational quenching of H2CO by molecular hydrogen: cross-sections, rates, pressure broadening
Comments:
Submitted: 2013-04-17, last modified: 2013-04-24
We compute the rotational quenching rates of the first 81 rotational levels of ortho- and para-H2CO in collision with ortho- and para-H2, for a temperature range of 10-300 K. We make use of the quantum close-coupling and coupled-states scattering methods combined with the high accuracy potential energy surface of Troscompt et al. (2009a). Rates are significantly different from the scaled rates of H2CO in collision with He; consequently, critical densities are noticeably lower. We compare a full close- coupling computation of pressure broadening cross sections with experimental data and show that our results are compatible with the low temperature measurements of Mengel & De Lucia (2000), for a spin temperature of H2 around 50 K.
[18]  oai:arXiv.org:1211.0514  [pdf] - 1157471
Water ice deuteration: a tracer of the chemical history of protostars
Comments: Accepted for publication to A&A. 24 pages, 10 figures
Submitted: 2012-11-02, last modified: 2012-12-13
Context. Millimetric observations have measured high degrees of molecular deuteration in several species seen around low-mass protostars. The Herschel Space Telescope, launched in 2009, is now providing new measures of the deuterium fractionation of water, the main constituent of interstellar ices. Aims. We aim at theoretically studying the formation and the deuteration of water, which is believed to be formed on interstellar grain surfaces in molecular clouds. Methods. We used our gas-grain astrochemical model GRAINOBLE, which considers the multilayer formation of interstellar ices. We varied several input parameters to study their impact on water deuteration. We included the treatment of ortho- and para-states of key species, including H2, which affects the deuterium fractionation of all molecules. The model also includes relevant laboratory and theoretical works on the water formation and deuteration on grain surfaces. In particular, we computed the transmission probabilities of surface reactions using the Eckart model, and we considered ice photodissociation following molecular dynamics simulations. Results. The use of a multilayer approach allowed us to study the influence of various parameters on the abundance and the deuteration of water. Deuteration of water is found to be very sensitive to the ortho-to-para ratio of H2 and to the total density, but it also depends on the gas/grain temperatures and the visual extinction of the cloud. Since the deuteration is very sensitive to the physical conditions, the comparison with sub-millimetric observation towards the low-mass protostar IRAS 16293 allows us to suggest that water ice is formed together with CO2 in molecular clouds with limited density, whilst formaldehyde and methanol are mainly formed in a later phase, where the condensation becomes denser and colder.
[19]  oai:arXiv.org:1210.7654  [pdf] - 1152536
Radiation thermo-chemical models of protoplanetary discs. IV Modelling CO ro-vibrational emission from Herbig Ae discs
Comments: 22 pages, 27 figures, accepted for publication in A&A
Submitted: 2012-10-29
The carbon monoxide rovibrational emission from discs around Herbig Ae stars and T Tauri stars with strong ultraviolet emissions suggests that fluorescence pumping from the ground X1 Sigma+ to the electronic A1 Pi state of CO should be taken into account in disc models. We implemented a CO model molecule that includes up to 50 rotational levels within nine vibrational levels for the ground and A excited states in the radiative photochemical code ProDiMo. We took CO collisions with hydrogen molecules, hydrogen atoms, helium, and electrons into account. We estimated the missing collision rates using standard scaling laws and discussed their limitations. UV fluorescence and IR pumping impact on the population of ro-vibrational v > 1 levels. The v = 1 rotational levels are populated at rotational temperatures between the radiation temperature around 4.6 micron and the gas kinetic temperature. The UV pumping efficiency increases with decreasing disc mass. The consequence is that the vibrational temperatures, which measure the relative populations between the vibrational levels, are higher than the disc gas kinetic temperatures (suprathermal population). Rotational temperatures from fundamental transitions derived using optically thick 12CO lines do not reflect the gas kinetic temperature. CO pure rotational levels with energies lower than 1000 K are populated in LTE but are sensitive to a number of vibrational levels included in the model. The 12CO pure rotational lines are highly optically thick for transition from levels up to Eupper=2000 K. (abridged)
[20]  oai:arXiv.org:1201.1785  [pdf] - 1092809
A study of deuterated water in the low-mass protostar IRAS16293-2422
Comments: 12 pages, 11 figures
Submitted: 2012-01-09, last modified: 2012-02-16
The HDO/H2O ratio is a powerful diagnostic to understand the evolution of water from the first stages of star formation to the formation of planets and comets. Our aim is to determine precisely the abundance distribution of HDO towards the low-mass protostar IRAS16293-2422 and learn more about the water formation mechanisms by determining the HDO/H2O abundance ratio. A spectral survey of the source IRAS16293-2422 was carried out in the framework of the CHESS Herschel Key program with the HIFI instrument, allowing the detection of numerous HDO lines. Other transitions have been observed previously with ground-based telescopes. The spherical Monte Carlo radiative transfer code RATRAN was used to reproduce the observed line profiles of HDO by assuming an abundance jump. To determine the H2O abundance throughout the envelope, a similar study was made of the H2-18O observed lines, as the H2O main isotope lines are contaminated by the outflows. We derive an inner HDO abundance of 1.7e-7 and an outer HDO abundance of 8e-11. To reproduce the HDO absorption lines, it is necessary to add an absorbing layer in front of the envelope. It may correspond to a water-rich layer created by the photodesorption of the ices at the edges of the molecular cloud. The HDO/H2O ratio is ~1.4-5.8% in the hot corino whereas it is ~0.2-2.2% in the outer envelope. It is estimated at ~4.8% in the added absorbing layer. Although it is clearly higher than the cosmic D/H abundance, the HDO/H2O ratio remains lower than the D/H ratio derived for other deuterated molecules observed in the same source. The similarity of the ratios derived in the hot corino and in the added absorbing layer suggests that water formed before the gravitational collapse of the protostar, contrary to formaldehyde and methanol, which formed later once the CO molecules had depleted on the grains.
[21]  oai:arXiv.org:1201.5887  [pdf] - 1093310
A KInetic Database for Astrochemistry (KIDA)
Comments: Accepted for publication in ApJS
Submitted: 2012-01-27
We present a novel chemical database for gas-phase astrochemistry. Named the KInetic Database for Astrochemistry (KIDA), this database consists of gas-phase reactions with rate coefficients and uncertainties that will be vetted to the greatest extent possible. Submissions of measured and calculated rate coefficients are welcome, and will be studied by experts before inclusion into the database. Besides providing kinetic information for the interstellar medium, KIDA is planned to contain such data for planetary atmospheres and for circumstellar envelopes. Each year, a subset of the reactions in the database (kida.uva) will be provided as a network for the simulation of the chemistry of dense interstellar clouds with temperatures between 10 K and 300 K. We also provide a code, named Nahoon, to study the time-dependent gas-phase chemistry of 0D and 1D interstellar sources.
[22]  oai:arXiv.org:1110.3948  [pdf] - 1084933
First detection of Hydrogen Chloride towards protostellar shocks
Comments: ApJ, in press
Submitted: 2011-10-18
We present the first detection of hydrogen chlorine in a protostellar shock, by observing the fundamental transition at 626 GHz with the Herschel HIFI spectrometer. We detected two of the three hyperfine lines, from which we derived a line opacity < 1. Using a non-LTE LVG code, we constrained the HCl column density, temperature and density of the emitting gas. The hypothesis that the emission originates in the molecular cloud is ruled out, as it would imply a too dense gas. Conversely, assuming that the emission originates in the 10"-15" size shocked gas previously observed at the IRAM PdB interferometer, we obtain: N(HCl)=0.7-2 x 10(13) cm-2, temperature > 15 K and density > 3 x 10(5) cm-3}. Combining with the Herschel HIFI CO(5-4) observations allows to further constrain the gas density and temperature, 10(5)-10(6) cm-3 and 120-250 K, as well as the HCl column density, 2 x 10(13) cm-2, and, finally, abundance: 3-6 x 10(-9). The estimated HCl abundance is consistent with that previously observed in low- and high- mass protostars. This puzzling result in the L1157-B1 shock, where species from volatile and refractory grains components are enhanced, suggests either that HCl is not the main reservoir of chlorine in the gas phase, against previous chemical models predictions, or that the elemental chlorine abundance is low in L1157-B1. Astrochemical modelling suggests that HCl is in fact formed in the gas phase, at low temperatures, prior to the occurance of the shock, and that the latter does not enhance its abundance.
[23]  oai:arXiv.org:1104.2715  [pdf] - 1053362
State-to-State Differential and Relative Integral Cross Sections for Rotationally Inelastic Scattering of H2O by Hydrogen
Comments: 46 pages
Submitted: 2011-04-14
State-to-state differential cross sections (DCSs) for rotationally inelastic scattering of H2O by H2 have been measured at 71.2 meV (574 cm-1) and 44.8 meV (361 cm-1) collision energy using crossed molecular beams combined with velocity map imaging. A molecular beam containing variable compositions of the (J = 0, 1, 2) rotational states of hydrogen collides with a molecular beam of argon seeded with water vapor that is cooled by supersonic expansion to its lowest para or ortho rotational levels (JKaKc= 000 and 101, respectively). Angular speed distributions of fully specified rotationally excited final states are obtained using velocity map imaging. Relative integral cross sections are obtained by integrating the DCSs taken with the same experimental conditions. Experimental state-specific DCSs are compared with predictions from fully quantum scattering calculations on the most complete H2O-H2 potential energy surface. Comparison of relative total cross sections and state-specific DCSs show excellent agreement with theory in almost all details
[24]  oai:arXiv.org:1012.5785  [pdf] - 1042797
Collisional excitation of doubly deuterated ammonia NHD2 by para-H2
Comments: 5 pages, 2 figures
Submitted: 2010-12-28
Collisional de-excitation rates of partially deuterated molecules are different from the fully hydrogenated species because of lowering of symmetry. We compute the collisional (de)excitation rates of ND2H by ground state para-H2, extending the previous results for He- lium. We describe the changes in the potential energy surface of NH3- H2 involved by the pres- ence of two deuterium nuclei. Cross sections are calculated within the full close-coupling ap- proach and augmented with coupled-state calculations. Collisional rate coefficients are given between 5 and 35 K, a range of temperatures which is relevant to cold interstellar conditions. We find that the collisional rates of ND2H by H2 are about one order of magnitude higher than those obtained with Helium as perturber. These results are essential to radiative transfer modelling and will allow to interpret the millimeter and submillimeter detections of ND2H with better constraints than previously.
[25]  oai:arXiv.org:1007.4410  [pdf] - 1033907
Ortho-to-para ratio of interstellar heavy water
Comments: 5 pages, 5 figures, accepted the A&A HIFI Special Issue as a letter
Submitted: 2010-07-26
Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.
[26]  oai:arXiv.org:0906.4468  [pdf] - 1002676
On the robustness of the ammonia thermometer
Comments: Accepted for publication in the MNRAS
Submitted: 2009-06-24
Ammonia inversion lines are often used as probes of the physical conditions in the dense ISM. The excitation temperature between the first two para metastable (rotational) levels is an excellent probe of the gas kinetic temperature. However, the calibration of this ammonia thermometer depends on the accuracy of the collisional rates with H2. Here we present new collisional rates for ortho-NH3 and para-NH3 colliding with para-H2 (J=0) and we investigate the effects of these new rates on the excitation of ammonia. Scattering calculations employ a new, high accuracy, potential energy surface computed at the coupled-cluster CCSD(T) level with a basis set extrapolation procedure. Rates are obtained for all transitions involving ammonia levels with J <= 3 and for kinetic temperatures in the range 5-100 K. We find that the calibration curve of the ammonia thermometer -- which relates the observed excitation temperature between the first two para metastable levels to the gas kinetic temperature -- does not change significantly when these new rates are used. Thus, the calibration of ammonia thermometer appears to be robust. Effects of the new rates on the excitation temperature of inversion and rotation-inversion transitions are also found to be small.
[27]  oai:arXiv.org:0708.0345  [pdf] - 3664
Quasi-classical rate coefficient calculations for the rotational (de)excitation of H2O by H2
Comments: 8 pages, 2 figures, 1 table (the online material (4 tables) can be obtained upon request to afaure@obs.ujf-grenoble.fr)
Submitted: 2007-08-02
The interpretation of water line emission from existing observations and future HIFI/Herschel data requires a detailed knowledge of collisional rate coefficients. Among all relevant collisional mechanisms, the rotational (de)excitation of H2O by H2 molecules is the process of most interest in interstellar space. To determine rate coefficients for rotational de-excitation among the lowest 45 para and 45 ortho rotational levels of H2O colliding with both para and ortho-H2 in the temperature range 20-2000 K. Rate coefficients are calculated on a recent high-accuracy H2O-H2 potential energy surface using quasi-classical trajectory calculations. Trajectories are sampled by a canonical Monte-Carlo procedure. H2 molecules are assumed to be rotationally thermalized at the kinetic temperature. By comparison with quantum calculations available for low lying levels, classical rates are found to be accurate within a factor of 1-3 for the dominant transitions, that is those with rates larger than a few 10^{-12}cm^{3}s^{-1}. Large velocity gradient modelling shows that the new rates have a significant impact on emission line fluxes and that they should be adopted in any detailed population model of water in warm and hot environments.
[28]  oai:arXiv.org:physics/0611258  [pdf] - 119196
Rotational Excitation of HC_3N by H_2 and He at low temperatures
Comments:
Submitted: 2006-11-27
Rates for rotational excitation of HC3N by collisions with He atoms and H2 molecules are computed for kinetic temperatures in the range 5-20K and 5-100K, respectively. These rates are obtained from extensive quantum and quasi-classical calculations using new accurate potential energy surfaces (PES).
[29]  oai:arXiv.org:astro-ph/0601384  [pdf] - 79242
Improved low-temperature rate constants for rotational excitation of CO by H_2
Comments: 8 pages, including 2 "online material" pages, 2 figures
Submitted: 2006-01-18
Cross sections for the rotational (de)excitation of CO by ground state para- and ortho-H_2 are obtained using quantum scattering calculations for collision energies between 1 and 520 cm^{-1}. A new CO-H_2 potential energy surface is employed and its quality is assessed by comparison with explicitly correlated CCSD(T)-R12 calculations. Rate constants for rotational levels of CO up to 5 and temperatures in the range 5-70 K are deduced. The new potential is found to have a strong influence on the resonance structure of the cross sections at very low collision energies. As a result, the present rates at 10 K differ by up to 50% with those obtained by \citet{flower01} on a previous, less accurate, potential energy surface.