GY 91, commonly categorized as a Class I young stellar object, is notable for disk dust substructures that have been hypothesized to trace early planet formation. Using the Atacama Large Millimeter/submillimeter Array 12 m and Atacama Compact Array, we present new Band 7 dust continuum and molecular line observations of GY 91 at an angular resolution of similar to 0.'' 3 (40 au). We report detections of CS J = 6-5, N2H+ J = 3-2, C18O J = 3-2, H2CS JKa,Kc=81,7-71,6 , H2CO JKa,Kc=40,4-30,3 , and H2CO JKa,Kc=42,3-32,2 , as well as a tentative detection of 13C18O J = 3-2. We observe azimuthal asymmetry in CS and H2CS emission, as well as radially structured H2CO 40,4-30,3 emission outside the dust continuum. C18O and H2CO 40,4-30,3 show significant cloud contamination, while CS and N2H+ are good tracers of Keplerian rotation originating from the disk. Envelope emission does not appear to contribute significantly either to the continuum or molecular line observations. GY 91's chemical properties appear in large part to resemble those of Class II disks, although observations of additional molecular probes should be obtained for a fuller comparison. With CS, we estimated a dynamical stellar mass of 0.58 M circle dot, which is higher than previous estimates from stellar evolutionary models (0.25 M circle dot). Using both radiative transfer modeling of the dust continuum and comparison of the C18O and N2H+ fluxes to literature thermochemical models, we estimate a disk mass of similar to 0.01 M circle dot.

[Abridged] We aim to study the radial and vertical extents of 12CO gas, millimeter dust thermal emission and optical/NIR scattered light by dust in disks. We analyze a sample of 14 highly inclined protoplanetary disks. We present ALMA high angular resolution band 7 (0.9 mm) continuum images and 12CO (3-2) gas moment maps as well as HST and VLT/SPHERE scattered light images. The majority of disks in our sample (11 out of 14) follow Rgas > Rdust,micron > Rdust,mm. The other 3 disks appear more extended in millimeter continuum than in scattered light. Highly inclined disks tend to appear less radially extended in CO gas line emission than in millimeter dust continuum compared to less inclined disks. This results from optical depth effects and/or radial drift. The known correlation between disk size and millimeter continuum and line fluxes are confirmed in our sample with highly inclined disks significantly fainter than disks seen at lower inclination for a given disk radius. We found that this correlation is significantly tightened once fluxes are corrected for the disk inclination, consistent with the disks being optically thick at millimeter wavelengths. Regarding the vertical extent defined as the apparent emitting height, most disks in our sample follow Hgas > Hdust, mm. This strengthens our previous findings that the millimeter dust is highly decoupled from the gas and forms a layer in the disk midplane due to vertical settling. Most disks appear more vertically extended in gas than in scattered light, suggesting that the micron-sized dust is not fully coupled to the gas. We also estimated dynamical masses using PV diagrams for the first time for most of the objects in our sample. We found an anti-correlation between the dynamical mass and the aspect ratio, emphasizing the dominant role of gravity in setting the disk vertical extent, but no correlation with the disk radius.

Turbulence is expected to transport angular momentum and drive mass accretion in protoplanetary disks. One way to directly measure turbulent motion in disks is through molecular line broadening. DM Tau is one of only a few disks with claimed detection of nonthermal line broadening of 0.25cs-0.33cs, where cs is the sound speed. Using the radiative transfer code mcfost within a Bayesian inference framework that evaluates over five million disk models to efficiently sample the parameter space, we fit high-resolution ( 0.'' 15 , 28 m s-1) 12CO J = 3-2 observations of DM Tau from the exoALMA Large Program. This approach enables us to simultaneously constrain the disk structure and kinematics, revealing a significant nonthermal contribution to the line width of similar to 0.4cs, inconsistent with purely thermal motions. Using the CO-based disk structure as a starting point, we reproduce the CS J = 7-6 emission well, demonstrating that the CS (which is more sensitive to nonthermal motions than CO) agrees with the turbulence inferred from the CO fit. Establishing a well-constrained background disk model further allows us to identify residual structures in the moment maps that deviate from the expected emission, revealing localized perturbations that may trace forming planets. This framework provides a powerful general approach for extracting disk structure and nonthermal broadening directly from molecular line data and can be applied to other disks with high-quality observations.

Context. The evolution of protoplanetary disks, especially in the early stages of planetary formation, as dust grows, is the cornerstone of the birth of planets. The mechanisms involved in the growth of sub-micrometric dust grains into planetesimals within a very short time frame are a challenging field of study, while the initial conditions remain relatively undefined. Aims. One of the main challenges is to unambiguously identify the dust properties within the disk, and our goal is to break this barrier by investigating the light scattered by dust particles lying on the protoplanetary disk surface from many recent promising observations. Methods. In this study, we used a set of 30 polarized light images composed of new VLT/SPHERE observations to examine the light scattered by dust grains. For each ring-shaped system, we used the new DRAGyS tool to estimate the disk geometry using the substructures visible on the surface and to extract the limb-brightening-corrected scattering phase function, which encodes the dust grains’ physical properties. Finally, we compared our results with the AggScatVIR database of numerical scattering phase functions of nonspherical dust. Results. We combined our measurement of disk geometry to estimate an average disk flaring of about 1.357. We note some general trends of dust populations in our results. First, we recovered the two categories of scattering phase functions based on their shape, as determined in previous studies. Category I is monotonically decreasing and can be explained by fractal organic aggregates with small monomers of 100nm, or compact aggregates with medium porosity and big monomers of 400nm. Category II is defined by a bellshaped scattering phase function and can be explained by sub-micrometric irregular grains or compact aggregates with low porosity. This statistical study offers general trends about dust populations, but the degeneracy is too strong to apply this method to a unique disk analysis. Conclusions. The extracted scattering phase function from the protoplanetary disk surface is a promising measurement for constraining the properties of dust within the disk. However, their interpretation should be combined with other observations. Multiwavelength investigation, i.e., combining SPFs with disk color measurements, is a promising option.

In the past few years, ALMA unveiled a variety of substructures (rings, spirals, crescents) in the continuum emission of most protoplanetary disks imaged at high spatial resolution. While the majority of disks presents axisymmetric ring-like structures in the dust brightness distribution, some sources display asymmetric morphologies (blobs, crescents) that have been often associated to vortices and/or mechanisms generated by the presence of one or more embedded planets. In this brief research report we present the analysis of the arc structure observed in the dust continuum emission of the disk around HD 163296, using high resolution ( similar to 8 au) matched continuum data from ALMA at four wavelengths. We characterize in detail the arc structures and present a kinematic signature observed in the CS(3-2) emission at the same location. Our results indicate that the crescent is caused by differential dust trapping in a local pressure maxima, for which plausible mechanisms can be the presence of a vortex or trapping in a Lagrangian point of the planet-star system.

The exoALMA large program provided an unprecedented view of the morphologies and kinematics of 15 circumstellar disks, offering a biased but homogenous and well-characterized sample for population-level analysis. Continuum observations revealed numerous dust substructures, known to be potential signatures of embedded planets. We analyze the observed dust morphologies with the simulation-based inference tool DBNets2.0, assuming these are due to embedded planets at fixed locations, to infer the system's properties. We estimate the putative planet mass, the disk alpha viscosity, scale height, and dust Stokes number that would reproduce 19 substructures in 13 of the 15 exoALMA disks. We compare our results with literature estimates derived with different methods, and find good agreement in most cases. We further explore the implications of the inferred disk properties for accretion, showing that for the Herbig stars in our sample, the implied viscous accretion timescales are too long to account for their observed stellar accretion rates. Regarding planet migration, our results favor inward migration, with only three putative planets expected to migrate outward. Finally, we check for correlations of the inferred disk and planet properties with the disks' gas-to-dust mass ratio, nonaxisymmetry index, and masses of the gas, dust, and host stars, finding no remarkable trends.

The sulfur chemistry in protoplanetary disks influences the properties of nascent planets, including potential habitability. Although the inventory of sulfur molecules in disks has gradually increased over the last decade, CS is still the most commonly observed sulfur-bearing species, and it is expected to be the dominant gas-phase sulfur carrier beyond the water snow line. Despite this, few dedicated multiline observations exist, and thus, the typical disk CS chemistry is not well constrained. Moreover, it is unclear how that chemistry—and in turn, the bulk volatile sulfur reservoir—varies with stellar and disk properties. Here, we present the largest survey of CS to date, combining both new and archival observations from the Atacama Large Millimeter/submillimeter Array, Submillimeter Array, and Northern Extended Millimeter Array of 12 planet-forming disks, covering a range of stellar spectral types and dust morphologies. Using these data, we derived disk-integrated CS gas excitation conditions in each source. Overall, CS chemistry appears similar across our sample with rotational temperatures of ≈10–40 K and column densities between 10 ^12 and 10 ^13 cm ^−2 . CS column densities do not show strong trends with most source properties, which broadly suggests that CS chemistry is not highly sensitive to disk structure or stellar characteristics. We do, however, identify a positive correlation between stellar X-ray luminosity and CS column density, which indicates that the dominant CS formation pathway is likely via ion-neutral reactions in the upper disk layers, where X-ray-enhanced S ^+ and C ^+ drive abundant CS production. Thus, using CS as a tracer of gas-phase sulfur abundance requires a nuanced approach that accounts for its emitting region and dependence on X-ray luminosity.

Vertical gas flows-such as winds and meridional circulations-are natural outcomes of protoplanetary disk processes and play a critical role in the earliest stages of planet formation. We analyze the vertical gas motions in 14 disks, as part of the exoALMA Large Program, focusing on the 12CO J = 3-2 and 13CO J = 3-2 emission lines. Using discminer to model the Keplerian velocity field, we extract line-of-sight velocity residuals and measure the radial and vertical components of the gas motion. Vertical motions are detected in most disks. Two types of patterns emerge: (1) oscillatory up/down flows, likely linked to instabilities; and (2) transitions from downward to upward motions that we interpret as the bases of disk winds. In most cases, the velocity amplitudes are of a few tens of m s-1. However, two disks, MWC758 and CQ Tau, show two spiral velocity features in their residual maps, red- and blueshifted, which we interpret as vertical velocities reaching up to similar to 350 m s-1 (similar to 0.7cs), consistent with the gas motion in eccentric disks. Fast upward motions (up to 500 m s-1; similar to 1.8cs) are also detected in the outer disk of MWC758. Synthetic observations from (magneto)hydrodynamic simulations validate the reliability of our method. Although strong molecular winds appear to be relatively rare in 12CO and 13CO, our study shows that, when traced by deep high-spectral-resolution line data, protoplanetary disks exhibit ubiquitous vertical flows. However, their overall velocity structure is highly complex, preventing the identification of a coherent dominant physical mechanism driving the vertical motions across all disks, thus requiring further theoretical investigation.

The exoALMA Large Program has revealed a wealth of substructures in the dust and molecular line emission of several protoplanetary discs, suggesting that planet formation may unfold within highly dynamic environments. Using synthetic observations of planet-disk interactions and disk instabilities, we demonstrate how the origin of these substructures can be investigated through a tomographic study of molecular lines, extending the scope of the analysis beyond line-centroid kinematics alone. Our results indicate that with only a few hours of ALMA integration at moderate angular resolution ( 0.'' 15-0.'' 30 ), it is possible to identify the key signatures driven by planets more massive than 0.1% of the stellar mass. These signatures manifest not only as deviations from Keplerian motion but also as localized line broadening, enabling accurate constraints on the orbital radius and azimuthal location of the planets. We further show that a diagnostic based on line skewness in spectrally resolved observations can help distinguish between planetary and instability-driven signatures, owing to the distinct degrees of velocity coherence associated with each mechanism. Finally, we apply this tomographic analysis to exoALMA CO line data for the disks of HD 135344B and MWC 758. In HD 135344B, we identify strongly localized velocity and line-width perturbations, suggesting the possibility of three massive planets embedded in the disk: one at R = 95 au, exterior to the continuum substructures, and two within dust gaps at R = 41 au and R = 73 au. For MWC 758, the dominance of vertical-velocity spirals over localized signatures is consistent with predictions from models of moderate disk eccentricities or warps, potentially induced by a substellar companion in the inner regions of the system.

Aims. Edge-on discs offer a unique opportunity to probe radial and vertical dust and gas distributions in the protoplanetary phase. This study aims to investigate the distribution of micron-sized dust particles in the Flying Saucer in Rho Ophiuchi by leveraging the unique observational conditions of a bright infrared background that enables the edge-on disc to be seen in both silhouette and scattered light at specific wavelengths. Methods. We used NIRSpec IFU observations from the JWST Edge-on Disc Ice program (JEDIce) of the Flying Saucer serendipitously observed against a Polycyclic Aromatic Hydrocarbons-emitting background to constrain the dust distribution and grain sizes via radiative transfer modelling. Results. The observation of the Flying Saucer in silhouette at 3.29 μm reveals that the midplane radial extent of small dust grains is ∼235 au, i.e. larger than the large-grain disc extent previously determined to be 190 au from millimetre data. The scattered light observed in emission probes micron-sized icy grains at large vertical distances above the midplane. The vertical extent of the disc silhouette is similar at visible, near-IR, and mid-IR wavelengths, corroborating the conclusion that dust settling is inefficient for grains as large as tens of microns, both vertically and radially.

Protoplanetary disks are the birthplaces of planetary systems, and deviations from Keplerian rotation imprinted in disk gas kinematics serve as key tracers of physical processes and the presence of protoplanets within disks. Using the CO (J = 3-2) data from the exoALMA Large Program, encompassing 15 disks, we constructed two-dimensional (2D) maps of centroid velocity, line width, and peak intensity, and extracted non-Keplerian deviations by subtracting smooth Keplerian models. This Letter provides the first systematic and uniform overview of 2D gas substructures across the entire exoALMA sample. We find that all targets exhibit large-scale deviations from smooth Keplerian disks, displaying a variety of morphologies, including spiral-like structures, arc- or ring-like features, and patterns indicative of variations in the emitting surface height. Nonaxisymmetric spiral-arm features are detected or suggested in five disks (CQ Tau, MWC 758, HD 135344B, HD 34282, and SY Cha), and are preferentially found in Herbig Ae/Fe systems. In contrast, some other sources (J1852, PDS 66, and V4046 Sgr), despite exhibiting noticeable deviations, appear to be dynamically quieter. This 2D atlas suggests that kinematic substructures are ubiquitous in large (greater than or similar to 100 au) protoplanetary disks with ages of a few million years, based on the observations obtained with sufficient sensitivity at a moderate-to-high spatial resolution of similar to 20 au and high-velocity resolution of similar to 0.1 km s-1.

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

Context. A considerable proportion of young stars belong to multiple star systems. Constraining the planet formation processes in multiple stellar systems is then key to understanding the global exoplanet population. Aims. This study focuses on investigating the dust reservoir within the triple system V892 Tau. Our objective is to establish constraints on the properties and characteristics of the dust present in the system's circumbinary ring. Methods. Based on archival ALMA and VLA data from 0.9 mm to 9.8 mm, we present a multi-wavelength analysis of the ring of V892 Tau. We first studied the spatial variation of the spectral index, before employing 3D full-radiative-transfer calculations to constrain the ring's geometry and the radial dependence of the dust-grain properties. Results. Spectral indices are consistent with non-dust emission in the vicinity of the central binary, and with dust emission in the ring likely remaining optically thick up to 3.0 mm. Our radiative transfer analysis supports these interpretations, yielding a model that reproduces the observed intensities within the 1 sigma uncertainties across all wavelengths. The resulting dust-surface density and temperature profiles both decrease with increasing radius, and are in agreement with values reported in the literature. Maximum grain sizes are constrained to 0.2 cm, with a size distribution power-law index -3.5. These results imply that the dust-grain fragmentation velocity does not exceed 8 m s(-1). Conclusions. Whilst our results suggest dust trapping at the cavity edge, they also suggest that tidal perturbations triggered by the central binary limit grain growth within the ring. This highlights the need to further constrain planet formation efficiency in multiple stellar systems, a goal that may be advanced by applying the methodology of this work to a wider sample of systems.

The bulk motion of the gas in protoplanetary disks around newborn stars is nearly Keplerian. By leveraging the high angular and spectral resolution of the Atacama Large Millimeter/submillimeter Array (ALMA), we can detect small-scale velocity perturbations in molecular line observations caused by local gas pressure variations in the disk, possibly induced by embedded protoplanets. This Letter presents the azimuthally averaged rotational velocity and its deviations from Keplerian rotation ( δυ _ϕ ) for the exoALMA sample, as measured in the ^12 CO J = 3–2 and ^13 CO J = 3–2 emission lines. The rotation signatures show evidence for vertically stratified disks, in which ^13 CO rotates faster than ^12 CO due to a distinct thermal gas pressure gradient at their emitting heights. We find δυ _ϕ substructures in the sample on both small (∼10 au) and large (∼100 au) radial scales, reaching deviations up to 15% from background Keplerian velocity in the most extreme cases. More than 75% of the rings and 80% of the gaps in the dust continuum emission resolved in δυ _ϕ are colocated with gas pressure maxima and minima, respectively. Additionally, gas pressure substructures are observed far beyond the dust continuum emission. For the first time, we determined the gas pressure derivative at the midplane from observations, and found it to align well with the dust substructures within the given uncertainties. Based on our findings, we conclude that gas pressure variations are likely the dominant mechanism for ring and gap formation in the dust continuum.

The key planet-formation processes in protoplanetary disks remain an active matter of research. One promising mechanism to radially and azimuthally trap millimeter-emitting dust grains, enabling them to concentrate and grow into planetesimals, is anticyclonic vortices. While dust observations have revealed crescent structures in several disks, observations of their kinematic signatures are still lacking. Studying the gas dynamics is, however, essential to confirm the presence of a vortex and understand its dust trapping properties. In this work, we make use of the high-resolution and sensitivity observations conducted by the exoALMA large program to search for such signatures in the ^12 CO and ^13 CO molecular line emission of four disks with azimuthal dust asymmetries: HD 135344B, HD 143006, HD 34282, and MWC 758. To assess the vortex features, we constructed an analytical vortex model and performed hydrodynamical simulations. For the latter, we assumed two scenarios: a vortex triggered at the edge of a dead zone and of a gap created by a massive embedded planet. These models reveal a complex kinematical morphology of the vortex. When compared to the data, we find that none of the sources show a distinctive vortex signature around the dust crescents in the kinematics. HD 135344B exhibits a prominent feature similar to the predictions from the simulations, thus making this the most promising target for sensitive follow-up studies at higher resolution and in particular with less abundant molecules at higher resolution and sensitivity to trace closer to the disk midplane.

We analyze the ^12 CO J = 3–2 data cubes of the disks in the exoALMA program. 13/15 disks reveal a variety of kinematic substructures in individual channels: large-scale arcs or spiral arms, localized velocity kinks, and/or multiple faints arcs that appear like filamentary structures on the disk surface. We find kinematic signatures that are consistent with planet wakes in six disks: AA Tau, SY Cha, J1842, J1615, LkCa 15, and HD 143006. Comparison with hydrodynamical and radiative transfer simulations suggests planets with orbital radii between 80 and 310 au and masses between 1 and 5 M _Jup . Additional kinematic substructures limit our ability to place tight constraints on the planet masses. When the inclination is favorable to separate the upper and lower surfaces (near 45°, i.e., in 7/15 disks), we always detect the vertical CO snowline and find that the ^12 CO freeze-out is partial in the disk midplane, with a depletion factor of ≈10 ^−3 –10 ^−2 compared to the warm molecular layer. In these same seven disks, we also systematically detect evidence of CO desorption in the outer regions.

Aims. The James Webb Space Telescope now enables the spectral study of ices with unprecedented sensitivity and angular resolution. Water ice plays a crucial role in the growth of grains and in planetary formation but its spatial distribution in protoplanetary disks is poorly constrained. To aid the interpretation of future observations, we study here for the first time how the water ice band depends on the observer's perspective and the location where spectra are measured within protoplanetary disks. Methods. Based on a standard protoplanetary disk model around a T Tauri star, we used the radiative transfer code MCFOST to extract water-ice spectra and to measure the depth and central wavelength of the water-ice band at different locations in the disk. Results. Even in the context of a spatially homogeneous ice mixture, the observed properties of water-ice bands depend on the inclination of the system as well as on the location in the disk from which the spectra are extracted. In particular, the wavelength of the band minimum can change by up to 0.17 mu m, comparable to the difference expected between amorphous and crystalline ices, for instance. This phenomenon stems from a balance between absorption and scattering and must be taken into account in detailed modeling of spatially resolved infrared spectroscopy of ices, including CO and CO2.

The early stages of planet formation, involving dust grain growth and planetesimals formation, remain shrouded in mystery. The analysis of the Scattering Phase Function (SPF) measured in disks surrounding young stars holds great potential for revealing crucial information about dust grain properties. Given the increasing number of high-quality datasets available, an efficient method to extract the SPF is required. DRAGyS is a tool designed for the quick and comprehensive analysis of ring-shaped protoplanetary disks. It directly estimates the disk geometry and extracts the total and polarized SPF from scattered light images, without requiring any radiative transfer modeling, a limitation of previous efforts. Key disk parameters (inclination, position angle, aspect ratio) are obtained by fitting ellipses to the disk intensity peaks from the ring surface, assuming the disks are circular. We validated the method using simulated disk images and then applied it to archival polarized-intensity images of nine images for six protoplanetary disks. DRAGyS provides a method to correct for the effect of limb brightening on the SPF. DRAGyS recovers well the injected geometry and the SPF from synthetic images where the parameters are known. When compared to previously published results extracted from images without considering limb brightening, DRAGyS yields similar results for the inclination, position angle, and SPF. We show that the effect of limb brightening on the SPF is significant, with consequences for the inference of dust properties. DRAGyS takes advantage of a fast and purely geometrical approach to estimate ringed-disk geometries. It allows the efficient extraction of SPF either globally or by sectors, allowing it to deal with disk asymmetries. By bypassing the need for a full modeling of the disk geometry before SPF extraction, DRAGyS is well suited to study large samples of disk images.

Context. Current models of binary systems often depend on simplified approximations of the radiation field, which are unlikely to accurately capture the complexities of asymmetric environments.Aims. We investigate the dynamical and chemical implications of a 3D asymmetric radiation field that accounts for the optical properties of sub-structures present in a protoplanetary disk, as well as the inclusion of a secondary radiation source in binary systems.Methods. We conducted a series of 3D-SPH hydrodynamical simulations using Phantom, coupled with the 3D Monte Carlo radiative transfer code Mcfost, to compute disk temperatures on-the-fly. We explored different binary-disk orientations (0 degrees and 30 degrees) for an eccentric binary, along with a constant dust-to-gas ratio and dust as a mixture prescription. We also simulated an outburst event as an example of a drastic increase in luminosity.Results. Heating from the secondary star inflates the outer disk, increasing the aspect ratio facing the companion by about 25% in inclined configurations compared to 10% in coplanar ones. Dust settling in the mid-plane enhances extinction along the disk plane, making the coplanar configuration cooler than the inclined one on the side of the disk facing the companion. Additional heating causes a shift in the snow line for species with freeze-out temperatures below 50 K, depending on the disk-binary inclination and binary phase. During outbursts, the aspect ratio doubles on the star-facing side and increases by 50% on the opposite side in inclined cases. The snow line shift would impact all the species considered in the outburst case.Conclusions. Protoplanetary disk heating in binary systems depends on stellar properties, the binary phase, and disk local and global characteristics. This results in temperature asymmetries, especially during secondary star outbursts, leading to variations in aspect ratio and snow lines that can affect chemistry and planet formation.

We report the serendipitous detection of the SO J _N = 6 _5 –5 _4 (219.949 GHz) rotational transition in archival Atacama Large Millimeter/submillimeter Array observations of the spiral hosting protoplanetary disks around CQ Tau (with ≈4.9 σ significance) and MWC 758 (with ≈3.4 σ significance). In the former, the SO emission comes in the shape of a ring, arises from the edge of the continuum cavity, and is qualitatively consistent, at the currently available spectral resolution, with being in Keplerian rotation. In the latter, instead, while arising primarily from inside the continuum cavity, the SO emission also extends to the continuum ring(s), and its morphology and kinematics are less clear. We put these sources in the context of the other protoplanetary disks where SO detections have been previously reported in the literature and discuss the possible origins of SO in terms of (thermal) desorption or formation in the gas-phase. We argue that these processes might be fostered by dynamical perturbations caused by unseen embedded massive companions, shadows, or late-time infall, thus suggesting a possible link between perturbed dynamics and SO emission in (these) protoplanetary disks. If confirmed, our interpretation would imply that chemical evolution timescales could be significantly shorter in these systems than is commonly assumed, indicating that dynamical perturbations might influence the composition of newborn (proto)planets by altering the volatile makeup of their formation environment.