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

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

(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

HM Lup is a young M-type star that accretes material from a circumstellar disk through a magnetosphere. Our aim is to study the inner disk structure of HM Lup and to characterize its variability. We used spectroscopic data from HST/STIS, X-Shooter, and ESPRESSO taken in the framework of the ULLYSES and PENELLOPE programs, together with photometric data from TESS and AAVSO. The 2021 TESS light curve shows variability typical for young stellar objects of the "accretion burster" type. The spectra cover the temporal evolution of the main burst in the 2021 TESS light curve. We compared the strength and morphology of emission lines from different species and ionization stages. We determined the mass accretion rate from selected emission lines and from the UV continuum excess emission at different epochs, and we examined its relation to the photometric light curves. The emission lines in the optical spectrum of HM Lup delineate a temperature stratification along the accretion flow. While the wings of the H I and He I lines originate near the star, the lines of species such as Na I, Mg I, Ca I, Ca II, Fe I, and Fe II are formed in an outer and colder region. The shape and periodicity of the 2019 and 2021 TESS light curves, when qualitatively compared to predictions from magnetohydrodynamic models, suggest that HM Lup was in a regime of unstable ordered accretion during the 2021 TESS observation due to an increase in the accretion rate. Although HM Lup is not an extreme accretor, it shows enhanced emission in the metallic species during this high accretion state that is produced by a density enhancement in the outer part of the accretion flow., Comment: 15 pages, 14 figures. Accepted for publication in A&A

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

Early dust grain growth in protostellar envelopes infalling on young discs has been suggested in recent studies, supporting the hypothesis that dust particles start to agglomerate already during the Class 0/I phase of young stellar objects (YSOs). If this early evolution were confirmed, it would impact the usually assumed initial conditions of planet formation, where only particles with sizes $\lesssim 0.25 \mu$m are usually considered for protostellar envelopes. We aim to determine the maximum grain size of the dust population in the envelope of the Class 0/I protostar L1527 IRS, located in the Taurus star-forming region (140 pc). We use Atacama Large millimetre/sub-millimetre Array (ALMA) and Atacama Compact Array (ACA) archival data and present new observations, in an effort to both enhance the signal-to-noise ratio of the faint extended continuum emission and properly account for the compact emission from the inner disc. Using observations performed in four wavelength bands and extending the spatial range of previous studies, we aim to place tight constraints on the spectral ($\alpha$) and dust emissivity ($\beta$) indices in the envelope of L1527 IRS. We find a rather flat $\alpha \sim$ 3.0 profile in the range 50-2000 au. Accounting for the envelope temperature profile, we derive values for the dust emissivity index, 0.9 < $\beta$ < 1.6, and reveal a tentative, positive outward gradient. This could be interpreted as a distribution of mainly ISM-like grains at 2000 au, gradually progressing to (sub-)millimetre-sized dust grains in the inner envelope, where at R=300 au, $\beta$ = 1.1 +/- 0.1. Our study supports a variation of the dust properties in the envelope of L1527 IRS. We discuss how this can be the result of in-situ grain growth, dust differential collapse from the parent core, or upward transport of disc large grains., Comment: Accepted for publication in A&A. Contains 18 pages, 21 figures, 5 tables Replacement on Nov 22 to change title number of FAUST series from "X" to "IX."

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

Elias 2-27 is a young star that hosts an extended, bright and inclined disk of dust and gas. The inclination and extreme flaring of the disk make Elias 2-27 an ideal target to study the vertical distribution of molecules, particularly CN. We directly trace the emission of CN in Elias 2-27 and compare it to previously published CO isotopologue data. CN $N = 3-2$ emission is analyzed in two different transitions $J = 7/2 - 5/2$ and $J = 5/2 - 3/2$, for which we detect two hyperfine group transitions. The vertical location of CN emission is traced directly from the channel maps, following geometrical methods that have been previously used to analyze the CO emission of Elias 2-27. Analytical models are used to parametrize the vertical profile of each molecule and study the extent of each tracer, additionally we compute radial profiles of column density and optical depth. We show that the vertical location of CN and CO isotopologues in Elias 2-27 is layered and consistent with predictions from thermochemical models. A north/south asymmetry in the radial extent of CN is detected and we find that the CN emission is mostly optically thin and constrained vertically to a thin slab at $z/r \sim$0.5. A column density of 10$^{14}$\,cm$^{-2}$ is measured in the inner disk which for the north side decreases to 10$^{12}$\,cm$^{-2}$ and for the south side to 10$^{13}$\,cm$^{-2}$ in the outer regions. In Elias 2-27, CN traces a vertically elevated region above the midplane, very similar to that traced by $^{12}$CO. The inferred CN properties are consistent with thermo-chemical disk models, in which CN formation is initiated by the reaction of N with UV-pumped H$_2$. The observed north/south asymmetry may be caused by either ongoing infall or by a warped inner disk. This study highlights the importance of tracing the vertical location of various molecules to constrain the disk physical conditions., Comment: Accepted for publication in A&A, 19 pages, 14 figures

(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

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

The evolution of young stars and disks is driven by the interplay of several processes, notably accretion and ejection of material. Critical to correctly describe the conditions of planet formation, these processes are best probed spectroscopically. About five-hundred orbits of the Hubble Space Telescope (HST) are being devoted in 2020-2022 to the ULLYSES public survey of about 70 low-mass (M<2Msun) young (age<10 Myr) stars at UV wavelengths. Here we present the PENELLOPE Large Program that is being carried out at the ESO Very Large Telescope (VLT) to acquire, contemporaneous to HST, optical ESPRESSO/UVES high-resolution spectra to investigate the kinematics of the emitting gas, and UV-to-NIR X-Shooter medium-resolution flux-calibrated spectra to provide the fundamental parameters that HST data alone cannot provide, such as extinction and stellar properties. The data obtained by PENELLOPE have no proprietary time, and the fully reduced spectra are made available to the whole community. Here, we describe the data and the first scientific analysis of the accretion properties for the sample of thirteen targets located in the Orion OB1 association and in the sigma-Orionis cluster, observed in Nov-Dec 2020. We find that the accretion rates are in line with those observed previously in similarly young star-forming regions, with a variability on a timescale of days of <3. The comparison of the fits to the continuum excess emission obtained with a slab model on the X-Shooter spectra and the HST/STIS spectra shows a shortcoming in the X-Shooter estimates of <10%, well within the assumed uncertainty. Its origin can be either a wrong UV extinction curve or due to the simplicity of this modelling, and will be investigated in the course of the PENELLOPE program. The combined ULLYSES and PENELLOPE data will be key for a better understanding of the accretion/ejection mechanisms in young stars., Comment: Accepted for publication on Astronomy & Astrophysics. 15 pages + appendix, language edited version

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

Intermediate-mass young stellar objects (YSOs) provide a link to understand how feedback from shocks and UV radiation scales from low to high-mass star forming regions. Aims: Our aim is to analyze excitation of CO and H$_2$O in deeply-embedded intermediate-mass YSOs and compare with low-mass and high-mass YSOs. Methods: Herschel/PACS spectral maps are analyzed for 6 YSOs with bolometric luminosities of $L_\mathrm{bol}\sim10^2 - 10^3$ $L_\odot$. The maps cover spatial scales of $\sim 10^4$ AU in several CO and H$_2$O lines located in the $\sim55-210$ $\mu$m range. Results: Rotational diagrams of CO show two temperature components at $T_\mathrm{rot}\sim320$ K and $T_\mathrm{rot}\sim700-800$ K, comparable to low- and high-mass protostars probed at similar spatial scales. The diagrams for H$_2$O show a single component at $T_\mathrm{rot}\sim130$ K, as seen in low-mass protostars, and about $100$ K lower than in high-mass protostars. Since the uncertainties in $T_\mathrm{rot}$ are of the same order as the difference between the intermediate and high-mass protostars, we cannot conclude whether the change in rotational temperature occurs at a specific luminosity, or whether the change is more gradual from low- to high-mass YSOs. Conclusions: Molecular excitation in intermediate-mass protostars is comparable to the central $10^{3}$ AU of low-mass protostars and consistent within the uncertainties with the high-mass protostars probed at $3\cdot10^{3}$ AU scales, suggesting similar shock conditions in all those sources., Comment: Accepted to Astronomy & Astrophysics. 4 pages, 5 figures, 3 tables