Search

The MIRI-MRS instrument onboard JWST allows for probing the molecular gas composition at mid-IR wavelengths at unprecedented resolution and sensitivity. It is important to study these features in low-mass embedded protostellar systems since the formation of planets is thought to start in this phase. We present JWST/MIRI-MRS data of 18 low-mass protostellar systems in the JOYS program, focusing on gas-phase molecular lines in spectra extracted from the central protostellar positions. Besides H2, the most commonly detected molecules are H2O, CO2, CO, and OH. Other molecules such as 13CO2, C2H2, 13CCH, HCN, C4H2, CH4, and SO2 are detected only toward at most three of the sources. The JOYS data also yield the surprising detection of SiO gas toward two sources (BHR71-IRS1, L1448-mm) and for the first time CS and NH3 at mid-IR wavelengths toward a low-mass protostar (B1-c). The temperatures derived for the majority of the molecules are 100-300 K, much lower than what is typically derived toward more evolved Class II sources (>500 K). Toward three sources (e.g., TMC1-W), hot (~1000 K) H2O is detected, indicative of the presence of hot molecular gas in the embedded disks, but such warm emission from other molecules is absent. The agreement in abundance ratios with respect to H2O between ice and gas point toward ice sublimation in a hot core for a few sources (e.g., B1-c) whereas their disagreement and velocity offsets hint at high-temperature (shocked) conditions toward other sources (e.g., L1448-mm, BHR71-IRS1). The typical temperatures of the gas-phase molecules of 100-300 K are consistent with both ice sublimation in hot cores as well as high-temperature gas phase chemistry. Molecular features originating from the inner embedded disks are not commonly detected, likely because they are too extincted even at mid-IR wavelengths by small not-settled dust grains in upper layers of the disk., Comment: 25 pages, 11 figures, accepted for publication in A&A

The carbon isotope ratio is a powerful tool for studying the evolution of stellar systems. Recent detections of CO isotopologues in disks and exoplanet atmospheres pointed towards significant fractionation in these systems. In order to understand the evolution of this quantity, it is crucial to trace the isotope abundance from stellar nurseries to planetary systems. During the protostellar stage the multiple vibrational modes of CO$_2$ and CO ice provide a unique opportunity to examine the carbon isotope ratio in the solid state. Now with the sensitivity of the \textit{James Webb Space Telescope}, these absorption features have become accessible at high S/N in Solar-mass systems. We quantify the $^{12}$CO$_2$/$^{13}$CO$_2$ and the $^{12}$CO/$^{13}$CO isotope ratios in 17 class 0/I low mass protostars from the $^{12}$CO$_2$ combination modes (2.70 ${\mu}$m and 2.77 ${\mu}$m), the $^{12}$CO$_2$ stretching mode (4.27 ${\mu}$m), the $^{13}$CO$_2$ stretching mode (4.39 ${\mu}$m), the $^{12}$CO$_2$ bending mode (15.2 ${\mu}$m), the $^{12}$CO stretching mode (4.67 ${\mu}$m) and the $^{13}$CO stretching mode (4.78 ${\mu}$m) using JWST observations. We also report a detection of the $^{12}$CO overtone mode at 2.35 ${\mu}$m. The $^{12}$CO$_2$/$^{13}$CO$_2$ ratios are in agreement and we find mean ratios of 85 $\pm$ 23, 76 $\pm$ 12 and 97 $\pm$ 17 for the 2.70 ${\mu}$m, 4.27 ${\mu}$m and the 15.2 ${\mu}$m bands, respectively. The main source of uncertainty stem from the error on the band strengths. The $^{12}$CO/$^{13}$CO ratios derived from the 4.67 ${\mu}$m bands are consistent, albeit elevated with respect to the $^{12}$CO$_2$/$^{13}$CO$_2$ ratios and we find a mean ratio of 165 $\pm$ 52. These findings indicate that ices leave the pre-stellar stage with elevated carbon isotope ratios relative to the interstellar medium and that fractionation becomes significant during the later stages.

We use the James Webb Space Telescope (JWST) and its Mid-Infrared Instrument (MIRI) (5-28 um), to study the embedded HH 211 flow. We map a 0.95'x0.22' region, covering the full extent of the blue-shifted lobe, the central protostellar region, and a small portion of the red-shifted lobe. The jet driving source is not detected even at the longest mid-IR wavelengths. The overall morphology of the flow consists of a highly collimated jet, mostly molecular (H2, HD) with an inner atomic ([FeI], [FeII], [SI], [NiII]) structure. The jet shocks the ambient medium, producing several large bow-shocks, rich in forbidden atomic and molecular lines, and is driving an H2 molecular outflow, mostly traced by low-J, v=0 transitions. Moreover, 0-0 S(1) uncollimated emission is also detected down to 2"-3" (~650-1000 au) from the source, tracing a cold (T=200-400 K), less dense and poorly collimated molecular wind. The atomic jet ([FeII] at 26 um) is detected down to ~130 au from source, whereas the lack of H2 emission close to the source is likely due to the large visual extinction. Dust continuum-emission is detected at the terminal bow-shocks, and in the blue- and red-shifted jet, being likely dust lifted from the disk. The jet shows an onion-like structure, with layers of different size, velocity, temperature, and chemical composition. Moreover, moving from the inner jet to the outer bow-shocks, different physical, kinematic and excitation conditions for both molecular and atomic gas are observed. The jet mass-flux rate, momentum, and momentum flux of the warm H2 component are up to one order of magnitude higher than those inferred from the atomic jet component. Our findings indicate that the warm H2 component is the primary mover of the outflow, namely it is the most significant dynamical component of the jet, in contrast to jets from more evolved YSOs, where the atomic component is dominant., Comment: Paper accepted in A&A for publication

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

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

Complex organic molecules (COMs) detected in the gas phase are thought to be mostly formed on icy grains, but no unambiguous detection of icy COMs larger than CH3OH has been reported so far. Exploring this matter in more detail has become possible with the JWST the critical 5-10 $\mu$m range. In the JOYS+ program, more than 30 protostars are being observed with the MIRI/MRS. This study explores the COMs ice signatures in the low and high-mass protostar, IRAS 2A and IRAS 23385, respectively. We fit continuum and silicate subtracted observational data with IR laboratory ice spectra. We use the ENIIGMA fitting tool to find the best fit between the lab data and the observations and to performs statistical analysis of the solutions. We report the best fits for the spectral ranges between 6.8 and 8.6 $\mu$m in IRAS 2A and IRAS 23385, originating from simple molecules, COMs, and negative ions. The strongest feature in this range (7.7 $\mu$m) is dominated by CH4 and has contributions of SO2 and OCN-. Our results indicate that the 7.2 and 7.4 $\mu$m bands are mostly dominated by HCOO-. We find statistically robust detections of COMs based on multiple bands, most notably CH3CHO, CH3CH2OH, and CH3OCHO. The likely detection of CH3COOH is also reported. The ice column density ratios between CH3CH2OH and CH3CHO of IRAS 2A and IRAS 23385, suggests that these COMs are formed on icy grains. Finally, the derived ice abundances for IRAS 2A correlate well with those in comet 67P/GC within a factor of 5. Based on the MIRI/MRS data, we conclude that COMs are present in interstellar ices, thus providing additional proof for a solid-state origin of these species in star-forming regions. The good correlation between the ice abundances in comet 67P and IRAS 2A is in line with the idea that cometary COMs can be inherited from the early protostellar phases., Comment: Accepted for publication in A&A

It is a challenge to study the formation process of massive stars: their formation time is short, they are few, often deeply embedded, and at relatively large distances. Our strategy is to study the outcome of the star formation process and to look for signatures remnant of the formation. We have access to a unique sample of (massive) pre-main-sequence (PMS) stars in the giant HII region M17, showing a photosphere and circumstellar disk. The aim is to determine the variability properties of the hot gaseous disks to understand the physical origin of the emission lines and identify dominant physical processes in these disks. We have obtained multiple-epoch (4-5 epochs) VLT/X-shooter spectra of six young stars in M17 covering about a decade. Using stacked spectra we update the spectral classification and identify circumstellar features. With the temporal variance method (TVS) we determine the extent and amplitude of the spectral line variations. The double-peaked emission lines in the PMS stars with gaseous disks are used to determine peak-to-peak velocities, V/R-ratios and the radial velocity of the systems. We identify many disk features, under which a new detection of CO bandhead and CI emission. In three of the stars we detect spectral variability, mainly in lines originating in the circumstellar disk, in a velocity range up to 320 km/s. In two PMS stars the ratio between the blue and red peaks shows a correlation with the peak-to-peak velocity, possibly explained by a spiral-arm structure in the disk. The PMS stars with variability are at similar positions in the HRD but show significant differences in disk lines and variability. The extent and timescale of the variability differs for each star and per line (sets). We find indications for an accretion flow, slow disk winds and/or disk structures in the hot gaseous inner disk as the cause of the variability in these PMS stars., Comment: 27 pages, 24 figures, accepted for publication in Astronomy and Astrophysics, abstract abbreviated

(abridged) We study and compare the warm (>100 K) and cold (<100 K) material toward the high-mass star-forming region IRAS 23385+6053 (IRAS 23385 hereafter) combining high angular resolution observations in the mid-infrared (MIR) with the JWST Observations of Young protoStars (JOYS) project and with the NOEMA at mm wavelengths at angular resolutions of 0.2"-1". The spatial morphology of atomic and molecular species is investigated by line integrated intensity maps. The temperature and column density of different gas components is estimated using H2 transitions (warm and hot component) and a series of CH3CN transitions as well as 3 mm continuum emission (cold component). Toward the central dense core in IRAS 23385 the material consists of relatively cold gas and dust (~50 K), while multiple outflows create heated and/or shocked H2 and show enhanced temperatures (~400 K) along the outflow structures. An energetic outflow with enhanced emission knots of [Fe II] and [Ni II] hints at J-type shocks, while two other outflows have enhanced emission of only H2 and [S I] caused by C-type shocks. The latter two outflows are also more prominent in molecular line emission at mm wavelengths (e.g., SiO, SO, H2CO, and CH3OH). Even higher angular resolution data are needed to unambiguously identify the outflow driving sources given the clustered nature of IRAS 23385. While most of the forbidden fine structure transitions are blueshifted, [Ne II] and [Ne III] peak at the source velocity toward the MIR source A/mmA2 suggesting that the emission is originating from closer to the protostar., Comment: 15 pages, 7 figures, accepted for publication in A&A

Complex organic molecules (COMs) have been observed to be abundant in the gas phase toward protostars. Deep line surveys have been carried out only for a limited number of well-known high-mass star forming regions using the Atacama Large Millimeter/submillimeter Array (ALMA), which has unprecedented resolution and sensitivity. Statistical studies on oxygen-bearing COMs (O-COMs) in high-mass protostars using ALMA are still lacking. With the recent CoCCoA survey, we are able to determine the column density ratios of six O-COMs with respect to methanol (CH$_3$OH) in a sample of 14 high-mass protostellar sources to investigate their origin through ice and/or gas-phase chemistry. The selected species are: acetaldehyde (CH$_3$CHO), ethanol (C$_2$H$_5$OH), dimethyl ether (DME, CH$_3$OCH$_3$), methyl formate (MF, CH$_3$OCHO), glycolaldehyde (GA, CH$_2$OHCHO), and ethylene glycol (EG, (CH$_2$OH)$_2$). DME and MF have the highest and most constant ratios within one order of magnitude, while the other four species have lower ratios and exhibit larger scatter by 1-2 orders of magnitude. We compare the O-COM ratios of high-mass CoCCoA sources with those of 5 low-mass protostars available from the literature, along with the results from experiments and simulations. We find that the O-COM ratios with respect to methanol are on the same level in both the high- and low-mass samples, which suggests that these species are mainly formed in similar environments during star formation, probably in ice mantles on dust grains during early pre-stellar stages. Current simulations and experiments can reproduce most observational trends with a few exceptions, and hypotheses exist to explain the differences between observations and simulations/experiments, such as the involvement of gas-phase chemistry and different emitting areas of molecules., Comment: Accepted for publication in Astronomy & Astrophysics. 17 pages and 7 figures in the main body, 15 pages and 8 figures in the appendix

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

Aims: The JWST program JOYS (JWST Observations of Young protoStars) aims at characterizing the physical and chemical properties of young high- and low-mass star-forming regions, in particular the unique mid-infrared diagnostics of the warmer gas and solid-state components. We present early results from the high-mass star formation region IRAS23385+6053. Methods: The JOYS program uses the MIRI MRS with its IFU to investigate a sample of high- and low-mass star-forming protostellar systems. Results: The 5 to 28mum MIRI spectrum of IRAS23385+6053 shows a plethora of features. While the general spectrum is typical for an embedded protostar, we see many atomic and molecular gas lines boosted by the higher spectral resolution and sensitivity compared to previous space missions. Furthermore, ice and dust absorption features are also present. Here, we focus on the continuum emission, outflow tracers like the H2, [FeII] and [NeII] lines as well as the potential accretion tracer Humphreys alpha HI(7--6). The short-wavelength MIRI data resolve two continuum sources A and B, where mid-infrared source A is associated with the main mm continuum peak. The combination of mid-infrared and mm data reveals a young cluster in its making. Combining the mid-infrared outflow tracer H2, [FeII] and [NeII] with mm SiO data shows a complex interplay of at least three molecular outflows driven by protostars in the forming cluster. Furthermore, the Humphreys alpha line is detected at a 3-4sigma level towards the mid-infrared sources A and B. Following Rigliaco et al. (2015), one can roughly estimate accretion luminosities and corresponding accretion rates between ~2.6x10^-6 and ~0.9x10^-4 M_sun/yr. This is discussed in the context of the observed outflow rates. Conclusions: The analysis of the MIRI MRS observations for this young high-mass star-forming region reveals connected outflow and accretion signatures., Comment: 12 pages, 9 figures, accepted for Astronomy & Astrophysics, the paper is also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

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

The deuteration of molecules forming in the ices such as methanol (CH$_3$OH) is sensitive to the physical conditions during their formation in dense cold clouds and can be probed through observations of deuterated methanol in hot cores. Observations with ALMA containing transitions of CH$_3$OH, CH$_2$DOH, CHD$_2$OH, $^{13}$CH$_3$OH, and CH$_3^{18}$OH are investigated. The column densities of CH$_2$DOH, CHD$_2$OH, and CH$_3$OH are determined for all sources, where the column density of CH$_3$OH is derived from optically thin $^{13}$C and $^{18}$O isotopologues. Consequently, the D/H ratio of methanol is derived taking statistical effects into account. Singly deuterated methanol (CH$_2$DOH) is detected toward 25 of the 99 sources in our sample of the high-mass protostars. Including upper limits, the $\rm (D/H)_{CH_3OH}$ ratio inferred from $N_\mathrm{CH_2DOH}/N_\mathrm{CH_3OH}$ was derived for 38 of the 99 sources and varies between $\sim10^{-3}-10^{-2}$. Including other high-mass hot cores from the literature, the mean methanol D/H ratio is $1.1\pm0.7\times10^{-3}$. This is more than one order of magnitude lower than what is seen for low-mass protostellar systems ($2.2\pm1.2\times10^{-2}$). Doubly deuterated methanol (CHD$_2$OH) is detected toward 11 of the 99 sources. Including upper limits for 15 sources, the $\rm (D/H)_{CH_2DOH}$ ratios derived from $N_\mathrm{CHD_2OH}/N_\mathrm{CH_2DOH}$ are more than two orders of magnitude higher than $\rm (D/H)_{CH_3OH}$ with an average of $2.0\pm0.8\times10^{-1}$ which is similar to what is found for low-mass sources. Comparison with literature GRAINOBLE models suggests that the high-mass prestellar phases are either warm ($>20$ K) or live shorter than the free-fall timescale. In contrast, for low-mass protostars, both a low temperature of $<15$ K and a prestellar phase timescale longer than the free-fall timescale are necessary., Comment: 17 pages, 9 figures, accepted to 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

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

[Abridged] Planet formation is expected to be enhanced around snowlines in protoplanetary disks, in particular around the water snowline. However, the close proximity of the water snowline to the host star and water in the Earth's atmosphere makes a direct detection of the water snowline in protoplanetary disks challenging. Following earlier work on protostellar envelopes, the aim of this research is to investigate the validity of HCO$^+$ and H$^{13}$CO$^+$, as tracers of the water snowline in protoplanetary disks, as HCO$^+$ is destroyed by gas-phase water. Two small chemical networks are used to predict the HCO$^+$ abundance in a typical Herbig Ae disk. Subsequently, the corresponding emission profiles are modelled for H$^{13}$CO$^+$ and HCO$^+$ $J=2-1$, which provides the best balance between brightness and optical depth effects of the continuum emission. The HCO$^+$ abundance jumps by two orders of magnitude just outside the water snowline at 4.5 AU. We find that the emission of H$^{13}$CO$^+$ and HCO$^+$ is ring-shaped due to three effects: destruction of HCO$^+$ by gas-phase water, continuum optical depth, and molecular excitation effects. The presence of gas-phase water causes an additional drop of only $\sim$13% and 24% in the center of the disk, for H$^{13}$CO$^+$ and HCO$^+$, respectively. For the much more luminous outbursting source V883Ori, our models predict that the effect of dust and excitation are not limiting if the snowline is located outside $\sim$40 AU. Our analysis of ALMA observations of HCO$^+$ $J=3-2$ is consistent with the water snowline located around 100 AU. The HCO$^+$ abundance drops steeply around the water snowline, but dust and excitation can conceal the drop in HCO$^+$ emission due to the water snowline. Therefore, locating the water snowline with HCO$^+$ in Herbig disks is very difficult, but it is possible for outbursting sources like V883Ori., Comment: Accepted for publication in Astronomy and Astrophysics. 13 pages, 11 figures + appendices

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

The process of massive star ($M\geq8~M_\odot$) formation is still poorly understood. Observations of massive young stellar objects (MYSOs) are challenging due to their rarity, short formation timescale, large distances, and high circumstellar extinction. Here, we present the results of a spectroscopic analysis of a population of MYSOs in the Large Magellanic Cloud (LMC). We took advantage of the spectral resolution and wavelength coverage of X-shooter (300-2500 nm), mounted on the European Southern Observatory Very Large Telescope, to detect characteristic spectral features in a dozen MYSO candidates near 30 Doradus, the largest starburst region in the Local Group hosting the most massive stars known. The X-shooter spectra are strongly contaminated by nebular emission. We used a scaling method to subtract the nebular contamination from our objects. We detect H$\alpha,\beta$, [O I] 630.0 nm, Ca II infrared triplet, [Fe II] 1643.5 nm, fluorescent Fe II 1687.8 nm, H$_2$ 2121.8 nm, Br$\gamma$, and CO bandhead emission in the spectra of multiple candidates. This leads to the spectroscopic confirmation of 10 candidates as bona fide MYSOs. We compare our observations with photometric observations from the literature and find all MYSOs to have a strong near-infrared excess. We compute lower limits to the brightness and luminosity of the MYSO candidates, confirming the near-infrared excess and the massive nature of the objects. No clear correlation is seen between the Br$\gamma$ luminosity and metallicity. Combining our sample with other LMC samples results in a combined detection rate of disk features such as fluorescent Fe II and CO bandheads which is consistent with the Galactic rate (40\%). Most of our MYSOs show outflow features., Comment: Accepted to A&A, 26 pages, 19 figures, 4 tables