Source author record

Valentine Wakelam

Valentine Wakelam appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2022arXiv

AB Aur, a Rosetta stone for studies of planet formation (II): H$_2$S detection and sulfur budget

The sulfur abundance is poorly known in most environments. Yet, deriving the sulfur abundance is key to understanding the evolution of the chemistry from molecular clouds to planetary atmospheres. We present observations of H$_2$S 110-101 at 168.763 GHz toward the Herbig Ae star AB Aur. We aim to study the abundance of sulfuretted species toward AB Aur and to constrain how different species and phases contribute to the sulfur budget. We present new NOrthern Extended Millimeter Array (NOEMA) interferometric observations of the continuum and H$_2$S 110-101 line at 168.763 GHz toward AB Aur. We derived radial and azimuthal profiles and used them to compare the geometrical distribution of different species in the disk. Assuming local thermodynamical equilibrium (LTE), we derived column density and abundance maps for H$_2$S, and we further used Nautilus to produce a more detailed model of the chemical abundances at different heights over the mid-plane at a distance of r=200 au. We have resolved H$_2$S emission in the AB Aur protoplanetary disk. The emission comes from a ring extending from 0.67 (109 au) to 1.69 (275 au). Under simple assumptions, we derived an abundance of (3.1$\pm$0.8)$\times$10$\rm ^{-10}$ with respect to H nuclei, which we compare with Nautilus models to deepen our understanding of the sulfur chemistry in protoplanetary disks. Chemical models indicate that H$-2$S is an important sulfur carrier in the solid and gas phase. We also find an important transition at a height of 12 au, where the sulfur budget moves from being dominated by ice species to being dominated by gas species. Studying sulfuretted species in detail in the different phases of the interstellar medium is key to solving the issue.

preprint2022arXiv

An Experimental and Theoretical Investigation of the Gas-Phase C(3P) + N2O Reaction. Low Temperature Rate Constants and Astrochemical Implications

The reaction between atomic carbon in its ground electronic state, C(3P), and nitrous oxide, N2O, has been studied below room temperature due to its potential importance for astrochemistry, with both species considered to be present at high abundance levels in a range of interstellar environments. On the experimental side, we measured rate constants for this reaction over the 50-296 K range using a continuous supersonic flow reactor. C(3P) atoms were generated by the pulsed photolysis of carbon tetrabromide at 266 nm and were detected by pulsed laser induced fluorescence at 115.8 nm. Additional measurements allowing the major product channels to be elucidated were also performed. On the theoretical side, statistical rate theory was used to calculate low temperature rate constants. These calculations employed the results of new electronic structure calculations of the 3A" potential energy surface of CNNO and provided a basis to extrapolate the measured rate constants to lower temperatures and pressures. The rate constant was found to increase monotonically as the temperature falls, reaching a value of k(C(3P)+N2O)(50 K) = (7.9 +- 0.8) x 10-11 cm3 s-1 at 50 K. As current astrochemical models do not include the C + N2O reaction, we tested the influence of this process on interstellar N2O and other related species using a gas-grain model of dense interstellar clouds. These simulations predict that N2O abundances decrease significantly at intermediate times (10^3 - 10^5 years) when gas-phase C(3P) abundances are high.

preprint2022arXiv

Investigating the hot molecular core, G10.47+0.03: A pit of nitrogen-bearing complex organic molecules

Recent observations have shown that Nitrogen-bearing complex organic species are present in large quantities in star-forming regions. Thus, investigating the N-bearing species in a hot molecular core, such as G10.47+0.03, is crucial to understanding the molecular complexity in star-forming regions. They also allow us to investigate the chemical and physical processes that determine the many phases during the structural and chemical evolution of the source in star-forming regions. The aim of this study is to investigate the spatial distribution and the chemical evolution states of N-bearing complex organic molecules in the hot core G10.47+0.03. We used the ALMA archival data of the hot molecular core G10.47+0.03. The extracted spectra were analyzed assuming LTE. Furthermore, robust methods such as MCMC and rotational diagram methods are implemented for molecules for which multiple transitions were identified to constrain the temperature and column density. Finally, we used the Nautilus gas-grain code to simulate the nitrogen chemistry in the hot molecular core. We carried out both 0D and 1D simulations of the source and compared with observational results. We report various transitions of nitrogen-bearing species (NH2CN, HC3N, HC5N, C2H3CN, C2H5CN, and H2NCH2CN) together with some of their isotopologues and isomers. Besides this, we also report the identification of CH3CCH and one of its isotopologues. The emissions originating from vinyl cyanide, ethyl cyanide, cyanoacetylene, and cyanamide are compact, which could be explained by our astrochemical modeling. Our 0D model shows that the chemistry of certain N-bearing molecules can be very sensitive to initial local conditions such as density or dust temperature. In our 1D model, simulated higher abundances of species such as HCN, HC3N, and HC5N toward the inner shells of the source confirm the observational findings.

preprint2022arXiv

The Si+SO$_2$ collision and an extended network of neutral-neutral reactions between silicon and sulphur bearing species

The Si+SO$_2$ reaction is investigated to verify its impact on the abundances of molecules with astrochemical interest, such as SiS, SiO, SO and others. According to our results Si($^3$P) and SO$_2$ react barrierlessly yielding only the monoxides SO and SiO as products. No favourable pathway has been found leading to other products, and this reaction should not contribute to SiS abundance. Furthermore, it is predicted that SiS is stable in collisions with O$_2$, and that S($^3$P)+SiO$_2$ and O($^3$P)+OSiS will also produce SO+SiO. Using these results and gathering further experimental and computational data from the literature, we provide an extended network of neutral-neutral reactions involving Si- and S-bearing molecules. The effects of these reactions were examined in a protostellar shock model, using the Nautilus gas-grain code. This consisted in simulating the physicochemical conditions of a shocked gas evolving from $i.$ primeval cold core, $ii.$ the shock region itself, $iii.$ and finally the gas bulk conditions after the passage of the shock. Emphasising on the cloud ages and including systematically these chemical reactions, we found that [SiS/H$_2$] can be of the order of $\sim$ 10$^{-8}$ in shocks that evolves from clouds of $t=1\times 10^6$ yr, whose values are mostly affected by the SiS+O $\longrightarrow$SiO+S reaction. Perspectives on further models along with observations are discussed in the context of sources harbouring molecular outflows.

preprint2022arXiv

Thermal desorption of interstellar ices. A review on the controlling parameters and their implications fromsnowlines to chemical complexity

The evolution of star-forming regions and their thermal balance are strongly influenced by their chemical composition, that, in turn, is determined by the physico-chemical processes that govern the transition between the gas phase and the solid state, specifically icy dust grains (e.g., particles adsorption and desorption). Gas-grain and grain-gas transitions as well as formation and sublimation of interstellar ices are thus essential elements of understanding astrophysical observations of cold environments (e.g., pre-stellar cores) where unexpected amounts of a large variety of chemical species have been observed in the gas phase. Adsorbed atoms and molecules also undergo chemical reactions which are not efficient in the gas phase. Therefore, the parameterization of the physical properties of atoms and molecules interacting with dust grain particles is clearly a key aspect to interpret astronomical observations and to build realistic and predictive astrochemical models. In this consensus evaluation, we focus on parameters controlling the thermal desorption of ices and how these determine pathways towards molecular complexity and define the location of snowlines, which ultimately influence the planet formation process. We review different crucial aspects of desorption parameters both from a theoretical and experimental point of view. We critically assess the desorption parameters commonly used in the astrochemical community for astrophysical relevant species and provide tables with recommended values. In addition, we show that a non-trivial determination of the pre-exponential factor nu using the Transition State Theory can affect the binding energy value. Finally, we conclude this work by discussing the limitations of theoretical and experimental approaches currently used to determine the desorption properties with suggestions for future improvements.

preprint2020arXiv

Constraints of the formation and abundances of methyl carbamate, a glycine isomer, in hot corinos

Methyl carbamate CH$_3$OC(O)NH$_2$ is an isomer of glycine. Quantum chemical analyses show that methyl carbamate is more stable isomer than glycine. Because of this, there could be a higher chance for methyl carbamte to exist in the interstellar medium as compared to glycine. Despite immense searches, till now glycine has not been detected in the ISM, therefore it is worthwhile to search its isomer methyl carbamate. In this paper, we present the constraints of methyl carbamate formation under the interstellar conditions. Large complex organic molecules are favorably produced in hot-corino environments of low mass protostars. We for the first time carried out astrochemical modeling focusing on the formation of methyl carbamate in physical conditions similar to hot-corino objects. Consequently, we examined ALMA archival data for existing spectral line observations toward hot corinos NGC1333 IRAS 4A2 and IRAS 16293B. Within the common spectral range towards these sources, we found three features are possibly related to the spectral transitions of methyl carbamate and consequently estimate the upper limit of column densities. Results of chemical modeling are consistent with the observational upper limit of estimated column density/abundance toward the sources. This may hint the validation of the proposed formation mechanism. Future observations using telescope like ngVLA may confirm the presence of MC toward the hot corinos.

preprint2016arXiv

Detection of CH3SH in protostar IRAS 16293-2422

The nature of the main sulphur reservoir in star forming regions is a long standing mystery. The observed abundance of sulphur-bearing species in dense clouds is only about 0.1 per cent of the same quantity in diffuse clouds. Therefore, the main sulphur species in star forming regions of the interstellar medium are still unknown. IRAS 16293-2422 is one of the regions where production of S-bearing species is favourable due to its conditions which allows the evaporation of ice mantles. We carried out observations in the 3 mm band towards the solar type protostar IRAS 16293-2422 with the IRAM 30m telescope. We observed a single frequency setup with the EMIR heterodyne 3 mm receiver with an Lower Inner (LI) tuning frequency of 89.98 GHz. Several lines of the complex sulphur species CH3SH were detected. Observed abundances are compared with simulations using the NAUTILUS gas-grain chemical model. Modelling results suggest that CH3SH has the constant abundance of 4e-9 (compared to H2) for radii lower than 200 AU and is mostly formed on the surfaces. Detection of CH3SH indicates that there may be several new families of S-bearing molecules (which could form starting from CH3SH) which have not been detected or looked for yet.

preprint2016arXiv

Gas and grain chemical composition in cold cores as predicted by the Nautilus 3-phase model

We present an extended version of the 2-phase gas-grain code NAUTILUS to the 3-phase modelling of gas and grain chemistry of cold cores. In this model, both the mantle and the surface are considered as chemically active. We also take into account the competition among reaction, diffusion and evaporation. The model predictions are confronted to ice observations in the envelope of low-mass and massive young stellar objects as well as toward background stars. Modelled gas-phase abundances are compared to species observed toward TMC-1 (CP) and L134N dark clouds. We find that our model successfully reproduces the observed ice species. It is found that the reaction-diffusion competition strongly enhances reactions with barriers and more specifically reactions with H2, which is abundant on grains. This finding highlights the importance to have a good approach to determine the abundance of H2 on grains. Consequently, it is found that the major N-bearing species on grains go from NH3 to N2 and HCN when the reaction-diffusion competition is accounted. In the gas-phase and before few 10^5 yrs, we find that the 3-phase model does not have a strong impact on the observed species compared to the 2-phase model. After this time, the computed abundances dramatically decrease due to the strong accretion on dust, which is not counterbalanced by the desorption less efficient than in the 2-phase model. This strongly constrains the chemical-age of cold cores to be of the order of few 10^5 yrs.

preprint2016arXiv

Methylacetylene (CH3CCH) and propene (C3H6) formation in cold dense clouds: a case of dust grain chemistry

We present an extensive review of gas phase reactions producing methylacetylene and propene showing that these relatively abundant unsaturated hydrocarbons cannot be synthesized through gas-phase reactions. We explain the formation of propene and methylacetylene through surface hydrogenation of C3 depleted onto interstellar ices, C3 being a very abundant species in the gas phase.

preprint2016arXiv

Survey Observations of a Possible Glycine Precursor, Methanimine (CH$_2$NH)

We conducted survey observations of a glycine precursor, methanimine or methylenimine (CH$_2$NH), with the NRO 45 m telescope and the SMT telescope towards 12 high-mass and two low-mass star-forming regions in order to increase number of CH$_2$NH sources and to better understand the characteristics of CH2NH sources. As a result of our survey, CH$_2$NH was detected in eight sources, including four new sources. The estimated fractional abundances were ~10$^8$ in Orion KL and G10.47+0.03, while they were ~10$^9$ towards the other sources. Our hydrogen recombination line and past studies suggest that CH$_2$NH-rich sources have less evolved HII regions. The less destruction rates by UV flux from the central star would be contributed to the high CH$_2$NH abundances towards CH$_2$NH-rich sources. Our gas-grain chemical simulations suggest that CH$_2$NH is mostly formed in the gas-phase by neutral-neutral reactions rather than the product of thermal evaporation from the dust surfaces.

preprint2016arXiv

Temperature Dependent Product Yields for the Spin Forbidden Singlet Channel of the C(3P) + C2H2 Reaction

The atomic hydrogen formation channels of the C + C2H2 reaction have been investigated using a continuous supersonic flow reactor over the 52 K 296 K temperature range. H-atoms were detected directly at 121.567 nm by vacuum ultraviolet laser induced fluorescence. Absolute H-atom yields were determined by comparison with the H-atom signal generated by the C + C2H4 reaction. The product yields agree with earlier crossed beam experiments employing universal detection methods. Incorporating these branching ratios in a gas-grain model of dense interstellar clouds increases the cC3H abundance. This reaction is a minor source of C3 containing molecules in the present simulations.

preprint2016arXiv

The C(3P) + NH3 reaction in interstellar chemistry: I. Investigation of the product formation channels

The product formation channels of ground state carbon atoms, C(3P), reacting with ammonia, NH3, have been investigated using two complementary experiments and electronic structure calculations. Reaction products are detected in a gas flow tube experiment (330 K, 4 Torr) using tunable VUV photoionization coupled with time of flight mass spectrometry. Temporal profiles of the species formed and photoionization spectra are used to identify primary products of the C + NH3 reaction. In addition, H-atom formation is monitored by VUV laser induced fluorescence from room temperature to 50 K in a supersonic gas flow generated by the Laval nozzle technique. Electronic structure calculations are performed to derive intermediates, transition states and complexes formed along the reaction coordinate. The combination of photoionization and laser induced fluorescence experiments supported by theoretical calculations indicate that in the temperature and pressure range investigated, the H + H2CN production channel represents 100% of the product yield for this reaction. Kinetics measurements of the title reaction down to 50 K and the effect of the new rate constants on interstellar nitrogen hydride abundances using a model of dense interstellar clouds are reported in paper II.

preprint2016arXiv

The C(3P) + NH3 reaction in interstellar chemistry: II. Low temperature rate constants and modeling of NH, NH2 and NH3 abundances in dense interstellar clouds

A continuous supersonic flow reactor has been used to measure rate constants for the C + NH3 reaction over the temperature range 50 to 296 K. C atoms were created by the pulsed laser photolysis of CBr4. The kinetics of the title reaction were followed directly by vacuum ultra-violet laser induced fluorescence (VUV LIF) of C loss and through H formation. The experiments show unambiguously that the reaction is rapid at 296 K, becoming faster at lower temperatures, reaching a value of 1.8 10-10 cm3 molecule-1 s-1 at 50 K. As this reaction is not currently included in astrochemical networks, its influence on interstellar nitrogen hydride abundances is tested through a dense cloud model including gas-grain interactions. In particular, the effect of the ortho-to-para ratio of H2 which plays a crucial role in interstellar NH3 synthesis is examined.

preprint2015arXiv

Chemistry in Disks. IX. Observations and modeling of HCO+ and DCO+ in DM Tau

We present resolved Plateau de Bure Array observations of DM Tau in lines of HCO+ (3-2), (1-0) and DCO+ (3-2). A power-law fitting approach allowed a derivation of column densities of these two molecules. A chemical inner hole of ~50 AU was found in both HCO+ and DCO+ with DCO+ emission extending to only 450 AU. An isotopic ratio of R_D = N(DCO+) / N(HCO+) was found to range from 0.1 at 50 AU and 0.2 at 450 AU. Chemical modeling allowed an exploration of the sensitivity of these molecular abundances to physical parameters out with temperature, finding that X-rays were the domination ionization source in the HCO+ molecular region and that R_D also is sensitive to the CO depletion. The ionization fraction, assuming a steady state system, was found to be x(e-) ~ 10$^{-7}$. Modeling suggests that HCO+ is the dominant charged molecule in the disk but its contribution to ionization fraction is dwarfed by atmoic ions such as C+, S+ and H+.

preprint2015arXiv

The Deuterium Fractionation Timescale in Dense Cloud Cores: A Parameter Space Exploration

The deuterium fraction [N$_2$D$^+$]/[N$_2$H$^+$], may provide information about the ages of dense, cold gas structures, important to compare with dynamical models of cloud core formation and evolution. Here we introduce a complete chemical network with species containing up to three atoms, with the exception of the Oxygen chemistry, where reactions involving H$_3$O$^+$ and its deuterated forms have been added, significantly improving the consistency with comprehensive chemical networks. Deuterium chemistry and spin states of H$_2$ and H$_3^+$ isotopologues are included in this primarily gas-phase chemical model. We investigate dependence of deuterium chemistry on model parameters: density ($n_{\rm H}$), temperature, cosmic ray ionization rate, and gas-phase depletion factor of heavy elements ($f_{\rm D}$). We also explore the effects of time-dependent freeze-out of gas-phase species and dynamical evolution of density at various rates relative to free-fall collapse. For a broad range of model parameters, the timescales to reach large values of $D_{\rm frac}^{\rm N_2H^+} \gtrsim 0.1$, observed in some low- and high-mass starless cores, are relatively long compared to the local free-fall timescale. These conclusions are unaffected by introducing time-dependent freeze-out and considering models with evolving density, unless the initial $f_{\rm D} \gtrsim$ 10. For fiducial model parameters, achieving $D_{\rm frac}^{\rm N_2H^+} \gtrsim 0.1$ requires collapse to be proceeding at rates at least several times slower than that of free-fall collapse, perhaps indicating a dynamically important role for magnetic fields in the support of starless cores and thus the regulation of star formation.

preprint2015arXiv

The interstellar chemistry of H2C3O isomers

We present the detection of two H2C3O isomers, propynal and cyclopropenone, toward various starless cores and molecular clouds, together with upper limits for the third isomer propadienone. We review the processes controlling the abundances of H2C3O isomers in interstellar media showing that the reactions involved are gas-phase ones. We show that the abundances of these species are controlled by kinetic rather than thermodynamic effects.

preprint2014arXiv

Grain Surface Reactions in Molecular Clouds: The Effect of Cosmic Rays and Quantum Tunneling

Grain-surface reactions play an essential role in interstellar chemistry, since dust grain catalyses reactions at its surface allowing for the formation of molecules. We used a chemical model in which both gas-phase and grain-surface reactions occur and studied particularly the diffusion mechanisms on the surface of the grains. Surface reactions can occur via thermal hopping when species cross over a potential barrier or via quantum tunneling when species cross through this barrier. We show that the thermal diffusion (hopping) can be much more efficient after a cosmic ray particle collides with a dust grain, heating it to a peak temperature of 70K. We present here the results of numerical simulations after including the quantum tunneling mechanism for species H, H2 and O and considering the effect of cosmic ray particle collision on the surface reactions. As a consequence, the gas-phase and grain-surface abundances are affected and we show that more complex molecules can be formed in molecular clouds.

preprint2014arXiv

Physical and chemical structure of planet-forming disks probed by millimeter observations and modeling

Protoplanetary disks composed of dust and gas are ubiquitous around young stars and are commonly recognized as nurseries of planetary systems. Their lifetime, appearance, and structure are determined by an interplay between stellar radiation, gravity, thermal pressure, magnetic field, gas viscosity, turbulence, and rotation. Molecules and dust serve as major heating and cooling agents in disks. Dust grains dominate the disk opacities, reprocess most of the stellar radiation, and shield molecules from ionizing UV/X-ray photons. Disks also dynamically evolve by building up planetary systems which drastically change their gas and dust density structures. Over the past decade significant progress has been achieved in our understanding of disk chemical composition thanks to the upgrade or advent of new millimeter/Infrared facilities (SMA, PdBI, CARMA, Herschel, e-VLA, ALMA). Some major breakthroughs in our comprehension of the disk physics and chemistry have been done since PPV. This review will present and discuss the impact of such improvements on our understanding of the disk physical structure and chemical composition.

preprint2014arXiv

The fast C(3P) + CH3OH reaction as an efficient loss process for gas-phase interstellar methanol

Rate constants for the C(3P) + CH3OH reaction have been measured in a continuous supersonic flow reactor over the range 50 K to 296 K. C(3P) was created by the in-situ pulsed laser photolysis of CBr4, a multiphoton process which also produced some C(1D), allowing us to investigate simultaneously the low temperature kinetics of the C(1D) + CH3OH reaction. C(1D) atoms were followed by an indirect chemiluminescent tracer method in the presence of excess CH3OH. C(3P) atoms were detected by the same chemiluminescence technique and also by direct vacuum ultra-violet laser induced fluorescence (VUV LIF). Secondary measurements of product H(2S) atom formation have been undertaken allowing absolute H atom yields to be obtained by comparison with a suitable reference reaction. In parallel, statistical calculations have been performed based on ab-initio calculations of the complexes, adducts and transition states (TSs) relevant to the title reaction. By comparison with the experimental H atom yields, the preferred reaction pathways could be determined, placing important constraints on the statistical calculations. The experimental and theoretical work are in excellent agreement, predicting a negative temperature dependence of the rate constant increasing from 2.2 x 10-11 cm3 molecule-1 s-1 at 296 K to 20.0 x 10-11 cm3 molecule-1 s-1 at 50 K. CH3 and HCO are found to be the major products under our experimental conditions. As this reaction is not considered in current astrochemical networks, its influence on interstellar methanol abundances is tested using a model of dense interstellar clouds.

preprint2014arXiv

The interstellar gas-phase chemistry of HCN and HNC

We review the reactions involving HCN and HNC in dark molecular clouds to elucidate new chemical sources and sinks of these isomers. We find that the most important reactions for the HCN-HNC system are Dissociative Recombination (DR) reactions of HCNH+ (HCNH+ + e-), the ionic CN + H3+, HCN + C+, HCN and HNC reactions with H+/He+/H3+/H3O+/HCO+, the N + CH2 reaction and two new reactions: H + CCN and C + HNC. We test the effect of the new rate constants and branching ratios on the predictions of gas-grain chemical models for dark cloud conditions. The rapid C + HNC reaction keeps the HCN/HNC ratio significantly above one as long as the carbon atom abundance remains high. However, the reaction of HCN with H3+ followed by DR of HCNH+ acts to isomerize HCN into HNC when carbon atoms and CO are depleted leading to a HCN/HNC ratio close to or slightly greater than 1. This agrees well with observations in TMC-1 and L134N taking into consideration the overestimation of HNC abundances through the use of the same rotational excitation rate constants for HNC as for HCN in many radiative transfer models.

preprint2013arXiv

Astrochemistry: Synthesis and Modelling

We discuss models that astrochemists have developed to study the chemical composition of the interstellar medium. These models aim at computing the evolution of the chemical composition of a mixture of gas and dust under as- trophysical conditions. These conditions, as well as the geometry and the physical dynamics, have to be adapted to the objects being studied because different classes of objects have very different characteristics (temperatures, densities, UV radia- tion fields, geometry, history etc); e.g., proto-planetary disks do not have the same characteristics as protostellar envelopes. Chemical models are being improved continually thanks to comparisons with observations but also thanks to laboratory and theoretical work in which the individual processes are studied.

preprint2013arXiv

Chemistry of Dark Clouds: Databases, Networks, and Models

Chemical models have been developed over the years by astrophysicists to study the pro- cesses at play in the various environments of the interstellar medium (ISM) that define the chemical composition of the gas and the dust. These qualitative aspects of the model predictions have been improved from a chemical point of view thanks to many recent developments of the experimental technics and theoretical methods that aim at studying the individual reactions in conditions as close to the ISM conditions as possible and characterize the rate constants of their efficiency. These models have also been more and more associated with dynamical evolution of the ISM physical conditions (for star forming regions for instance) since the chemical composition is far from steady-state in such regions. In this paper, we try to assess the state of the art concerning the chemical modeling of dark clouds, the initial step for the formation of stars and disks.

preprint2013arXiv

Survival of interstellar molecules to prestellar dense core collapse and early phases of disk formation

An outstanding question of astrobiology is the link between the chemical composition of planets, comets, and other solar system bodies and the molecules formed in the interstellar medium. Understanding the chemical and physical evolution of the matter leading to the formation of protoplanetary disks is an important step for this. We provide some new clues to this long-standing problem using three-dimensional chemical simulations of the early phases of disk formation: we interfaced the full gas-grain chemical model Nautilus with the radiation-magnetohydrodynamic model RAMSES, for different configurations and intensities of magnetic field. Our results show that the chemical content (gas and ices) is globally conserved during the collapsing process, from the parent molecular cloud to the young disk surrounding the first Larson core. A qualitative comparison with cometary composition suggests that comets are constituted of different phases, some molecules being direct tracers of interstellar chemistry, while others, including complex molecules, seem to have been formed in disks, where higher densities and temperatures allow for an active grain surface chemistry. The latter phase, and its connection with the formation of the first Larson core, remains to be modeled.

preprint2013arXiv

The gas-phase chemistry of carbon chains in dark cloud chemical models

We review the reactions between carbon chain molecules and radicals, namely Cn, CnH, CnH2, C2n+1O, CnN, HC2n+1N, with C, N and O atoms. Rate constants and branching ratios for these processes have been re-evaluated using experimental and theoretical literature data. In total 8 new species have been introduced, 41 new reactions have been proposed and 122 rate coefficients from kida.uva.2011 (Wakelam et al. 2012) have been modified. We test the effect of the new rate constants and branching ratios on the predictions of gas-grain chemical models for dark cloud conditions using two different C/O elemental ratios. We show that the new rate constants produce large differences in the predicted abundances of carbon chains since the formation of long chains is less effective. The general agreement between the model predictions and observed abundances in the dark cloud TMC-1 (CP) is improved by the new network and we find that C/O ratios of 0.7 and 0.95 both produce a similar agreement for different times. The general agreement for L134N (N) is not significantly changed. The current work specifically highlights the importance of O + CnH and N + CnH reactions. As there are very few experimental or theoretical data for the rate constants of these reactions we highlight the need for experimental studies of the O + CnH and N + CnH reactions, particularly at low temperature.

preprint2013arXiv

Water in Protoplanetary Disks: Deuteration and Turbulent Mixing

We investigate water and deuterated water chemistry in turbulent protoplanetary disks. Chemical rate equations are solved with the diffusion term, mimicking turbulent mixing in vertical direction. Water near the midplane is transported to the disk atmosphere by turbulence and destroyed by photoreactions to produce atomic oxygen, while the atomic oxygen is transported to the midplane and reforms water and/or other molecules. We find that this cycle significantly decreases column densities of water ice at r < 30 AU, where dust temperatures are too high to reform water ice effectively. The radial extent of such region depends on the desorption energy of atomic hydrogen. Our model indicates that water ice could be deficient even outside the sublimation radius. Outside this radius, the cycle decreases the D/H ratio of water ice from 2x10^-2, which is set by the collapsing core model, to 10^-4-10^-2 in 10^6 yr, without significantly decreasing the water ice column density. The resultant D/H ratios depend on the strength of mixing and the radial distance from the central star. Our finding suggests that the D/H ratio of cometary water (10^-3-10^-4) could be established (i.e. cometary water could be formed) in the solar nebula, even if the D/H ratio of water ice delivered to the disk was very high (10^-2).

preprint2012arXiv

Chemistry in the First Hydrostatic Core Stage By Adopting Three-Dimensional Radiation Hydrodynamic Simulations

We investigate molecular evolution from a molecular cloud core to a first hydrostatic core in three spatial dimensions. We perform a radiation hydrodynamic simulation in order to trace fluid parcels, in which molecular evolution is investigated, using a gas-phase and grain-surface chemical reaction network. We derive spatial distributions of molecular abundances and column densities in the core harboring the first core. We find that the total of gas and ice abundances of many species in a cold era (10 K) remain unaltered until the temperature reaches ~500 K. The gas abundances in the warm envelope and the outer layer of the first core (T < 500 K) are mainly determined via the sublimation of ice-mantle species. Above 500 K, the abundant molecules, such as H2CO, start to be destroyed, and simple molecules, such as CO, H2O and N2 are reformed. On the other hand, some molecules are effectively formed at high temperature; carbon-chains, such as C2H2 and cyanopolyynes, are formed at the temperature of >700 K. We also find that large organic molecules, such as CH3OH and HCOOCH3, are associated with the first core (r < 10 AU). Although the abundances of these molecules in the first core stage are comparable or less than in the protostellar stage (hot corino), reflecting the lower luminosity of the central object, their column densities in our model are comparable to the observed values toward the prototypical hot corino, IRAS 16293-2422. We propose that these large organic molecules can be good tracers of the first cores.

preprint2012arXiv

Elemental nitrogen partitioning in dense interstellar clouds

Many chemical models of dense interstellar clouds predict that the majority of gas-phase elemental nitrogen should be present as N2, with an abundance approximately five orders of magnitude less than that of hydrogen. As a homonuclear diatomic molecule, N2 is difficult to detect spectroscopically through infrared or millimetre-wavelength transitions so its abundance is often inferred indirectly through its reaction product N2H+. Two main formation mechanisms each involving two radical-radical reactions are the source of N2 in such environments. Here we report measurements of the low temperature rate constants for one of these processes, the N + CN reaction down to 56 K. The effect of the measured rate constants for this reaction and those recently determined for two other reactions implicated in N2 formation are tested using a gas-grain model employing a critically evaluated chemical network. We show that the amount of interstellar nitrogen present as N2 depends on the competition between its gas-phase formation and the depletion of atomic nitrogen onto grains. As the reactions controlling N2 formation are inefficient, we argue that N2 does not represent the main reservoir species for interstellar nitrogen. Instead, elevated abundances of more labile forms of nitrogen such as NH3 should be present on interstellar ices, promoting the eventual formation of nitrogen-bearing organic molecules.

preprint2012arXiv

From Prestellar to Protostellar Cores II. Time Dependence and Deuterium Fractionation

We investigate the molecular evolution and D/H abundance ratios that develop as star formation proceeds from a dense-cloud core to a protostellar core, by solving a gas-grain reaction network applied to a 1-D radiative hydrodynamic model with infalling fluid parcels. Spatial distributions of gas and ice-mantle species are calculated at the first-core stage, and at times after the birth of a protostar. Gas-phase methanol and methane are more abundant than CO at radii $r\lesssim 100$ AU in the first-core stage, but gradually decrease with time, while abundances of larger organic species increase. The warm-up phase, when complex organic molecules are efficiently formed, is longer-lived for those fluid parcels in-falling at later stages. The formation of unsaturated carbon chains (warm carbon-chain chemistry) is also more effective in later stages; C$^+$, which reacts with CH$_4$ to form carbon chains, increases in abundance as the envelope density decreases. The large organic molecules and carbon chains are strongly deuterated, mainly due to high D/H ratios in the parent molecules, determined in the cold phase. We also extend our model to simulate simply the chemistry in circumstellar disks, by suspending the 1-D infall of a fluid parcel at constant disk radii. The species CH$_3$OCH$_3$ and HCOOCH$_3$ increase in abundance in $10^4-10^5$ yr at the fixed warm temperature; both also have high D/H ratios.

preprint2012arXiv

Review of OCS gas-phase reactions in dark cloud chemical models

The association reaction S + CO {\to} OCS + hnu has been identified as being particularly important for the prediction of gas-phase OCS abundances by chemical models of dark clouds. We performed detailed ab-initio calculations for this process in addition to undertaking an extensive review of the neutral-neutral reactions involving this species which might be important in such environments. The rate constant for this association reaction was estimated to be several orders of magnitude smaller than the one present in current astrochemical databases. The new rate for this reaction and the introduction of other processes, notably OH + CS {\to} OCS + H and C + OCS {\to} CO + CS, dramatically changes the OCS gas-phase abundance predicted by chemical models for dark clouds. The disagreement with observations in TMC-1 (CP) and L134N (N), suggests that OCS may be formed on grain surfaces as is the case for methanol. The observation of solid OCS on interstellar ices supports this hypothesis.

preprint2011arXiv

CID: Chemistry in disks VI.sulfur-bearing molecules in the protoplanetary disks surrounding LkCa15, MWC480, DM Tau, and GO Tau

We study the content in S-bearing molecules of protoplanetary disks around low-mass stars. We used the new IRAM 30-m receiver EMIR to perform simultaneous observations of the $1_{10}-1_{01}$ line of H$_2$S at 168.8 GHz and $2_{23}-1_{12}$ line of SO at 99.3 GHz. We compared the observational results with predictions coming from the astrochemical code NAUTILUS, which has been adapted to protoplanetary disks. The data were analyzed together with existing CS J=3-2 observations. We fail to detect the SO and H$_2$S lines, although CS is detected in LkCa15, DM\,Tau, and GO\,Tau but not in MWC\,480. However, our new upper limits are significantly better than previous ones and allow us to put some interesting constraints on the sulfur chemistry. Our best modeling of disks is obtained for a C/O ratio of 1.2, starting from initial cloud conditions of H density of $2\times 10^5$ cm$^{-3}$ and age of $10^6$ yr. The results agree with the CS data and are compatible with the SO upper limits, but fail to reproduce the H$_2$S upper limits. The predicted H$_2$S column densities are too high by at least one order of magnitude. H$_2$S may remain locked onto grain surfaces and react with other species, thereby preventing the desorption of H$_2$S.

preprint2006arXiv

Chemical sensitivity to the ratio of the cosmic-ray ionization rates of He and H2 in dense clouds

Aim: To determine whether or not gas-phase chemical models with homogeneous and time-independent physical conditions explain the many observed molecular abundances in astrophysical sources, it is crucial to estimate the uncertainties in the calculated abundances and compare them with the observed abundances and their uncertainties. Non linear amplification of the error and bifurcation may limit the applicability of chemical models. Here we study such effects on dense cloud chemistry. Method: Using a previously studied approach to uncertainties based on the representation of rate coefficient errors as log normal distributions, we attempted to apply our approach using as input a variety of different elemental abundances from those studied previously. In this approach, all rate coefficients are varied randomly within their log normal (Gaussian) distribution, and the time-dependent chemistry calculated anew many times so as to obtain good statistics for the uncertainties in the calculated abundances. Results: Starting with so-called ``high-metal'' elemental abundances, we found bimodal rather than Gaussian like distributions for the abundances of many species and traced these strange distributions to an extreme sensitivity of the system to changes in the ratio of the cosmic ray ionization rate zeta\_He for He and that for molecular hydrogen zeta\_H2. The sensitivity can be so extreme as to cause a region of bistability, which was subsequently found to be more extensive for another choice of elemental abundances. To the best of our knowledge, the bistable solutions found in this way are the same as found previously by other authors, but it is best to think of the ratio zeta\_He/zeta\_H2 as a control parameter perpendicular to the ''standard'' control parameter zeta/n\_H.