Source author record

Cecilia Ceccarelli

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

47works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

47 published item(s)

preprint2025arXiv

FAUST XXX: Dust enhancement in the young binary L1551 IRS 5

Young binary stars with discs provide unique laboratories to study the earliest stages of planet formation in star-forming environments. The detection of substructure in discs around Class I protostars challenges current models of disc evolution, suggesting that planets may form earlier than previously expected ($<1$ Myr). In the context of the ALMA Large Program FAUST, we present observations of the circumbinary disc (CBD) around the young binary system L1551 IRS 5. The CBD exhibits two prominent over-densities in the continuum emission at the edge of the cavity, with the Northern over-density being about 20% brighter than the Southern one. By analysing the disc morphology and kinematics of L1551 IRS 5, we delineate dynamical constraints on the binary's orbital parameters. Additionally, we present 3D hydrodynamical models of the CBD to predict both the dust and the gas surface densities. Then, we compare the resulting synthetic observations with ALMA observations of the continuum emission at 1.3 mm and the C$^{18}$O line emission. Our analysis suggests that the density enhancements observed with ALMA in L1551 IRS 5 can be caused by interactions between the binary stars and the CBD, leading to dust concentration within the disc. We conclude that the observed over-density corresponds to a location where could potentially grow under favourable conditions.

preprint2025arXiv

Gas-phase formation routes of dimethyl sulfide in the interstellar medium

Context: Dimethyl sulfide (DMS; CH$_3$SCH$_3$) is an organosulfur compound that has been suggested as a potential biosignature in exoplanetary atmospheres. In addition to its tentative detections toward the sub-Neptune planet K2-18b, DMS has been detected in the coma of the 67/P comet and toward the galactic center molecular cloud G+0.693-0.027. However, its formation routes have not been characterized yet. Aims: In this work, we have investigated three gas-phase reactions (CH$_3$SH + CH$_3$OH$_2^+$, CH$_3$OH + CH$_3$SH$_2^+$, and the CH$_3$ + CH$_3$S radiative association), aiming at characterizing DMS formation routes in shocked molecular clouds and star-forming regions. Methods: We have performed dedicated quantum and kinetics calculations to evaluate the reaction rate coefficients as a function of temperature to be included in astrochemical models. Results: Among the investigated processes, the reaction between methanethiol (CH$_3$SH) and protonated methanol (CH$_3$OH$_2^+$)(possibly followed by a gentle proton transfer to ammonia) is a compelling candidate to explain the formation of DMS in the galactic center molecular cloud G+0.693-0.027. The CH$_3$ + CH$_3$S radiative association does not seem to be a very efficient process, with the exclusion of cold clouds, provided that the thiomethoxy radical (CH$_3$S) is available. This work does not deal directly with the possible formation of DMS in the atmosphere of exoplanets. However, it clearly indicates that there are efficient abiotic formation routes of this interesting species.

preprint2022arXiv

A train of shocks at 3000 au scale? Exploring the clash of an expanding bubble into the NGC 1333 IRAS 4 region. SOLIS XIV

There is evidence that the star formation process is linked to the intricate net of filaments in molecular clouds, which may be also due to gas compression from external triggers. We studied the southern region of the Perseus NGC 1333 molecular cloud, known to be heavily shaped by similar external triggers, to shed light on the process that perturbed the filament where the Class 0 IRAS4 protostars lie. We use new IRAM-NOEMA observations of SiO and CH3OH, both known to trace violent events as shocks, toward IRAS 4A as part of the Large Program Seeds Of Life in Space (SOLIS). We detected three parallel elongated ($>$6000 au) structures, called fingers, with narrow line profiles (~1.5 $km s^{-1}$) peaked at the cloud systemic velocity, tracing gas with high density (5-20 $10^5 cm^{-3}$) and high temperature (80-160 K). They are chemically different, with the northern finger traced by both SiO and CH3OH ([CH3OH]/[SiO]~160-300), while the other two only by SiO ([CH3OH]/[SiO]$<$ 40). Among various possibilities, a train of three shocks, distanced by $>$5000 yr, would be consistent with the observations if a substantial fraction of silicon, frozen onto the grain mantles, is released by the shocks.We suggest that the shock train is due to an expanding gas bubble, coming behind NGC 1333 from the southwest and clashing against the filament, where IRAS 4A lies. Finally, we propose a solution to the two-decades long debate on the nature and origin of the widespread narrow SiO emission observed in the south part of NGC 1333, namely that it is due to unresolved trains of shocks.

preprint2022arXiv

CH$_3$CN deuteration in the SVS13-A Class I hot-corino. SOLIS XV

We studied the line emission from CH3CN and its deuterated isotopologue CH$_2$DCN towards the prototypical Class I object SVS13-A, where the deuteration of a large number of species has already been reported. Our goal is to measure the CH$_3$CN deuteration in a Class I protostar, for the first time, in order to constrain the CH$_3$CN formation pathways and the chemical evolution from the early prestellar core and Class 0 to the evolved Class I stages. We imaged CH2DCN towards SVS13-A using the IRAM NOEMA interferometer at 3mm in the context of the Large Program SOLIS (with a spatial resolution of 1.8"x1.2"). The NOEMA images have been complemented by the CH$_3$CN and CH$_2$DCN spectra collected by the IRAM-30m Large Program ASAI, that provided an unbiased spectral survey at 3mm, 2mm, and 1.3mm. The observed line emission has been analysed using LTE and non-LTE LVG approaches. The NOEMA/SOLIS images of CH2DCN show that this species emits in an unresolved area centered towards the SVS13-A continuum emission peak, suggesting that methyl cyanide and its isotopologues are associated with the hot corino of SVS13-A, previously imaged via other iCOMs. In addition, we detected 41 and 11 ASAI transitions of CH$_3$CN and CH2DCN, respectively, which cover upper level energies (Eup) from 13 to 442 K and from 18 K to 200 K, respectively. The derived [CH2DCN]/[CH3CN] ratio is $\sim$9\%. This value is consistent with those measured towards prestellar cores and a factor 2-3 higher than those measured in Class 0 protostars. Contrarily to what expected for other molecular species, the CH3CN deuteration does not show a decrease in SVS13-A with respect to measurements in younger prestellar cores and Class 0 protostars. Finally, we discuss why our new results suggest that CH3CN was likely synthesised via gas-phase reactions and frozen onto the dust grain mantles during the cold prestellar phase.

preprint2022arXiv

Chemical exploration of Galactic cold cores

A solar-type system starts from an initial molecular core that acquires organic complexity as it evolves. The so-called prestellar cores that can be studied are rare, which has hampered our understanding of how organic chemistry sets in and grows. Aims. We selected the best prestellar core targets from the cold core catalogue that represent a diversity in terms of their environment to explore their chemical complexity: 1390 (in the compressed shell of Lambda Ori), 869 (in the MBM12 cloud), and 4149 (in the California nebula). We obtained a spectral survey with the IRAM 30 m telescope in order to explore the molecular complexity of the cores. We carried out a radiative transfer analysis of the detected transitions in order to place some constraints on the physical conditions of the cores and on the molecular column densities. We also used the molecular ions in the survey to estimate the cosmic-ray ionisation rate and the S/H initial elemental abundance using a gas-phase chemical model to reproduce their abundances. We found large differences in the molecular complexity (deuteration, complex organic molecules, sulphur, carbon chains, and ions) and compared their chemical properties with a cold core and two prestellar cores. The chemical diversity we found in the three cores seems to be correlated with their chemical evolution: two of them are prestellar (1390 and 4149), and one is in an earlier stage (869). The influence of the environment is likely limited because cold cores are strongly shielded from their surroundings. The high extinction prevents interstellar UV radiation from penetrating deeply into the cores. Higher spatial resolution observations of the cores are therefore needed to constrain the physical structure of the cores, as well as a larger-scale distribution of molecular ions to understand the influence of the environment on their molecular complexity.

preprint2022arXiv

FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417

We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed.

preprint2022arXiv

Non-Energetic Formation of Ethanol via CCH Reaction with Interstellar H2O Ices. A Computational Chemistry Study

Ethanol (CH$_3$CH$_2$OH) is a relatively common molecule, often found in star forming regions. Recent studies suggest that it could be a parent molecule of several so-called interstellar complex organic molecules (iCOMs). Yet, the formation route of this species remains debated. In the present work, we study the formation of ethanol through the reaction of CCH with one H$_2$O molecule belonging to the ice, as a test case to investigate the viability of chemical reactions based on a "radical + ice component" scheme as an alternative mechanism for the synthesis of iCOMs, beyond the usual radical-radical coupling. This has been done by means of DFT calculations adopting two clusters of 18 and 33 water molecules as ice models. Results indicate that CH$_3$CH$_2$OH can potentially be formed by this proposed reaction mechanism. The reaction of CCH with H$_2$O on the water ice clusters can be barrierless (thanks to the help of boundary icy water molecules acting as proton transfer assistants) leading to the formation of vinyl alcohol precursors (H$_2$CCOH and CHCHOH). Subsequent hydrogenation of vinyl alcohol yielding ethanol is the only step presenting a low activation energy barrier. We finally discuss the astrophysical implications of these findings.

preprint2022arXiv

Seeds of Life in Space (SOLIS). XIII. Nitrogen fractionation towards the protocluster OMC-2 FIR4

Isotopic fractionation is an important tool to investigate the chemical history of our Solar System (SS). In particular, the isotopic fraction of nitrogen (14N/15N) is lower in comets and other pristine SS bodies with respect to the value measured for the protosolar nebula, suggesting a local chemical enrichment of 15N during the SS formation. Therefore, interferometric studies of nitrogen fractionation in SS precursors are imperative for us to obtain clues about our astrochemical origins. In this work, we investigated the variation of the 14N/15N ratio in one of the closest analogues of the environment in which the SS was born: the protocluster OMC-2 FIR4. We present the first comparison at high angular resolution between HCN and N2H+ using interferometric data. We analysed observations of the HCN isotopologues H13CN and HC15N in the OMC-2 FIR4 protocluster, specifically the transitions H13CN (1-0) and HC15N (1-0), from NOEMA within the context of the IRAM Seeds Of Life In Space Large Program. We combined our results with analysis of archival data obtained with ALMA of N2H+ and its 15N isotopologues. Our results show a small regional variation in the ratio for HCN from ~250 to 500. The ratios in the central regions of FIR4, where the candidate protostars are located, are largely consistent (~300). They also show little variation from the part of the protocluster known to harbour a high cosmic-ray ionisation rate, to the portion with lower rate. We also found a small variation in the ratio of N2H+ across different regions from ~200 to ~400. These results suggest that local changes in the physical parameters occurring on small linear scales probed by our observations do not seem to affect the 14N/15N ratio in either HCN or N2H+ and hence that this is independent of the molecule used. Moreover, the high level of irradiation due to cosmic rays does not affect the N-fractionation either.

preprint2022arXiv

The two hot corinos of the SVS13-A protostellar binary system: counterposed siblings

We present ALMA high-angular resolution ($\sim$ 50 au) observations of the Class I binary system SVS13-A. We report images of SVS13-A in numerous interstellar complex organic molecules: CH$_{\rm 3}$OH, $^{13}$CH$_{\rm 3}$OH, CH$_{\rm 3}$CHO, CH$_{\rm 3}$OCH$_{\rm 3}$, and NH$_{\rm 2}$CHO. Two hot corinos at different velocities are imaged in VLA4A (V$_{sys}$= +7.7 km s$^{-1}$) and VLA4B (V$_{sys}$= +8.5 km s$^{-1}$). From a non-LTE analysis of methanol lines we derive a gas density of 3 $\times$ 10$^8$ cm$^{-3}$, and gas temperatures of 140 K and 170 K for VLA4A and VLA4B, respectively. For the other species the column densities are derived from a LTE analysis. Formamide, which is the only N-bearing species detected in our observations, is more prominent around VLA4A, while dimethyl ether, methanol and acetaldehyde are associated with both VLA4A and VLA4B. We derive in the two hot corinos abundance ratios of $\sim$ 1 for CH$_{\rm 3}$OH, $^{13}$CH$_{\rm 3}$OH, and CH$_{\rm 3}$OCH$_{\rm 3}$, $\sim$ 2 for CH$_{\rm 3}$CHO, and $\sim$ 4 for NH$_{\rm 2}$CHO. The present dataset supports a chemical segregation between the different species inside the binary system. The emerging picture is that of an onion-like structure of the two SVS13-A hot corinos, caused by the different binding energies of the species, also supported by ad hoc quantum chemistry calculations. In addition, the comparison between molecular and dust maps suggests that the interstellar complex organic molecules emission originates from slow shocks produced by accretion streamers impacting the VLA4A and VLA4B disks and enriching the gas-phase component.

preprint2022arXiv

Theoretical distribution of the ammonia binding energy at interstellar icy grains: a new computational framework

The binding energies (BE) of molecules on the interstellar grains are crucial in the chemical evolution of the interstellar medium (ISM). Both temperature programmed desorption (TPD) laboratory experiments and quantum chemistry computations have often provided, so far, only single values of the BE for each molecule. This is a severe limitation, as the ices enveloping the grain mantles are structurally amorphous, giving rise to a manifold of possible adsorption sites, each with different BEs. However, the ice amorphous nature prevents the knowledge of structural details, hindering the development of a common accepted atomistic icy model. In this work, we propose a computational framework that closely mimics the formation of the interstellar grain mantle through a water by water accretion. On that grain, an unbiased random (but well reproducible) positioning of the studied molecule is then carried out. Here we present the test case of NH$_3$, an ubiquitous species in the molecular ISM. We provide the BE distribution computed by a hierarchy approach, using the semiempirical xTB-GFN2 as low-level method to describe the whole icy cluster combined with the B97D3 DFT functional as high-level method on the local zone of the NH$_3$ interaction. The final ZPE corrected BE is computed at ONIOM(DLPNO-CCSD(T)//B97D3:xTB-GFN2) level, ensuring the best cost/accuracy ratio. The main peak of the predicted NH$_3$ BE distribution is in agreement with experimental TPD and literature computed data. A second broad peak at very low BE values is also present, never detected before. It may provide the solution to a long-standing puzzle about the presence of gaseous NH$_3$ observed also in cold ISM objects.

preprint2022arXiv

Tracking the ice mantle history in the Solar-type Protostars of NGC 1333 IRAS 4

To understand the origin of the diversity observed in exoplanetary systems, it is crucial to characterize the early stages of their formation, represented by Solar-type protostars. Likely, the gaseous chemical content of these objects directly depends on the composition of the dust grain mantles formed before the collapse. Directly retrieving the ice mantle composition is challenging, but it can be done indirectly by observing the major components, such as NH3 and CH3OH at cm wavelengths, once they are released into the gas-phase during the warm protostellar stage. We observed several CH3OH and NH3 lines toward three Class 0 protostars in NGC1333 (IRAS 4A1, IRAS 4A2, and IRAS 4B), at high angular resolution (1"; ~300 au) with the VLA interferometer at 24-26 GHz. Using a non-LTE LVG analysis, we derived a similar NH3/CH3OH abundance ratio in the three protostars (<0.5, 0.015-0.5, and 0.003-0.3 for IRAS 4A1, 4A2, and 4B, respectively). Hence, we infer they were born from pre-collapse material with similar physical conditions. Comparing the observed abundance ratios with astrochemical model predictions, we constrained the dust temperature at the time of the mantle formation to be ~17 K, which coincides with the average temperature of the southern NGC 1333 diffuse cloud. We suggest that a brutal event started the collapse that eventually formed IRAS 4A1, 4A2 and 4B, which,therefore, did not experience the usual pre-stellar core phase. This event could be the clash of a bubble with NGC 1333 south, that has previously been evoked in the literature.

preprint2020arXiv

A 3mm chemical exploration of small organics in Class I YSOs

There is mounting evidence that the composition and structure of planetary systems are intimately linked to their birth environments. During the past decade, several spectral surveys probed the chemistry of the earliest stages of star formation and of late planet-forming disks. However, very little is known about the chemistry of intermediate protostellar stages, i.e. Class I Young Stellar Objects (YSOs), where planet formation may have already begun. We present here the first results of a 3mm spectral survey performed with the IRAM-30m telescope to investigate the chemistry of a sample of seven Class I YSOs located in the Taurus star-forming region. These sources were selected to embrace the wide diversity identified for low-mass protostellar envelope and disk systems. We present detections and upper limits of thirteen small ($N_{\rm atoms}\leq3$) C, N, O, and S carriers - namely CO, HCO$^+$, HCN, HNC, CN, N$_2$H$^+$, C$_2$H, CS, SO, HCS$^+$, C$_2$S, SO$_2$, OCS - and some of their D, $^{13}$C, $^{15}$N, $^{18}$O, $^{17}$O, and $^{34}$S isotopologues. Together, these species provide constraints on gas-phase C/N/O ratios, D- and $^{15}$N-fractionation, source temperature and UV exposure, as well as the overall S-chemistry. We find substantial evidence of chemical differentiation among our source sample, some of which can be traced back to Class I physical parameters, such as the disk-to-envelope mass ratio (proxy for Class I evolutionary stage), the source luminosity, and the UV-field strength. Overall, these first results allow us to start investigating the astrochemistry of Class I objects, however, interferometric observations are needed to differentiate envelope versus disk chemistry.

preprint2020arXiv

ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT) V: Sample, overview, and demography of disk molecular emission

We present an overview of the ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT), a campaign devoted to the characterization of the molecular emission from partly embedded, young stars. The project aims at better understanding the gaseous products that are delivered to planets by means of high-resolution maps of assorted lines probing disks at the time of the planet formation (less than 1 Myr). Nine different molecules are surveyed by our observations of six Class I/flat-spectrum sources. A series of accompanying articles analyze specific targets and molecules. Here we describe the sample and provide a general overview of the results, focusing on the spatial distribution, column densities, and abundance ratios of H$_2$CO, CS, and CN. The results of this work are a first step toward the characterization of the disk chemical evolution that need to be complemented by further observations of less exceptional disks and customized thermo-chemical modeling.

preprint2020arXiv

Chemical desorption versus energy dissipation: insights from ab-initio molecular dynamics of HCO formation

Molecular clouds are the cold regions of the Milky Way where stars form. They are enriched by rather complex molecules. Many of these molecules are believed to be synthesized on the icy surfaces of the interstellar submicron-sized dust grains that permeate the Galaxy. At 10 K thermal desorption is ineffcient and, therefore, why these molecules are found in the cold gas has tantalized astronomers for years. The assumption of the current models, called chemical desorption, is that the molecule formation energy released by the chemical reaction at the grain surface is partially absorbed by the grain and the remaining one causes the ejection of the newly formed molecule into the gas. Here we report an accurate ab-initio molecular dynamics simulations aimed to study the fate of the energy released by the first reaction of the H addition chain on CO, CO + H $\rightarrow$ HCO, occurring on a crystalline ice surface model. We show that about 90% of the HCO formation energy is injected towards the ice in the first picosecond, leaving HCO with an energy content (10-15 kJ mol$^{-1}$) more than a factor two lower than its adsorption energy (30 kJ mol$^{-1}$). As a result, in agreement with laboratory experiments, we conclude that chemical desorption is ineffcient for this specific system, namely H + CO on crystalline ice. We suspect this behavior to be quite general when dealing with hydrogen bonds, which are responsible of both the cohesive energy of the ice mantle and the interaction with adsorbates, as the HCO radical, even though ad hoc simulations are needed to draw specific conclusions on other systems.

preprint2020arXiv

H$_2$ mass-velocity relationship from 3D numerical simulations of jet-driven molecular outflows

Previous numerical studies have shown that in protostellar outflows, the mass-velocity distribution $m(v)$ can be well described by a broken power law $\propto v^{- γ}$. On the other hand, recent observations of a sample of outflows show that the CO intensity-velocity distribution, closely related to $m(v)$, follows an exponential law $\propto \exp(-v/v_0)$. In the present work, we revisit the physical origin of the mass-velocity relationship $m(v)$ in jet-driven protostellar outflows. We investigate the respective contributions of the different regions of the outflow, from the swept-up ambient gas to the jet. We performed 3D numerical simulations of a protostellar jet propagating into a molecular cloud using the hydrodynamical code Yguazu-a. The code takes into account atomic and ionic species and was modified to include the H$_2$ gas. We find that by excluding the jet contribution, $m(v)$ is satisfyingly fitted with a single exponential law, with $v_0$ well in the range of observational values. The jet contribution results in additional components in the mass-velocity relationship. This empirical mass-velocity relationship is found to be valid locally in the outflow. The exponent $v_0$ is almost constant in time and for a given level of mixing between the ambient medium and the jet material. In general, $v_0$ displays only a weak spatial dependence. A simple modeling of the L1157 outflow successfully reproduces the various components of the observed CO intensity-velocity relationship. Our simulations indicate that these components trace the outflow cavity of swept-up gas and the material entrained along the jet, respectively. The CO intensity-velocity exponential law is naturally explained by the jet-driven outflow model. The entrained material plays an important role in shaping the mass-velocity profile.

preprint2020arXiv

Hot Corinos Chemical Diversity: Myth or Reality?

After almost 20 years of hunting, only about a dozen hot corinos, hot regions enriched in interstellar complex organic molecules (iCOMs), are known. Of them, many are binary systems with the two components showing drastically different molecular spectra. Two obvious questions arise. Why are hot corinos so difficult to find and why do their binary components seem chemically different? The answer to both questions could be a high dust opacity that would hide the molecular lines. To test this hypothesis, we observed methanol lines at centimeter wavelengths, where dust opacity is negligible, using the Very Large Array interferometer. We targeted the NGC 1333 IRAS 4A binary system, for which one of the two components, 4A1, has a spectrum deprived of iCOMs lines when observed at millimeter wavelengths, while the other component, 4A2, is very rich in iCOMs. We found that centimeter methanol lines are similarly bright toward 4A1 and 4A2. Their non-LTE analysis indicates gas density and temperature ($\geq2\times10^6$ cm$^{-3}$ and 100--190 K), methanol column density ($\sim10^{19}$ cm$^{-2}$) and extent ($\sim$35 au in radius) similar in 4A1 and 4A2, proving that both are hot corinos. Furthermore, the comparison with previous methanol line millimeter observations allows us to estimate the optical depth of the dust in front of 4A1 and 4A2, respectively. The obtained values explain the absence of iCOMs line emission toward 4A1 at millimeter wavelengths and indicate that the abundances toward 4A2 are underestimated by $\sim$30\%. Therefore, centimeter observations are crucial for the correct study of hot corinos, their census, and their molecular abundances.

preprint2020arXiv

Seeds of Life in Space (SOLIS). VI. Chemical evolution of sulfuretted species along the outflows driven by the low-mass protostellar binary NGC1333-IRAS4A

Context. Low-mass protostars drive powerful molecular outflows that can be observed with mm and sub-mm telescopes. Various sulfuretted species are known to be bright in shocks and could be used to infer the physical and chemical conditions throughout the observed outflows. Aims. The evolution of sulfur chemistry is studied along the outflows driven by the NGC1333-IRAS4A protobinary system located in the Perseus cloud to constrain the physical and chemical processes at work in shocks. Methods. We observed various transitions from OCS, CS, SO, and SO$_2$ towards NGC1333-IRAS4A in the 1.3, 2, and 3mm bands using the IRAM NOEMA array and we interpreted the observations through the use of the Paris-Durham shock model. Results. The targeted species clearly show different spatial emission along the two outflows driven by IRAS4A. OCS is brighter on small and large scales along the south outflow driven by IRAS4A1, whereas SO$_2$ is detected rather along the outflow driven by IRAS4A2 that is extended along the north east - south west (NE-SW) direction. Column density ratio maps estimated from a rotational diagram analysis allowed us to confirm a clear gradient of the OCS/SO$_2$ column density ratio between the IRAS4A1 and IRAS4A2 outflows. SO is detected at extremely high radial velocity up to 25 km/s relative to the source velocity, clearly allowing us to distinguish the two outflows on small scales. Conclusions. The observed chemical differentiation between the two outflows of the IRAS4A system could be explained by a different chemical history. The outflow driven by IRAS4A1 is likely younger and more enriched in species initially formed in interstellar ices, such as OCS, and recently sputtered into the shock gas. In contrast, the longer and likely older outflow triggered by IRAS4A2 is more enriched in species that have a gas phase origin, such as SO$_2$.

preprint2016arXiv

A Search for O_2 in CO-depleted Molecular Cloud Cores with Herschel

The general lack of molecular oxygen in molecular clouds is an outstanding problem in astrochemistry. Extensive searches with SWAS, Odin and Herschel have only produced two detections; upper limits to the O_2 abundance in the remaining sources observed are about 1000 times lower than predicted by chemical models. Previous atomic oxygen observations and inferences from observations of other molecules indicated that high abundances of O atoms might be present in dense cores exhibiting large amounts of CO depletion. Theoretical arguments concerning the oxygen gas-grain interaction in cold dense cores suggested that, if O atoms could survive in the gas after most of the rest of the heavy molecular material has frozen out on to dust, then O_2 could be formed efficiently in the gas. Using Herschel HIFI we searched a small sample of four depletion cores - L1544, L694-2, L429, Oph D - for emission in the low excitation O_2 N_J=3_3-1_2 line at 487.249 GHz. Molecular oxygen was not detected and we derive upper limits to its abundance in the range N(O_2)/N(H_2) = (0.6 - 1.6)x10^{-7}. We discuss the absence of O_2 in the light of recent laboratory and observational studies.

preprint2016arXiv

Abundance of HOCO+ and CO2 in the outer layers of the L1544 prestellar core

The L1544 prestellar core has been observed as part of the ASAI IRAM Large Program at 3 mm. These observations led to the detection of many complex molecules. In this Letter, we report the detection of two lines, at 85.5 GHz (4,0,4-3,0,3) and 106.9 GHz (5,0,5-4,0,4), respectively, of the protonated carbon dioxide ion, HOCO+. We also report the tentative detection of the line at 100.4 GHz (5,0,5-4,0,4) of DOCO+. The non-LTE analysis of the detected lines shows that the HOCO+ emission originates in the external layer where non-thermal desorption of other species has previously been observed. Its abundance is (5 +/- 2) e-11. Modelling of the chemistry involved in the formation and destruction of HOCO+ provides a gaseous CO2 abundance of 2e-7 (with respect to H2) with an upper limit of 2e-6.

preprint2016arXiv

Discovery of a Hot Corino in the Bok Globule B335

We report the first evidence of a hot corino in a Bok globule. This is based on the ALMA observations in the 1.2 mm band toward the low-mass Class 0 protostar IRAS 19347+0727 in B335. Saturated complex organic molecules (COMs), CH$_3$CHO, HCOOCH$_3$, and NH$_2$CHO, are detected in a compact region within a few 10 au around the protostar. Additionally, CH$_3$OCH$_3$, C$_2$H$_5$OH, C$_2$H$_5$CN, and CH$_3$COCH$_3$ are tentatively detected. Carbon-chain related molecules, CCH and c-C$_3$H$_2$, are also found in this source, whose distributions are extended over a few 100 au scale. On the other hand, sulfur-bearing molecules CS, SO, and SO$_2$, have both compact and extended components. Fractional abundances of the COMs relative to H$_2$ are found to be comparable to those in known hot-corino sources. Though the COMs lines are as broad as 5-8 km s$^{-1}$, they do not show obvious rotation motion in the present observation. Thus, the COMs mainly exist in a structure whose distribution is much smaller than the synthesized beam (0."58 x 0."52).

preprint2016arXiv

Geometric and Kinematic Structure of the Outflow/Envelope System of L1527 Revealed by Subarcsecond-resolution Observation of CS

Subarcsecond-resolution images of the rotational line emissions of CS and c-C$_3$H$_2$ obtained toward the low-mass protostar IRAS 04368$+$2557 in L1527 with the Atacama Large Millimeter/submillimeter Array are investigated to constrain the orientation of the outflow/envelope system. The distribution of CS consists of an envelope component extending from north to south and a faint butterfly-shaped outflow component. The kinematic structure of the envelope is well reproduced by a simple ballistic model of an infalling rotating envelope. Although the envelope has a nearly edge-on configuration, the inclination angle of the rotation axis from the plane of the sky is found to be 5$^\circ$, where we find that the western side of the envelope faces the observer. This configuration is opposite to the direction of the large-scale ($\sim$ 10$^4$ AU) outflow suggested previously from the $^{12}$CO ($J$=3$-$2) observation, and to the morphology of infrared reflection near the protostar ($\sim$ 200 AU). The latter discrepancy could originate from high extinction by the outflow cavity of the western side, these discrepancies or may indicate that the outflow axis is not parallel to the rotation axis of the envelope. Position-velocity diagrams show the accelerated outflow cavity wall, and its kinematic structure in the 2000 AU scale is explained by a standard parabolic model with the inclination angle derived from the analysis of the envelope. The different orientation of the outflow between the small and large scale implies a possibility of precession of the outflow axis. The shape and the velocity of the outflow in the vicinity of the protostar are compared with those of other protostars.

preprint2016arXiv

H2S in the L1157-B1 Bow shock

Sulfur-bearing molecules are highly reactive in the gas phase of the ISM. However, the form in which most of the sulfur is locked onto interstellar dust grains is unknown. By taking advantage of the short time-scales of shocks in young molecular outflows, one could track back the main form of sulfur in the ices. In this paper, six transitions of H$_2$S and its isotopologues in the L1157-B1 bowshock have been detected using data from the Herschel-CHESS survey and the IRAM-30m ASAI large program. These detections are used to calculate the properties of H$_2$S gas in L1157-B1 through use of a rotation diagram and to explore the possible carriers of sulfur on the grains. The isotopologue detections allow the first calculation of the H$_2$S deuteration fraction in an outflow from a low mass protostar. The fractional abundance of H$_2$S in the region is found to be 6.0$\times$10$^{-7}$ and the deuteration fraction is 2$\times$10$^{-2}$. In order to investigate the form of sulfur on the grains, a chemical model is run with four different networks, each with different branching ratios for the freeze out of sulfur bearing species into molecules such as OCS and H$_2$S. It is found that the model best fits the data when at least half of each sulfur bearing species hydrogenates when freezing. We therefore conclude that a significant fraction of sulfur in L1157-B1 is likely to be locked in H$_2$S on the grains.

preprint2016arXiv

Infalling-Rotating Motion and Associated Chemical Change in the Envelope of IRAS 16293-2422 Source A Studied with ALMA

We have analyzed rotational spectral line emission of OCS, CH3OH, HCOOCH3, and H2CS observed toward the low-mass Class 0 protostellar source IRAS 16293-2422 Source A at a sub-arcsecond resolution (~0".6 x 0".5) with ALMA. Significant chemical differentiation is found at a 50 AU scale. The OCS line is found to well trace the infalling-rotating envelope in this source. On the other hand, the CH3OH and HCOOCH3 distributions are found to be concentrated around the inner part of the infalling-rotating envelope. With a simple ballistic model of the infalling-rotating envelope, the radius of the centrifugal barrier (a half of the centrifugal radius) and the protostellar mass are evaluated from the OCS data to be from 40 to 60 AU and from 0.5 to 1.0 Msun, respectively, assuming the inclination angle of the envelope/disk structure to be 60 degrees (90 degrees for the edge-on configuration). Although the protostellar mass is correlated with the inclination angle, the radius of the centrifugal barrier is not. This is the first indication of the centrifugal barrier of the infalling-rotating envelope in a hot corino source. CH3OH and HCOOCH3 may be liberated from ice mantles due to weak accretion shocks around the centrifugal barrier, and/or due to protostellar heating. The H2CS emission seems to come from the disk component inside the centrifugal barrier in addition to the envelope component. The centrifugal barrier plays a central role not only in the formation of a rotationally-supported disk but also in the chemical evolution from the envelope to the protoplanetary disk.

preprint2016arXiv

Subarcsecond Analysis of Infalling-Rotating Envelope around the Class I Protostar IRAS 04365+2535

Sub-arcsecond images of the rotational line emission of CS and SO have been obtained toward the Class I protostar IRAS 04365$+$2535 in TMC-1A with ALMA. A compact component around the protostar is clearly detected in the CS and SO emission. The velocity structure of the compact component of CS reveals infalling-rotating motion conserving the angular momentum. It is well explained by a ballistic model of an infalling-rotating envelope with the radius of the centrifugal barrier (a half of the centrifugal radius) of 50 AU, although the distribution of the infalling gas is asymmetric around the protostar. The distribution of SO is mostly concentrated around the radius of the centrifugal barrier of the simple model. Thus a drastic change in chemical composition of the gas infalling onto the protostar is found to occur at a 50 AU scale probably due to accretion shocks, demonstrating that the infalling material is significantly processed before being delivered into the disk.

preprint2016arXiv

The (impossible?) formation of acetaldehyde on the grain surfaces: insights from quantum chemical calculations

Complex Organic Molecules (COMs) have been detected in the interstellar medium (ISM). However, it is not clear whether their synthesis occurs on the icy surfaces of interstellar grains or via a series of gas-phase reactions. As a test case of the COMs synthesis in the ISM, we present new quantum chemical calculations on the formation of acetaldehyde (CH3CHO) from the coupling of the HCO and CH3 radicals, both in gas phase and on water ice surfaces. The binding energies of HCO and CH3 on the amorphous water ice were also computed (2333 and 734 K, respectively). Results indicate that, in gas phase, the products could be either CH3CHO, CH4 + CO, or CH3OCH, depending on the relative orientation of the two radicals. However, on the amorphous water ice, only the CH4 + CO product is possible due to the geometrical constraints imposed by the water ice surface. Therefore, acetaldehyde cannot be synthesized by the CH3 + HCO coupling on the icy grains. We discuss the implications of these results and other cases, such as ethylene glycol and dimethyl ether, in which similar situations can occur, suggesting that formation of these molecules on the grain surfaces might be unlikely.

preprint2015arXiv

Abundance Anomaly of the $^{13}$C Isotopic Species of c-C$_3$H$_2$ in the Low-Mass Star Formation Region L1527

The rotational spectral lines of c-C$_3$H$_2$ and two kinds of the $^{13}$C isotopic species, c-$^{13}$CCCH$_2$ ($C_{2v}$ symmetry) and c-CC$^{13}$CH$_2$ ($C_s$ symmetry) have been observed in the 1-3 mm band toward the low-mass star-forming region L1527. We have detected 7, 3, and 6 lines of c-C$_3$H$_2$, c-$^{13}$CCCH$_2$ , and c-CC$^{13}$CH$_2$, respectively, with the Nobeyama 45 m telescope, and 34, 6, and 13 lines, respectively, with the IRAM 30 m telescope, where 7, 2, and 2 transitions, respectively, are observed with the both telescopes. With these data, we have evaluated the column densities of the normal and $^{13}$C isotopic species. The [c-C$_3$H$_2$]/[c-$^{13}$CCCH$_2$] ratio is determined to be $310\pm80$, while the [c-C$_3$H$_2$]/[c-CC$^{13}$CH$_2$] ratio is determined to be $61\pm11$. The [c-C$_3$H$_2$]/[c-$^{13}$CCCH$_2$] and [c-C$_3$H$_2$]/[c-CC$^{13}$CH$_2$] ratios expected from the elemental $^{12}$C/$^{13}$C ratio are 60-70 and 30-35, respectively, where the latter takes into account the statistical factor of 2 for the two equivalent carbon atoms in c-C$_3$H$_2$. Hence, this observation further confirms the dilution of the $^{13}$C species in carbon-chain molecules and their related molecules, which are thought to originate from the dilution of $^{13}$C$^+$ in the gas-phase C$^+$ due to the isotope exchange reaction: $\mathrm{^{13}C^++CO\rightarrow{}^{13}CO+C^+}$. Moreover, the abundances of the two $^{13}$C isotopic species are different from each other. The ratio of c-$\mathrm{^{13}CCCH_2}$ species relative to c-$\mathrm{CC^{13}CH_2}$ is determined to be $0.20\pm0.05$. If $^{13}$C were randomly substituted for the three carbon atoms, the [c-$\mathrm{^{13}CCCH_2}$]/[c-$\mathrm{CC^{13}CH_2}$] ratio would be 0.5. Hence, the observed ratio indicates that c-$\mathrm{CC^{13}CH_2}$ exists more favorably. Possible origins of the different abundances are discussed.

preprint2015arXiv

Constraining the abundances of complex organics in the inner regions of solar-type protostars

The high abundances of Complex Organic Molecules (COMs) with respect to methanol, the most abundant COM, detected towards low-mass protostars, tend to be underpredicted by astrochemical models. This discrepancy might come from the large beam of the single-dish telescopes, encompassing several components of the studied protostar, commonly used to detect COMs. To address this issue, we have carried out multi-line observations of methanol and several COMs towards the two low-mass protostars NGC1333-IRAS2A and -IRAS4A with the Plateau de Bure interferometer at an angular resolution of 2 arcsec, resulting in the first multi-line detection of the O-bearing species glycolaldehyde and ethanol and of the N-bearing species ethyl cyanide towards low-mass protostars other than IRAS 16293. The high number of detected transitions from COMs (more than 40 methanol transitions for instance) allowed us to accurately derive the source size of their emission and the COMs column densities. The COMs abundances with respect to methanol derived towards IRAS2A and IRAS4A are slightly, but not substantitally, lower than those derived from previous single-dish observations. The COMs abundance ratios do not vary significantly with the protostellar luminosity, over five orders of magnitude, implying that low-mass hot corinos are quite chemically rich as high-mass hot cores. Astrochemical models still underpredict the abundances of key COMs, such as methyl formate or di-methyl ether, suggesting that our understanding of their formation remains incomplete.

preprint2015arXiv

Detectability of deuterated water in prestellar cores

Water is an important molecule in the chemical and thermal balance of dense molecular gas, but knowing its history through-out the various stages of the star formation is a fundamental problem. Its molecular deuteration provides us with a crucial clue to its formation history. H$_2$O has recently been detected for the first time towards the prestellar core L1544 with the Herschel Space Observatory with a high spectral resolution (HIFI instrument). Prestellar cores provide the original reservoir of material from which future planetary systems are built, but few observational constraints exist on the formation of water and none on its deuteration before the collapse starts and a protostar forms at the centre. We report on new APEX observations of the ground state 1$_{0,1}$-0$_{0,0}$ HDO transition at 464 GHz towards the prestellar core L1544. The line is undetected, and we present an extensive study of the conditions for its detectability in cold and dense cloud cores. The water and deuterated water abundances have been estimated using an advanced chemical model simplified for the limited number of reactions or processes that are active in cold regions (< 15 K). We use the LIME radiative transfer code to compute the expected intensity and profile of both H$_2$O and HDO lines and compare them with the observations. We present several ad hoc profiles that best-fit the observations and compare the profiles with results from an astrochemical modelling, coupling gas phase and grain surface chemistry. Our comparison between observations, radiative transfer, and chemical modelling shows the limits of detectability for singly deuterated water, through the ground-state transitions 1$_{0,1}$-0$_{0,0}$ and 1$_{1,1}$-0$_{0,0}$ at 464.9 and 893.6 GHz, respectively, with both single-dish telescope and interferometric observations.

preprint2015arXiv

Formation of complex organic molecules in cold objects: the role of gas phase reactions

While astrochemical models are successful in reproducing many of the observed interstellar species, they have been struggling to explain the observed abundances of complex organic molecules. Current models tend to privilege grain surface over gas phase chemistry in their formation. One key assumption of those models is that radicals trapped in the grain mantles gain mobility and react on lukewarm (>30 K) dust grains. Thus, the recent detections of methyl formate (MF) and dimethyl ether (DME) in cold objects represent a challenge and may clarify the respective role of grain surface and gas phase chemistry. We propose here a new model to form DME and MF with gas phase reactions in cold environments, where DME is the precursor of MF via an efficient reaction overlooked by previous models. Furthermore, methoxy, a precursor of DME, is also synthetized in the gas phase from methanol, which is desorbed by a non-thermal process from the ices. Our new model reproduces fairy well the observations towards L1544. It also explains, in a natural way, the observed correlation between DME and MF. We conclude that gas phase reactions are major actors in the formation of MF, DME and methoxy in cold gas. This challenges the exclusive role of grain-surface chemistry and favours a combined grain-gas chemistry.

preprint2014arXiv

Cosmic ray induced ionisation of a molecular cloud shocked by the W28 supernova remnant

Cosmic rays are an essential ingredient in the evolution of the interstellar medium, as they dominate the ionisation of the dense molecular gas, where stars and planets form. However, since they are efficiently scattered by the galactic magnetic fields, many questions remain open, such as where exactly they are accelerated, what is their original energy spectrum, and how they propagate into molecular clouds. In this work we present new observations and discuss in detail a method that allows us to measure the cosmic ray ionisation rate towards the molecular clouds close to the W28 supernova remnant. To perform these measurements, we use CO, HCO$^+$, and DCO$^+$ millimetre line observations and compare them with the predictions of radiative transfer and chemical models away from thermodynamical equilibrium. The CO observations allow us to constrain the density, temperature, and column density towards each observed position, while the DCO$^+$/HCO$^+$ abundance ratios provide us with constraints on the electron fraction and, consequently, on the cosmic ray ionisation rate. Towards positions located close to the supernova remnant, we find cosmic ray ionisation rates much larger ($\gtrsim100$) than those in standard galactic clouds. Conversely, towards one position situated at a larger distance, we derive a standard cosmic ray ionisation rate. Overall, these observations support the hypothesis that the $γ$ rays observed in the region have a hadronic origin. In addition, based on CR diffusion estimates, we find that the ionisation of the gas is likely due to $0.1 -1$ GeV cosmic rays. Finally, these observations are also in agreement with the global picture of cosmic ray diffusion, in which the low-energy tail of the cosmic ray population diffuses at smaller distances than the high-energy counterpart.

preprint2014arXiv

Deuterium Fractionation: the Ariadne's Thread from the Pre-collapse Phase to Meteorites and Comets today

The Solar System formed about 4.6 billion years ago from a condensation of matter inside a molecular cloud. Trying to reconstruct what happened is the goal of this chapter. For that, we put together our understanding of Galactic objects that will eventually form new suns and planetary systems, with our knowledge on comets, meteorites and small bodies of the Solar System today. Our specific tool is the molecular deuteration, namely the amount of deuterium with respect to hydrogen in molecules. This is the Ariadne's thread that helps us to find the way out from a labyrinth of possible histories of our Solar System. The chapter reviews the observations and theories of the deuterium fractionation in pre-stellar cores, protostars, protoplanetary disks, comets, interplanetary dust particles and meteorites and links them together trying to build up a coherent picture of the history of the Solar System formation. We emphasise the interdisciplinary nature of the chapter, which gathers together researchers from different communities with the common goal of understanding the Solar System history.

preprint2014arXiv

The origin of complex organic molecules in prestellar cores

Complex organic molecules (COMs) have been detected in a variety of environments, including cold prestellar cores. Given the low temperature of these objects, these last detections challenge existing models. We report here new observations towards the prestellar core L1544. They are based on an unbiased spectral survey of the 3mm band at the IRAM-30m telescope, as part of the Large Program ASAI. The observations allow us to provide the full census of the oxygen bearing COMs in this source. We detected tricarbon monoxide, methanol, acetaldehyde, formic acid, ketene, and propyne with abundances varying from 5e-11 to 6e-9. The non-LTE analysis of the methanol lines shows that they are likely emitted at the border of the core, at a radius of ~8000 AU where T~10 K and nH2~2e4 cm-3. Previous works have shown that water vapour is enhanced in the same region because of the photodesorption of water ices. We propose that a non-thermal desorption mechanism is also responsible for the observed emission of methanol and COMs from the same layer. The desorbed oxygen and a tiny amount of desorbed methanol and ethene are enough to reproduce the abundances of tricarbon monoxide, methanol, acetaldehyde and ketene measured in L1544. These new findings open the possibility that COMs in prestellar cores originate in a similar outer layer rather than in the dense inner cores, as previously assumed, and that their formation is driven by the non-thermally desorbed species.

preprint2013arXiv

Water deuterium fractionation in the inner regions of two solar type protostars

The [HDO]/[H2O] ratio is a crucial parameter for probing the history of water formation. So far, it has been measured for only three solar type protostars and yielded different results, possibly pointing to a substantially different history in their formation. In the present work, we report new interferometric observations of the HDO 4 2,2 - 4 2,3 line for two solar type protostars, IRAS2A and IRAS4A, located in the NGC1333 region. In both sources, the detected HDO emission originates from a central compact unresolved region. Comparison with previously published interferometric observations of the H218$O 3 1,3 - 2 2,0 line shows that the HDO and H$_2$O lines mostly come from the same region. A non-LTE LVG analysis of the HDO and H218$O line emissions, combined with published observations, provides a [HDO]/[H2O] ratio of 0.3 - 8 % in IRAS2A and 0.5 - 3 % in IRAS4A. First, the water fractionation is lower than that of other molecules such as formaldehyde and methanol in the same sources. Second, it is similar to that measured in the solar type protostar prototype, IRAS16293-2422, and, surprisingly enough, larger than that measured in NGC1333 IRAS4B. {The comparison of the measured values towards IRAS2A and IRAS4A with the predictions of our gas-grain model GRAINOBLE gives similar conclusions to those for IRAS 16293, arguing that these protostars {share} a similar chemical history, although they are located in different clouds.

preprint2012arXiv

A study of deuterated water in the low-mass protostar IRAS16293-2422

The HDO/H2O ratio is a powerful diagnostic to understand the evolution of water from the first stages of star formation to the formation of planets and comets. Our aim is to determine precisely the abundance distribution of HDO towards the low-mass protostar IRAS16293-2422 and learn more about the water formation mechanisms by determining the HDO/H2O abundance ratio. A spectral survey of the source IRAS16293-2422 was carried out in the framework of the CHESS Herschel Key program with the HIFI instrument, allowing the detection of numerous HDO lines. Other transitions have been observed previously with ground-based telescopes. The spherical Monte Carlo radiative transfer code RATRAN was used to reproduce the observed line profiles of HDO by assuming an abundance jump. To determine the H2O abundance throughout the envelope, a similar study was made of the H2-18O observed lines, as the H2O main isotope lines are contaminated by the outflows. We derive an inner HDO abundance of 1.7e-7 and an outer HDO abundance of 8e-11. To reproduce the HDO absorption lines, it is necessary to add an absorbing layer in front of the envelope. It may correspond to a water-rich layer created by the photodesorption of the ices at the edges of the molecular cloud. The HDO/H2O ratio is ~1.4-5.8% in the hot corino whereas it is ~0.2-2.2% in the outer envelope. It is estimated at ~4.8% in the added absorbing layer. Although it is clearly higher than the cosmic D/H abundance, the HDO/H2O ratio remains lower than the D/H ratio derived for other deuterated molecules observed in the same source. The similarity of the ratios derived in the hot corino and in the added absorbing layer suggests that water formed before the gravitational collapse of the protostar, contrary to formaldehyde and methanol, which formed later once the CO molecules had depleted on the grains.

preprint2012arXiv

Formaldehyde and methanol deuteration in protostars: fossiles from a past fast high density pre-collapse phase

Extremely high deuteration of several molecules have been observed around low mass protostars since a decade. Among them, formaldehyde and methanol present particularly high deuteration, with observations of abundant doubly and triply deuterated forms. Both species are thought to be mainly formed on interstellar grains during the low temperature and dense pre-collapse phase by H and D atom additions on the iced CO. We present here a theoretical study of the formaldehyde and methanol deuteration obtained with our gas-grain model, GRAINOBLE. This model takes into account the multilayer nature of the mantle and explores the robustness of the results against the uncertainties of poorly constrained chemical and surface model parameters. The comparison of the model predictions with the observations leads to two major results: i) the observed high deuteration is obtained during the last phase of the pre-collapse stage, when the density reaches 5 10^6 cm^-3, and this phase is fast, lasting only several thousands years. ii) D and H abstraction and substitution reactions are crucial in making up the observed deuteration ratios; This work shows the power of chemical composition as a tool to reconstruct the past history of protostars.

preprint2012arXiv

Herschel CHESS discovery of the fossil cloud that gave birth to the Trapezium and Orion KL

As part of the Herschel Guaranteed Time Key Programme CHESS, we present the discovery of a diffuse gas component in the foreground of the intermediate-mass protostar OMC-2 FIR 4, located in the Orion A region. Making use of the full HIFI spectrum of OMC-2 FIR 4 obtained in CHESS, we detected several ground-state lines from OH+, H2O+, HF, and CH+, all of them seen in absorption against the dust continuum emission of the protostar's envelope. The lines peak at a velocity of 9 km/s, which is blue-shifted by 2 km/s with respect to the systemic velocity of OMC-2 FIR 4 (V = 11.4 km/s). We derived column densities for each species, as well as an upper limit to the column density of the undetected H3O+. In order to model and characterise the foreground cloud, we used the Meudon PDR code to run a homogeneous grid of models that spans a reasonable range of densities, visual extinctions, cosmic ray ionisation rates and far-ultraviolet (FUV) radiation fields. The results of our modelling indicate that the foreground cloud is composed of predominantly neutral diffuse gas (n(H) = 100 cm-3) and is heavily irradiated by an external source of FUV that most likely arises from the nearby Trapezium OB association. The cloud is 6 pc thick and bears many similarities with the so-called C+ interface between Orion-KL and the Trapezium cluster, 2 pc south of OMC-2 FIR 4. We conclude that the foreground cloud we detected is an extension of the C+ interface seen in the direction of Orion KL, and interpret it to be the remains of the parental cloud of OMC-1, which extends from OMC-1 up to OMC-2.

preprint2012arXiv

Hydrogen Fluoride in High-Mass Star-forming Regions

Hydrogen fluoride has been established to be an excellent tracer of molecular hydrogen in diffuse clouds. In denser environments, however, the HF abundance has been shown to be approximately two orders of magnitude lower. We present Herschel/HIFI observations of HF J=1-0 toward two high-mass star formation sites, NGC6334 I and AFGL 2591. In NGC6334 I the HF line is seen in absorption in foreground clouds and the source itself, while in AFGL 2591 HF is partially in emission. We find an HF abundance with respect to H2 of 1.5e-8 in the diffuse foreground clouds, whereas in the denser parts of NGC6334 I, we derive a lower limit on the HF abundance of 5e-10. Lower HF abundances in dense clouds are most likely caused by freeze out of HF molecules onto dust grains in high-density gas. In AFGL 2591, the view of the hot core is obstructed by absorption in the massive outflow, in which HF is also very abundant 3.6e-8) due to the desorption by sputtering. These observations provide further evidence that the chemistry of interstellar fluorine is controlled by freeze out onto gas grains.

preprint2012arXiv

Our astrochemical heritage

Our Sun and planetary system were born about 4.5 billion years ago. How did this happen and what is our heritage from these early times? This review tries to address these questions from an astrochemical point of view. On the one hand, we have some crucial information from meteorites, comets and other small bodies of the Solar System. On the other hand, we have the results of studies on the formation process of Sun-like stars in our Galaxy. These results tell us that Sun-like stars form in dense regions of molecular clouds and that three major steps are involved before the planet formation period. They are represented by the pre-stellar core, protostellar envelope and protoplanetary disk phases. Simultaneously with the evolution from one phase to the other, the chemical composition gains increasing complexity. In this review, we first present the information on the chemical composition of meteorites, comets and other small bodies of the Solar System, which is potentially linked to the first phases of the Solar System's formation. Then we describe the observed chemical composition in the pre-stellar core, protostellar envelope and protoplanetary disk phases, including the processes that lead to them. Finally, we draw together pieces from the different objects and phases to understand whether and how much we inherited chemically from the time of the Sun's birth.

preprint2012arXiv

Rotating Disks and Non-Kinematic Double Peaks

Double-peaked line profiles are commonly considered a hallmark of rotating disks, with the distance between the peaks a measure of the rotation velocity. However, double-peaks can arise also from radiative transfer effects in optically thick non-rotating sources. Utilizing exact solutions of the line transfer problem we present a detailed study of line emission from geometrically thin Keplerian disks. We derive the conditions for emergence of kinematic double peaks in optically thin and thick disks, and find that it is generally impossible to disentangle the effects of kinematics and line opacity in observed double-peaked profiles. Unless supplemented by additional information, a double-peaked profile alone is not a reliable indicator of a rotating disk. In certain circumstances, triple and quadruple profiles might be better indicators of rotation in optically thick disks.

preprint2012arXiv

The CHESS survey of the L1157-B1 shock: the dissociative jet shock as revealed by Herschel--PACS

Outflows generated by protostars heavily affect the kinematics and chemistry of the hosting molecular cloud through strong shocks that enhance the abundance of some molecules. L1157 is the prototype of chemically active outflows, and a strong shock, called B1, is taking place in its blue lobe between the precessing jet and the hosting cloud. We present the Herschel-PACS 55--210 micron spectra of the L1157-B1 shock, showing emission lines from CO, H2O, OH, and [OI]. The spatial resolution of the PACS spectrometer allows us to map the warm gas traced by far-infrared (FIR) lines with unprecedented detail. The rotational diagram of the high-Jup CO lines indicates high-excitation conditions (Tex ~ 210 +/- 10 K). We used a radiative transfer code to model the hot CO gas emission observed with PACS and in the CO (13-12) and (10-9) lines measured by Herschel-HIFI. We derive 200<Tkin<800 K and n>10^5 cm-3. The CO emission comes from a region of about 7 arcsec located at the rear of the bow shock where the [OI] and OH emission also originate. Comparison with shock models shows that the bright [OI] and OH emissions trace a dissociative J-type shock, which is also supported by a previous detection of [FeII] at the same position. The inferred mass-flux is consistent with the "reverse" shock where the jet is impacting on the L1157-B1 bow shock. The same shock may contribute significantly to the high-Jup CO emission.

preprint2012arXiv

Water ice deuteration: a tracer of the chemical history of protostars

Context. Millimetric observations have measured high degrees of molecular deuteration in several species seen around low-mass protostars. The Herschel Space Telescope, launched in 2009, is now providing new measures of the deuterium fractionation of water, the main constituent of interstellar ices. Aims. We aim at theoretically studying the formation and the deuteration of water, which is believed to be formed on interstellar grain surfaces in molecular clouds. Methods. We used our gas-grain astrochemical model GRAINOBLE, which considers the multilayer formation of interstellar ices. We varied several input parameters to study their impact on water deuteration. We included the treatment of ortho- and para-states of key species, including H2, which affects the deuterium fractionation of all molecules. The model also includes relevant laboratory and theoretical works on the water formation and deuteration on grain surfaces. In particular, we computed the transmission probabilities of surface reactions using the Eckart model, and we considered ice photodissociation following molecular dynamics simulations. Results. The use of a multilayer approach allowed us to study the influence of various parameters on the abundance and the deuteration of water. Deuteration of water is found to be very sensitive to the ortho-to-para ratio of H2 and to the total density, but it also depends on the gas/grain temperatures and the visual extinction of the cloud. Since the deuteration is very sensitive to the physical conditions, the comparison with sub-millimetric observation towards the low-mass protostar IRAS 16293 allows us to suggest that water ice is formed together with CO2 in molecular clouds with limited density, whilst formaldehyde and methanol are mainly formed in a later phase, where the condensation becomes denser and colder.

preprint2011arXiv

Multilayer modeling of porous grain surface chemistry I. The GRAINOBLE model

Mantles of iced water, mixed with CO, H2CO, and CH3OH are formed during the so called prestellar core phase. In addition, radicals are also thought to be formed on the grain surfaces, and to react to form complex organic molecules later on, during the warm-up phase of the protostellar evolution. We aim to study the formation of the grain mantles during the prestellar core phase and the abundance of H2CO, CH3OH, and radicals trapped in them. We have developed a macrosopic statistic multilayer model that follows the formation of grain mantles with time and that includes two effects that may increase the number of radicals trapped in the mantles: i) at each time of the mantle formation, only the surface layer is chemically active rather than the entire bulk, and ii) the porous structure of grains allows to trap reactive particles. The model considers a network of H, O and CO forming neutral species such as water, CO, formaldehyde, and methanol, plus several radicals. We run a large grid of models to study the impact of the mantle multilayer nature and grain porous structure. In addition, we explored the influence of the uncertainty of other key parameters on the mantle composition. Our model predicts relatively large abundances of radicals. In addition, the multilayer approach makes it possible to follow the chemical differentiation within the grain mantle, showing that the mantles are far from being uniform. For example, methanol is mostly present in the outer layers of the mantles whereas CO and other reactive species are trapped in the inner layers. The overall mantle composition depends on the density and age of the prestellar core, and on some microscopic parameters. Comparison with observations allows us to constrain the value of few parameters and provide some indications on the physical conditions during the formation of the ices.

preprint2010arXiv

The solar type protostar IRAS16293-2422: new constraints on the physical structure

Context: The low mass protostar IRAS16293-2422 is a prototype Class 0 source with respect to the studies of the chemical structure during the initial phases of life of Solar type stars. Aims: In order to derive an accurate chemical structure, a precise determination of the source physical structure is required. The scope of the present work is the derivation of the structure of IRAS16293-2422. Methods: We have re-analyzed all available continuum data (single dish and interferometric, from millimeter to MIR) to derive accurate density and dust temperature profiles. Using ISO observations of water, we have also reconstructed the gas temperature profile. Results: Our analysis shows that the envelope surrounding IRAS16293-2422 is well described by the Shu "inside-out" collapsing envelope model or a single power-law density profile with index equal to 1.8. In contrast to some previous studies, our analysis does not show evidence of a large (>/- 800 AU in diameter) cavity. Conclusions: Although IRAS16293-2422 is a multiple system composed by two or three objects, our reconstruction will be useful to derive the chemical structure of the large cold envelope surrounding these objects and the warm component, treated here as a single source, from single-dish observations of molecular emission.

preprint2010arXiv

TIMASSS : The IRAS16293-2422 Millimeter And Submillimeter Spectral Survey: Tentative Detection of Deuterated Methyl Formate (DCOOCH3)

High deuterium fractionation is observed in various types of environment such as prestellar cores, hot cores and hot corinos. It has proven to be an efficient probe to study the physical and chemical conditions of these environments. The study of the deuteration of different molecules helps us to understand their formation. This is especially interesting for complex molecules such as methanol and bigger molecules for which it may allow to differentiate between gas-phase and solid-state formation pathways. Methanol exhibits a high deuterium fractionation in hot corinos. Since CH3OH is thought to be a precursor of methyl formate we expect that deuterated methyl formate is produced in such environments. We have searched for the singly-deuterated isotopologue of methyl formate, DCOOCH3, in IRAS 16293-2422, a hot corino well-known for its high degree of methanol deuteration. We have used the IRAM/JCMT unbiased spectral survey of IRAS 16293-2422 which allows us to search for the DCOOCH3 rotational transitions within the survey spectral range (80-280 GHz, 328-366 GHz). The expected emission of deuterated methyl formate is modelled at LTE and compared with the observations.} We have tentatively detected DCOOCH3 in the protostar IRAS 16293-2422. We assign eight lines detected in the IRAM survey to DCOOCH3. Three of these lines are affected by blending problems and one line is affected by calibration uncertainties, nevertheless the LTE emission model is compatible with the observations. A simple LTE modelling of the two cores in IRAS 16293-2422, based on a previous interferometric study of HCOOCH3, allows us to estimate the amount of DCOOCH3 in IRAS 16293-2422. Adopting an excitation temperature of 100 K and a source size of 2\arcsec and 1\farcs5 for the A and B cores, respectively, we find that N(A,DCOOCH3) = N(B,DCOOCH3) ~ 6.10^14 /cm2. The derived deuterium fractionation is ~ 15%, consistent with values for other deuterated species in this source and much greater than that expected from the deuterium cosmic abundance. DCOOCH3, if its tentative detection is confirmed, should now be considered in theoretical models that study complex molecule formation and their deuteration mechanisms. Experimental work is also needed to investigate the different chemical routes leading to the formation of deuterated methyl formate.

preprint2008arXiv

Chemical modeling of L183 (= L134N) : an estimate of the ortho/para H2 ratio

Context. The high degree of deuteration observed in some prestellar cores depends on the ortho-to-para H2 ratio through the H3+ fractionation. Aims. We want to constrain the ortho/para H2 ratio across the L183 prestellar core. This is mandatory to correctly describe the deuter- ation amplification phenomenon in depleted cores such as L183 and to relate the total (ortho+para) H2D+ abundance to the sole ortho-H2D+ column density measurement. Methods. To constrain this ortho/para H2 ratio and derive its profile, we make use of the N2D+ /N2H+ ratio and of the ortho-H2D+ observations performed across the prestellar core. We use two simple chemical models limited to an almost totally depleted core description. New dissociative recombination and trihydrogen cation-dihydrogen reaction rates (including all isotopologues) are presented in this paper and included in our models. Results. We estimate the H2D+ ortho/para ratio in the L183 cloud, and constrain the H2 ortho/para ratio : we show that it is varying across the prestellar core by at least an order of magnitude being still very high (~0.1) in most of the cloud. Our time-dependent model indicates that the prestellar core is presumably older than 1.5-2 x 10^5 years but that it may not be much older. We also show that it has reached its present density only recently and that its contraction from a uniform density cloud can be constrained. Conclusions. A proper understanding of deuteration chemistry cannot be attained without taking into account the whole ortho/para family of molecular hydrogen and trihydrogen cation isotopologues as their relations are of utmost importance in the global scheme. Tracing the ortho/para H2 ratio should also give useful constraints on the dynamical evolution of prestellar cores.

preprint2006arXiv

Hot corinos in NGC1333-IRAS4B and IRAS2A

(Abridged) We aim to enlarge the number of known hot corinos and carry out a first comparative study with hot cores. The ultimate goal is to understand whether complex organic molecules form in the gas phase or on grain surfaces, and what the possible key parameters are. We observed millimeter rotational transitions of HCOOH, HCOOCH3, CH3OCH3, CH3CN, and C2H5CN in a sample of low-mass protostars with the IRAM-30m. Using the rotational diagram method coupled with the information about the sources' structure, we calculate the abundances of the observed molecules. To interpret these abundances, we review the proposed formation processes of the above molecules. We report the detection of HCOOCH3 and/or CH3CN towards NGC1333-IRAS4B and NGC1333-IRAS2A. We find that abundance ratios of O-bearing molecules to methanol or formaldehyde in hot corinos are comparable and about unity, and are relatively (depending on how the ratios are determined) higher than those in hot cores and in Galactic center clouds. So far, complex organic molecules were detected in all the hot corinos where they were searched for, suggesting that it is a common phase for low-mass protostars. While some evidence points to grain-surface synthesis (either in the cold or warm-up phase) of these molecules (in particular for HCOOH and HCOOCH3), the present data do not allow us to disregard gas-phase formation. More observational, laboratory, and theoretical studies are required to improve our understanding of hot corinos.

preprint2002arXiv

The Degree of CO Depletion in Pre-stellar Cores

We present new results on CO depletion in a sample of nearby pre-stellar cores, based on observations of the millimeter C17O and C18O lines and the 1.3 mm dust emission with the IRAM 30m telescope. In most cases, the distribution of CO is much flatter than that of the dust, whereas other tracers, like N2H+, still probe the latter. In the centre of these objects, we estimate CO to be underabundant by a factor 4-15 depending on the cores. The CO underabundance is more pronounced in the central regions and appears to decrease with increasing distance from the core centre. This underabundance is most likely due to the freezing out of CO onto the dust grains in the cold, dense parts of the cores. We find evidence for an increase of the CO depletion degree with the core density.