The problem of understanding how a newly forming solar system manages to lose angular momentum from the stage of collapsing core to the final main sequence phase is a long-standing one. This is especially pertinent regarding the evolution of the accretion disks. Both theoretical studies and recent observational surveys have questioned the paradigm of viscosity-driven disk evolution, while at the same time alternative scenarios considering the action of magnetically driven winds have gained credit. In recent years, observers have struggled to provide constraints to establish which mechanism is principally responsible, and one of the most interesting findings proved to be the detection and characterisation of rotational motions in the wind components, coupled with the rotation signatures in the associated disk. Such results allowed the estimate of the angular momentum transported by the winds and confirmed the validity of the magneto-hydrodynamic models for the acceleration of outflows and simultaneous disk braking. Importantly, the rotation estimates proved to allow the derivation of the region of origin of the winds in the disk (‘footpoint’), thus enabling a direct evaluation of the feedback of the outflows on the disk in the planet formation region. Here we present an overview on the observational efforts in this direction, from the first detection of jet rotation with the Hubble Space Telescope, to the current state of play, in which ALMA is providing the most stringent constraints to both wind models and related disk evolution. Finally, we mention the future observational directions to further explore this intriguing topic.

Jets and outflows are ubiquitous phenomena associated with the formation of young stellar objects (YSOs). They play a crucial role in removing angular momentum from the accreting system and in regulating star-formation efficiency. Theoretical studies and observations with ALMA and VLA have shown that jets and winds may have a crucial role in promoting dust growth in the envelope-disc system and in shaping the physical and chemical properties of the surrounding environment. Despite these significant advances, many fundamental questions remain unanswered regarding the acceleration, collimation, and chemical impact of jets and outflows from YSOs. The SKA-project will overcome the limitations of current mm/cm-facilities by enabling high-angular resolution and high-sensitivity cm-observations, crucial for probing jets/outflows near YSOs. Radio recombination lines, combined with proper motions, offer a unique opportunity to study the 3D-kinematics of jets. Non-thermal linearly polarised synchrotron emission will allow measuring magnetic field strength and morphology at unprecedented scales of a few au. Observations of dust emission in outflow cavities will allow studying how dust grows and is eventually transported from the disc to the envelope and back. Finally, the SKA-project will allow exploring the dust composition and chemical enrichment in shocks, where sputtering/shattering of grains cause the release of their mantles and refractory cores in the gas-phase. Complementary to ALMA's detection of simple and complex organic molecules, the SKAO will probe, for the first time, long carbon chains/rings, several Cl-, Al-, Mg-, and other metal-bearing species (missed by current sub-mm facilities).

Protoplanetary disks represent a crucial stage in the evolution of Young Stellar Objects towards the formation of fully formed planetary systems. While substantial progress has been made in the last decades in the characterization of the dust and molecular gas in these systems, the ionized component remains poorly understood. Ionized gas traces important processes such as photoevaporation, accretion, disk winds, and jets, and therefore is key to studying disk dynamics, evolution, and ultimately planet formation. In this paper, we investigate the capabilities of the forthcoming SKA telescope to probe this component in protoplanetary disks within nearby star forming regions. We present state-of-the-art simulations of photoevaporative, magneto-thermal, and magnetohydrodynamic winds, and generate theoretical predictions and synthetic SKAO observations to assess its potential in detecting and characterizing free-free emission and Hydrogen recombination lines. Finally, we discuss synergies with complementary facilities and how they will provide a comprehensive, multi-scale view of disk winds and offer critical insights on the mechanisms driving disk evolution and the onset of planet formation.

We analyze the H _2 emission observed in the HH 46 Class I system as part of PROtostellar JEts Cradle Tested with JWST (PROJECT-J), to investigate the origin and excitation of the warm molecular outflow. We used NIRSpec and MIRI spectral maps (1.6–27.9 μ m) to trace the structure and physical conditions of the outflow. By fitting the H _2 rotational diagrams with a multitemperature gas model, we derived key physical parameters, including temperature, extinction, column densities, and the ortho-to-para ratio. This information is combined with a detailed kinematical analysis and comparison with irradiated shock models. We find no evidence of H _2 temperature or velocity stratification from the axis to the edge of the outflow, as would be expected in MHD disk-wind models and as observed in other outflows. Instead, the observations suggest that the H _2 emission arises from shock interactions of jet bow shocks and/or wide-angle winds with the ambient medium and cavity walls. NIRSpec emission and velocity maps reveal expanding molecular shells, likely driven by the less luminous source in the binary system. We infer an accretion rate of ≲10 ^−9 M _⊙ yr ^−1 for the secondary source, approximately one order of magnitude lower than that of the primary. The H _2 emission is consistent with excitation by low-velocity (∼10 km s ^−1 ) J-type shocks, irradiated by an external UV field that may originate from strong dissociative shocks driven by the atomic jet. Future JWST observations will further constrain the evolution of the expanding shell and the mechanisms driving the outflow.

Context. Characterising stellar jet asymmetries is key to setting robust constraints on jet launching models and improving our understanding of the underlying mechanisms behind jet launching. Aims. We aim to characterise the asymmetric properties of the bipolar jet coming from the Classical T Tauri Star Th 28. Methods. We combined data from integral field spectroscopy with VLT/MUSE and high-resolution spectra from VLT/X-shooter to map the optical emission line ratios in both jet lobes. We carried out a diagnostic analysis of these ratios to compare the density, electron temperature, and ionisation fraction within both lobes. The mass accretion rate was derived from the emission lines at the source and compared with the mass outflow rate derived for both lobes, using the estimated densities and measured [O I]lambda 6300 and [S II]lambda 6731 luminosities. Results. The blue-shifted jet exhibits a significantly higher electron temperature and moderately higher ionisation fraction than the red-shifted jet. In contrast to previous studies, we also estimated higher densities, denoted as n(H), in the blue-shifted jet by a factor of similar to 2. These asymmetries are traced to within 1 '' (160 au) of the source in the line ratio maps. We find <(M)over dot>(acc) = 2.4 x 10(-7) M-circle dot yr(-1), with an estimated obscuration factor of similar to 54 due to grey scattering around the star. Estimated values of <(M)over dot>(out) range between 0.66 and 13.7 x 10(-9) M-circle dot yr(-1) in the blue-shifted jet and 5-9 x 10(-9) M-circle dot yr(-1) in the red-shifted jet. Conclusions. The emission line maps and diagnostic results suggest that the jet asymmetries originate close to the source and are likely to be intrinsic to the jet. Furthermore, the combined dataset offers access to a broad array of accretion tracers. In turn, this enables a more accurate estimation of the mass accretion rate, revealing a value of Macc that is higher by a factor >350 than would otherwise be determined.

How have planets formed in the Solar System? And what chemical composition they inherited from their natal environment? Is the chemical composition passed unaltered from the earliest stages of the formation of the Sun to its disk and then to the planets which assembled in the disk? Or does it reflects chemical processes occurring in the disk and/or during the planet formation process? And what was the role of comets in the delivery of volatiles and prebiotic compounds to early Earth? A viable way to answer these questions is to observe protoplanetary disks around young Sun-like stars and compare their chemical composition with that of the early Solar System, which is imprinted in comets. The impacting images recently obtained by millimetre arrays of antennas such as ALMA provided the first observational evidence of ongoing planet formation in 0.1-1 million years old disks, through rings and gaps in their dust and gas distribution. The chemical composition of the forming planets and small bodies clearly depends on the location and timescale for their formation and is intimately connected to the spatial distribution and abundance of the various molecular species in the disk. The chemical characterisation of disks is therefore crucial. This field, however, is still in its infancy, because of the small sizes of disks (~100 au) and to the low gas-phase abundance of molecules (abundances with respect to H2 down to 10-12), which requires an unprecedented combination of angular resolution and sensitivity. I will show the first pioneering results obtained as part of the ALMA chemical survey of protoplanetary disks in the Taurus star forming region (ALMA-DOT program). Thanks to the ALMA images at ~20 au resolution, we recovered the radial distribution and abundance of diatomic molecules (CO and CN), S-bearing molecules (CS, SO, SO2, H2CS), as well as simple organics (H2CO and CH3OH) which are key for the formation of prebiotic compounds. Enhanced H2CO emission in the cold outer disk, outside the CO snowline, suggests that organic molecules may be efficiently formed in disks on the icy mantles of dust grain. This could be the dawn of ice chemistry in the disk, producing ices rich of complex organic molecules (COMs) which could be incorporated by the bodies forming in the outer disk region, such as comets. The next step is the comparison of the molecules radial distribution and abundance in disks with the chemical composition of comets, which are the leftover building blocks of giant planet cores and other planetary bodies. The first pioneering results in this direction have been obtained thanks to the ESA’s Rosetta mission, which allowed obtaining in situ measurements of the COMs abundance on the comet 67P/Churyumov-Gerasimenko. The comparison with three protostellar solar analogs observed on Solar System scales has shown comparable COMs abundance, implying that the volatile composition of comets and planetesimals may be partially inherited from the protostellar stage. The advent of new mission, devoted to sample return such as AMBITION will allow us to do a step ahead in this direction.

We present the first results of the JWST program PROJECT-J (PROtostellar JEts Cradle Tested with JWST ), designed to study the Class I source HH46 IRS and its outflow through NIRSpec and MIRI spectroscopy (1.66 to 28 micron). The data provide line-images (~ 6.6" in length with NIRSpec, and up to 20" with MIRI) revealing unprecedented details within the jet, the molecular outflow and the cavity. We detect, for the first time, the red-shifted jet within ~ 90 au from the source. Dozens of shock-excited forbidden lines are observed, including highly ionized species such as [Ne III] 15.5 micron, suggesting that the gas is excited by high velocity (> 80 km/s) shocks in a relatively high density medium. Images of H2 lines at different excitations outline a complex molecular flow, where a bright cavity, molecular shells, and a jet-driven bow-shock interact with and are shaped by the ambient conditions. Additional NIRCam 2 micron images resolve the HH46 IRS ~ 110 au binary system and suggest that the large asymmetries observed between the jet and the H2 wide angle emission could be due to two separate outflows being driven by the two sources. The spectra of the unresolved binary show deep ice bands and plenty of gaseous lines in absorption, likely originating in a cold envelope or disk. In conclusion, JWST has unraveled for the first time the origin of the HH46 IRS complex outflow demonstrating its capability to investigate embedded regions around young stars, which remain elusive even at near-IR wavelengths.

Planet-forming disks are not isolated systems. Their interaction with the surrounding medium affects their mass budget and chemical content. In the context of the ALMA-DOT program, we obtained high-resolution maps of assorted lines from six disks that are still partly embedded in their natal envelope. In this work, we examine the SO and SO$_2$ emission that is detected from four sources: DG Tau, HL Tau, IRAS 04302+2247, and T Tau. The comparison with CO, HCO$^+$, and CS maps reveals that the SO and SO$_2$ emission originates at the intersection between extended streamers and the planet-forming disk. Two targets, DG Tau and HL Tau, offers clear cases of inflowing material inducing an accretion shock on the disk material. The measured rotational temperatures and radial velocities are consistent with this view. In contrast to younger Class 0 sources, these shocks are confined to the specific disk region impacted by the streamer. In HL Tau, the known accreting streamer induces a shock in the disk outskirt, and the released SO and SO$_2$ molecules spiral toward the star in a few hundreds years. These results suggest that shocks induced by late accreting material may be common in the disks of young star-forming regions with possible consequences on the chemical composition and mass content of the disk. They also highlight the importance of SO and SO$_2$ line observations to probe accretion shocks from a larger sample.

SHARK-NIR is an instrument which provides direct imaging, coronagraphic imaging, dual band imaging and low resolution spectroscopy in Y, J and H bands, taking advantage of the outstanding performance of the Large Binocular Telescope AO systems. Binocular observations will be provided used in combination with SHARK-VIS (operating in V band) and LMIRCam of LBTI (operating from K to M bands), in a way to exploit coronagraphic simultaneous observations in three different wavelengths. A wide variety of coronagraphic techniques have been implemented in SHARK-NIR, ranging from conventional ones such as the Gaussian Lyot, to others quite robust to misalignments such as the Shaped Pupil, to eventually techniques more demanding in term of stability during the observation, as the Four Quadrant; the latter is giving in theory and simulations outstanding contrast, and it is supported in term of stability by the SHARK-NIR internal fast tip-tilt loop and local NCPA correction, which should ensure the necessary stability allowing this technique to operate at its best. The main science case is of course exoplanets search and characterization and young stellar systems, jets and disks characterization, although the LBT AO extreme performance, allowing to reach excellent correction even at very faint magnitudes, may open to science previously difficult to be achieved, as for example AGN and QSO morphological studies. The institutes participating to the SHARK-NIR consortium which designed and built the instrument are Istituto Nazionale di Astro Fisica (INAF, Italy), the Max Planck Institute for Astronomy (MPIA, Heidelberg, Germany) and University of Arizona/Steward Observatory (UoA/SO, Tucson, Az, USA). We report here about the SHARK-NIR status, that should achieve first light at LBT before the end of 2022.

The evolution of young stars and disks is driven by the interplay of several processes, notably the accretion and ejection of material. These processes, critical to correctly describe the conditions of planet formation, are best probed spectroscopically. Between 2020 and 2022, about 500orbits of the Hubble Space Telescope (HST) are being devoted in to the ULLYSES public survey of about 70 low-mass (M⋆ ≤ 2 M⊙) young (age < 10 Myr) stars at UV wavelengths. Here, we present the PENELLOPE Large Program carried out with the ESO Very Large Telescope (VLT) with the aim of acquiring, contemporaneously to the HST, optical ESPRESSO/UVES high-resolution spectra for the purpose of investigating the kinematics of the emitting gas, along with 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 being made available to the whole community. Here, we describe the data and the first scientific analysis of the accretion properties for the sample of 13 targets located in the Orion OB1 association and in the σ-Orionis cluster, observed in November–December 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 (≲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%, which is well within the assumed uncertainty. Its origin can be either due to an erroneous UV extinction curve or to the simplicity of the modeling and, thus, this question will form the basis of the investigation undertaken over the course of the PENELLOPE program. The combined ULLYSES and PENELLOPE data will be key in attaining a better understanding of the accretion and ejection mechanisms in young stars.

Context. The role of bipolar jets in the formation of stars, and in particular how they are launched, is still not well understood. Aims. We probe the protostellar jet launching mechanism using high-resolution observations of the near-infrared (IR) [Fe II ] λ 1.53,1.64 μ m emission lines. Methods. We consider the case of the bipolar jet from Classical T Tauri star, DO Tau, and investigate the jet morphology and kinematics close to the star (within 140 au) using AO-assisted IFU observations from GEMINI/NIFS. Results. We find that the brighter, blueshifted jet is collimated very quickly after it is launched. This early collimation requires the presence of magnetic fields. We confirm velocity asymmetries between the two lobes of the bipolar jet, and also confirm no time variability in the asymmetry over a 20-year interval. This sustained asymmetry is in accordance with recent simulations of magnetised disc winds. We examine the data for signatures of jet rotation. We report an upper limit on differences in radial velocity of 6.3 and 8.7 km s −1 for the blue- and redshifted jets, respectively. Interpreting this as an upper limit on jet rotation implies that any steady, axisymmetric magneto-centrifugal model of jet launching is constrained to a launch radius in the disc plane of r 0 < 0.5 and 0.3 au for the blue- and redshifted jets, respectively. This supports an X-wind or narrow disc-wind model. However, the result pertains only to the observed high-velocity [Fe II ] emission, and does not rule out a wider flow launched from a wider radius. We report the detection of small-amplitude jet axis wiggling in both lobes. We rule out orbital motion of the jet source as the cause. Precession can better account for the observations but requires double the precession angle, and a different phase for the counter-jet. Such non-solid body precession could arise from an inclined massive Jupiter companion, or a warping instability induced by launching a magnetic disc wind. Conclusions. Overall, our observations are consistent with an origin of the DO Tau jets from the inner regions of the disc.

The rotating disk of gas and dust observed around young stars is the birthplace of planets. The study of its properties is fundamental to understand the formation and evolution of nascent worlds. Polarimetric data obtained with the Atacama Large Millimeter/Submillimeter Array (ALMA) are a new powerful tool to investigate the properties of the dust and magnetic structure in these systems.

Context. Th 28 is a Classical T Tauri star in the Lupus 3 cloud that drives an extended bipolar jet. Previous studies of the inner jet identified signatures of rotation around the outflow axis, a key result for theories of jet launching. Thus this is an important source in which to investigate the poorly understood jet launching mechanism. Aims. In this study we investigate the morphology and kinematics of the Th 28 micro-jets, with the aim of characterising their structure and outflow activity, using optical integral-field spectroscopy observations obtained with VLT/MUSE. Methods. We use spectro-imaging and position–velocity maps to investigate the kinematic and morphological features of the jet and to obtain a catalogue of emission lines in which the jet is visible. A Lucy-Richardson deconvolution procedure is used to differentiate the structure of the inner micro-jet region in selected emission lines. Spatial profiles extracted perpendicular to the jet axis are fitted to investigate the jet width, opening angle, and the evolution of the jet axis. Results. We confirm the previously identified knot HHW2 within the red-shifted jet and identify three additional knots in each lobe for the first time. We also find [O III]λ5007 emission from the blue-shifted micro-jet, including the knot closest to the star. Proper motions for the innermost knots on each side are estimated to be 0′′.35 yr−1 and 0′′.47 yr−1 for the red- and blue-shifted jets, respectively. Based on this we show that new knots are ejected on an approximate timescale of 10–15 yr. Gaussian fitting to the jet axis centroids shows a point-symmetric wiggle within the inner portion of both micro-jets, indicating precession of the jet. We use the jet shape to measure a precession period of 8 yr, with a half-opening angle β < 0.6°. This precession may provide an alternative explanation for the rotation signatures previously reported. Conclusions. We find that these parameters are compatible with precession due to a brown dwarf companion orbiting at a separation of ≤0.3 au. Further observations with higher spatial resolution may help to clarify the source of this precession.

We present near-infrared [Fe ii] images of four Class 0/I jets (HH 1/2, HH 34, HH 111, HH 46/47) observed with the Hubble Space Telescope (HST) Wide Field Camera 3. The unprecedented angular resolution allows us to measure proper motions, jet widths and trajectories, and extinction along the jets. In all cases, we detect the counterjet, which was barely visible or invisible at shorter wavelengths. We measure tangential velocities of a few hundred kilometers per second, consistent with previous HST measurements over 10 years ago. We measure the jet width as close as a few tens of au from the star, revealing high collimations of about 2° for HH 1, HH 34, and HH 111 and about 8° for HH 46, all of which are preserved up to large distances. For HH 34, we find evidence of a larger initial opening angle of about 7°. Measurement of knot positions reveals deviations in trajectory of both the jet and counterjet of all sources. Analysis of asymmetries in the inner knot positions for HH 111 suggests the presence of a low mass stellar companion at separation 20–30 au. Finally, we find extinction values of 15–20 mag near the source, which gradually decrease moving downstream along the jet. These observations have allowed us to study the counterjet at unprecedentedly high angular resolution and will be a valuable reference for planning future JWST mid-infrared observations that will peer even closer into the jet engine.

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

Context. Star formation (SF) is a multi-scale process in which the mode of fragmentation of the collapsing clump on scales of 0.1–1 pc determines the mass reservoir and affects the accretion process of the individual protostars on scales of 10–100 au. Aims. We want to investigate the nearby (located at 1.63 ± 0.05 kpc) high-mass star-forming region IRAS 21078+5211 at linear scales from ~1 pc down to ~10 au. Methods. We combine the data of two recent programs: the NOrthern Extended Millimeter Array large project CORE and the Protostellar Outflows at the EarliesT Stages (POETS) survey. The former provides images of the 1 mm dust continuum and molecular line emissions with a linear resolution of ≈600 au covering a field of view up to ≈0.5 pc. The latter targets the ionized gas and 22 GHz water masers, mapping linear scales from a few 103 au down to a few astronomical units. Results. In IRAS 21078+5211, a highly fragmented cluster (size ~0.1 pc) of molecular cores is observed, located at the density peak of an elongated (size ~1 pc) molecular cloud. A small (≈1 km s−1 per 0.1 pc) LSR velocity (VLSR) gradient is detected across the major axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of ≈10−4 M⊙ yr−1. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. The masses of the young stellar objects (YSOs) inside these three cores are estimated in the range 1–6 M⊙. We reveal an extended (size ~0.1 pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by the compact (size ≲1000 au) radio jet discovered in the POETS survey, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (≈14 km s−1 over 500 au) VLSR gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6 ± 2.0 M⊙. The water masers observed in the POETS survey emerge within 100–300 au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be ≈ 16 au and the corresponding launching radius and terminal velocity to be ≈ 2.2 au and ≈ 200 km s−1, respectively. Conclusions. Complementing high-angular resolution, centimeter and millimeter interferometric observations in thermal tracers with Very Long Baseline Interferometry of molecular masers, is invaluable in studying high-mass SF. The combination of these twodatasets allows us to connect the events that we see at large scales, as clump fragmentation and mass flows, with the physical processes identified at small scales, specifically, accretion and ejection in disk-jet systems.

We present an overview of the ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT), a campaign devoted to the characterization of the molecular emission from partly embedded young stars. The project is aimed at attaining a better understanding of the gaseous products delivered to planets by means of high-resolution maps of the assorted lines probing disks at the time of planet formation (≲1 Myr). Nine different molecules are surveyed through our observations of six Class I/flat-spectrum sources. As part of a series of articles analyzing specific targets and molecules, in this work we describe the sample and provide a general overview of the results, focusing specifically on the spatial distribution, column densities, and abundance ratios of H2CO, CS, and CN. In these embedded sources, the 12CO emission is dominated by envelope and outflow emission while the CS and, especially, the H2CO are good tracers of the gaseous disk structure. The spatial distribution and brightness of the o-H2CO 31,2−21,1 and CS 5−4 lines are very similar to each other and across all targets. The CN 2−1 line emission is fainter and distributed over radii larger than the dust continuum. The H2CO and CS emission is always dimmed in the inner ~50 au. While the suppression by the dusty disk and absorption by the line-of-sight material significantly contributes to this inner depression, an actual decrease in the column density is plausible in most cases, making the observed ring-like morphology realistic. We also found that the gaseous disk extent, when traced by H2CO (150−390 au), is always 60% larger than the dust disk. This systematic discrepancy may, in principle, be explained by the different optical depth of continuum and line emission without invoking any dust radial drift. Finally, the o-H2CS 71,6−61,5 and CH3OH 50,5−40,4 line emission are detected in two disks and one disk, respectively, while the HDO is never detected. The H2CO column densities are 12−50 times larger than those inferred for Class II sources while they are in line with those of other Class 0/I. The CS column densities are lower than those of H2CO, which is an opposite trend with regard to Class II objects. We also inferred abundance ratios between the various molecular species finding, among others, a H2CS/H2CO ratio that is systematically lower than unity (0.4−0.7 in HL Tau, 0.1 − 0.2 in IRAS 04302+2247, and <0.4 in all other sources), as well as a CH3OH/H2CO ratio (<0.7 in HL Tau and 0.5−0.7 in IRAS 04302+2247) that is lower than the only available estimate in a protoplanetary disks (1.3 in TW Hya) and between one and two orders of magnitude lower than those of the hot corinos around Class 0 protostars. These results are a first step toward the characterization of the disk’s chemical evolution, which ought to be complemented by subsequent observations of less exceptional disks and customized thermo-chemical modeling.