-
AB Aur, a Rosetta stone for studies of planet formation (III): continuum observations at 2 and 7 mm
Authors:
Pablo Rivière-Marichalar,
Enrique Macías,
Clément Baruteau,
Asunción Fuente,
Roberto Neri,
Álvaro Ribas,
Gisela Esplugues,
David Navarro-Almaida,
Mayra Osorio,
Guillem Anglada
Abstract:
Observational constraints on dust properties in protoplanetary disks are key to better understanding disks' evolution. We continue our exploration of the protoplanetary disk around AB Aur by characterizing its dust properties. We present ALMA observations at 2.2 mm and VLA observations at 6.8 mm. Together with previous ALMA and NOEMA observations at 0.87 and 1.1 mm, these observations are used to…
▽ More
Observational constraints on dust properties in protoplanetary disks are key to better understanding disks' evolution. We continue our exploration of the protoplanetary disk around AB Aur by characterizing its dust properties. We present ALMA observations at 2.2 mm and VLA observations at 6.8 mm. Together with previous ALMA and NOEMA observations at 0.87 and 1.1 mm, these observations are used to compute global spectral index profiles as well as spectral index maps. On the interpretation side, we present the results of a simple isothermal slab model to help constrain dust properties along the ring of continuum emission. We also present results of dust radiative transfer calculations based on a disc-planet hydrodynamical simulation to explain how the azimuthal contrast ratio of the ring emission varies with millimeter wavelength. The spectral energy distribution and the radial profiles of the spectral index indicate that the radiation from the compact source towards the center is not dominated by dust thermal emission, but most likely by free-free emission originating in the radio jet: it constitutes 93% of the emission at 6.8 mm, and 37% at 0.87 mm. The protoplanetary disk has a typical spectral index of 2.3. We estimate a dust disk mass of 8$\rm \times 10^{-5}$ M$\rm _{\odot}$ which, assuming a mean gas-to-dust ratio of 40, gives a total disk mass of 3.2$\rm \times 10^{-3}$ M$\rm _{\odot}$. The azimuthal contrast ratio of the ring outside the millimeter cavity is smaller at 2.2 mm than at 1.1 mm, in agreement with previous findings. The VLA image shows several knots of $5σ$ emission all along the ring, which are consistent with the ring emission being nearly axisymmetric at that wavelength. The decrease in the azimuthal contrast ratio of the ring emission from 0.87 to 6.8 mm can be explained by a dust-losing decaying vortex at the outer edge of a planet gap.
△ Less
Submitted 20 December, 2023;
originally announced December 2023.
-
Evolution of Chemistry in the envelope of Hot Corinos (ECHOS). I. Extremely young sulphur chemistry in the isolated Class 0 object B335
Authors:
G. Esplugues,
M. Rodríguez-Baras,
D. San Andrés,
D. Navarro-Almaida,
A. Fuente,
P. Rivière-Marichalar,
Á. Sánchez-Monge,
M. N. Drozdovskaya,
S. Spezzano,
P. Caselli
Abstract:
Within the project Evolution of Chemistry in the envelope of HOt corinoS (ECHOS), we present a study of sulphur chemistry in the envelope of the Class 0 source B335 through observations in the spectral range 7, 3, and 2 mm. We have modelled observations assuming LTE and LVG approximation. We have also used the code Nautilus to study the time evolution of sulphur species. We have detected 20 sulphu…
▽ More
Within the project Evolution of Chemistry in the envelope of HOt corinoS (ECHOS), we present a study of sulphur chemistry in the envelope of the Class 0 source B335 through observations in the spectral range 7, 3, and 2 mm. We have modelled observations assuming LTE and LVG approximation. We have also used the code Nautilus to study the time evolution of sulphur species. We have detected 20 sulphur species with a total gas-phase S abundance similar to that found in the envelopes of other Class 0 objects, but with significant differences in the abundances between sulphur carbon chains and sulphur molecules containing oxygen and nitrogen. Our results highlight the nature of B335 as a source especially rich in sulphur carbon chains unlike other Class 0 sources. The low presence or absence of some molecules, such as SO and SO+, suggests a chemistry not particularly influenced by shocks. We, however, detect a large presence of HCS+ that, together with the low rotational temperatures obtained for all the S species (<15 K), reveals the moderate or low density of the envelope of B335. We also find that observations are better reproduced by models with a sulphur depletion factor of 10 with respect to the sulphur cosmic elemental abundance. The comparison between our model and observational results for B335 reveals an age of 10$^4$$<$t$<$10$^5$ yr, which highlights the particularly early evolutionary stage of this source. B335 presents a different chemistry compared to other young protostars that have formed in dense molecular clouds, which could be the result of accretion of surrounding material from the diffuse cloud onto the protostellar envelope of B335. In addition, the analysis of the SO2/C2S, SO/CS, and HCS+/CS ratios within a sample of prestellar cores and Class 0 objects show that they could be used as good chemical evolutionary indicators of the prestellar to protostellar transition.
△ Less
Submitted 4 September, 2023;
originally announced September 2023.
-
Gas phase Elemental abundances in Molecular cloudS (GEMS). IX. Deuterated compounds of H2S in starless cores
Authors:
Marina Rodríguez-Baras,
Gisela Esplugues,
Asunción Fuente,
Silvia Spezzano,
Paola Caselli,
Jean-Christophe Loison,
Evelyne Roueff,
David Navarro-Almaida,
Rafael Bachiller,
Rafael Martín-Doménech,
Izaskun Jiménez-Serra,
Leire Beitia-Antero,
Romane Le Gal
Abstract:
H2S is thought to be the main sulphur reservoir in the ice, being therefore a key molecule to understand sulphur chemistry in the star formation process and to solve the missing sulphur problem. The H2S deuterium fraction can be used to constrain its formation pathways. We investigate for the first time the H2S deuteration in a large sample of starless cores (SC). We use observations of the GEMS I…
▽ More
H2S is thought to be the main sulphur reservoir in the ice, being therefore a key molecule to understand sulphur chemistry in the star formation process and to solve the missing sulphur problem. The H2S deuterium fraction can be used to constrain its formation pathways. We investigate for the first time the H2S deuteration in a large sample of starless cores (SC). We use observations of the GEMS IRAM 30m Large Program and complementary IRAM 30m observations. We consider a sample of 19 SC in Taurus, Perseus, and Orion, detecting HDS in 10 and D2S in five. The H2S single and double deuterium fractions are analysed with regard to their relation with the cloud physical parameters, their comparison with other interstellar sources, and their comparison with deuterium fractions in early stage star-forming sources of c-C3H2, H2CS, H2O, H2CO, and CH3OH. We obtain a range of X(HDS)/X(H2S)~0.025-0.2 and X(D2S)/X(HDS)~0.05-0.3. H2S single deuteration shows an inverse relation with the cloud kinetic temperature. H2S deuteration values in SC are similar to those observed in Class 0. Comparison with other molecules in other sources reveals a general trend of decreasing deuteration with increasing temperature. In SC and Class 0 objects H2CS and H2CO present higher deuteration fractions than c-C3H2, H2S, H2O, and CH3OH. H2O shows single and double deuteration values one order of magnitude lower than those of H2S and CH3OH. Differences between c-C3H2, H2CS and H2CO deuterium fractions and those of H2S, H2O, and CH3OH are related to deuteration processes produced in gas or solid phases, respectively. We interpret the differences between H2S and CH3OH deuterations and that of H2O as a consequence of differences on the formation routes in the solid phase, particularly in terms of the different occurrence of the D-H and H-D substitution reactions in the ice, together with the chemical desorption processes.
△ Less
Submitted 1 September, 2023;
originally announced September 2023.
-
Gas phase Elemental abundances in Molecular cloudS (GEMS) VIII. Unlocking the CS chemistry: the CH + S$\rightarrow$ CS + H and C$_2$ + S$\rightarrow$ CS + C reactions
Authors:
Carlos M. R. Rocha,
Octavio Roncero,
Niyazi Bulut,
Piotr Zuchowski,
David Navarro-Almaida,
Asuncion Fuente,
Valentine Wakelam,
Jean-Christophe Loison,
Evelyne Roueff,
Javier R. Goicoechea,
Gisela Esplugues,
Leire Beitia-Antero,
Paola Caselli,
Valerio Lattanzi,
Jaime Pineda,
Romane Le Gal,
Marina Rodriguez-Baras,
Pablo Riviere-Marichalar
Abstract:
We revise the rates of reactions CH + S -> CS + H and C_2 + S -> CS + C, important CS formation routes in dark and diffuse warm gas. We performed ab initio calculations to characterize the main features of all the electronic states correlating to the open shell reactants. For CH+S we have calculated the full potential energy surfaces for the lowest doublet states and the reaction rate constant wit…
▽ More
We revise the rates of reactions CH + S -> CS + H and C_2 + S -> CS + C, important CS formation routes in dark and diffuse warm gas. We performed ab initio calculations to characterize the main features of all the electronic states correlating to the open shell reactants. For CH+S we have calculated the full potential energy surfaces for the lowest doublet states and the reaction rate constant with a quasi-classical method. For C_2+S, the reaction can only take place through the three lower triplet states, which all present deep insertion wells. A detailed study of the long-range interactions for these triplet states allowed to apply a statistic adiabatic method to determine the rate constants. This study of the CH + S reaction shows that its rate is nearly independent on the temperature in a range of 10-500 K with an almost constant value of 5.5 10^{-11} cm^3/s at temperatures above 100~K. This is a factor \sim 2-3 lower than the value obtained with the capture model. The rate of the reaction C_2 + S depends on the temperature taking values close to 2.0 10^{-10} cm^3/s at low temperatures and increasing to 5. 10^{-10} cm^3/s for temperatures higher than 200~K. Our modeling provides a rate higher than the one currently used by factor of \sim 2. These reactions were selected for involving open-shell species with many degenerate electronic states, and the results obtained in the present detailed calculations provide values which differ a factor of \sim 2-3 from the simpler classical capture method. We have updated the sulphur network with these new rates and compare our results in the prototypical case of TMC1 (CP). We find a reasonable agreement between model predictions and observations with a sulphur depletion factor of 20 relative to the sulphur cosmic abundance, but it is not possible to fit all sulphur-bearing molecules better than a factor of 10 at the same chemical time.
△ Less
Submitted 1 July, 2023;
originally announced July 2023.
-
Linking the dust and chemical evolution: Taurus and Perseus -- New collisional rates for HCN, HNC, and their C, N, and H isotopologues
Authors:
D. Navarro-Almaida,
C. T. Bop,
F. Lique,
G. Esplugues,
M. Rodríguez-Baras,
C. Kramer,
C. E. Romero,
A. Fuente,
P. Caselli,
P. Riviére-Marichalar,
J. M. Kirk,
A. Chacón-Tanarro,
E. Roueff,
T. Mroczkowski,
T. Bhandarkar,
M. Devlin,
S. Dicker,
I. Lowe,
B. Mason,
C. L. Sarazin,
J. Sievers
Abstract:
HCN, HNC, and their isotopologues are ubiquitous molecules that can serve as chemical thermometers and evolutionary tracers to characterize star-forming regions. Despite their importance in carrying information that is vital to studies of the chemistry and evolution of star-forming regions, the collision rates of some of these molecules have not been available for rigorous studies in the past. We…
▽ More
HCN, HNC, and their isotopologues are ubiquitous molecules that can serve as chemical thermometers and evolutionary tracers to characterize star-forming regions. Despite their importance in carrying information that is vital to studies of the chemistry and evolution of star-forming regions, the collision rates of some of these molecules have not been available for rigorous studies in the past. We perform an up-to-date gas and dust chemical characterization of two different star-forming regions, TMC 1-C and NGC 1333-C7, using new collisional rates of HCN, HNC, and their isotopologues. We investigated the possible effects of the environment and stellar feedback in their chemistry and their evolution. With millimeter observations, we derived their column densities, the C and N isotopic fractions, the isomeric ratios, and the deuterium fractionation. The continuum data at 3 mm and 850 $μ$m allowed us to compute the emissivity spectral index and look for grain growth as an evolutionary tracer. The H$^{13}$CN/HN$^{13}$C ratio is anticorrelated with the deuterium fraction of HCN, thus it can readily serve as a proxy for the temperature. The spectral index $(β\sim 1.34-2.09)$ shows a tentative anticorrelation with the H$^{13}$CN/HN$^{13}$C ratio, suggesting grain growth in the evolved, hotter, and less deuterated sources. Unlike TMC 1-C, the south-to-north gradient in dust temperature and spectral index observed in NGC 1333-C7 suggests feedback from the main NGC 1333 cloud. With this up-to-date characterization of two star-forming regions, we found that the chemistry and the physical properties are tightly related. The dust temperature, deuterium fraction, and the spectral index are complementary evolutionary tracers. The large-scale environmental factors may dominate the chemistry and evolution in clustered star-forming regions.
△ Less
Submitted 15 December, 2022;
originally announced December 2022.
-
Gas phase Elemental abundances in Molecular cloudS (GEMS) VII. Sulfur elemental abundance
Authors:
A. Fuente,
P. Rivière-Marichalar,
L. Beitia-Antero,
P. Caselli,
V. Wakelam,
G. Esplugues,
M. Rodríguez-Baras,
D. Navarro-Almaida,
M. Gerin,
C. Kramer,
R. Bachiller,
J. R. Goicoechea,
I. Jiménez-Serra,
J. C. Loison,
A. Ivlev,
R. Martín-Doménech,
S. Spezzano,
O. Roncero,
G. Muñoz-Caro,
S. Cazaux,
N. Marcelino
Abstract:
Gas phase Elemental abundances in molecular CloudS (GEMS) is an IRAM 30m large program aimed at determining the elemental abundances of carbon (C), oxygen (O), nitrogen (N), and sulfur (S) in a selected set of prototypical star-forming filaments. In particular, the elemental abundance of S remains uncertain by several orders of magnitude and its determination is one of the most challenging goals o…
▽ More
Gas phase Elemental abundances in molecular CloudS (GEMS) is an IRAM 30m large program aimed at determining the elemental abundances of carbon (C), oxygen (O), nitrogen (N), and sulfur (S) in a selected set of prototypical star-forming filaments. In particular, the elemental abundance of S remains uncertain by several orders of magnitude and its determination is one of the most challenging goals of this program. We have carried out an extensive chemical modeling of the fractional abundances of CO, HCO$^+$, HCN, HNC, CS, SO, H$_2$S, OCS, and HCS$^+$ to determine the sulfur depletion toward the 244 positions in the GEMS database. These positions sample visual extinctions from A$_V$ $\sim$ 3 mag to $>$50 mag, molecular hydrogen densities ranging from a few 10$^3$~cm$^{-3}$ to 3$\times$10$^6$~cm$^{-3}$, and T$_k$ $\sim$ 10$-$35 K. Most of the positions in Taurus and Perseus are best fitted assuming early-time chemistry, t=0.1 Myr, $ζ_{H_2}$$\sim$ (0.5$-$1)$\times$10$^{-16}$ s$^{-1}$, and [S/H]$\sim$1.5$\times$10$^{-6}$. On the contrary, most of the positions in Orion are fitted with t=1~Myr and $ζ_{H_2}$$\sim$ 10$^{-17}$ s$^{-1}$. Moreover, $\sim$40% of the positions in Orion are best fitted assuming the undepleted sulfur abundance, [S/H]$\sim$1.5$\times$10$^{-5}$. Our results suggest that sulfur depletion depends on the environment. While the abundances of sulfur-bearing species are consistent with undepleted sulfur in Orion, a depletion factor of $\sim$20 is required to explain those observed in Taurus and Perseus. We propose that differences in the grain charge distribution in the envelopes of the studied clouds might explain these variations. The shocks associated with past and ongoing star formation could also contribute to enhance [S/H] in Orion.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
AB Aur, a Rosetta stone for studies of planet formation (II): H$_2$S detection and sulfur budget
Authors:
Pablo Rivière-Marichalar,
Asunción Fuente,
Gisela Esplugues,
Valentine Wakelam,
Romane le Gal,
Clément Baruteau,
Álvaro Ribas,
Enrique Macías,
Roberto Neri,
David Navarro-Almaida
Abstract:
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 ph…
▽ More
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.
△ Less
Submitted 14 July, 2022;
originally announced July 2022.
-
Gas phase Elemental abundances in Molecular cloudS (GEMS) VI. A sulphur journey across star-forming regions: study of thioformaldehyde emission
Authors:
G. Esplugues,
A. Fuente,
D. Navarro-Almaida,
M. Rodriguez-Baras,
L. Majumdar,
P. Caselli,
V. Wakelam,
E. Roueff,
R. Bachiller,
S. Spezzano,
P. Riviere-Marichalar,
R. Martin-Domenech,
G. M. Muñoz Caro
Abstract:
In the context of the IRAM 30m Large Program GEMS, we present a study of thioformaldehyde in several starless cores located in star-forming filaments of Taurus, Perseus, and Orion. We investigate the influence of the environmental conditions on the abundances of these molecules in the cores, and the effect of time evolution. We have modelled the observed lines of H2CS, HDCS, and D2CS using the rad…
▽ More
In the context of the IRAM 30m Large Program GEMS, we present a study of thioformaldehyde in several starless cores located in star-forming filaments of Taurus, Perseus, and Orion. We investigate the influence of the environmental conditions on the abundances of these molecules in the cores, and the effect of time evolution. We have modelled the observed lines of H2CS, HDCS, and D2CS using the radiative transfer code RADEX. We have also used the chemical code Nautilus to model the evolution of these species depending on the characteristics of the starless cores. We derive column densities and abundances for all the cores. We also derive deuterium fractionation ratios, Dfrac, to determine and compare the evolutionary stage between different parts of each star-forming region. Our results indicate that the north region of the B213 filament in Taurus is more evolved than the south, while the north-eastern part of Perseus presents an earlier evolutionary stage than the south-western zone. Model results also show that Dfrac decreases with the cosmic-ray ionisation rate, while it increases with density and with the degree of sulphur depletion. In particular, we only reproduce the observations when the initial sulphur abundance in the starless cores is at least one order of magnitude lower than the solar elemental sulphur abundance. The progressive increase in HDCS/H2CS and D2CS/H2CS with time makes these ratios powerful tools for deriving the chemical evolutionary stage of starless cores. However, they cannot be used to derive the temperature of these regions, since both ratios present a similar evolution at two different temperature ranges (7-11 K and 15-19 K). Regarding chemistry, (deuterated) thioformaldehyde is mainly formed through gas-phase reactions (double-replacement and neutral-neutral displacement reactions), while surface chemistry plays an important role as a destruction mechanism.
△ Less
Submitted 6 April, 2022;
originally announced April 2022.
-
H2CS deuteration maps towards the pre-stellar core L1544
Authors:
S. Spezzano,
O. Sipilä,
P. Caselli,
S. S. Jensen,
S. Czakli,
L. Bizzocchi,
J. Chantzos,
G. Esplugues,
A. Fuente,
F. Eisenhauer
Abstract:
Deuteration is a crucial tool to understand the complexity of interstellar chemical processes, especially when they involve the interplay of gas-phase and grain-surface chemistry. In the case of multiple deuteration, comparing observation with the results of chemical modelling is particularly effective to study how molecules are inherited in the different stages within the process of star and plan…
▽ More
Deuteration is a crucial tool to understand the complexity of interstellar chemical processes, especially when they involve the interplay of gas-phase and grain-surface chemistry. In the case of multiple deuteration, comparing observation with the results of chemical modelling is particularly effective to study how molecules are inherited in the different stages within the process of star and planet formation. We aim to study the the D/ H ratio in H2CS across the prototypical pre-stellar core L1544. This study allows us to test current gas-dust chemical models involving sulfur in dense cores. We present here single-dish observations of H2CS, HDCS and D2CS with the IRAM 30m telescope. We analyse their column densities and distributions, and compare these observations with gas-grain chemical models. The deuteration maps of H2CS in L1544 are compared with the deuteration maps of methanol, H2CO, N2H+ and HCO+ towards the same source. Furthermore, the single and double deuteration of H2CS towards the dust peak of L1544 is compared with H2CO and c-C3H2. The difference between the deuteration of these molecules in L1544 is discussed and compared with the prediction of chemical models. The maximum deuterium fractionation for the first deuteration of H2CS is N(HDCS)/N(H2CS)$\sim$30$\%$ and is located towards the north-east at a distance of about 10000 AU from the dust peak. While for c-C3H2 the first and second deuteration have a similar efficiency, for H2CS and H2CO the second deuteration is more efficient, leading to D2CX/HDCX$\sim$100$\%$ (with X= O or S). Our results imply that the large deuteration of H2CO and H2CS observed in protostellar cores as well as in comets is likely inherited from the pre-stellar phase. However, the comparison with state-of-the-art chemical models suggests that the reaction network for the formation of the doubly deuterated H2CS and H2CO it is not complete yet.
△ Less
Submitted 15 March, 2022;
originally announced March 2022.
-
Gas phase Elemental abundances in Molecular cloudS (GEMS). IV. Observational results and statistical trends
Authors:
M. Rodríguez-Baras,
A. Fuente,
P. Riviére-Marichalar,
D. Navarro-Almaida,
P. Caselli,
M. Gerin,
C. Kramer,
E. Roueff,
V. Wakelam,
G. Esplugues,
S. García-Burillo,
R. Le Gal,
S. Spezzano,
T. Alonso-Albi,
R. Bachiller,
S. Cazaux,
B. Commercon,
J. R. Goicoechea,
J. C. Loison,
S. P. Treviño-Morales,
O. Roncero,
I. Jiménez-Serra,
J. Laas,
A. Hacar,
J. Kirk
, et al. (11 additional authors not shown)
Abstract:
Gas phase Elemental abundances in Molecular CloudS (GEMS) is an IRAM 30m Large Program designed to estimate the S, C, N, and O depletions and gas ionization degree, X(e-), in a set of star-forming filaments of Taurus, Perseus and Orion. Our immediate goal is to build up a complete database of molecular abundances that can serve as an observational basis for estimating X(e-) and the C, O, N, and S…
▽ More
Gas phase Elemental abundances in Molecular CloudS (GEMS) is an IRAM 30m Large Program designed to estimate the S, C, N, and O depletions and gas ionization degree, X(e-), in a set of star-forming filaments of Taurus, Perseus and Orion. Our immediate goal is to build up a complete database of molecular abundances that can serve as an observational basis for estimating X(e-) and the C, O, N, and S depletions through chemical modeling. We observed and derived the abundances of 14 species (13CO, C18O, HCO+, H13CO+, HC18O+, HCN, H13CN, HNC, HCS+, CS, SO, 34SO, H2S, and OCS) in 244 positions, covering the AV 3 to 100 mag, n(H2) a few 10$^{3}$ to 10$^6$ cm$^{-3}$, and Tk 10 to 30 K ranges in these clouds, avoiding protostars, HII regions, and outflows. A statistical analysis is carried out to identify general trends between different species and with physical parameters. Relations between molecules reveal strong linear correlations which define three different families: (1) 13CO and C18O; (2) H13CO+, HC18O+, H13CN, and HNC; and (3) the S-bearing molecules. The abundances of the CO isotopologs increase with the gas kinetic temperature until TK 15 K. For higher temperatures, the abundance remains constant with a scatter of a factor of 3. The abundances of H13CO+, HC18O+, H13CN, and HNC are well correlated with each other, and all of them decrease with molecular hydrogen density, following the law n(H2)$^{-0.8\pm0.2}$. The abundances of S-bearing species also decrease with n(H2) at a rate of (S-bearing/H)gas n(H2)$^{-0.6\pm0.1}$. The abundances of molecules belonging to groups 2 and 3 do not present any clear trend with gas temperature. At scales of molecular clouds, the C18O abundance is the quantity that better correlates with the cloud mass. We discuss the utility of the 13CO/C18O, HCO+/H13CO+, and H13CO+/H13CN abundance ratios as chemical diagnostics of star formation in external galaxies.
△ Less
Submitted 25 February, 2021;
originally announced February 2021.
-
Gas-phase Elemental abundances in Molecular cloudS (GEMS) III. Unlocking the CS chemistry: the CS+O reaction
Authors:
N. Bulut,
O. Roncero,
A. Aguado,
J. -C. Loison,
D. Navarro-Almaida,
V. Wakelam,
A. Fuente,
E. Roueff,
R. Le Gal,
P. Caselli,
M. Gerinm K. M. Hickson,
S. Spezzano,
P. Riviere-Marichalar,
T. Alonso-Albi,
R. Bachiller,
I. Jimenez-Serra,
C. Kramer,
B. Tercero,
M. Rodriguez-Baras,
S. Garcia-Burillo,
J. R. Goicoechea,
S. P. Treviño-Morales,
G. Esplugues,
S. Cazaux,
B. Commercon
, et al. (14 additional authors not shown)
Abstract:
CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the el…
▽ More
CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the elemental abundance of sulfur. The CS+O -> CO + S reaction has been proposed as a relevant CS destruction mechanism at low temperatures, and could explain the discrepancy between models and observations. Its reaction rate has been experimentally measured at temperatures of 150-400 K, but the extrapolation to lower temperatures is doubtful. Here we calculate the CS+O reaction rate at temperatures <150 K which are prevailing in the interstellar medium. We performed ab initio calculations to obtain the three lowest PES of the CS+O system. These PESs are used to study the reaction dynamics, using several methods to eventually calculate the CS+O thermal reaction rates. We compare the results of our theoretical calculations for 150-400 K with those obtained in the laboratory. Our detailed theoretical study on the CS+O reaction, which is in agreement with the experimental data obtained at 150-400 K, demonstrates the reliability of our approach. After a careful analysis at lower temperatures, we find that the rate constant at 10 K is negligible, which is consistent with the extrapolation of experimental data using the Arrhenius expression. We use the updated chemical network to model the sulfur chemistry in TMC1 based on molecular abundances determined from GEMS project observations. In our model, we take into account the expected decrease of the cosmic ray ionization rate along the cloud. The abundance of CS is still overestimated when assuming the cosmic value for the sulfur abundance.
△ Less
Submitted 18 December, 2020;
originally announced December 2020.
-
Dust temperature and time-dependent effects in the chemistry of photodissociation regions
Authors:
Gisela Esplugues,
Stephanie Cazaux,
Paola Caselli,
Seyit Hocuk,
Marco Spaans
Abstract:
When studying the chemistry of PDRs, time dependence becomes important as visual extinction increases, since certain chemical timescales are comparable to the cloud lifetime. Dust temperature is also a key factor, since it significantly influences gas temperature and mobility on dust grains, determining the chemistry occurring on grain surfaces. We present a study of the dust temperature impact an…
▽ More
When studying the chemistry of PDRs, time dependence becomes important as visual extinction increases, since certain chemical timescales are comparable to the cloud lifetime. Dust temperature is also a key factor, since it significantly influences gas temperature and mobility on dust grains, determining the chemistry occurring on grain surfaces. We present a study of the dust temperature impact and time effects on the chemistry of different PDRs, using an updated version of the Meijerink PDR code and combining it with the time-dependent code Nahoon. We find the largest temperature effects in the inner regions of high $G$$_{\mathrm{0}}$ PDRs, where high dust temperatures favour the formation of simple oxygen-bearing molecules (especially that of O$_2$), while the formation of complex organic molecules is much more efficient at low dust temperatures. We also find that time-dependent effects strongly depend on the PDR type, since long timescales promote the destruction of oxygen-bearing molecules in the inner parts of low $G$$_{\mathrm{0}}$ PDRs, while favouring their formation and that of carbon-bearing molecules in high $G$$_{\mathrm{0}}$ PDRs. From the chemical evolution, we also conclude that, in dense PDRs, CO$_2$ is a late-forming ice compared to water ice, and confirm a layered ice structure on dust grains, with H$_2$O in lower layers than CO$_2$. Regarding steady state, the PDR edge reaches chemical equilibrium at early times ($\lesssim$10$^5$ yr). This time is even shorter ($<$10$^4$ yr) for high $G$$_{\mathrm{0}}$ PDRs. By contrast, inner regions reach equilibrium much later, especially low $G$$_{\mathrm{0}}$ PDRs, where steady state is reached at $\sim$10$^6$-10$^7$ yr.
△ Less
Submitted 6 April, 2019;
originally announced April 2019.
-
The first frost in the Pipe Nebula
Authors:
Miwa Goto,
J. D. Bailey,
Seyit Hocuk,
Paola Caselli,
Gisela B. Esplugues,
Stephanie Cazaux,
Marco Spaans
Abstract:
Spectroscopic studies of ices in nearby star-forming regions indicate that ice mantles form on dust grains in two distinct steps, starting with polar ice formation (H2O rich) and switching to apolar ice (CO rich). We test how well the picture applies to more diffuse and quiescent clouds where the formation of the first layers of ice mantles can be witnessed. Medium-resolution near-infrared spectra…
▽ More
Spectroscopic studies of ices in nearby star-forming regions indicate that ice mantles form on dust grains in two distinct steps, starting with polar ice formation (H2O rich) and switching to apolar ice (CO rich). We test how well the picture applies to more diffuse and quiescent clouds where the formation of the first layers of ice mantles can be witnessed. Medium-resolution near-infrared spectra are obtained toward background field stars behind the Pipe Nebula. The water ice absorption is positively detected at 3.0 micron in seven lines of sight out of 21 sources for which observed spectra are successfully reduced. The peak optical depth of the water ice is significantly lower than those in Taurus with the same visual extinction. The source with the highest water-ice optical depth shows CO ice absorption at 4.7 micron as well. The fractional abundance of CO ice with respect to water ice is 16+7-6 %, and about half as much as the values typically seen in low-mass star-forming regions. A small fractional abundance of CO ice is consistent with some of the existing simulations. Observations of CO2 ice in the early diffuse phase of a cloud play a decisive role in understanding the switching mechanism between polar and apolar ice formation.
△ Less
Submitted 16 October, 2017;
originally announced October 2017.
-
Parameterizing the interstellar dust temperature
Authors:
S. Hocuk,
L. Szucs,
P. Caselli,
S. Cazaux,
M. Spaans,
G. B. Esplugues
Abstract:
The temperature of interstellar dust particles is of great importance to astronomers. It plays a crucial role in the thermodynamics of interstellar clouds, because of the gas-dust collisional coupling. It is also a key parameter in astrochemical studies that governs the rate at which molecules form on dust. In 3D (magneto)hydrodynamic simulations often a simple expression for the dust temperature…
▽ More
The temperature of interstellar dust particles is of great importance to astronomers. It plays a crucial role in the thermodynamics of interstellar clouds, because of the gas-dust collisional coupling. It is also a key parameter in astrochemical studies that governs the rate at which molecules form on dust. In 3D (magneto)hydrodynamic simulations often a simple expression for the dust temperature is adopted, because of computational constraints, while astrochemical modelers tend to keep the dust temperature constant over a large range of parameter space. Our aim is to provide an easy-to-use parametric expression for the dust temperature as a function of visual extinction ($A_{\rm V}$) and to shed light on the critical dependencies of the dust temperature on the grain composition. We obtain an expression for the dust temperature by semi-analytically solving the dust thermal balance for different types of grains and compare to a collection of recent observational measurements. We also explore the effect of ices on the dust temperature. Our results show that a mixed carbonaceous-silicate type dust with a high carbon volume fraction matches the observations best. We find that ice formation allows the dust to be warmer by up to 15% at high optical depths ($A_{\rm V}> 20$ mag) in the interstellar medium. Our parametric expression for the dust temperature is presented as $T_{\rm d} = \left[ 11 + 5.7\times \tanh\bigl( 0.61 - \log_{10}(A_{\rm V})\bigr) \right] \, χ_{\rm uv}^{1/5.9}$, where $χ_{\rm uv}$ is in units of the Draine (1978) UV field
△ Less
Submitted 16 April, 2017; v1 submitted 10 April, 2017;
originally announced April 2017.
-
Surface chemistry in photodissociation regions
Authors:
G. B. Esplugues,
S. Cazaux,
R. Meijerink,
M. Spaans,
P. Caselli
Abstract:
The presence of dust can strongly affect the chemical composition of the interstellar medium. We model the chemistry in photodissociation regions (PDRs) using both gas-phase and dust-phase chemical reactions. Our aim is to determine the chemical compositions of the interstellar medium (gas/dust/ice) in regions with distinct (molecular) gas densities that are exposed to radiation fields with differ…
▽ More
The presence of dust can strongly affect the chemical composition of the interstellar medium. We model the chemistry in photodissociation regions (PDRs) using both gas-phase and dust-phase chemical reactions. Our aim is to determine the chemical compositions of the interstellar medium (gas/dust/ice) in regions with distinct (molecular) gas densities that are exposed to radiation fields with different intensities. We have significantly improved the Meijerink PDR code by including 3050 new gas-phase chemical reactions and also by implementing surface chemistry. In particular, we have included 117 chemical reactions occurring on grain surfaces covering different processes, such as adsorption, thermal desorption, chemical desorption, two-body reactions, photo processes, and cosmic-ray processes on dust grains. We obtain abundances for different gas and solid species as a function of visual extinction, depending on the density and radiation field. We also analyse the rates of the formation of CO2 and H2O ices in different environments. In addition, we study how chemistry is affected by the presence/absence of ice mantles (bare dust or icy dust) and the impact of considering different desorption probabilities. The type of substrate (bare dust or icy dust) and the probability of desorption can significantly alter the chemistry occurring on grain surfaces, leading to differences of several orders of magnitude in the abundances of gas-phase species, such as CO, H2CO, and CH3OH. The type of substrate, together with the density and intensity of the radiation field, also determine the threshold extinction to form ices of CO2 and H2O. We also conclude that H2CO and CH3OH are mainly released into the gas phase of low, far-ultraviolet illuminated PDRs through chemical desorption upon two-body surface reactions, rather than through photodesorption.
△ Less
Submitted 22 April, 2016;
originally announced April 2016.
-
Modelling the sulphur chemistry evolution in Orion KL
Authors:
G. B. Esplugues,
S. Viti,
J. R. Goicoechea,
J. Cernicharo
Abstract:
We study the sulphur chemistry evolution in the Orion KL along the gas and grain phases of the cloud. We investigate the processes that dominate the sulphur chemistry and to determine how physical and chemical parameters, such as the final star mass and the initial elemental abundances, influence the evolution of the hot core and of the surrounding outflows and shocked gas (the plateau). We indepe…
▽ More
We study the sulphur chemistry evolution in the Orion KL along the gas and grain phases of the cloud. We investigate the processes that dominate the sulphur chemistry and to determine how physical and chemical parameters, such as the final star mass and the initial elemental abundances, influence the evolution of the hot core and of the surrounding outflows and shocked gas (the plateau). We independently modelled the chemistry evolution of both components using the time-dependent gas-grain model UCL_CHEM and considering two different phase calculations. Phase I starts with the collapsing cloud and the depletion of atoms and molecules onto grain surfaces. Phase II starts when a central protostar is formed and the evaporation from grains takes place. We show how the gas density, the gas depletion efficiency, the initial sulphur abundance, the shocked gas temperature and the different chemical paths on the grains leading to different reservoirs of sulphur on the mantles affect sulphur-bearing molecules at different evolutionary stages. We also compare the predicted column densities with those inferred from observations of the species SO, SO2, CS, OCS, H2S and H2CS. The models that reproduce the observations of the largest number of sulphur-bearing species are those with an initial sulphur abundance of 0.1 times the sulphur solar abundance and a density of at least n_H=5x10^6 cm^-3 in the shocked gas region. We conclude that most of the sulphur atoms were ionised during Phase I, consistent with an inhomogeneous and clumpy region where the UV interstellar radiation penetrates leading to sulphur ionisation. We also conclude that the main sulphur reservoir on the ice mantles was H2S. In addition, we deduce that a chemical transition currently takes place in the shocked gas, where SO and SO2 gas-phase formation reactions change from being dominated by O2 to being dominated by OH.
△ Less
Submitted 9 June, 2014;
originally announced June 2014.
-
Herschel observations of EXtra-Ordinary Sources: Analysis of the HIFI 1.2 THz Wide Spectral Survey Toward Orion KL I. Methods
Authors:
Nathan R. Crockett,
Edwin A. Bergin,
Justin L. Neill,
Cécile Favre,
Peter Schilke,
Dariusz C. Lis,
Tom A. Bell,
Geoffrey Blake,
José Cernicharo,
Martin Emprechtinger,
Gisela B. Esplugues,
Harshal Gupta,
Maria Kleshcheva,
Steven Lord,
Nuria Marcelino,
Brett A. McGuire,
John Pearson,
Thomas G. Phillips,
Rene Plume,
Floris van der Tak,
Belén Tercero,
Shanshan Yu
Abstract:
We present a comprehensive analysis of a broad band spectral line survey of the Orion Kleinmann-Low nebula (Orion KL), one of the most chemically rich regions in the Galaxy, using the HIFI instrument on board the Herschel Space Observatory. This survey spans a frequency range from 480 to 1907 GHz at a resolution of 1.1 MHz. These observations thus encompass the largest spectral coverage ever obtai…
▽ More
We present a comprehensive analysis of a broad band spectral line survey of the Orion Kleinmann-Low nebula (Orion KL), one of the most chemically rich regions in the Galaxy, using the HIFI instrument on board the Herschel Space Observatory. This survey spans a frequency range from 480 to 1907 GHz at a resolution of 1.1 MHz. These observations thus encompass the largest spectral coverage ever obtained toward this high-mass star-forming region in the sub-mm with high spectral resolution, and include frequencies $>$ 1 THz where the Earth's atmosphere prevents observations from the ground. In all, we detect emission from 39 molecules (79 isotopologues). Combining this dataset with ground based mm spectroscopy obtained with the IRAM 30 m telescope, we model the molecular emission from the mm to the far-IR using the XCLASS program which assumes local thermodynamic equilibrium (LTE). Several molecules are also modeled with the MADEX non-LTE code. Because of the wide frequency coverage, our models are constrained by transitions over an unprecedented range in excitation energy. A reduced $χ^{2}$ analysis indicates that models for most species reproduce the observed emission well. In particular, most complex organics are well fit by LTE implying gas densities are high ($>$10$^6$ cm$^{-3}$) and excitation temperatures and column densities are well constrained. Molecular abundances are computed using H$_{2}$ column densities also derived from the HIFI survey. The distribution of rotation temperatures, $T_{\rm rot}$, for molecules detected toward the hot core is significantly wider than the compact ridge, plateau, and extended ridge $T_{\rm rot}$ distributions, indicating the hot core has the most complex thermal structure.
△ Less
Submitted 15 May, 2014; v1 submitted 9 May, 2014;
originally announced May 2014.
-
Extended warm gas in Orion KL as probed by methyl cyanide
Authors:
T. A. Bell,
J. Cernicharo,
S. Viti,
N. Marcelino,
Aina Palau,
G. B. Esplugues,
B. Tercero
Abstract:
In order to study the temperature distribution of the extended gas within the Orion Kleinmann-Low nebula, we have mapped the emission by methyl cyanide (CH3CN) in its J=6_K-5_K, J=12_K-11_K, J=13_K-12_K, and J=14_K-13_K transitions at an average angular resolution of ~10 arcsec (22 arcsec for the 6_K-5_K lines), as part of a new 2D line survey of this region using the IRAM 30m telescope. These ful…
▽ More
In order to study the temperature distribution of the extended gas within the Orion Kleinmann-Low nebula, we have mapped the emission by methyl cyanide (CH3CN) in its J=6_K-5_K, J=12_K-11_K, J=13_K-12_K, and J=14_K-13_K transitions at an average angular resolution of ~10 arcsec (22 arcsec for the 6_K-5_K lines), as part of a new 2D line survey of this region using the IRAM 30m telescope. These fully sampled maps show extended emission from warm gas to the northeast of IRc2 and the distinct kinematic signatures of the hot core and compact ridge source components. We have constructed population diagrams for the four sets of K-ladder emission lines at each position in the maps and have derived rotational excitation temperatures and total beam-averaged column densities from the fitted slopes. In addition, we have fitted LVG model spectra to the observations to determine best-fit physical parameters at each map position, yielding the distribution of kinetic temperatures across the region. The resulting temperature maps reveal a region of hot (T > 350 K) material surrounding the northeastern edge of the hot core, whereas the column density distribution is more uniform and peaks near the position of IRc2. We attribute this region of hot gas to shock heating caused by the impact of outflowing material from active star formation in the region, as indicated by the presence of broad CH3CN lines. This scenario is consistent with predictions from C-shock chemical models that suggest that gas-phase methyl cyanide survives in the post-shock gas and can be somewhat enhanced due to sputtering of grain mantles in the passing shock front.
△ Less
Submitted 12 February, 2014;
originally announced February 2014.
-
A combined IRAM and Herschel/HIFI study of cyano(di)acetylene in Orion KL: tentative detection of DC3N
Authors:
G. B. Esplugues,
J. Cernicharo,
S. Viti,
J. R. Goicoechea,
B. Tercero,
N. Marcelino,
Aina Palau,
T. A. Bell,
E. A. Bergin,
N. R. Crockett,
S. Wang
Abstract:
We present a study of cyanoacetylene (HC3N) and cyanodiacetylene (HC5N) in Orion KL, through observations from two line surveys performed with the IRAM 30m telescope and the HIFI instrument on board the Herschel telescope. The frequency ranges covered are 80-280 GHz and 480-1906 GHz. We model the observed lines of HC3N, HC5N, their isotopologues (including DC3N), and vibrational modes, using a non…
▽ More
We present a study of cyanoacetylene (HC3N) and cyanodiacetylene (HC5N) in Orion KL, through observations from two line surveys performed with the IRAM 30m telescope and the HIFI instrument on board the Herschel telescope. The frequency ranges covered are 80-280 GHz and 480-1906 GHz. We model the observed lines of HC3N, HC5N, their isotopologues (including DC3N), and vibrational modes, using a non-LTE radiative transfer code. To investigate the chemical origin of HC3N and DC3N in Orion KL, we use a time-dependent chemical model. We detect 40 lines of the ground state of HC3N and 68 lines of its 13C isotopologues. We also detect 297 lines of six vibrational modes of this molecule (nu_7, 2nu_7, 3nu_7, nu_6, nu_5, and nu_6+nu_7) and 35 rotational lines of the ground state of HC5N. We report the first tentative detection of DC3N in a giant molecular cloud with a DC3N/HC3N abundance ratio of 0.015. We provide column densities and isotopic and molecular abundances. We also perform a 2x2" map around Orion IRc2 and we present maps of HC3N lines and maps of lines of the HC3N vibrational modes nu_6 and nu_7. In addition, a comparison of our results for HC3N with those in other clouds allows us to derive correlations between the column density, the FWHM, the mass, and the luminosity of the clouds. The high column densities of HC3N obtained in the hot core, make this molecule an excellent tracer of hot and dense gas. In addition, the large frequency range covered reveals the need to consider a temperature and density gradient in the hot core in order to obtain better line fits. The high D/H ratio (comparable to that obtained in cold clouds) that we derive suggests a deuterium enrichment. Our chemical models indicate that the possible deuterated HC3N present in Orion KL is formed during the gas-phase. This fact provides new hints concerning the processes leading to deuteration.
△ Less
Submitted 2 September, 2013;
originally announced September 2013.
-
A line confusion-limited millimeter survey of Orion KL. III. Sulfur oxide species
Authors:
G. B. Esplugues,
B. Tercero,
J. Cernicharo,
J. R. Goicoechea,
Aina Palau,
N. Marcelino,
T. A. Bell
Abstract:
We present a study of the sulfur-bearing species detected in a line confusion-limited survey towards Orion KL performed with the IRAM 30m telescope in the range 80-281 GHz. The study is part of an analysis of the line survey divided into families of molecules. Our aim is to derive accurate physical conditions and molecular abundances in the different components of Orion KL from observed SO and SO2…
▽ More
We present a study of the sulfur-bearing species detected in a line confusion-limited survey towards Orion KL performed with the IRAM 30m telescope in the range 80-281 GHz. The study is part of an analysis of the line survey divided into families of molecules. Our aim is to derive accurate physical conditions and molecular abundances in the different components of Orion KL from observed SO and SO2 lines. First we assumed LTE conditions obtain rotational temperatures. We then used a radiative transfer model, assuming either LVG or LTE excitation to derive column densities of these molecules in the different components of Orion KL. We have detected 68 lines of SO, 34SO, 33SO, and S18O and 653 lines of SO2, 34SO2, 33SO2, SO18O and SO2 v2=1. We provide column densities for all of them and also upper limits for the column densities of S17O, 36SO, 34S18O, SO17O and 34SO2 v2=1 and for several undetected sulfur-bearing species. In addition, we present 2'x2' maps around Orion IRc2 of SO2 transitions with energies from 19 to 131 K and also maps with four transitions of SO, 34SO and 34SO2. We observe an elongation of the gas along the NE-SW direction. An unexpected emission peak appears at 20.5 km/s in most lines of SO and SO2. A study of the spatial distribution of this emission feature shows that it is a new component ~5" in diameter, which lies ~4" west of IRc2. We suggest the emission from this feature is related to shocks associated to the BN object. The highest column densities for SO and SO2 are found in the high-velocity plateau (a region dominated by shocks) and in the hot core. These values are up to three orders of magnitude higher than the results for the ridge components. We also find high column densities for their isotopologues in both components. Therefore, we conclude that SO and SO2 are good tracers, not only of regions affected by shocks, but also of regions with warm dense gas.
△ Less
Submitted 23 June, 2013;
originally announced June 2013.
