Recognition: 1 theorem link
· Lean TheoremPRODIGE -- envelope to disk with NOEMA VIII. Sulfur oxides trace a shock caused by a streamer in the inner envelope of a protostar
Pith reviewed 2026-05-15 08:35 UTC · model grok-4.3
The pith
Sulfur dioxide peaks mark a low-velocity shock where a streamer collides with the inner envelope of protostar Per-emb 50.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the weaker SO2 peak at approximately 400 au is blueshifted with respect to the velocity field predicted by an infalling-rotating envelope, placing it at the location where a streamer impacts the inner envelope and drives a low-velocity shock; the brighter inner peak at 180 au follows envelope kinematics and may trace disk-envelope interface shocks, while both sites display SO2/SO ratios consistent with shocks of 3-4 km/s.
What carries the argument
Comparison of observed SO2 velocities against an infalling-rotating envelope kinematic model, combined with Bayesian selection of kinematic components and LTE/non-LTE derivation of temperature and density, to isolate the streamer-induced shock signature.
If this is right
- Streamers can drive observable shocks that change sulfur chemistry at the envelope-disk boundary.
- SO2 and SO abundance ratios serve as diagnostics for low-velocity shocks in embedded protostellar systems.
- Both the inner envelope and the disk can be chemically and thermally altered by streamer impact before material reaches the disk.
- Multiple kinematic components traced by SO2 can be separated to reveal distinct shock locations.
Where Pith is reading between the lines
- Similar sulfur enhancements in other protostars may trace unseen streamers rather than internal envelope processes.
- The shock heating could affect the initial conditions for disk chemistry once material lands on the disk.
- Higher-resolution observations of additional tracers could map the full streamer trajectory and its landing point.
Load-bearing premise
The kinematic model of a pure infalling-rotating envelope correctly predicts the expected velocity field without significant additional perturbations from the streamer.
What would settle it
A high-resolution position-velocity diagram in which the 400 au peak shows no blueshift relative to the envelope model, or abundance ratios inconsistent with 3-4 km/s shocks.
read the original abstract
(Abridged) Recently, streamers have been observed causing shocks at the outer edge of protoplanetary disks. The study of sulfur-bearing species can help us to understand the physical and chemical changes caused by infalling streamers toward their landing positions. We study the physical properties traced by SO$_2$ and SO toward the Class I protostar Per-emb 50, which is possibly related to the streamer infalling toward its disk. We present new NOEMA A-array observations as part of the large program "Protostars and Disks: Global Evolution" (PRODIGE). We analyzed the morphology of SO$_2$ and SO, and complement our interpretations with additional H_$2$CO and CO data from the same program. We compared the SO$_2$ and SO morphology with an infalling-rotating model. We applied Bayesian model selection to the brightest SO$_2$ line to disentangle the different kinematic components traced by this molecule. We used Local Thermodynamic Equilibrium (LTE) and non-LTE analyses to determine the temperature and density of the SO$_2$ emission. There are two separate peaks of SO$_2$ emission offset toward the southwest of Per-emb 50, one brighter (peak 1) at about 180 au from the protostar, and a weaker one (peak 2) at about 400 au. Peak 2 is blueshifted with respect to an infalling-rotating envelope. We propose that this peak is caused by the shock between the inner envelope and the streamer. Peak 1 is consistent with the expected envelope motion, and could thus be caused by shocks at the disk-envelope interface, but potential streamer influence cannot be neglected. Both peaks show abundance ratios consistent with a low velocity shock ($\sim 3-4$ \kms) when compared with shock models. Streamers can affect the physical and chemical structure of both disks and envelopes, suggesting that streamers can play an important role in shaping both structures in the embedded stages of star formation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports NOEMA A-array observations of SO2 and SO emission toward the Class I protostar Per-emb 50. Two distinct SO2 peaks are identified southwest of the source: a brighter peak 1 at ~180 au consistent with envelope motion and a weaker peak 2 at ~400 au that is blueshifted relative to an infalling-rotating envelope model. The authors interpret peak 2 as tracing a low-velocity (~3-4 km/s) shock induced by an infalling streamer, supported by Bayesian model selection on the brightest SO2 line to separate kinematic components, LTE and non-LTE analyses for temperature and density, and abundance ratios compared to existing shock-chemistry models. Additional H2CO and CO data are mentioned to complement the interpretation.
Significance. If the streamer-shock interpretation is confirmed, the result provides direct observational evidence that infalling streamers can drive shocks and alter chemistry within the inner envelopes of embedded protostars, extending prior work on disk-edge shocks to smaller scales. Strengths include the application of Bayesian model selection to disentangle velocity components and the joint LTE/non-LTE characterization of physical conditions; these elements make the physical conclusions more robust than morphology alone.
major comments (2)
- [Kinematic comparison section] Kinematic comparison section: the infalling-rotating envelope model is used to establish the expected velocity field against which the blueshift of peak 2 is measured, yet no formal uncertainties on the fitted parameters (centrifugal radius, specific angular momentum, inclination) or residual velocity maps at the ~400 au location of peak 2 are presented. Without these, it is impossible to determine whether the reported blueshift exceeds plausible model systematics or projection effects at that radius.
- [Shock-chemistry comparison] Shock-chemistry comparison: the inference that both peaks are consistent with ~3-4 km/s shocks rests on abundance ratios matched to a discrete grid of existing models. The manuscript does not report the goodness-of-fit metric or likelihood surface across the full velocity grid, leaving open whether other velocities or non-shock processes (e.g., enhanced UV or cosmic-ray ionization) could reproduce the observed ratios equally well.
minor comments (2)
- [Abstract] The abstract states that H2CO and CO data 'complement our interpretations' but does not specify which lines or how they constrain the shock versus envelope scenario; a brief sentence clarifying their quantitative role would improve clarity.
- [Figures] Figure captions for the SO2 integrated-intensity and velocity maps should explicitly state the velocity integration ranges and the contour levels used to define peak 1 and peak 2.
Simulated Author's Rebuttal
We thank the referee for their thorough review and positive recommendation. We address the major comments point-by-point below and have revised the manuscript to incorporate additional details on the kinematic modeling and shock chemistry analysis.
read point-by-point responses
-
Referee: Kinematic comparison section: the infalling-rotating envelope model is used to establish the expected velocity field against which the blueshift of peak 2 is measured, yet no formal uncertainties on the fitted parameters (centrifugal radius, specific angular momentum, inclination) or residual velocity maps at the ~400 au location of peak 2 are presented. Without these, it is impossible to determine whether the reported blueshift exceeds plausible model systematics or projection effects at that radius.
Authors: We agree with the referee that formal uncertainties and residual maps are important for robust interpretation. The original analysis used a Bayesian approach for model selection, but parameter uncertainties were not explicitly tabulated. In the revised manuscript, we will add a table or text reporting the posterior uncertainties on the centrifugal radius, specific angular momentum, and inclination from the fit. We will also include a residual velocity map for the SO2 emission, highlighting the ~400 au region to show that the observed blueshift of peak 2 is statistically significant compared to the model residuals and projection effects. revision: yes
-
Referee: Shock-chemistry comparison: the inference that both peaks are consistent with ~3-4 km/s shocks rests on abundance ratios matched to a discrete grid of existing models. The manuscript does not report the goodness-of-fit metric or likelihood surface across the full velocity grid, leaving open whether other velocities or non-shock processes (e.g., enhanced UV or cosmic-ray ionization) could reproduce the observed ratios equally well.
Authors: We thank the referee for pointing this out. While the abundance ratios were compared to the grid, no explicit fit metric was provided. In the revision, we will report the reduced chi-squared values for the best-matching models at 3-4 km/s and discuss the values for other velocities in the grid. We will also note that the kinematic evidence (blueshift and morphology) supports the shock interpretation over alternatives like UV or cosmic-ray ionization, which would not produce the observed velocity offsets. A full likelihood surface is not feasible with the current discrete models but we will expand the discussion. revision: yes
Circularity Check
No circularity: central claim is observational interpretation against standard external kinematic and shock-chemistry models
full rationale
The paper observes two SO2 peaks, compares their morphology and velocity to a pre-existing infalling-rotating envelope prescription, applies standard Bayesian component separation and LTE/non-LTE analysis, and interprets the blueshift of peak 2 as a streamer-induced shock on the basis of abundance ratios matching external low-velocity shock models. None of the load-bearing steps (model comparison, Bayesian selection, abundance-to-shock inference) reduce by construction to quantities fitted or defined inside the paper itself; the kinematic reference frame is imported rather than self-derived, and no parameter is renamed as a prediction. This is the normal non-circular case for an observational interpretation paper.
Axiom & Free-Parameter Ledger
free parameters (2)
- SO2 column density and excitation temperature
- Shock velocity grid point (~3-4 km/s)
axioms (2)
- domain assumption The velocity field of the envelope follows the standard infalling-rotating model without additional perturbations from the streamer.
- domain assumption Sulfur oxide abundance ratios are set primarily by shock heating and not by other processes such as UV irradiation or cosmic-ray ionization.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.