Search

The gas-phase abundance of carbon, x_C = C/H, and its depletion factors are essential parameters for understanding the gas and solid compositions that are ultimately incorporated into planets. The majority of protoplanetary disks are born in clusters and, as a result, are exposed to external FUV radiation. These FUV photons potentially affect the disk's evolution, chemical composition, and line excitation. We present the first detection of the [CI]609um fine-structure line of neutral carbon (CI), achieved with ALMA, toward one of these disks, d203-506, in the Orion Nebula Cluster. We also report the detection of CI forbidden and permitted lines (from electronically excited states up to 10 eV) observed with JWST in the IR. These lines trace the irradiated outer disk and photo-evaporative wind. Contrary to the common belief that these IR lines are C+ recombination lines, we find that they are dominated by FUV-pumping of CI followed by fluorescence cascades. They trace the transition from atomic to molecular gas, and their intensities scale with G0. The lack of outstanding IR OI fluorescent emission, however, implies a sharper attenuation of external FUV radiation with E > 12 eV (~Lyman-beta). This is related to a lower effective FUV dust absorption cross section compared to that of interstellar grains, implying a more prominent role for FUV shielding by the CI photoionization continuum. The [CI]609um intensity is proportional to N(CI) and can be used to infer x_C. We derive x_C ~ 1.4E-4. This implies that there is no major depletion of volatile carbon compared to x_C measured in the natal cloud, hinting at a young disk. We also show that external FUV radiation impacts the outer disk and wind by vertically shifting the water freeze-out depth, which results in less efficient grain growth and settling. This shift leads to nearly solar gas-phase C/O abundance ratios in these irradiated layers., Comment: Accepted for publication in A&A Letters. 14 pages including Appendices

A rich inventory of complex organic molecules (COMs) has been observed in high abundances in the gas phase toward Class 0 protostars. These molecules are suggested to be formed in ices and sublimate in the warm inner envelope close to the protostar. However, only the most abundant COM, methanol (CH3OH), has been firmly detected in ices before the era of James Webb Space Telescope (JWST). Now it is possible to detect the interstellar ices of other COMs and constrain their ice column densities quantitatively. We aim to determine the column densities of several oxygen-bearing COMs (O-COMs) in both gas and ice for two low-mass protostellar sources, NGC 1333 IRAS 2A and B1-c, as case studies in our JWST Observations of Young protoStars (JOYS+) program. By comparing the column density ratios w.r.t. CH3OH between both phases measured in the same sources, we can probe into the evolution of COMs from ice to gas in the early stages of star formation. We are able to fit the fingerprints range of COM ices between 6.8 and 8.8 um in the JWST/MIRI-MRS spectra of B1-c using similar components as recently used for IRAS 2A. We claim detection of CH4, OCN-, HCOO-, HCOOH, CH3CHO, C2H5OH, CH3OCH3, CH3OCHO, and CH3COCH3 in B1-c, and upper limits are estimated for SO2, CH3COOH, and CH3CN. The comparison of O-COM ratios w.r.t CH3OH between ice and gas shows two different cases. 1) the column density ratios of CH3OCHO and CH3OCH3 match well between the two phases, which may be attributed to a direct inheritance from ice to gas or strong chemical links with CH3OH. 2) the ice ratios of CH3CHO and C2H5OH w.r.t. CH3OH are higher than the gas ratios by 1-2 orders of magnitudes. This difference can be explained by the gas-phase reprocessing following sublimation, or different spatial distributions of COMs in the envelope., Comment: 42 pages (22 main text, 20 appendix); 27 figures (12 in main text, 15 in appendix); 5 tables (2 in main text, 3 in appendix) Accepted for publication in A&A

Very low-mass stars (those <0.3 solar masses) host orbiting terrestrial planets more frequently than other types of stars, but the compositions of those planets are largely unknown. We use mid-infrared spectroscopy with the James Webb Space Telescope to investigate the chemical composition of the planet-forming disk around ISO-ChaI 147, a 0.11 solar-mass star. The inner disk has a carbon-rich chemistry: we identify emission from 13 carbon-bearing molecules including ethane and benzene. We derive large column densities of hydrocarbons indicating that we probe deep into the disk. The high carbon to oxygen ratio we infer indicates radial transport of material within the disk, which we predict would affect the bulk composition of any planets forming in the disk., Comment: Published, 36 pages, 8 figures

Column density ratios of complex organic molecules are generally constant across protostellar systems with some low-level scatter. However, the scatter in formamide (NH$_2$CHO) to methanol (CH$_3$OH) column density ratio is one of the highest. This larger scatter is sometimes interpreted as evidence of gas-phase formation of NH$_2$CHO. In this work we propose an alternative interpretation in which this scatter is produced by differences in the snowline locations related to differences in binding energies of these species and the small-scale structure of the envelope and the disk system. We also include CH$_3$CN in our work as a control molecule which has a similar binding energy to CH$_3$OH. We use radiative transfer models to calculate the emission from these species in protostellar systems with and without disks. The abundances of these species are parameterized in our models. Then we fit the calculated emission lines to find the column densities as done in real observations. We find a correction factor of ~10 to be multiplied by gas-phase $N_{NH_2CHO}/N_{CH_3OH}$ to give the true abundance ratio of these two species in the ices. We find that models with different physical parameters produce a scatter in $N_{NH_2CHO}/N_{CH_3OH}$, comparable with that of observations. The scatter in $N_{NH_2CHO}/N_{CH_3OH}$ is larger than that of $N_{CH_3CN}/N_{CH_3OH}$ in models consistent with the observations. We show that the scatter in $N_{NH_2CHO}/N_{CH_3OH}$ will be lower if we correct for the difference in sublimation temperatures of these two species in observations of ~40 protostellar systems with ALMA. The scatter in $N_{NH_2CHO}/N_{CH_3OH}$ can be partially explained by the difference in their binding energies. We conclude that gas-phase chemistry routes for NH$_2$CHO are not necessary to explain the larger scatter of $N_{NH_2CHO}/N_{CH_3OH}$ compared with other ratios., Comment: Accepted for publication in A&A

Context. Two protoplanets have recently been discovered within the PDS 70 protoplanetary disk. JWST/NIRCam offers a unique opportunity to characterize them and their birth environment at wavelengths difficult to access from the ground. Aims. We aim to image the circumstellar environment of PDS 70 at 1.87 $\mu$m and 4.83 $\mu$m, assess the presence of Pa-$\alpha$ emission due to accretion onto the protoplanets, and probe any IR excess indicative of heated circumplanetary material. Methods. We obtain non-coronagraphic JWST/NIRCam images of PDS 70 within the MINDS (MIRI mid-INfrared Disk Survey) program. We leverage the Vortex Image Processing (VIP) package for data reduction, and develop dedicated routines for optimal stellar PSF subtraction, unbiased imaging of the disk, and protoplanet flux measurement in this type of dataset. A radiative transfer model of the disk is used to disentangle the contributions from the disk and the protoplanets. Results. We re-detect both protoplanets and identify extended emission after subtracting a disk model, including a large-scale spiral-like feature. We interpret its signal in the direct vicinity of planet c as tracing the accretion stream feeding its circumplanetary disk, while the outer part of the feature may rather reflect asymmetric illumination of the outer disk. We also report a bright signal consistent with a previously proposed protoplanet candidate enshrouded in dust, near the 1:2:4 mean-motion resonance with planets b and c. The 1.87 $\mu$m flux of planet b is consistent with atmospheric model predictions, but not that of planet c. We discuss potential origins for this discrepancy, including significant Pa-$\alpha$ line emission. The 4.83 $\mu$m fluxes of planets b and c suggest enshrouding dust or heated CO emission from their circumplanetary environment., Comment: 6+11 pages, 3+10 figures (text+appendix). Accepted for publication in A&A

Protostellar disks are known to accrete, however, the exact mechanism that extracts the angular momentum and drives accretion in the low-ionization "dead" region of the disk is under debate. In recent years, magneto-hydrodynamic (MHD) disk winds have become a popular solution. Yet, observations of these winds require both high spatial resolution (${\sim}10$s au) and high sensitivity, which has resulted in only a handful of MHD disk wind candidates so far. In this work we present high angular resolution (${\sim}30$ au) ALMA observations of the emblematic L1448-mm protostellar system and find suggestive evidence for an MHD disk wind. The disk seen in dust continuum (${\sim}0.9$ mm) has a radius of ${\sim}23$ au. Rotating infall signatures in H$^{13}$CO$^+$ indicate a central mass of $0.4\pm 0.1$ M$_\odot$ and a centrifugal radius similar to the dust disk radius. Above the disk, we unveil rotation signatures in the outflow traced by H$^{13}$CN, CH$_3$OH, and SO lines and find a kinematical structure consistent with theoretical predictions for MHD disk winds. This is the first detection of an MHD disk wind candidate in H$^{13}$CN and CH$_3$OH. The wind launching region estimated from cold MHD wind theory extends out to the disk edge. The magnetic lever arm parameter would be $\lambda_{\phi} \simeq 1.7$, in line with recent non-ideal MHD disk models. The estimated mass-loss rate is ${\sim}4$ times the protostellar accretion rate ($\dot{M}_{\rm acc} \simeq 2 \times 10^{-6} M_{\odot}/yr$) and suggests that the rotating wind could carry enough angular momentum to drive disk accretion., Comment: Accepted for publication in A&A

Terrestrial and sub-Neptune planets are expected to form in the inner ($<10~$AU) regions of protoplanetary disks. Water plays a key role in their formation, although it is yet unclear whether water molecules are formed in-situ or transported from the outer disk. So far Spitzer Space Telescope observations have only provided water luminosity upper limits for dust-depleted inner disks, similar to PDS 70, the first system with direct confirmation of protoplanet presence. Here we report JWST observations of PDS 70, a benchmark target to search for water in a disk hosting a large ($\sim54~$AU) planet-carved gap separating an inner and outer disk. Our findings show water in the inner disk of PDS 70. This implies that potential terrestrial planets forming therein have access to a water reservoir. The column densities of water vapour suggest in-situ formation via a reaction sequence involving O, H$_2$, and/or OH, and survival through water self-shielding. This is also supported by the presence of CO$_2$ emission, another molecule sensitive to UV photodissociation. Dust shielding, and replenishment of both gas and small dust from the outer disk, may also play a role in sustaining the water reservoir. Our observations also reveal a strong variability of the mid-infrared spectral energy distribution, pointing to a change of inner disk geometry., Comment: To appear in Nature on 24 July 2023. 21 pages, 10 figures; includes extended data. Part of the JWST MINDS Guaranteed Time Observations program's science enabling products. Spectra downloadable on Zenodo at https://zenodo.org/record/7991022

Early results from the JWST-MIRI guaranteed time programs on protostars (JOYS) and disks (MINDS) are presented. Thanks to the increased sensitivity, spectral and spatial resolution of the MIRI spectrometer, the chemical inventory of the planet-forming zones in disks can be investigated with unprecedented detail across stellar mass range and age. Here data are presented for five disks, four around low-mass stars and one around a very young high-mass star. The mid-infrared spectra show some similarities but also significant diversity: some sources are rich in CO2, others in H2O or C2H2. In one disk around a very low-mass star, booming C2H2 emission provides evidence for a ``soot'' line at which carbon grains are eroded and sublimated, leading to a rich hydrocarbon chemistry in which even di-acetylene (C4H2) and benzene (C6H6) are detected (Tabone et al. 2023). Together, the data point to an active inner disk gas-phase chemistry that is closely linked to the physical structure (temperature, snowlines, presence of cavities and dust traps) of the entire disk and which may result in varying CO2/H2O abundances and high C/O ratios >1 in some cases. Ultimately, this diversity in disk chemistry will also be reflected in the diversity of the chemical composition of exoplanets., Comment: 17 pages, 8 figures. Author's version of paper submitted to Faraday Discussions January 18 2023, Accepted March 16 2023

Constraining the strength of gas turbulence in protoplanetary discs is an open problem that has relevant implications for the physics of gas accretion and planet formation. In this work, we gauge the amount of turbulence in 6 of the discs observed in the DSHARP programme by indirectly measuring the vertical distribution of their dust component. We employ the differences in the gap contrasts observed along the major and the minor axes due to projection effects, and build a radiative transfer model to reproduce these features for different values of the dust scale heights. We find that (a) the scale heights that yield a better agreement with data are generally low ($\lesssim 4$ AU at a radial distance of $100$ AU), and in almost all cases we are only able to place upper limits on their exact values; these conclusions imply (assuming an average Stokes number of $\approx10^{-2}$) low turbulence levels of $\alpha_{\rm SS}\lesssim10^{-3}-10^{-4}$; (b) for the 9 other systems we considered out of the DSHARP sample, our method yields no significant constraints on the disc vertical structure; we conclude that this is because these discs have either a low inclination or gaps that are not deep enough. Based on our analysis we provide an empirical criterion to assess whether a given disc is suitable to measure the vertical scale height., Comment: Accepted for publication in MNRAS. 13 pages + appendix, 12 figures

We show that the distribution of observed accretion rates is a powerful diagnostic of protoplanetary disc physics. Accretion due to turbulent ("viscous") transport of angular momentum results in a fundamentally different distribution of accretion rates than accretion driven by magnetised disc winds. We find that a homogeneous sample of $\gtrsim$300 observed accretion rates would be sufficient to distinguish between these two mechanisms of disc accretion at high confidence, even for pessimistic assumptions. Current samples of T Tauri star accretion rates are not this large, and also suffer from significant inhomogeneity, so both viscous and wind-driven models are broadly consistent with the existing observations. If accretion is viscous, the observed accretion rates require low rates of disc photoevaporation ($\lesssim$$10^{-9}$M$_{\odot}$yr$^{-1}$). Uniform, homogeneous surveys of stellar accretion rates can therefore provide a clear answer to the long-standing question of how protoplanetary discs accrete., Comment: 10 pages, 8 figures. Accepted for publication in MNRAS

Earth is deficient in carbon and nitrogen by up to ${\sim}4$ orders of magnitude compared with the Sun. Destruction of (carbon- and nitrogen-rich) refractory organics in the high-temperature planet forming regions could explain this deficiency. Assuming a refractory cometary composition for these grains, their destruction enhances nitrogen-containing oxygen-poor molecules in the hot gas ($\gtrsim 300$K) after the initial formation and sublimation of these molecules from oxygen-rich ices in the warm gas (${\sim}150$K). Using observations of $37$ high-mass protostars with ALMA, we find that oxygen-containing molecules (CH$_3$OH and HNCO) systematically show no enhancement in their hot component. In contrast, nitrogen-containing, oxygen-poor molecules (CH$_3$CN and C$_2$H$_3$CN) systematically show an enhancement of a factor ${\sim} 5$ in their hot component, pointing to additional production of these molecules in the hot gas. Assuming only thermal excitation conditions, we interpret these results as a signature of destruction of refractory organics, consistent with the cometary composition. This destruction implies a higher C/O and N/O in the hot gas than the warm gas, while, the exact values of these ratios depend on the fraction of grains that are effectively destroyed. This fraction can be found by future chemical models that constrain C/O and N/O from the abundances of minor carbon, nitrogen and oxygen carriers presented here., Comment: Accepted for publication in ApJ Letters

Carbon is an essential element for life but how much can be delivered to young planets is still an open question. The chemical characterization of planet-forming disks is a crucial step in our understanding of the diversity and habitability of exoplanets. Very low-mass stars ($<0.2~M_{\odot}$) are interesting targets because they host a rich population of terrestrial planets. Here we present the JWST detection of abundant hydrocarbons in the disk of a very low-mass star obtained as part of the MIRI mid-INfrared Disk Survey (MINDS). In addition to very strong and broad emission from C$_2$H$_2$ and its $^{13}$C$^{12}$CH$_2$ isotopologue, C$_4$H$_2$, benzene, and possibly CH$_4$ are identified, but water, PAH and silicate features are weak or absent. The lack of small silicate grains implies that we can look deep down into this disk. These detections testify to an active warm hydrocarbon chemistry with a high C/O ratio in the inner 0.1 au of this disk, perhaps due to destruction of carbonaceous grains. The exceptionally high C$_2$H$_2$/CO$_2$ and C$_2$H$_2$/H$_2$O column density ratios suggest that oxygen is locked up in icy pebbles and planetesimals outside the water iceline. This, in turn, will have significant consequences for the composition of forming exoplanets., Comment: version submitted to Nature Astronomy

[Abridged] Most well-resolved disks observed with ALMA show signs of dust traps. These dust traps set the chemical composition of the planet forming material in these disks, as the dust grains with their icy mantles are trapped at specific radii and could deplete the gas and dust of volatiles at smaller radii. In this work we analyse the first detection of nitric oxide (NO) in a protoplanetary disk. We aim to constrain the nitrogen chemistry and the gas-phase C/O ratio in the highly asymmetric dust trap in the Oph-IRS 48 disk. We use ALMA observations of NO, CN, C$_2$H, and related molecules and model the effect of the dust trap on the physical and chemical structure using the thermochemical code DALI. Furthermore, we explore how ice sublimation contributes to the observed emission lines. NO is only observed at the location of the dust trap but CN and C$_2$H are not detected in the Oph-IRS 48 disk. This results in an CN/NO column density ratio of $< 0.05$ and thus a low C/O ratio at the location of the dust trap. The main gas-phase formation pathways to NO through OH and NH in the fiducial model predict NO emission that is an order of magnitude lower than is observed. The gaseous NO column density can be increased by factors ranging from 2.8 to 10 when the H$_2$O and NH$_3$ gas abundances are significantly boosted by ice sublimation. However, these models are inconsistent with the upper limits on the H$_2$O and OH column densities derived from observations. We propose that the NO emission in the Oph-IRS 48 disk is closely related to the nitrogen containing ices sublimating in the dust trap. The non-detection of CN constrains the C/O ratio both inside and outside the dust trap to be $< 1$ if all nitrogen initially starts as N$_2$ and $\leq 0.6$, consistent with the Solar value, if (part of) the nitrogen initially starts as N or NH$_3$., Comment: Accepted for publication in Astronomy and Astrophysics

Complex organic molecules are only detected toward a fraction of high-mass protostars. The goal of this work is to investigate whether high-mass disks can explain the lack of methanol emission from some massive protostellar systems. We consider an envelope-only and an envelope-plus-disk model and use RADMC-3D to calculate the methanol emission. High and low millimeter (mm) opacity dust are considered for both models separately and the methanol abundance is parameterized. Viscous heating is included due to the high accretion rates of these objects in the disk. In contrast with low-mass protostars, the presence of a disk does not significantly affect the temperature structure and methanol emission. The shadowing effect of the disk is not as important for high-mass objects and the disk mid-plane is hot because of viscous heating, which is effective due to the high accretion rates. Consistent with observations of infrared absorption lines toward high-mass protostars, we find a vertical temperature inversion, i.e. higher temperatures in the disk mid-plane than the disk surface, at radii < 50au for the models with $L=10^4$ L$_{\odot}$ and large mm opacity dust as long as the envelope mass is >550 M$_{\odot}$. The large observed scatter in methanol emission from massive protostars can be mostly explained toward lower luminosity objects with the envelope-plus-disk models including low and high mm opacity dust. The methanol emission variation toward sources with high luminosities cannot be explained by models with or without a disk. However, the $L/M$ of these objects suggest that they could be associated with hypercompact/ultracompact HII regions. Therefore, the low methanol emission toward the high-luminosity sources can be explained by them hosting an HII region where methanol is absent., Comment: 25 pages, 24 figures, Accepted for publication in A&A

We aim to directly trace the vertical location of the emitting surface of multiple molecular tracers in protoplanetary disks. Our sample of disks includes Elias 2-27, WaOph 6 and the sources targeted by the MAPS ALMA Large Program. The set of molecules studied include CO isotopologues in various transitions, HCN, CN, H2CO, HCO+, C2H and c-C3H2. The vertical emitting region is determined directly from the channel maps, implementing accurate masking of the channel emission to recover the vertical location of the emission surface even at large radial distances from the star and for low-SNR lines. The vertical location of the emitting layer is obtained for 4-10 lines in each disk. IM Lup, HD163296 and MWC 480 12CO and 13CO show vertical modulations, which are coincident with dust gaps and kinematical perturbations. We also present estimates of the gas pressure scale height in the disks from the MAPS sample. Compared to physical-chemical models we find good agreement with the vertical location of CO isotopologues. In HD 163296 CN and HCN trace a similar intermediate layer, for the other disks, the UV flux tracers and the vertical profiles of HCN and C2H are lower than predicted in theoretical models. HCN and H2CO show a highly structured vertical profile, possibly indicative of different formation pathways. It is possible to trace the vertical locations of multiple molecular species that trace a wide variety of physical and chemical disk properties. The distribution of CO isotopologues are found at a wide range of vertical heights $z/r = $ 0.5-0.05. Other molecular lines are mostly found at $z/r \leq $0.15. The vertical layering of molecules is in agreement with theory in some systems, but not in all, therefore dedicated chemical-physical models are needed to further study and understand the emission surfaces., Comment: Accepted for publication in A&A. 29 pages, 28 figures

No statistical study of COMs toward a large sample of high-mass protostars with ALMA has been carried out so far. We aim to study six N-bearing species: CH$_3$CN, HNCO, NH$_2$CHO, C$_2$H$_5$CN, C$_2$H$_3$CN and CH$_3$NH$_2$ in a large sample of high-mass protostars. From the ALMAGAL survey, 37 of the most line-rich hot molecular cores are selected. Next, we fit their spectra and find column densities and excitation temperatures of the above N-bearing species, in addition to CH$_3$OH. We (tentatively) detect CH$_3$NH$_2$ in $\sim32%$ of the sources. We find three groups of species when comparing their excitation temperatures: hot (NH$_2$CHO; Tex > 250 K), warm (C$_2$H$_3$CN, HN$^{13}$CO and CH$_{3}^{13}$CN; 100 K < Tex < 250 K) and cold species (CH$_3$OH and CH$_3$NH$_2$; Tex < 100 K). This temperature segregation reflects the trend seen in their sublimation temperatures and validates the idea of onion-like structure of COMs around protostars. Moreover, the molecules studied here show constant column density ratios across low- and high-mass protostars with scatter less than a factor $\sim3$ around the mean. The constant column density ratios point to a common formation environment of COMs or their precursors, most likely in the pre-stellar ices. The scatter around the mean of the ratios, although small, varies depending on the species considered. This spread can either have a physical origin (source structure, line or dust optical depth) or a chemical one. Formamide is most prone to the physical effects as it is tracing the closest regions to the protostars, whereas such effects are small for other species. Assuming that all molecules form in the pre-stellar ices, the scatter variations could be explained by differences in lifetimes or physical conditions of the pre-stellar clouds. If the pre-stellar lifetimes are the main factor, they should be similar for low- and high-mass protostars., Comment: 49 pages, 46 figures, Accepted for publication in A&A

Some protostellar systems show little or no millimetre line emission of complex organics. This can be interpreted as a low abundance of these molecules, alternatively they could be present in the system but are not seen in the gas. The goal is to investigate the latter hypothesis for methanol. We will attempt to answer the question: Does the presence of a disk and optically thick dust reduce methanol emission even if methanol is abundant in the ices and gas? Using the radiative transfer code RADMC-3D, methanol emission lines from an envelope-only model and an envelope-plus-disk model are calculated and compared with each other and the observations. Methanol emission from the envelope-only model is always stronger than from the envelope-plus-disk model by at least a factor 2 as long as the disk radius is larger than 30 au (for L=8 L$_{\odot}$). In most cases, this is due to lower temperatures and, hence, the smaller amount of warm methanol inside the snow surface of the envelope-plus-disk model. The intensities drop by more than an order of magnitude for models including high mm opacity dust grains and disk radii of at least 50 au (for L=8 L$_{\odot}$) due to continuum over-subtraction. The line intensities from the envelope-only models overproduce the observations of protostars with lower methanol emission even with large dust optical depth effects. The envelope-plus-disk models can explain the bulk of the observations. However, they can only reproduce the observations of sources with high luminosities and low methanol emission when dust optical depth effects especially continuum over-subtraction in the disk becomes effective. Therefore, both the effects of disk and dust optical depth should be considered to explain the observations. In conclusion: Absence of methanol emission does not imply absence of methanol molecules in either gas or ice., Comment: Accepted to A&A

Complex organic molecules (COMs) are often observed toward embedded Class 0 and I protostars. However, not all Class 0 and I protostars exhibit COMs emission. In this work, variations in methanol (CH$_3$OH) emission are studied to test if absence of CH$_3$OH emission can be linked to source properties. Combining both new and archival observations with ALMA and sources from the literature, a sample of 184 low-mass and high-mass protostars is investigated. The warm (T > 100 K) gaseous CH$_3$OH mass, $M_{\rm CH_3OH}$, is determined for each source using primarily optically thin isotopologues. On average, Class I protostellar systems seem to have less warm $M_{\rm CH_3OH}$ ($<10^{-10}$ M$_\odot$) than younger Class 0 sources ($\sim10^{-7}$ M$_\odot$). High-mass sources in our sample show higher warm $M_{\rm CH_3OH}$ up to $10^{-7}-10^{-3}$ M$_\odot$. To take into account the effect of the source's overall mass on $M_{\rm CH_3OH}$, a normalized CH$_3$OH mass is defined as $M_{\rm CH_3OH}/M_{\rm dust,0}$, where $M_{\rm dust,0}$ is the cold + warm dust mass within a fixed radius. Excluding upper limits, a simple power-law fit to the normalized warm CH$_3$OH masses results in $M_{\rm CH_3OH}/M_{\rm dust,0}\propto L_{\rm bol}^{0.70\pm0.05}$. This is in good agreement with a simple hot core toy model which predicts that the normalized $M_{\rm CH_3OH}$ increases with $L_{\rm bol}^{0.75}$ due to the snowline moving outward. Sources for which the size of the disk is equivalent or smaller than the estimated 100 K radius agree well with the best-fit power-law model, whereas sources with significantly larger disks show up to two orders of magnitude lower normalized warm CH$_3$OH masses. Based on the latter results, we suggest that source structure such as a disk can result in colder gas and thus fewer COMs in the gas phase. Additionally, optically thick dust can hide the emission of COMs., Comment: 22 pages, 6 figures, accepted for publication in A&A

(Abridged) Recent surveys of young star formation regions have shown that the average Class II object does not have enough dust mass to make the cores of giant planets. Younger Class 0/I objects have enough dust in their embedded disk, which begs the questions: can the first steps of planet formation occur in these younger systems? The first step is building the first planetesimals, generally believed to be the product of the streaming instability. Hence the question can be restated: are the physical conditions of embedded disks conducive to the growth of the streaming instability? Here we model the collapse of a `dusty' proto-stellar cloud to show that if there is sufficient drift between the falling gas and dust, regions of the embedded disk can become sufficiently enhanced in dust to drive the streaming instability. We include four models, three with different dust grain sizes and one with a different initial cloud angular momentum to test a variety of collapse trajectories. We find a `sweet spot' for planetesimal formation for grain sizes of a few 10s of micron since they fall sufficiently fast relative to the gas to build a high dust-to-gas ratio along the disk midplane, but have slow enough radial drift speeds in the embedded disk to maintain the high dust-to-gas ratio. Unlike the gas, which is held in hydrostatic equilibrium for a time due to gas pressure, the dust can begin collapsing from all radii at a much earlier time. The streaming instability can produce at least between 7-35 M$_\oplus$ of planetesimals in the Class 0/I phase of our smooth embedded disks, depending on the size of the falling dust grains. This first generation of planetesimals could represent the first step in planet formation, and occurs earlier in the lifetime of the young star than is traditionally thought., Comment: Accepted for publication in A&A, 22 pages, 19 figures

As material from an infalling protostellar envelope hits the forming disk, an accretion shock may develop which could (partially) alter the envelope material entering the disk. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) indicate that emission originating from warm SO and SO$_2$ might be good tracers of such accretion shocks. The goal of this work is to test under what shock conditions the abundances of gas-phase SO and SO$_2$ increase in an accretion shock at the disk-envelope interface. Detailed shock models including gas dynamics are computed using the Paris-Durham shock code for non-magnetized J-type accretion shocks in typical inner envelope conditions. The effect of pre-shock density, shock velocity, and strength of the ultraviolet (UV) radiation field on the abundance of warm SO and SO$_2$ is explored. Warm gas-phase chemistry is efficient in forming SO under most J-type shock conditions considered. In lower-velocity (~3 km/s) shocks, the abundance of SO is increased through subsequent reactions starting from thermally desorbed CH$_4$ toward H$_2$CO and finally SO. In higher velocity (>4 km/s) shocks, both SO and SO$_2$ are formed through reactions of OH and atomic S. The strength of the UV radiation field is crucial for SO and in particular SO$_2$ formation through the photodissociation of H$_2$O. Thermal desorption of SO and SO$_2$ ice is only relevant in high-velocity (>5 km/s) shocks at high densities ($>10^7$ cm$^{-3}$). Warm emission from SO and SO$_2$ is a possible tracer of accretion shocks at the disk-envelope interface as long as a local UV field is present. Additional observations with ALMA at high-angular resolution could provide further constraints. Moreover, the James Webb Space Telescope will give access to other possible slow, dense shock tracers such as H$_2$, H$_2$O, and [S I] 25$\mu$m., Comment: Accepted to A&A, 18 pages, 17 figures

The chemical inventory of planets is determined by the physical and chemical processes that govern the early phases of star formation. The aim is to investigate N-bearing complex organic molecules towards two Class 0 protostars (B1-c and S68N) at millimetre wavelengths with ALMA. Next, the results of the detected N-bearing species are compared with those of O-bearing species for the same and other sources. ALMA observations in Band 6 ($\sim$ 1 mm) and Band 5 ($\sim$ 2 mm) are studied at $\sim$ 0.5" resolution, complemented by Band 3 ($\sim$ 3 mm) data in a $\sim$ 2.5" beam. NH2CHO, C2H5CN, HNCO, HN13CO, DNCO, CH3CN, CH2DCN, and CHD2CN are identified towards the investigated sources. Their abundances relative to CH3OH and HNCO are similar for the two sources, with column densities that are typically an order of magnitude lower than those of O-bearing species. The largest variations, of an order of magnitude, are seen for NH2CHO abundance ratios with respect to HNCO and CH3OH and do not correlate with the protostellar luminosity. In addition, within uncertainties, the N-bearing species have similar excitation temperatures to those of O-bearing species ($\sim$ 100 $\sim$ 300 K). The similarity of most abundances with respect to HNCO, including those of CH2DCN and CHD2CN, hints at a shared chemical history, especially the high D/H ratio in cold regions prior to star formation. However, some of the variations in abundances may reflect the sensitivity of the chemistry to local conditions such as temperature (e.g. NH2CHO), while others may arise from differences in the emitting areas of the molecules linked to their different binding energies in the ice. The two sources discussed here add to the small number of sources with such a detailed chemical analysis on Solar System scales. Future JWST data will allow a direct comparison between the ice and gas abundances of N-bearing species., Comment: Accepted to A&A, 41 pages, 37 figures

As a part of the CALYPSO large programme, we constrain the properties of protostellar jets and outflows in a sample of 21 Class 0 protostars with internal luminosities, Lint, from 0.035 to 47 Lsun. We analyse high angular resolution (~0.5"-1") IRAM PdBI observations in CO (2-1), SO ($5_6-4_5$), and SiO (5-4). CO (2-1), which probes outflowing gas, is detected in all the sources (for the first time in SerpS-MM22 and SerpS-MM18b). Collimated high-velocity jets in SiO (5-4) are detected in 67% of the sources (for the first time in IRAS4B2, IRAS4B1, L1448-NB, SerpS-MM18a), and 77% of these also show jet/outflow emission in SO ($5_6-4_5$). In 5 sources (24% of the sample) SO ($5_6-4_5$) probes the inner envelope and/or the disk. The CALYPSO survey shows that the outflow phenomenon is ubiquitous and that the detection rate of high-velocity jets increases with protostellar accretion, with at least 80% of the sources with Lint>1 Lsun driving a jet. The protostellar flows exhibit an onion-like structure, where the SiO jet (opening angle ~10$^o$) is nested into a wider angle SO (~15$^o$) and CO (~25$^o$) outflow. On scales >300 au the SiO jets are less collimated than atomic jets from Class II sources (~3$^o$). Velocity asymmetry between the two jet lobes are detected in one third of the sources, similarly to Class II atomic jets, suggesting that the same launching mechanism is at work. Most of the jets are SiO rich (SiO/H2 from >2.4e-7 to >5e-6), which indicates efficient release of >1%-10% of silicon in gas phase likely in dust-free winds, launched from inside the dust sublimation radius. The mass-loss rates (from ~7e-8 to ~3e-6 Msun/yr) are larger than what was measured for Class II jets. Similarly to Class II sources, the mass-loss rates are ~1%-50% of the mass accretion rates suggesting that the correlation between ejection and accretion in young stars holds from 1e4 yr up to a few Myr., Comment: Accepted for publication by A&A, 40 pages, 7 figures + appendixes

Complex organic molecules (COMs) are thought to form on icy dust grains in the earliest phase of star formation. The evolution of these COMs from the youngest Class 0/I protostellar phases toward the more evolved Class II phase is still not fully understood. Since planet formation seems to start early, and mature disks are too cold for characteristic COM emission lines, studying the inventory of COMs on solar system scales in the Class 0/I stage is relevant. ALMA Band 3 (3 mm) and Band 6 (1 mm) observations are obtained of seven Class 0 protostars in the Perseus and Serpens star-forming regions. By modeling the inner protostellar region using 'LTE' models, the excitation temperature and column densities are determined for several O-bearing COMs. B1-c, B1-bS, and Serpens S68N show COM emission, i.e, three out of the seven sources. No clear correlation seems to exist between the occurrence of COMs and source luminosity. The abundances of several COMs with respect to CH3OH are remarkably similar for the three COM-rich sources, and to IRAS 16293-2422B and HH 212. For other COMs the abundances differ by up to an order of magnitude, indicating that local source conditions are case determining. B1-c hosts a cold ($T_{ex}\approx60$ K), more extended component of COM emission with a column density of typically a few % of the warm/hot ($T_{ex}\sim 200$ K), central component. A D/H ratio of 1-3 % is derived based on the CH2DOH/CH3OH ratio suggesting a temperature of $\sim$15~K during the formation of methanol. This ratio is consistent with other low-mass protostars. Future mid-infrared facilities such as JWST/MIRI will be essential to directly observe COM ices. Combining this with a larger sample of COM-rich sources with ALMA will allow for directly linking ice and gas-phase abundances in order to constrain the routes that produce and maintain chemical complexity during the star formation process., Comment: Accepted to A&A, 44 pages, 17 figures

Large millimeter interferometers are revealing a growing number of rotating outflows, which are suggested to trace magneto-centrifugal disk winds (MHD DWs). However, their impact on disk accretion is not yet well quantified. Here we identify systematic biases in retrieving the true launch zone, magnetic lever arm, and angular momentum flux of an MHD DW from apparent rotation signatures. Synthetic position-velocity cuts are constructed from self-similar MHD DWs over a broad range of parameters, and three different methods are applied for estimating the specific angular momentum. We find that the launch radius inferred using the well-known relation from Anderson et al. (2006) can markedly differ from the true outermost launch radius $r_{out}$ of the DW. The "double-peak separation" and "flow width" methods provide only a strict lower limit to $r_{out}$. This bias is independent of angular resolution and can reach a factor ten. In contrast, the "rotation curve" method gives a good estimate of $r_{out}$ when the flow is well resolved, and an upper limit otherwise. The magnetic lever arm is always underestimated. Only comparison with synthetic predictions can take into account properly all observational effects. As an application, we present a comparison with ALMA observations of HH212 at resolutions from 250 au to 16 au, which represents the most stringent observational test of MHD DW to date. This comparison confirms our predicted biases for the double-peak separation method, and the large $r_{out}\sim40~$au and small magnetic lever arm first suggested by Tabone et al. (2017). We also derive the first accurate analytical expression for the fraction of disk angular momentum extracted by an MHD disk wind of given radial extent, magnetic lever arm, and mass flux. Application to HH212 confirms that MHD DWs are serious candidates for the steady angular momentum extraction process in young disks., Comment: Accepted for publication in A&A

Recent ALMA observations suggest that the highest velocity part of molecular protostellar jets are launched from the dust-sublimation regions of the accretion disks (<0.3 au). However, formation and survival of molecules in inner protostellar disk winds, in the presence of a harsh FUV radiation field and the absence of dust, remain unexplored. We aim at determining if simple molecules can be synthesized and spared in fast and collimated dust-free disk winds or if a fraction of dust is necessary to explain the observed molecular abundances. This work is based on the Paris-Durham shock code designed to model irradiated environments. Fundamental properties of the dust-free chemistry are investigated from single point models. A laminar 1D disk wind model is then built using a parametric flow geometry. This model includes time-dependent chemistry and the attenuation of the radiation field by gas-phase photoprocesses. We show that a small fraction of H2 (< 1e-2), primarily formed through the H- route, can efficiently initiate molecule synthesis such as CO and SiO above TK ~ 800 K. The attenuation of the radiation field by atomic species (eg. C, Si, S) proceeds through continuum self-shielding. This process ensures efficient formation of CO, OH, SiO, H2O through neutral-neutral reactions, and the survival of these molecules. Class 0 dust-free winds with high mass-loss rates ($\dot{M}_w >$ 2e-6 Msun/yr) are predicted to be rich in molecules if warm (TK > 800 K). The molecular content of disk winds is very sensitive to the presence of dust and a mass-fraction of surviving dust as small as 1e-5 significantly increases the H2O and SiO abundances. Chemistry of high-velocity jets is a powerful tool to probe their content in dust and uncover their launching point. Models of internal shocks are required to fully exploit the current (sub-)millimeter observations and prepare future JWST observations., Comment: accepted for publication in A&A

Understanding the formation mechanisms of protoplanetary disks and multiple systems, and their pristine properties, is a key question for modern astrophysics. The properties of the youngest disks, embedded in rotating infalling protostellar envelopes, have largely remained unconstrained up to now. In the framework of the IRAM-PdBI CALYPSO survey, we have obtained sub-arcsecond observations of the dust continuum emission at 231 GHz and 94 GHz, for a sample of 16 solar-type Class 0 protostars. In an attempt to identify disk-like structures embedded at small scales in the protostellar envelopes, we model the dust continuum emission visibility profiles using both Plummer-like envelope models and envelope models including additional Gaussian disk-like components. Our analysis shows that in the CALYPSO sample, 11 of the 16 Class 0 protostars are better reproduced by models including a disk-like dust continuum component contributing to the flux at small scales, but less than 25% of these candidate protostellar disks are resolved at radii > 60 au. Including all available literature constraints on Class 0 disks at subarcsecond scales, we show that our results are representative: most (> 72% in a sample of 26 protostars) Class 0 protostellar disks are small and emerge only at radii < 60 au. Our multiplicity fraction at scales 100-5000 au is in global agreement with the multiplicity properties of Class I protostars at similar scales. We confront our observational constraints on the disk size distribution in Class 0 protostars to the typical disk properties from protostellar formation models. Because they reduce the centrifugal radius, and produce a disk size distribution peaking at radii <100 au during the main accretion phase, magnetized models of rotating protostellar collapse are favored by our observations., Comment: Accepted for publication in Astronomy \& Astrophysics. 24 pages of main text, 22 pages of additional content (under the form of three appendices) - revised version including language and typos corrections

The HH 212 protostellar system, in Orion B, has been mapped thanks to ALMA-Band 7 Cycle 1 and Cycle 4 observations of dueterated water (HDO) and acetaldehyde (CH$_3$CHO) emission with an angular resolution down to $\sim$0.15 arcsec (60 au). Many emission lines due to 14 CH$_3$CHO and 1 HDO transitions at high excitation ($E_{\rm u}$ between 163 K and 335 K) have been imaged in the inner $\sim$ 70 au region. The local thermal equilibrium analysis of the CH$_3$CHO emission leads to a temperature of 78$\pm$14 K and a column density of 7.6$\pm$3.2 $\times$ 10$^{15}$ cm$^{-2}$, which, when $N_{\rm H_2}$ of 10$^{24}$ cm$^{-2}$ is assumed, leads to an abundance of $X_{\rm CH_3CHO}$ $\simeq$ 8 $\times$ 10$^{-9}$. The large velocity gradient analysis of the HDO emission also places severe constraints on the volume density, n$_{\rm H_2}$ $\geq$ 10$^8$ cm$^{-3}$. The line profiles are 5--7 km s$^{-1}$ wide, and CH$_3$CHO and HDO both show a $\pm$ 2 km s$^{-1}$ velocity gradient over a size of $\sim$ 70 au (blue-shifted emission towards the north-west and red-shifted emission towards the south-east) along the disc equatorial plane, in agreement with what was found so far using other molecular tracers. The kinematics of CH$_3$CHO and HDO are consistent with the occurrence of a centrifugal barrier, that is, the infalling envelope-rotating disc ring, which is chemically enriched through low-velocity accretion shocks. The emission radius is $\sim$ 60 au, in good agreement with what was found before for another interstellar complex organic molecule such as NH$_2$CHO. We support a vertical structure for the centrifugal barrier, suggesting the occurrence of two outflowing, expanding, and rotating rings above and below (of about 40-45 au) the optically thick equatorial disc plane. It is tempting to speculate that these rings could probe the basis of a wind launched from this region., Comment: Astronomy & Astrophysics, in press

When interacting, the solar wind and the ionised atmosphere of a comet exchange energy and momentum. Our aim is to understand the influence of the average Parker spiral configuration of the solar wind magnetic field on this interaction. We compare the theoretical expectations of an analytical generalised gyromotion with Rosetta observations at comet 67P/Churyumov-Gerasimenko. A statistical approach allows one to overcome the lack of upstream solar wind measurement. We find that additionally to their acceleration along (for cometary pick-up ions) or against (for solar wind ions) the upstream electric field orientation and sense, the cometary pick-up ions are drifting towards the dawn side of the coma, while the solar wind ions are drifting towards the dusk side of the coma, independent of the heliocentric distance. The dynamics of the interaction is not taking place in a plane, as often assumed in previous works., Comment: 6 pages, 3 figures

We aim at analytically modelling the solar wind proton trajectories during their interaction with a partially ionised cometary atmosphere, not in terms of bulk properties of the flow but in terms of single particle dynamics. We first derive a generalised gyromotion, in which the electric field is reduced to its motional component. Steady-state is assumed, and simplified models of the cometary density and of the electron fluid are used to express the force experienced by individual solar wind protons during the interaction. A three-dimensional (3D) analytical expression of the gyration of two interacting plasma beams is obtained. Applying it to a comet case, the force on protons is always perpendicular to their velocity and has an amplitude proportional to 1/r^2. The solar wind deflection is obtained at any point in space. The resulting picture presents a caustic of intersecting trajectories, and a circular region is found that is completely free of particles. The particles do not lose any kinetic energy and this absence of deceleration, together with the solar wind deflection pattern and the presence of a solar wind ion cavity, is in good agreement with the general results of the Rosetta mission. The qualitative match between the model and the in situ data highlights how dominant the motional electric field is throughout most of the interaction region for the solar wind proton dynamics. The model provides a simple general kinetic description of how momentum is transferred between these two collisionless plasmas. It also shows the potential of this semi-analytical model for a systematic quantitative comparison to the data., Comment: 12 pages, 6 figures

The first 1000 km of the ion tail of comet 67P/Churyumov-Gerasimenko were explored by the European Rosetta spacecraft, 2.7 au away from the Sun. We characterised the dynamics of both the solar wind and the cometary ions on the night-side of the comet's atmosphere. We analysed in situ ion and magnetic field measurements and compared the data to a semi-analytical model. The cometary ions are observed flowing close to radially away from the nucleus during the entire excursion. The solar wind is deflected by its interaction with the new-born cometary ions. Two concentric regions appear, an inner region dominated by the expanding cometary ions and an outer region dominated by the solar wind particles. The single night-side excursion operated by Rosetta revealed that the near radial flow of the cometary ions can be explained by the combined action of three different electric field components, resulting from the ion motion, the electron pressure gradients, and the magnetic field draping. The observed solar wind deflection is governed mostly by the motional electric field., Comment: 11 pages, 7 figures

The molecular richness of fast protostellar jets within 20-100 au of their source, despite strong ultraviolet irradiation, remains a challenge for the models investigated so far. We aim to investigate the effect of interaction between a time-variable jet and a surrounding steady disk wind, to assess the possibility of jet chemical enrichement by the wind, and the characteristic signatures of such a configuration. We have constructed an analytic model of a jet bow shock driven into a surrounding slower disk wind in the thin shell approximation. The refilling of the post bow shock cavity from below by the disk wind is also studied. An extension of the model to the case of two or more successive internal working surfaces (IWS) is made. We then compared this analytic model with numerical simulations with and without a surrounding disk wind. We find that at early times (of order the variability period), jet bow shocks travel in refilled pristine disk-wind material, before interacting with the cocoon of older bow shocks. This opens the possibility of bow shock chemical enrichment (if the disk wind is molecular and dusty) and of probing the unperturbed disk wind structure near the jet base. Several distinctive signatures of the presence of a surrounding disk wind are identified, in the bow shock morphology and kinematics. Numerical simulations validate our analytical approach and further show that at large scale, the passage of many jet IWS inside a disk wind produces a stationary V-shaped cavity, closing down onto the axis at a finite distance from the source., Comment: Accepted for publication in A&A. Version uncorrected for language

We wish to constrain the possible contribution of a magnetohydrodynamic disk wind (DW) to the HH212 molecular jet. We mapped the flow base with ALMA Cycle 4 at 0.13"$\sim$ 60 au resolution and compared these observations with synthetic DW predictions. We identified, in SO/SO$_2$, a rotating flow that is wider and slower than the axial SiO jet. The broad outflow cavity seen in C$^{34}$S is not carved by a fast wide-angle wind but by this slower agent. Rotation signatures may be fitted by a DW of a moderate lever arm launched out to $\sim$ 40 au with SiO tracing dust-free streamlines from 0.05-0.3 au. Such a DW could limit the core-to-star efficiency to $\leq$50%., Comment: Accepted for publication in A&A letters, final language edited version

Context: Methanol is thought to be mainly formed during the prestellar phase and its deuterated form keeps memory of the conditions at that epoch. Thanks to the unique combination of high angular resolution and sensitivity provided by ALMA, we wish to measure methanol deuteration in the planet formation region around a Class 0 protostar and to understand its origin. Aims: We mapped both the $^{13}$CH$_3$OH and CH$_2$DOH distribution in the inner regions ($\sim$100 au) of the HH212 system in Orion B. To this end, we used ALMA Cycle 1 and Cycle 4 observations in Band 7 with angular resolution down to $\sim$0.15$"$. Results: We detected 6 lines of $^{13}$CH$_3$OH and 13 lines of CH$_2$DOH with upper level energies up to 438 K in temperature units. We derived a rotational temperature of (171 $\pm$ 52) K and column densities of 7$\times$10$^{16}$ cm$^{-2}$ ($^{13}$CH$_3$OH) and 1$\times$10$^{17}$ cm$^{-2}$ (CH$_2$DOH), respectively. Consequently, the D/H ratio is (2.4 $\pm$ 0.4)$\times$10$^{-2}$, a value lower by an order of magnitude with respect to what was previously measured using single dish telescopes toward protostars located in Perseus. Our findings are consistent with the higher dust temperatures in Orion B with respect to that derived for the Perseus cloud. The emission is tracing a rotating structure extending up to 45 au from the jet axis and elongated by 90 au along the jet axis. So far, the origin of the observed emission appears to be related with the accretion disk. Only higher spatial resolution measurements however, will be able to disentangle between different possible scenarios: disk wind, disk atmosphere, or accretion shocks., Comment: Astronomy & Astrophysics Letter, in press

Fast jets are thought to be a crucial ingredient of star formation because they might extract angular momentum from the disk and thus allow mass accretion onto the star. However, it is unclear whether jets are ubiquitous, and likewise, their contribution to mass and angular momentum extraction during protostar formation remains an open question. Our aim is to investigate the ejection process in the low-mass Class 0 protostar L1157. This source is associated with a spectacular bipolar outflow, and the recent detection of high-velocity SiO suggests the occurrence of a jet. Observations of CO 2-1 and SiO 5-4 at 0.8" resolution were obtained with the IRAM Plateau de Bure Interferometer as part of the CALYPSO large program. The jet and outflow structure were fit with a precession model. We derived the column density of CO and SiO, as well as the jet mass-loss rate and mechanical luminosity. High-velocity CO and SiO emission resolve for the first time the first 200 au of the outflow-driving molecular jet. The jet is strongly asymmetric, with the blue lobe 0.65 times slower than the red lobe. This suggests that the large-scale asymmetry of the outflow is directly linked to the jet velocity and that the asymmetry in the launching mechanism has been at work for the past 1800 yr. Velocity asymmetries are common in T Tauri stars, which suggests that the jet formation mechanism from Class 0 to Class II stages might be similar. Our model simultaneously fits the inner jet and the clumpy 0.2 pc scale outflow by assuming that the jet precesses counter-clockwise on a cone inclined by 73 degree to the line of sight with an opening angle of 8 degree on a period of 1640 yr. The estimated jet mass flux and mechanical luminosity are 7.7e-7 Msun/yr, and 0.9 Lsun, indicating that the jet could extract at least 25% of the gravitational energy released by the forming star., Comment: 6 pages, 3 figures, accepted for publication in A&A Letter