Showing 1–36 of 36 results for author: Vasyunin, A I

Deuterium fractionation in cold dense cores in the low-mass star forming region L1688

Authors: I. V. Petrashkevich, A. F. Punanova, P. Caselli, O. Sipilä, J. E. Pineda, R. K. Friesen, M. G. Korotaeva, A. I. Vasyunin

Abstract: In this work, we study deuterium fractionation in four starless cores in the low-mass star-forming region L1688 in the Ophiuchus molecular cloud. We study how the deuterium fraction ($R_D$) changes with environment, compare deuteration of ions and neutrals, core centre and its envelope, and attempt to reproduce the observed results with a gas-grain chemical model. We chose high and low gas density… ▽ More In this work, we study deuterium fractionation in four starless cores in the low-mass star-forming region L1688 in the Ophiuchus molecular cloud. We study how the deuterium fraction ($R_D$) changes with environment, compare deuteration of ions and neutrals, core centre and its envelope, and attempt to reproduce the observed results with a gas-grain chemical model. We chose high and low gas density tracers to study both core centre and the envelope. With the IRAM 30m antenna, we mapped N$_2$H$^+$(1-0), N$_2$D$^+$(1-0), H$^{13}$CO$^+$ (1-0) and (2-1), DCO$^+$(2-1), and $p$-NH$_2$D(1$_{11}$-1$_{01}$) towards the chosen cores. The missing $p$-NH$_3$ and N$_2$H$^+$(1-0) data were taken from the literature. To measure the molecular hydrogen column density, dust and gas temperature within the cores, we used the Herschel/SPIRE dust continuum emission data, the GAS survey data (ammonia), and the COMPLETE survey data to estimate the upper limit on CO depletion. We present the deuterium fraction maps for three species towards four starless cores. Deuterium fraction of the core envelopes traced by DCO$^+$/H$^{13}$CO$^+$ is one order of magnitude lower ($\sim$0.08) than that of the core central parts traced by the nitrogen-bearing species ($\sim$0.5). Deuterium fraction increases with the gas density as indicated by high deuterium fraction of high gas density tracers and low deuterium fraction of lower gas density tracers and by the decrease of $R_D$ with core radii, consistent with the predictions of the chemical model. Our model results show a good agreement with observations for $R_D$(N$_2$D$^+$/N$_2$H$^+$) and R$_D$(DCO$^+$/HCO$^+$) and underestimate the $R_D$(NH$_2$D/NH$_3$). △ Less

Submitted 10 January, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: 30 pages with 22 figures, accepted for publication in MNRAS

MSC Class: 70F15 (Primary) ACM Class: J.2

Nuclear spin ratios of deuterated ammonia in prestellar cores. LAsMA observations of H-MM1 and Oph D

Authors: Jorma Harju, Jaime E. Pineda, Olli Sipilä, Paola Caselli, Arnaud Belloche, Friedrich Wyrowski, Wiebke Riedel, Elena Redaelli, Anton I. Vasyunin

Abstract: We determine the ortho/para ratios of NH2D and NHD2 in two dense, starless cores, where their formation is supposed to be dominated by gas-phase reactions, which, in turn, is predicted to result in deviations from the statistical spin ratios. The Large APEX sub-Millimeter Array (LAsMA) multibeam receiver of the Atacama Pathfinder EXperiment (APEX) telescope was used to observe the prestellar cores… ▽ More We determine the ortho/para ratios of NH2D and NHD2 in two dense, starless cores, where their formation is supposed to be dominated by gas-phase reactions, which, in turn, is predicted to result in deviations from the statistical spin ratios. The Large APEX sub-Millimeter Array (LAsMA) multibeam receiver of the Atacama Pathfinder EXperiment (APEX) telescope was used to observe the prestellar cores H-MM1 and Oph D in Ophiuchus in the ground-state lines of ortho and para NH2D and NHD2. The fractional abundances of these molecules were derived employing 3D radiative transfer modelling, using different assumptions about the abundance profiles as functions of density. We also ran gas-grain chemistry models with different scenarios concerning proton or deuteron exchanges and chemical desorption from grains to find out if one of these models can reproduce the observed spin ratios. The observationally deduced ortho/para ratios of NH2D and NHD2 are in both cores within 10% of their statistical values 3 and 2, respectively, and taking 3-sigma limits, deviations from these of about 20% are allowed. Of the chemistry models tested here, the model that assumes proton hop (as opposed to full scrambling) in reactions contributing to ammonia formation, and a constant efficiency of chemical desorption, comes nearest to the observed abundances and spin ratios. The nuclear spin ratios derived here are in contrast with spin-state chemistry models that assume full scrambling in proton donation and hydrogen abstraction reactions leading to deuterated ammonia. The efficiency of chemical desorption influences strongly the predicted abundances of NH3, NH2D, and NHD2, but has a lesser effect on their ortho/para ratios. For these the proton exchange scenario in the gas is decisive. We suggest that this is because of rapid re-processing of ammonia and related cations by gas-phase ion-molecule reactions. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: accepted for publication in Astronomy & Astrophysics

Modelling Deuterated Isotopologues of Methanol toward the Pre-Stellar Core L1544

Authors: W. Riedel, O. Sipilä, E. Redaelli, P. Caselli, A. I. Vasyunin, F. Dulieu, N. Watanabe

Abstract: Aims. We aim to improve a previous model for the prediction of column densities and deuterium fractions of non- and singly deuterated methanol. Thereby, we try to identify crucial chemical and physical parameters, for which the study of deuteration could provide valuable additional constraints. Methods. We employed a gas-grain chemical code to devise a model that is in agreement with the observe… ▽ More Aims. We aim to improve a previous model for the prediction of column densities and deuterium fractions of non- and singly deuterated methanol. Thereby, we try to identify crucial chemical and physical parameters, for which the study of deuteration could provide valuable additional constraints. Methods. We employed a gas-grain chemical code to devise a model that is in agreement with the observed column density and deuterium fraction profiles of the innermost region of the pre-stellar core L1544. For that purpose, we developed a new treatment of reactive desorption, deriving an individual reactive desorption efficiency for every product species in a chemical reaction, that depends on the reaction enthalpy and type of underlying surface. Furthermore, we explored several options to promote the diffusion of hydrogen and deuterium atoms over the surface of interstellar dust grains, in order to increase methanol formation. Results. Our fiducial model employs diffusion by quantum tunneling of hydrogen and deuterium atoms, resulting in CH$_3$OH and CH$_2$DOH column densities that are approximately an order of magnitude lower than the observed values, which improves the results compared to the previous model by a factor 10. The $N$(CH$_2$DOH)/$N$(CH$_3$OH) ratio is reproduced within a factor of 1.2 for the centre and 1.8 for the position of the methanol peak. Given the large uncertainties that chemical models typically have, we consider our predictions to be in agreement with the observations. In general, we conclude that a diffusion process with a high diffusion rate needs to be employed to obtain methanol column densities that are in accordance with the observed values. Also, we find that the introduction of abstraction reactions into the methanol formation scheme suppresses deuteration, when used in combination with a high diffusion rate. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Journal ref: A&A 680, A87 (2023)

Impact of HAC evolution on the formation of small hydrocarbons in the Orion Bar and the Horsehead PDRs

Authors: M. S. Murga, A. I. Vasyunin, M. S. Kirsanova

Abstract: We study evolution of hydrogenated amorphous carbon (HAC) grains under harsh UV radiation in photo-dissociation regions (PDRs) near young massive stars. Our aim is to evaluate the impact of the HAC grains on formation of observed small hydrocarbons: C$_2$H, C$_2$H$_2$, C$_3$H$^+$, C$_3$H, C$_3$H$_2$, C$_4$H, in PDRs. We developed a microscopic model of the HAC grains based on available experimenta… ▽ More We study evolution of hydrogenated amorphous carbon (HAC) grains under harsh UV radiation in photo-dissociation regions (PDRs) near young massive stars. Our aim is to evaluate the impact of the HAC grains on formation of observed small hydrocarbons: C$_2$H, C$_2$H$_2$, C$_3$H$^+$, C$_3$H, C$_3$H$_2$, C$_4$H, in PDRs. We developed a microscopic model of the HAC grains based on available experimental results. The model includes processes of photo- and thermodesorption, accretion of hydrogen and carbon atoms and subsequent formation of carbonaceous mantle on dust surface. H$_2$, CH$_4$, C$_2$H$_2$, C$_2$H$_4$, C$_2$H$_6$, C$_3$H$_4$, C$_3$H$_6$, C$_3$H$_8$ are considered as the main fragments of the HAC photo-destruction. We simulated evolution of the HAC grains under the physical conditions of two PDRs, the Orion Bar and the Horsehead nebula. We estimated the production rates of the HAC fragments in gas phase chemical reactions and compared them with the production rates of fragments due to the HAC destruction. The latter rates may dominate under some conditions, namely, at A$_V$=0.1 in both PDRs. We coupled our model with the gas-grain chemical model MONACO and calculated abundances of observed small hydrocarbons. We conclude that the contribution of the HAC destruction fragments to chemistry is not enough to match the observed abundances, although it increases the abundances by several orders of magnitude in the Orion Bar at A$_V$=0.1. Additionally, we found that the process of carbonaceous mantle formation on dust surface can be an inhibitor for the formation of observed small hydrocarbons in PDRs. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Comments: Accepted to MNRAS

The complex organic molecular content in the L1517B starless core

Authors: Andrés Megías, Izaskun Jiménez-Serra, Jesús Martín-Pintado, Anton I. Vasyunin, Silvia Spezzano, Paola Caselli, Giuliana Cosentino, Serena Viti

Abstract: Recent observations of the pre-stellar core L1544 and the younger starless core L1498 have revealed that complex organic molecules (COMs) are enhanced in the gas phase toward their outer and intermediate-density shells. Our goal is to determine the level of chemical complexity toward the starless core L1517B, which seems younger than L1498, and compare it with the other two previously studied core… ▽ More Recent observations of the pre-stellar core L1544 and the younger starless core L1498 have revealed that complex organic molecules (COMs) are enhanced in the gas phase toward their outer and intermediate-density shells. Our goal is to determine the level of chemical complexity toward the starless core L1517B, which seems younger than L1498, and compare it with the other two previously studied cores to see if there is a chemical evolution within the cores. We have carried out 3 mm high-sensitivity observations toward two positions in the L1517B starless core: the core's centre and the position where the methanol emission peaks (at a distance of $\sim$5000 au from the core's centre). Our observations reveal that a lower number of COMs and COM precursors are detected in L1517B with respect to L1498 and L1544, and also show lower abundances. Besides methanol, we only detected CH$_3$O, H$_2$CCO, CH$_3$CHO, CH$_3$CN, CH$_3$NC, HCCCN, and HCCNC. Their measured abundances are $\sim$3 times larger toward the methanol peak than toward the core's centre, mimicking the behaviour found toward the more evolved cores L1544 and L1498. We propose that the differences in the chemical complexity observed between the three studied starless cores are a consequence of their evolution, with L1517B being the less evolved one, followed by L1498 and L1544. Chemical complexity in these cores seems to increase over time, with N-bearing molecules forming first and O-bearing COMs forming at a later stage as a result of the catastrophic depletion of CO. △ Less

Submitted 16 August, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: 18 pages, 13 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 519, Issue 2, February 2023, Pages 1601-1617

Physical and chemical vertical structure of magnetostatic accretion disks of young stars

Authors: Sergey A. Khaibrakhmanov, Alexander E. Dudorov, Anton I. Vasyunin, Mikhail Yu. Kiskin

Abstract: The vertical structure of the accretion disks of young stars with fossil large-scale magnetic field is studied. The equations of magnetostatic equilibrium of the disk are solved taking into account the stellar gravity, gas and magnetic pressure, turbulent heating, and heating by stellar radiation. The modelled physical structure of the disk is used to simulate its chemical structure, in particular… ▽ More The vertical structure of the accretion disks of young stars with fossil large-scale magnetic field is studied. The equations of magnetostatic equilibrium of the disk are solved taking into account the stellar gravity, gas and magnetic pressure, turbulent heating, and heating by stellar radiation. The modelled physical structure of the disk is used to simulate its chemical structure, in particular, to study the spatial distribution of CN molecules. The disk of the typical T Tauri star is considered. Simulations show that the temperature within the disk in the region $r<50$ au decreases with height and density profiles are steeper than in the isothermal case. Outside the `dead' zone, vertical profiles of the azimuthal component of the magnetic field are nonmonotonic, and the magnetic field strength maximum is reached within the disk. The magnetic pressure gradient can cause an increase in the disk thickness in comparison with the hydrostatic one. The CN molecule concentration is maximum near the photosphere and in the disk atmosphere where the magnetic field strength at chosen parameters is $\sim 0.01$ G. Measurements of the Zeeman splitting of CN lines in the submm range can be used to determine the magnetic field strength in these regions of accretion disks. △ Less

Submitted 23 October, 2022; originally announced October 2022.

Comments: 3 pages, 1 figure, published in Bulletin of the Lebedev Physics Institute

Journal ref: Bulletin of the Lebedev Physics Institute, 2021, vol. 48, no. 10, pp. 312-316

Mass ejection and time variability in protostellar outflows: Cep E. SOLIS XVI

Authors: A. de A. Schutzer, P. R. Rivera-Ortiz, B. Lefloch, A. Gusdorf, C. Favre, D. Segura-Cox, A. Lopez-Sepulcre, R. Neri, J. Ospina-Zamudio, M. De Simone, C. Codella, S. Viti, L. Podio, J. Pineda, R. O'Donoghue, C. Ceccarelli, P. Caselli, F. Alves, R. Bachiller, N. Balucani, E. Bianchi, L. Bizzocchi, S. Bottinelli, E. Caux, A. Chacón-Tanarro , et al. (24 additional authors not shown)

Abstract: Protostellar jets are an important agent of star formation feedback, tightly connected with the mass-accretion process. The history of jet formation and mass-ejection provides constraints on the mass accretion history and the nature of the driving source. We want to characterize the time-variability of the mass-ejection phenomena at work in the Class 0 protostellar phase, in order to better unders… ▽ More Protostellar jets are an important agent of star formation feedback, tightly connected with the mass-accretion process. The history of jet formation and mass-ejection provides constraints on the mass accretion history and the nature of the driving source. We want to characterize the time-variability of the mass-ejection phenomena at work in the Class 0 protostellar phase, in order to better understand the dynamics of the outflowing gas and bring more constraints on the origin of the jet chemical composition and the mass-accretion history. We have observed the emission of the CO 2-1 and SO N_J=5_4-4_3 rotational transitions with NOEMA, towards the intermediate-mass Class 0 protostellar system Cep E. The CO high-velocity jet emission reveals a central component associated with high-velocity molecular knots, also detected in SO, surrounded by a collimated layer of entrained gas. The gas layer appears to accelerate along the main axis over a length scale delta_0 ~700 au, while its diameter gradually increases up to several 1000au at 2000au from the protostar. The jet is fragmented into 18 knots of mass ~10^-3 Msun, unevenly distributed between the northern and southern lobes, with velocity variations up to 15 km/s close to the protostar, well below the jet terminal velocities. The knot interval distribution is approximately bimodal with a scale of ~50-80yr close to the protostar and ~150-200yr at larger distances >12". The mass-loss rates derived from knot masses are overall steady, with values of 2.7x10^-5 Msun/yr (8.9x10^-6 Msun/yr) in the northern (southern) lobe. The interaction of the ambient protostellar material with high-velocity knots drives the formation of a molecular layer around the jet, which accounts for the higher mass-loss rate in the north. The jet dynamics are well accounted for by a simple precession model with a period of 2000yr and a mass-ejection period of 55yr. △ Less

Submitted 18 March, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 13 pages, 9 figures, 3 table. Accepted in A&A

Journal ref: A&A 662, A104 (2022)

Dark cloud-type chemistry in PDRs with moderate UV field

Authors: Maria S. Kirsanova, Anna F. Punanova, Dmitry A. Semenov, Anton I. Vasyunin

Abstract: We present a study of emission lines of small hydrocarbons C$_2$H and $c$-C$_3$H$_2$, and COMs precursors H$_2$CO and CH$_3$OH in order to better understand the possible chemical link between the molecular abundances and UV radiation field in photodissociation regions (PDRs). We study two PDRs around extended and compact HII regions with $G \leq 50$~Habings in the S235 star-forming complex. We fin… ▽ More We present a study of emission lines of small hydrocarbons C$_2$H and $c$-C$_3$H$_2$, and COMs precursors H$_2$CO and CH$_3$OH in order to better understand the possible chemical link between the molecular abundances and UV radiation field in photodissociation regions (PDRs). We study two PDRs around extended and compact HII regions with $G \leq 50$~Habings in the S235 star-forming complex. We find the highest abundances of both hydrocarbons on the edges of molecular clumps, while $c$-C$_3$H$_2$ is also abundant in the low-density expanding PDR around compact HII region S235\,A. We see the highest methanol column density towards the positions with the UV~field $G\approx 20-30$~Habings and explain them by reactive desorption from the dust grains. The $N_{\rm C_2H}/N_{\rm CH_3OH}$ ratio is lower by a factor of few or the order of magnitude in comparison with the Horsehead and Orion Bar PDRs. The ratio is similar to the value observed in hot corinos in the Perseus cloud. We conclude that ion-molecular and grain surface chemical routes rule the molecular abundances in the PDRs, and the PDRs inherit molecular abundances from the previous dark stage of molecular cloud evolution in spite of massive stars already emitting in optics. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: accepted by MNRAS

The complex organic molecular content in the L1498 starless core

Authors: Izaskun Jimenez-Serra, Anton I. Vasyunin, Silvia Spezzano, Paola Caselli, Giuliana Cosentino, Serena Viti

Abstract: Observations carried out toward starless and pre-stellar cores have revealed that complex organic molecules are prevalent in these objects, but it is unclear what chemical processes are involved in their formation. Recently, it has been shown that complex organics are preferentially produced at an intermediate-density shell within the L1544 pre-stellar core at radial distances of ~4000 au with res… ▽ More Observations carried out toward starless and pre-stellar cores have revealed that complex organic molecules are prevalent in these objects, but it is unclear what chemical processes are involved in their formation. Recently, it has been shown that complex organics are preferentially produced at an intermediate-density shell within the L1544 pre-stellar core at radial distances of ~4000 au with respect to the core center. However, the spatial distribution of complex organics has only been inferred toward this core and it remains unknown whether these species present a similar behaviour in other cores. We report high-sensitivity observations carried out toward two positions in the L1498 pre-stellar core, the dust peak and a position located at a distance of ~11000 au from the center of the core where the emission of CH$_3$OH peaks. Similarly to L1544, our observations reveal that small O-bearing molecules and N-bearing species are enhanced by factors ~4-14 toward the outer shell of L1498. However, unlike L1544, large O-bearing organics such as CH3CHO, CH3OCH3 or CH3OCHO are not detected within our sensitivity limits. For N-bearing organics, these species are more abundant toward the outer shell of the L1498 pre-stellar core than toward the one in L1544. We propose that the differences observed between O-bearing and N-bearing species in L1498 and L1544 are due to the different physical structure of these cores, which in turn is a consequence of their evolutionary stage, with L1498 being younger than L1544. △ Less

Submitted 18 May, 2021; originally announced May 2021.

Comments: Accepted for Astrophysical Journal

A New Method for Simulating Photoprocesses in Astrochemical Models

Authors: Ella Mullikin, Hannah Anderson, Natalie O'Hern, Megan Farrah, Christopher R. Arumainayagam, Ewine F. van Dishoeck, Perry A. Gerakines, Anton I. Vasyunin, Liton Majumdar, Paola Caselli, Christopher N. Shingledecker

Abstract: We propose a new model for treating solid-phase photoprocesses in interstellar ice analogues. In this approach, photoionization and photoexcitation are included in more detail, and the production of electronically-excited (suprathermal) species is explicitly considered. In addition, we have included non-thermal, non-diffusive chemistry to account for the low-temperature characteristic of cold core… ▽ More We propose a new model for treating solid-phase photoprocesses in interstellar ice analogues. In this approach, photoionization and photoexcitation are included in more detail, and the production of electronically-excited (suprathermal) species is explicitly considered. In addition, we have included non-thermal, non-diffusive chemistry to account for the low-temperature characteristic of cold cores. As an initial test of our method, we have simulated two previous experimental studies involving the UV irradiation of pure solid O$_2$. In contrast to previous solid-state astrochemical model calculations which have used gas-phase photoabsorption cross-sections, we have employed solid-state cross-sections in our calculations. This method allows the model to be tested using well-constrained experiments rather than poorly constrained gas-phase abundances in ISM regions. Our results indicate that inclusion of non-thermal reactions and suprathermal species allows for reproduction of low-temperature solid-phase photoprocessing that simulate interstellar ices within cold ($\sim$ 10 K) dense cores such as TMC-1. △ Less

Submitted 5 February, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: ApJ, accepted: 15 pages, 3 figures

Water and Methanol Ice in L1544

Authors: Miwa Goto, A. I. Vasyunin, B. M. Giuliano, I. Jiménez-Serra, P. Caselli, C. G. Román-Zúñiga, J. Alves

Abstract: Methanol and complex organic molecules have been found in cold starless cores, where a standard warm-up scenario would not work because of the absence of heat sources. A recent chemical model attributed the presence of methanol and large organics to the efficient chemical desorption and a class of neutral-neutral reactions that proceed fast at low temperatures in the gas phase. The model calls for… ▽ More Methanol and complex organic molecules have been found in cold starless cores, where a standard warm-up scenario would not work because of the absence of heat sources. A recent chemical model attributed the presence of methanol and large organics to the efficient chemical desorption and a class of neutral-neutral reactions that proceed fast at low temperatures in the gas phase. The model calls for a high abundance of methanol ice at the edge of the CO freeze-out zone in cold cloud cores. We performed medium resolution spectroscopy toward 3 field stars behind the starless core L1544 at 3 $μ$m to constrain the methanol ice abundance and compare it with the model predictions. One of the field stars shows a methanol-ice abundance of 11% with respect to water ice. This is higher than the typical methanol abundance previously found in cold cloud cores (4%), but is 4.5 times smaller than predicted. The reason for the disagreement between the observations and the model calculations is not yet understood. △ Less

Submitted 20 December, 2020; originally announced December 2020.

Comments: 9 pages, 6 figures, accepted to A&A

Journal ref: A&A 651, A53 (2021)

Impact of PAH photodissociation on the formation of small hydrocarbons in the Orion Bar and the Horsehead PDRs

Authors: M. S. Murga, M. S. Kirsanova, A. I. Vasyunin, Ya. N. Pavlyuchenkov

Abstract: We study whether polycyclic aromatic hydrocarbons (PAHs) can be a weighty source of small hydrocarbons in photo-dissociation regions (PDRs). We modeled the evolution of 20 specific PAH molecules in terms of dehydrogenation and destruction of the carbon skeleton under the physical conditions of two well-studied PDRs, the Orion Bar and the Horsehead nebula which represent prototypical examples of PD… ▽ More We study whether polycyclic aromatic hydrocarbons (PAHs) can be a weighty source of small hydrocarbons in photo-dissociation regions (PDRs). We modeled the evolution of 20 specific PAH molecules in terms of dehydrogenation and destruction of the carbon skeleton under the physical conditions of two well-studied PDRs, the Orion Bar and the Horsehead nebula which represent prototypical examples of PDRs irradiated by "high" and "low" ultraviolet radiation field. PAHs are described as microcanonical systems. The acetylene molecule is considered as the main carbonaceous fragment of the PAH dissociation as it follows from laboratory experiments and theory. We estimated the rates of acetylene production in gas phase chemical reactions and compared them with the rates of the acetylene production through the PAH dissociation. It is found that the latter rates can be higher than the former rates in the Orion Bar at $A_{\rm V}<1$ and also at $A_{\rm V}>3.5$. In the Horsehead nebula, the chemical reactions provide more acetylene than the PAH dissociation. The produced acetylene participate in the reactions of the formation of small hydrocarbons (C$_2$H, C$_3$H, C$_3$H$^{+}$, C$_3$H$_2$, C$_4$H). Acetylene production via the PAH destruction may increase the abundances of small hydrocarbons produced in gas phase chemical reactions in the Orion Bar only at $A_{\rm V}>3.5$. In the Horsehead nebula, the contribution of PAHs to the abundances of the small hydrocarbons is negligible. We conclude that the PAHs are not a major source of small hydrocarbons in both PDRs except some locations in the Orion Bar. △ Less

Submitted 13 July, 2020; originally announced July 2020.

Comments: Accepted to MNRAS. 13 pages, 10 figures

Optical properties and dust temperatures in clumpy diffuse medium

Authors: Andrei B. Ostrovskii, S. Yu. Parfenov, A. I. Vasyunin, A. V. Ivlev, V. A. Sokolova

Abstract: In this study, we explore the impact of inhomogeneities in the spatial distribution of interstellar dust on spatial scales of $\le1$ au caused by ion shadowing forces on the optical properties of diffuse interstellar medium (ISM) as well as on the dust temperature. We show that recently proposed possibility that interstellar dust grains in the diffuse ISM are grouped in spherical cloudlets (clumps… ▽ More In this study, we explore the impact of inhomogeneities in the spatial distribution of interstellar dust on spatial scales of $\le1$ au caused by ion shadowing forces on the optical properties of diffuse interstellar medium (ISM) as well as on the dust temperature. We show that recently proposed possibility that interstellar dust grains in the diffuse ISM are grouped in spherical cloudlets (clumps) may significantly affect the observed optical properties of the diffuse ISM in comparison to that calculated under the commonly accepted assumption on the uniform dust/gas mixture if the size of clumps $\gtrsim0.1$ au. We found that opacity of an arbitrary region of diffuse ISM quickly decreases with growth of dusty clumps. We also studied the dependence of opacity and dust temperature inside the dusty clumps on their size. We show that the clumps larger than 0.1 au are opaque for far ultraviolet radiation. Dust temperature exhibit a gradient inside a clump, decreasing from the edge to the center by several degrees for a clump of a size of 0.1 au and larger. We argue that dust temperatures and high opacity within clumps larger than 0.1 au may facilitate somewhat more efficient synthesis of molecules on surfaces of interstellar grains in the diffuse ISM than it was anticipated previously. On the other hand, the presence of clumps with sizes below 0.1 au makes small or negligible influence on the optical properties of the diffuse ISM in comparison to the case with uniform dust/gas mixture. △ Less

Submitted 26 May, 2020; v1 submitted 21 May, 2020; originally announced May 2020.

Comments: 13 pages, 7 figures, accepted to MNRAS

Seeds of Life in Space (SOLIS).VII. Discovery of a cold dense methanol blob toward the L1521F VeLLO system

Authors: C. Favre, C. Vastel, I. Jimenez-Serra, D. Quénard, P. Caselli, C. Ceccarelli, A. Chacón-Tanarro, F. Fontani, J. Holdship, Y. Oya, A. Punanova, N. Sakai, S. Spezzano, S. Yamamoto, R. Neri, A. López-Sepulcre, F. Alves, R. Bachiller, N. Balucani, E. Bianchi, L. Bizzocchi, C. Codella, E. Caux, M. De Simone, J. Enrique Romero , et al. (18 additional authors not shown)

Abstract: The SOLIS (Seeds Of Life In Space) IRAM/NOEMA Large Program aims at studying a set of crucial complex organic molecules in a sample of sources, with well-known physical structure, covering the various phases of Solar-type star formation. One representative object of the transition from the prestellar core to the protostar phases has been observed toward the Very Low Luminosity Object (VeLLO) calle… ▽ More The SOLIS (Seeds Of Life In Space) IRAM/NOEMA Large Program aims at studying a set of crucial complex organic molecules in a sample of sources, with well-known physical structure, covering the various phases of Solar-type star formation. One representative object of the transition from the prestellar core to the protostar phases has been observed toward the Very Low Luminosity Object (VeLLO) called L1521F. This type of source is important to study to make the link between prestellar cores and Class 0 sources and also to constrain the chemical evolution during the process of star formation. Two frequency windows (81.6-82.6 GHz and 96.65-97.65 GHz) were used to observe the emission from several complex organics toward the L1521F VeLLO. Only 2 transitions of methanol (A+, E2) have been detected in the narrow window centered at 96.7 GHz (with an upper limit on E1) in a very compact emission blob (~7'' corresponding to ~1000au) toward the NE of the L1521F protostar. The CS 2-1 transition is also detected within the WideX bandwidth. Consistently, with what has been found in prestellar cores, the methanol emission appears ~1000au away from the dust peak. The location of the methanol blob coincides with one of the filaments previously reported in the literature. The Tex of the gas inferred from methanol is (10$\pm$2) K, while the H2 gas density (estimated from the detected CS 2-1 emission and previous CS 5-4 ALMA obs.) is a factor >25 higher than the density in the surrounding environment (n(H2) >10$^{7}$ cm$^{-3}$). From its compactness, low excitation temperature and high gas density, we suggest that the methanol emission detected with NOEMA is either a cold and dense shock-induced blob, recently formed ($\leq$ few hundred years) by infalling gas or a cold and dense fragment that may have just been formed as a result of the intense gas dynamics found within the L1521F VeLLO system. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 11 pages, 7 figures, 3 tables, accepted for publication by A&A

Journal ref: A&A 635, A189 (2020)

Efficient methanol production on the dark side of a prestellar core

Authors: Jorma Harju, Jaime E. Pineda, Anton I. Vasyunin, Paola Caselli, Stella S. R. Offner, Alyssa A. Goodman, Mika Juvela, Olli Sipilae, Alexandre Faure, Romane Le Gal, Pierre Hily-Blant, Joao Alves, Luca Bizzocchi, Andreas Burkert, Hope Chen, Rachel K. Friesen, Rolf Guesten, Philip C. Myers, Anna Punanova, Claire Rist, Erik Rosolowsky, Stephan Schlemmer, Yancy Shirley, Silvia Spezzano, Charlotte Vastel , et al. (1 additional authors not shown)

Abstract: We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particul… ▽ More We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity. △ Less

Submitted 4 May, 2020; v1 submitted 27 March, 2019; originally announced March 2019.

Comments: Accepted for publication in the Astrophysical Journal

O2 signature in thin and thick O2-H2O ices

Authors: B. Müller, B. M. Giuliano, L. Bizzocchi, A. I. Vasyunin, P. Caselli

Abstract: Aims. In this paper we investigate the detectability of the molecular oxygen in icy dust grain mantles towards astronomical objects. Methods. We present a systematic set of experiments with O2-H2O ice mixtures designed to disentangle how the molecular ratio affects the O2 signature in the mid- and near-infrared spectral regions. All the experiments were conducted in a closed-cycle helium cryostat… ▽ More Aims. In this paper we investigate the detectability of the molecular oxygen in icy dust grain mantles towards astronomical objects. Methods. We present a systematic set of experiments with O2-H2O ice mixtures designed to disentangle how the molecular ratio affects the O2 signature in the mid- and near-infrared spectral regions. All the experiments were conducted in a closed-cycle helium cryostat coupled to a Fourier transform infrared spectrometer. The ice mixtures comprise varying thicknesses from 8 $\times$ 10$^{-3}$ to 3 $μ$m. The absorption spectra of the O2-H2O mixtures are also compared to the one of pure water. In addition, the possibility to detect the O2 in icy bodies and in the interstellar medium is discussed. Results. We are able to see the O2 feature at 1551 cm$^{-1}$ even for the most diluted mixture of H2O : O2 = 9 : 1, comparable to a ratio of O2/H2O = 10 % which has already been detected in situ in the coma of the comet 67P/Churyumov-Gerasimenko. We provide an estimate for the detection of O2 with the future mission of the James Webb Space Telescope (JWST). △ Less

Submitted 25 October, 2018; originally announced October 2018.

Comments: 11 pages, 10 figures, article in press, to appear in A&A 2018

Journal ref: A&A 620, A46 (2018)

Physical properties and chemical composition of the cores in the California molecular cloud

Authors: Guo-Yin Zhang, Jin-Long Xu, A. I. Vasyunin, D. A. Semenov, Jun-Jie Wang, Sami Dib, Tie Liu, Sheng-Yuan Liu, Chuan-Peng Zhang, Xiao-Lan Liu, Ke Wang, Di Li, Zhong-Zu Wu, Jing-Hua Yuan, Da-Lei Li, Yang Gao

Abstract: We aim to reveal the physical properties and chemical composition of the cores in the California molecular cloud (CMC), so as to better understand the initial conditions of star formation. We made a high-resolution column density map (18.2") with Herschel data, and extracted a complete sample of the cores in the CMC with the \textsl{fellwalker} algorithm. We performed new single-pointing observati… ▽ More We aim to reveal the physical properties and chemical composition of the cores in the California molecular cloud (CMC), so as to better understand the initial conditions of star formation. We made a high-resolution column density map (18.2") with Herschel data, and extracted a complete sample of the cores in the CMC with the \textsl{fellwalker} algorithm. We performed new single-pointing observations of molecular lines near 90 GHz with the IRAM 30m telescope along the main filament of the CMC. In addition, we also performed a numerical modeling of chemical evolution for the cores under the physical conditions. We extracted 300 cores, of which 33 are protostellar and 267 are starless cores. About 51\% (137 of 267) of the starless cores are prestellar cores. Three cores have the potential to evolve into high-mass stars. The prestellar core mass function (CMF) can be well fit by a log-normal form. The high-mass end of the prestellar CMF shows a power-law form with an index $α=-0.9\pm 0.1$ that is shallower than that of the Galactic field stellar mass function. Combining the mass transformation efficiency ($\varepsilon$) from the prestellar core to the star of $15\pm 1\%$ and the core formation efficiency (CFE) of 5.5\%, we suggest an overall star formation efficiency of about 1\% in the CMC. In the single-pointing observations with the IRAM 30m telescope, we find that 6 cores show blue-skewed profile, while 4 cores show red-skewed profile. [$\rm {HCO}^{+}$]/[HNC] and [$\rm {HCO}^{+}$]/$\rm [N_{2}H^{+}]$ in protostellar cores are higher than those in prestellar cores; this can be used as chemical clocks. The best-fit chemical age of the cores with line observations is $\sim 5\times 10^4$~years. △ Less

Submitted 12 March, 2019; v1 submitted 22 October, 2018; originally announced October 2018.

Comments: Accepted by Astronomy & Astrophysics (A&A)

Journal ref: A&A 620, A163 (2018)

Seeds of Life in Space (SOLIS). III. Zooming into the methanol peak of the pre-stellar core L1544

Authors: Anna Punanova, Paola Caselli, Siyi Feng, Ana Chacón-Tanarro, Cecilia Ceccarelli, Roberto Neri, Francesco Fontani, Izaskun Jiménez-Serra, Charlotte Vastel, Luca Bizzocchi, Andy Pon, Anton I. Vasyunin, Silvia Spezzano, Pierre Hily-Blant, Leonardo Testi, Serena Viti, Satoshi Yamamoto, Felipe Alves, Rafael Bachiller, Nadia Balucani, Eleonora Bianchi, Sandrine Bottinelli, Emmanuel Caux, Rumpa Choudhury, Claudio Codella , et al. (19 additional authors not shown)

Abstract: Towards the pre-stellar core L1544, the methanol (CH$_3$OH) emission forms an asymmetric ring around the core centre, where CH$_3$OH is mostly in solid form, with a clear peak 4000~au to the north-east of the dust continuum peak. As part of the NOEMA Large Project SOLIS (Seeds of Life in Space), the CH$_3$OH peak has been spatially resolved to study its kinematics and physical structure and to inv… ▽ More Towards the pre-stellar core L1544, the methanol (CH$_3$OH) emission forms an asymmetric ring around the core centre, where CH$_3$OH is mostly in solid form, with a clear peak 4000~au to the north-east of the dust continuum peak. As part of the NOEMA Large Project SOLIS (Seeds of Life in Space), the CH$_3$OH peak has been spatially resolved to study its kinematics and physical structure and to investigate the cause behind the local enhancement. We find that methanol emission is distributed in a ridge parallel to the main axis of the dense core. The centroid velocity increases by about 0.2~km~s$^{-1}$ and the velocity dispersion increases from subsonic to transonic towards the central zone of the core, where the velocity field also shows complex structure. This could be indication of gentle accretion of material onto the core or interaction of two filaments, producing a slow shock. We measure the rotational temperature and show that methanol is in local thermodynamic equilibrium (LTE) only close to the dust peak, where it is significantly depleted. The CH$_3$OH column density, $N_{tot}({\rm CH_3OH})$, profile has been derived with non-LTE radiative transfer modelling and compared with chemical models of a static core. The measured $N_{tot}({\rm CH_3OH})$ profile is consistent with model predictions, but the total column densities are one order of magnitude lower than those predicted by models, suggesting that the efficiency of reactive desorption or atomic hydrogen tunnelling adopted in the model may be overestimated; or that an evolutionary model is needed to better reproduce methanol abundance. △ Less

Submitted 2 February, 2018; originally announced February 2018.

Comments: 19 pages, 18 figures. Accepted by ApJ

Seeds Of Life In Space (SOLIS): The organic composition diversity at 300--1000 au scale in Solar-type star forming regions

Authors: C. Ceccarelli, P. Caselli, F. Fontani, R. Neri, A. Lopez-Sepulcre, C. Codella, S. Feng, I. Jimenez-Serra, B. Lefloch, J. E. Pineda, C. Vastel, F. Alves, R. Bachiller, N. Balucani, E. Bianchi, L. Bizzocchi, S. Bottinelli, E. Caux, A. Chacon-Tanarro, R. Choudhury, A. Coutens, F. Dulieu, C. Favre, P. Hily-Blant, J. Holdship , et al. (21 additional authors not shown)

Abstract: 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 Pro… ▽ More 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. △ Less

Submitted 28 October, 2017; originally announced October 2017.

Comments: Accepted for publication on The Astrophysical Journal

Seeds of Life in Space (SOLIS) III. Formamide in protostellar shocks: evidence for gas-phase formation

Authors: C. Codella, C. Ceccarelli, P. Caselli, N. Balucani, V. Baroneınst, F. Fontani, B. Lefloch, L. Podio, S. Viti, S. Feng, R. Bachiller, E. Bianchi, F. Dulieu, I. Jiménez-Serra, J. Holdship, R. Neri, J. Pineda, A. Pon, I. Sims, S. Spezzano, A. I. Vasyunin, F. Alves, L. Bizzocchi, S. Bottinelli, E. Caux , et al. (25 additional authors not shown)

Abstract: 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 quest… ▽ More 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. △ Less

Submitted 15 August, 2017; originally announced August 2017.

Comments: 7 pages, 1 table, 5 figures, A&A Letters, in press

SOLIS II. Carbon-chain growth in the Solar-type protocluster OMC2-FIR4

Authors: F. Fontani, C. Ceccarelli, C. Favre, P. Caselli, R. Neri, I. R. Sims, C. Kahane, F. Alves, N. Balucani, E. Bianchi, E. Caux, A. Jaber Al-Edhari, A. Lopez-Sepulcre, J. E. Pineda, R. Bachiller, L. Bizzocchi, S. Bottinelli, A. Chacon-Tanarro, R. Choudhury, C. Codella, A. Coutens, F. Dulieu, S. Feng, A. Rimola, P. Hily-Blant , et al. (20 additional authors not shown)

Abstract: 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 galaxi… ▽ More 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 $ζ$ 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 $ζ$ 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. △ Less

Submitted 5 July, 2017; originally announced July 2017.

Comments: 12 pages, 8 figures, accepted for publication in A&A

Formation of Complex Molecules in Prestellar Cores: a Multilayer Approach

Authors: Anton I. Vasyunin, Paola Caselli, François Dulieu, Izaskun Jiménez-Serra

Abstract: We present the results of chemical modeling of complex organic molecules (COMs) under conditions typical for prestellar cores. We utilize an advanced gas-grain astrochemical model with updated gas-phase chemistry, with a multilayer approach to ice-surface chemistry and an up-to-date treatment of reactive desorption based on recent experiments of Minissale et al. (2016). With the chemical model, ra… ▽ More We present the results of chemical modeling of complex organic molecules (COMs) under conditions typical for prestellar cores. We utilize an advanced gas-grain astrochemical model with updated gas-phase chemistry, with a multilayer approach to ice-surface chemistry and an up-to-date treatment of reactive desorption based on recent experiments of Minissale et al. (2016). With the chemical model, radial profiles of molecules including COMs are calculated for the case of the prototypical prestellar core L1544 at the timescales when the modeled depletion factor of CO becomes equal to that observed. We find that COMs can be formed efficiently in L1544 up to the fractional abundances of 10(-10) wrt. total hydrogen nuclei. Abundances of many COMs such as CH3OCH3, HCOOCH3, and others peak at similar radial distances of ~2000-4000 AU. Gas-phase abundances of COMs depend on the efficiency of reactive desorption, which in turn depends on the composition of the outer monolayers of icy mantles. In prestellar cores, the outer monolayers of mantles likely include large fractions of CO and its hydrogenation products, which may increase the efficiency of reactive desorption according to Minissale et al. (2016), and makes the formation of COMs efficient under conditions typical for prestellar cores, although this assumption is yet to be confirmed experimentally. The hydroxyl radical (OH) appears to play an important role in gas-phase chemistry of COMs, which makes it deserving further detailed studies. △ Less

Submitted 12 May, 2017; originally announced May 2017.

Comments: 59 pages, 11 figures, 4 tables. Accepted to The Astrophysical Journal

The Spatial Distribution of Complex Organic Molecules in the L1544 Pre-stellar Core

Authors: Izaskun Jimenez-Serra, Anton I. Vasyunin, Paola Caselli, Nuria Marcelino, Nicolas Billot, Serena Viti, Leonardo Testi, Charlotte Vastel, Bertrand Lefloch, Rafael Bachiller

Abstract: The detection of complex organic molecules (COMs) toward cold sources such as pre-stellar cores (with T<10 K), has challenged our understanding of the formation processes of COMs in the interstellar medium. Recent modelling on COM chemistry at low temperatures has provided new insight into these processes predicting that COM formation depends strongly on parameters such as visual extinction and th… ▽ More The detection of complex organic molecules (COMs) toward cold sources such as pre-stellar cores (with T<10 K), has challenged our understanding of the formation processes of COMs in the interstellar medium. Recent modelling on COM chemistry at low temperatures has provided new insight into these processes predicting that COM formation depends strongly on parameters such as visual extinction and the level of CO freeze out. We report deep observations of COMs toward two positions in the L1544 pre-stellar core: the dense, highly-extinguished continuum peak with Av>=30 mag within the inner 2700 au; and a low-density shell with average Av~7.5-8 mag located at 4000 au from the core's center and bright in CH3OH. Our observations show that CH3O, CH3OCH3 and CH3CHO are more abundant (by factors ~2-10) toward the low-density shell than toward the continuum peak. Other COMs such as CH3OCHO, c-C3H2O, HCCCHO, CH2CHCN and HCCNC show slight enhancements (by factors <=3) but the associated uncertainties are large. This suggests that COMs are actively formed and already present in the low-density shells of pre-stellar cores. The modelling of the chemistry of O-bearing COMs in L1544 indicates that these species are enhanced in this shell because i) CO starts freezing out onto dust grains driving an active surface chemistry; ii) the visual extinction is sufficiently high to prevent the UV photo-dissociation of COMs by the external interstellar radiation field; and iii) the density is still moderate to prevent severe depletion of COMs onto grains. △ Less

Submitted 16 September, 2016; originally announced September 2016.

Comments: 16 pages, 3 figures, 2 tables. Accepted for publication in ApJ Letters

Water deuteration and ortho-to-para nuclear spin ratio of H2 in molecular clouds formed via the accumulation of HI gas

Authors: K. Furuya, Y. Aikawa, U. Hincelin, G. E. Hassel, E. A. Bergin, A. I. Vasyunin, Eric Herbst

Abstract: We investigate the water deuteration ratio and ortho-to-para nuclear spin ratio of H2 (OPR(H2)) during the formation and early evolution of a molecular cloud, following the scenario that accretion flows sweep and accumulate HI gas to form molecular clouds. We follow the physical evolution of post-shock materials using a one-dimensional shock model, with post-processing gas-ice chemistry simulation… ▽ More We investigate the water deuteration ratio and ortho-to-para nuclear spin ratio of H2 (OPR(H2)) during the formation and early evolution of a molecular cloud, following the scenario that accretion flows sweep and accumulate HI gas to form molecular clouds. We follow the physical evolution of post-shock materials using a one-dimensional shock model, with post-processing gas-ice chemistry simulations. This approach allows us to study the evolution of the OPR(H2) and water deuteration ratio without an arbitrary assumption concerning the initial molecular abundances, including the initial OPR(H2). When the conversion of hydrogen into H2 is almost complete, the OPR(H2) is already much smaller than the statistical value of three due to the spin conversion in the gas phase. As the gas accumulates, the OPR(H2) decreases in a non-equilibrium manner. We find that water ice can be deuterium-poor at the end of its main formation stage in the cloud, compared to water vapor observed in the vicinity of low-mass protostars where water ice is likely sublimated. If this is the case, the enrichment of deuterium in water should mostly occur at somewhat later evolutionary stages of star formation, i.e., cold prestellar/protostellar cores. The main mechanism to suppress water ice deuteration in the cloud is the cycle of photodissociation and reformation of water ice, which efficiently removes deuterium from water ice chemistry. The removal efficiency depends on the main formation pathway of water ice. The OPR(H2) plays a minor role in water ice deuteration at the main formation stage of water ice. △ Less

Submitted 23 November, 2015; v1 submitted 17 October, 2015; originally announced October 2015.

Comments: Minor changes after language edition, including typo corrections in Eqs. (A.11) and (B.5) (Accepted for publication in A&A)

Journal ref: A&A 584, A124 (2015)

The HNC/HCN Ratio in Star-Forming Regions

Authors: Dawn Graninger, Eric Herbst, Karin I. Oberg, Anton I. Vasyunin

Abstract: HNC and HCN, typically used as dense gas tracers in molecular clouds, are a pair of isomers that have great potential as a temperature probe because of temperature dependent, isomer-specific formation and destruction pathways. Previous observations of the HNC/HCN abundance ratio show that the ratio decreases with increasing temperature, something that standard astrochemical models cannot reproduce… ▽ More HNC and HCN, typically used as dense gas tracers in molecular clouds, are a pair of isomers that have great potential as a temperature probe because of temperature dependent, isomer-specific formation and destruction pathways. Previous observations of the HNC/HCN abundance ratio show that the ratio decreases with increasing temperature, something that standard astrochemical models cannot reproduce. We have undertaken a detailed parameter study on which environmental characteristics and chemical reactions affect the HNC/HCN ratio and can thus contribute to the observed dependence. Using existing gas and gas-grain models updated with new reactions and reaction barriers, we find that in static models the H + HNC gas-phase reaction regulates the HNC/HCN ratio under all conditions, except for very early times. We quantitively constrain the combinations of H abundance and H + HNC reaction barrier that can explain the observed HNC/HCN temperature dependence and discuss the implications in light of new quantum chemical calculations. In warm-up models, gas-grain chemistry contributes significantly to the predicted HNC/HCN ratio and understanding the dynamics of star formation is therefore key to model the HNC/HCN system. △ Less

Submitted 21 April, 2014; originally announced April 2014.

Comments: 14 pages, 11 figures, Accepted for publication in ApJ

Complex organic molecules in protoplanetary disks

Authors: Catherine Walsh, T. J. Millar, Hideko Nomura, Eric Herbst, Susanna L. Widicus Weaver, Yuri Aikawa, Jake C. Laas, Anton I. Vasyunin

Abstract: (Abridged) Protoplanetary disks are vital objects in star and planet formation, possessing all the material which may form a planetary system orbiting the new star. We investigate the synthesis of complex organic molecules (COMs) in disks to constrain the achievable chemical complexity and predict species and transitions which may be observable with ALMA. We have coupled a 2D model of a protoplane… ▽ More (Abridged) Protoplanetary disks are vital objects in star and planet formation, possessing all the material which may form a planetary system orbiting the new star. We investigate the synthesis of complex organic molecules (COMs) in disks to constrain the achievable chemical complexity and predict species and transitions which may be observable with ALMA. We have coupled a 2D model of a protoplanetary disk around a T Tauri star with a gas-grain chemical network including COMs. We compare compare synthesised line intensities and calculated column densities with observations and determine those COMs which may be observable in future. COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances 1e-6 - 1e-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, 1e-12 - 1e-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 several T Tauri star-disk systems. 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. 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. △ Less

Submitted 3 March, 2014; originally announced March 2014.

Comments: 34 pages, 12 figures, 7 tables

Journal ref: Walsh et al. 2014, A&A, 536, A33 (available here: http://www.aanda.org/articles/aa/abs/2014/03/aa22446-13/aa22446-13.html)

Organic species in Infrared Dark Clouds

Authors: T. Vasyunina, A. I. Vasyunin, E. Herbst, H. Linz, M. Voronkov, T. Britton, I. Zinchenko, F. Schuller

Abstract: It is currently assumed that infrared dark clouds (IRDCs) represent the earliest evolutionary stages of high-mass stars ($>$ 8 M$_{\odot}$). Submillimeter and millimeter-wave studies performed over the past 15 years show that IRDCs possess a broad variety of properties, and hence a wide range of problems and questions that can be tackled. In this paper, we report an investigation of the molecular… ▽ More It is currently assumed that infrared dark clouds (IRDCs) represent the earliest evolutionary stages of high-mass stars ($>$ 8 M$_{\odot}$). Submillimeter and millimeter-wave studies performed over the past 15 years show that IRDCs possess a broad variety of properties, and hence a wide range of problems and questions that can be tackled. In this paper, we report an investigation of the molecular composition and chemical processes in two groups of IRDCs. Using the Mopra, APEX, and IRAM radio telescopes over the last four years, we have collected molecular line data for CO, H$_2$CO, HNCO, CH$_3$CCH, CH$_3$OH, CH$_3$CHO, CH$_3$OCHO, and CH$_3$OCH$_3$. For all of these species we estimated molecular abundances. We then undertook chemical modeling studies, concentrating on the source IRDC028.34+0.06, and compared observed and modeled abundances. This comparison showed that to reproduce observed abundances of complex organic molecules (COMs), a 0-D gas-grain model with constant physical conditions is not sufficient. We achieved greater success with the use of a warm-up model, in which warm-up from 10 K to 30 K occurs following a cold phase. △ Less

Submitted 30 October, 2013; originally announced October 2013.

Comments: accepted for publication in ApJ

Reactive Desorption and Radiative Association as Possible Drivers of Complex Molecule Formation in the Cold Interstellar Medium

Authors: A. I. Vasyunin, Eric Herbst

Abstract: The recent discovery of terrestrial-type organic species such as methyl formate and dimethyl ether in the cold interstellar gas has proved that the formation of organic matter in the Galaxy begins at a much earlier stage of star formation than was thought before. This discovery represents a challenge for astrochemical modelers. The abundances of these molecules cannot be explained by the previousl… ▽ More The recent discovery of terrestrial-type organic species such as methyl formate and dimethyl ether in the cold interstellar gas has proved that the formation of organic matter in the Galaxy begins at a much earlier stage of star formation than was thought before. This discovery represents a challenge for astrochemical modelers. The abundances of these molecules cannot be explained by the previously developed "warm-up" scenario, in which organic molecules are formed via diffusive chemistry on surfaces of interstellar grains starting at 30 K, and then released to the gas at higher temperatures during later stages of star formation. In this article, we investigate an alternative scenario in which complex organic species are formed via a sequence of gas-phase reactions between precursor species formed on grain surfaces and then ejected into the gas via efficient reactive desorption, a process in which non-thermal desorption occurs as a result of conversion of the exothermicity of chemical reactions into the ejection of products from the surface. The proposed scenario leads to reasonable if somewhat mixed results at temperatures as low as 10 K and may be considered as a step towards the explanation of abundances of terrestrial-like organic species observed during the earliest stages of star formation. △ Less

Submitted 15 April, 2013; v1 submitted 28 March, 2013; originally announced March 2013.

Comments: 23 pages, 5 figures, accepted to ApJ

A Unified Monte Carlo Treatment of Gas-Grain Chemistry for Large Reaction Networks. II. A Multiphase Gas-Surface-Layered Bulk Model

Authors: A. I. Vasyunin, E. Herbst

Abstract: The observed gas-phase molecular inventory of hot cores is believed to be significantly impacted by the products of chemistry in interstellar ices. In this study, we report the construction of a full macroscopic Monte Carlo model of both the gas-phase chemistry and the chemistry occurring in the icy mantles of interstellar grains. Our model treats icy grain mantles in a layer-by-layer manner, whic… ▽ More The observed gas-phase molecular inventory of hot cores is believed to be significantly impacted by the products of chemistry in interstellar ices. In this study, we report the construction of a full macroscopic Monte Carlo model of both the gas-phase chemistry and the chemistry occurring in the icy mantles of interstellar grains. Our model treats icy grain mantles in a layer-by-layer manner, which incorporates laboratory data on ice desorption correctly. The ice treatment includes a distinction between a reactive ice surface and an inert bulk. The treatment also distinguishes between zeroth and first order desorption, and includes the entrapment of volatile species in more refractory ice mantles. We apply the model to the investigation of the chemistry in hot cores, in which a thick ice mantle built up during the previous cold phase of protostellar evolution undergoes surface reactions and is eventually evaporated. For the first time, the impact of a detailed multilayer approach to grain mantle formation on the warm-up chemistry is explored. The use of a multilayer ice structure has a mixed impact on the abundances of organic species formed during the warm-up phase. For example, the abundance of gaseous HCOOCH3 is lower in the multilayer model than in previous grain models that do not distinguish between layers (so-called "two phase" models). Other gaseous organic species formed in the warm-up phase are affected slightly. Finally, we find that the entrapment of volatile species in water ice can explain the two-jump behavior of H2CO previously found in observations of protostars. △ Less

Submitted 14 November, 2012; v1 submitted 13 November, 2012; originally announced November 2012.

Comments: 50 pages, 14 figures, accepted to ApJ

Chemical modeling of Infrared Dark Clouds: the Role of Surface Chemistry

Authors: T. Vasyunina, A. I. Vasyunin, E. Herbst, H. Linz

Abstract: We simulate the chemistry of infrared dark clouds (IRDCs) with a model in which the physical conditions are homogeneous and time-independent. The chemistry is solved as a function of time with three networks: one purely gas-phase, one that includes accretion and desorption, and one, the complete gas-grain network, that includes surface chemistry in addition. We compare our results with observed mo… ▽ More We simulate the chemistry of infrared dark clouds (IRDCs) with a model in which the physical conditions are homogeneous and time-independent. The chemistry is solved as a function of time with three networks: one purely gas-phase, one that includes accretion and desorption, and one, the complete gas-grain network, that includes surface chemistry in addition. We compare our results with observed molecular abundances for two representative IRDCs -- IRDC013.90-1 and IRDC321.73-1 -- using the molecular species N$_2$H$^+$, HC$_3$N, HNC, HCO$^+$, HCN, C$_2$H, NH$_3$ and CS. IRDC013.90-1 is a cold IRDC, with a temperature below 20 K, while IRDC321.73-1 is somewhat warmer, in the range 20 - 30 K. We find that the complete gas-grain model fits the data very well, but that the goodness-of-fit is not sharply peaked at a particular temperature. Surface processes are important for the explanation of the high gas-phase abundance of N$_2$H$^+$ in IRDC321.73-1. The general success of the 0-D model in reproducing single-dish observations of our limited sample of 8 species shows that it is probably sufficient for an explanation of this type of data. To build and justify more complicated models, including spatial temperature and density structure, contraction, and heating, we require high-resolution interferometric data. △ Less

Submitted 25 March, 2012; originally announced March 2012.

Comments: accepted for publication in ApJ

Discovery of interstellar mercapto radicals (SH) with the GREAT instrument on SOFIA

Authors: D. A. Neufeld, E. Falgarone, M. Gerin, B. Godard, E. Herbst, G. Pineau des Forêts, A. I. Vasyunin, R. Güsten, H. Wiesemeyer, O. Ricken

Abstract: We report the first detection of interstellar mercapto radicals, obtained along the sight-line to the submillimeter continuum source W49N. We have used the GREAT instrument on SOFIA to observe the 1383 GHz Doublet Pi 3/2 J = 5/2 - 3/2 lambda doublet in the upper sideband of the L1 receiver. The resultant spectrum reveals SH absorption in material local to W49N, as well as in foreground gas, unasso… ▽ More We report the first detection of interstellar mercapto radicals, obtained along the sight-line to the submillimeter continuum source W49N. We have used the GREAT instrument on SOFIA to observe the 1383 GHz Doublet Pi 3/2 J = 5/2 - 3/2 lambda doublet in the upper sideband of the L1 receiver. The resultant spectrum reveals SH absorption in material local to W49N, as well as in foreground gas, unassociated with W49N, that is located along the sight-line. For the foreground material at velocities in the range 37 - 44 km/s with respect to the local standard of rest, we infer a total SH column density ~ 2.6 E+12 cm-2, corresponding to an abundance of ~ 7 E-9 relative to H2, and yielding an SH/H2S abundance ratio ~ 0.13. The observed SH/H2S abundance ratio is much smaller than that predicted by standard models for the production of SH and H2S in turbulent dissipation regions and shocks, and suggests that the endothermic neutral-neutral reaction SH + H2 -> H2S + H must be enhanced along with the ion-neutral reactions believed to produce CH+ and SH+ in diffuse molecular clouds. △ Less

Submitted 14 February, 2012; originally announced February 2012.

Comments: Accepted for publication in Astronomy and Astrophysics (SOFIA/GREAT special issue)

New extended deuterium fractionation model: assessment at dense ISM conditions and sensitivity analysis

Authors: T. Albertsson, D. A. Semenov, A. I. Vasyunin, Th. Henning, E. Herbst

Abstract: Observations of deuterated species are useful in probing the temperature, ionization level, evolutionary stage, chemistry, and thermal history of astrophysical environments. The analysis of data from ALMA and other new telescopes requires an elaborate model of deuterium fractionation. This paper presents a publicly available chemical network with multi-deuterated species and an extended, up-to-dat… ▽ More Observations of deuterated species are useful in probing the temperature, ionization level, evolutionary stage, chemistry, and thermal history of astrophysical environments. The analysis of data from ALMA and other new telescopes requires an elaborate model of deuterium fractionation. This paper presents a publicly available chemical network with multi-deuterated species and an extended, up-to-date set of gas-phase and surface reactions. To test this network, we simulate deuterium fractionation in diverse interstellar sources. Two cases of initial abundances are considered: i) atomic except for H2 and HD, and ii) molecular from a prestellar core. We reproduce the observed D/H ratios of many deuterated molecules, and sort the species according to their sensitivity to temperature gradients and initial abundances. We find that many multiply-deuterated species produced at 10 K retain enhanced D/H ratios at temperatures $\la 100$ K. We study how recent updates to reaction rates affect calculated D/H ratios, and perform a detailed sensitivity analysis of the uncertainties of the gas-phase reaction rates in the network. We find that uncertainties are generally lower in dark cloud environments than in warm IRDCs and that uncertainties increase with the size of the molecule and number of D-atoms. A set of the most problematic reactions is presented. We list potentially observable deuterated species predicted to be abundant in low- and high-mass star-formation regions. △ Less

Submitted 4 June, 2013; v1 submitted 12 October, 2011; originally announced October 2011.

Comments: 30 pages, 8 figures, accepted for publication in ApJS

Journal ref: ApJS, 207 (2013), 27

Impact of grain evolution on the chemical structure of protoplanetary disks

Authors: A. I. Vasyunin, D. S. Wiebe, T. Birnstiel, S. Zhukovska, Th. Henning, C. P. Dullemond

Abstract: We study the impact of dust evolution in a protoplanetary disk around a T Tauri star on the disk chemical composition. For the first time we utilize a comprehensive model of dust evolution which includes growth, fragmentation and sedimentation. Specific attention is paid to the influence of grain evolution on the penetration of the UV field in the disk. A chemical model that includes a comprehensi… ▽ More We study the impact of dust evolution in a protoplanetary disk around a T Tauri star on the disk chemical composition. For the first time we utilize a comprehensive model of dust evolution which includes growth, fragmentation and sedimentation. Specific attention is paid to the influence of grain evolution on the penetration of the UV field in the disk. A chemical model that includes a comprehensive set of gas phase and grain surface chemical reactions is used to simulate the chemical structure of the disk. The main effect of the grain evolution on the disk chemical composition comes from sedimentation, and, to a lesser degree, from the reduction of the total grain surface area. The net effect of grain growth is suppressed by the fragmentation process which maintains a population of small grains, dominating the total grain surface area. We consider three models of dust properties. In model GS both growth and sedimentation are taken into account. In models A5 and A4 all grains are assumed to have the same size (10(-5) cm and 10(-4) cm, respectively) with constant gas-to-dust mass ratio of 100. Like in previous studies, the "three-layer" pattern (midplane, molecular layer, hot atmosphere) in the disk chemical structure is preserved in all models, but shifted closer to the midplane in models with increased grain size (GS and A4). Unlike other similar studies, we find that in models GS and A4 column densities of most gas-phase species are enhanced by 1-3 orders of magnitude relative to those in a model with pristine dust (A5), while column densities of their surface counterparts are decreased. We show that column densities of certain species, like C2H, HC(2n+1)N (n=0-3), H2O and some other molecules, as well as the C2H2/HCN abundance ratio which are accessible with Herschel and ALMA can be used as observational tracers of early stages of the grain evolution process in protoplanetary disks. △ Less

Submitted 19 November, 2010; originally announced November 2010.

Comments: 50 pages, 4 tables, 11 figures, accepted to the ApJ

A Unified Monte Carlo Treatment of Gas-Grain Chemistry for Large Reaction Networks. I. Testing Validity of Rate Equations in Molecular Clouds

Authors: A. I. Vasyunin, D. A. Semenov, D. S. Wiebe, Th. Henning

Abstract: In this study we demonstrate for the first time that the unified Monte Carlo approach can be applied to model gas-grain chemistry in large reaction networks. Specifically, we build a time-dependent gas-grain chemical model of the interstellar medium, involving about 6000 gas-phase and 200 grain surface reactions. This model is used to test the validity of the standard and modified rate equation… ▽ More In this study we demonstrate for the first time that the unified Monte Carlo approach can be applied to model gas-grain chemistry in large reaction networks. Specifically, we build a time-dependent gas-grain chemical model of the interstellar medium, involving about 6000 gas-phase and 200 grain surface reactions. This model is used to test the validity of the standard and modified rate equation methods in models of dense and translucent molecular clouds and to specify under which conditions the use of the stochastic approach is desirable. We found that at temperatures 25--30 K gas-phase abundances of H$_2$O, NH$_3$, CO and many other gas-phase and surface species in the stochastic model differ from those in the deterministic models by more than an order of magnitude, at least, when tunneling is accounted for and/or diffusion energies are 3x lower than the binding energies. In this case, surface reactions, involving light species, proceed faster than accretion of the same species. In contrast, in the model without tunneling and with high binding energies, when the typical timescale of a surface recombination is greater than the timescale of accretion onto the grain, we obtain almost perfect agreement between results of Monte Carlo and deterministic calculations in the same temperature range. At lower temperatures ($\sim10$ K) gaseous and, in particular, surface abundances of most important molecules are not much affected by stochastic processes. △ Less

Submitted 9 October, 2008; originally announced October 2008.

Comments: 33 pages, 9 figures, 1 table. Accepted for publication in ApJ

Journal ref: Astrophys.J.691:1459-1469,2009

Chemistry in Protoplanetary Disks: A Sensitivity Analysis

Authors: A. I. Vasyunin, D. Semenov, Th. Henning, V. Wakelam, E. Herbst, A. M. Sobolev

Abstract: We study how uncertainties in the rate coefficients of chemical reactions in the RATE06 database affect abundances and column densities of key molecules in protoplanetary disks. We randomly varied the gas-phase reaction rates within their uncertainty limits and calculated the time-dependent abundances and column densities using a gas-grain chemical model and a flaring steady-state disk model. We… ▽ More We study how uncertainties in the rate coefficients of chemical reactions in the RATE06 database affect abundances and column densities of key molecules in protoplanetary disks. We randomly varied the gas-phase reaction rates within their uncertainty limits and calculated the time-dependent abundances and column densities using a gas-grain chemical model and a flaring steady-state disk model. We find that key species can be separated into two distinct groups according to the sensitivity of their column densities to the rate uncertainties. The first group includes CO, C$^+$, H$_3^+$, H$_2$O, NH$_3$, N$_2$H$^+$, and HCNH$^+$. For these species, the column densities are not very sensitive to the rate uncertainties but the abundances in specific regions are. The second group includes CS, CO$_2$, HCO$^+$, H$_2$CO, C$_2$H, CN, HCN, HNC and other, more complex species, for which high abundances and abundance uncertainties co-exist in the same disk region, leading to larger scatters in the column densities. However, even for complex and heavy molecules, the dispersion in their column densities is not more than a factor of ~4. We perform a sensitivity analysis of the computed abundances to rate uncertainties and identify those reactions with the most problematic rate coefficients. We conclude that the rate coefficients of about a hundred of chemical reactions need to be determined more accurately in order to greatly improve the reliability of modern astrochemical models. This improvement should be an ultimate goal for future laboratory studies and theoretical investigations. △ Less

Submitted 20 September, 2007; originally announced September 2007.

Comments: 42 pages, 6 figures, 3 tables, accepted for publication in ApJ

On the Influence of Uncertainties in Chemical Reaction Rates on Results of the Astrochemical Modelling

Authors: A. I. Vasyunin, A. M. Sobolev, D. S. Wiebe, D. A. Semenov

Abstract: With the chemical reaction rate database UMIST95 (Millar et al. 1997) we analyze how uncertainties in rate constants of gas-phase chemical reactions influence the modelling of molecular abundances in the interstellar medium. Random variations are introduced into the rate constants to estimate the scatter in theoretical abundances. Calculations are performed for dark and translucent molecular clo… ▽ More With the chemical reaction rate database UMIST95 (Millar et al. 1997) we analyze how uncertainties in rate constants of gas-phase chemical reactions influence the modelling of molecular abundances in the interstellar medium. Random variations are introduced into the rate constants to estimate the scatter in theoretical abundances. Calculations are performed for dark and translucent molecular clouds where gas phase chemistry is adequate. Similar approach was used by Pineau des Forets & Roueff (2000) for the study of chemical bistability. All the species are divided into 6 sensitivity groups according to the value of the scatter in their model abundances computed with varied rate constants. It is shown that the distribution of species within these groups depends on the number of atoms in a molecule and on the adopted physical conditions. The simple method is suggested which allows to single out reactions that are most important for the evolution of a given species. △ Less

Submitted 19 November, 2003; originally announced November 2003.

Comments: 4 pages. To appear in the proceedings of the 4th Cologne-Bonn Zermatt Symposium

Journal ref: Astron.Lett. 30 (2004) 566-576