pith. sign in

arxiv: 2604.11904 · v2 · submitted 2026-04-13 · 🌌 astro-ph.SR · astro-ph.EP· astro-ph.GA

JOYS: Launching and destruction of dust in protostellar jets. The case of BHR71-IRS1 with JWST/MIRI

{\L}ukasz Tychoniec , Logan Francis , Maria Gabriela Navarro , Jakobus M. Vorster , Ewine F. van Dishoeck , Alessio Caratti o Garatti , Korash Assani , Valentin J. M. Le Gouellec
show 19 more authors
Beno\^it Tabone Pamela Klaassen Adriaan G. M. Janssen Kay Justtanont Daniel Harsono Pooneh Nazari Simon Reyes Katerina Slavicinska Caroline Gieser Tyler Bourke Yao-Lun Yang Brunella Nisini Teresa Giannini Henrik Beuther R. Devaraj Thomas P. Ray Nashanty G. C. Brunken Yuan Chen Martijn L. van Gelder
This is my paper

Pith reviewed 2026-05-10 15:09 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.EPastro-ph.GA
keywords protostellar jetsdust launchingClass 0 protostarJWST/MIRIshock processingmid-infrared continuumBHR71-IRS1refractory depletions
0
0 comments X

The pith

JWST observations show dust grains are launched in a Class 0 jet and partly survive shock processing.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The study uses JWST/MIRI spectral imaging of the BHR71-IRS1 outflow to map atomic fine-structure lines and construct PSF-subtracted mid-IR continuum maps. Line intensities are compared to shock models to derive declining shock velocities and densities along the jet, plus depletions in gas-phase iron and nickel relative to solar values. Continuum spectral energy distributions reveal both warm and cold dust components aligned with the jet knots and outflow cone. These results indicate that solids can be lifted from planet-forming disks into outflows at the earliest stages and resist complete vaporization in shocks. The findings matter because they connect disk solids directly to material carried away in jets during star formation.

Core claim

The ionized jet is spatially resolved with a unique set of refractory, volatile, and noble-gas lines concentrated in four bright knots that wiggle along the axis. PSF-subtracted continuum maps show extended mid-IR emission co-spatial with the jet bullets, with SED fits indicating a 200-400 K warm dust component and a 70-90 K cold component. Shock modeling of the mid-IR lines constrains a drop in velocity from 70 to 35 km/s and pre-shock density from above 10^5 to 4x10^4 cm^-3 with distance, while measurable depletions of Fe and Ni indicate that a substantial fraction of refractories remains locked in grains despite the shocks.

What carries the argument

PSF-subtracted mid-IR continuum maps from the MIRI-MRS cubes combined with shock models applied to the intensities of fine-structure atomic lines to quantify both physical conditions and elemental depletions.

If this is right

  • Dust can be launched from the disk into the jet already during the Class 0 phase.
  • A measurable fraction of refractory elements stays locked in grains after passing through the shocks.
  • Shock velocity and density decrease systematically with distance from the protostar.
  • Mid-IR fine-structure lines provide a direct probe of both gas-phase depletions and dust presence in outflows.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The survival of grains implies that some solids can be transported outward and potentially re-accreted or incorporated into later planet-forming material.
  • Similar mid-IR continuum signatures could be searched for in other Class 0 sources to test whether dust launching is common.
  • The observed wiggling knots and declining shock conditions offer a template for modeling how jets evolve and process material over hundreds of au.

Load-bearing premise

The extended mid-IR continuum after point-spread-function subtraction represents thermal emission from dust grains launched in the jet rather than scattered light or other non-thermal contributions.

What would settle it

Higher-resolution or polarimetric imaging that attributes the continuum to scattered light instead of thermal dust, or shock models that reproduce the observed line depletions with no surviving grains at all.

Figures

Figures reproduced from arXiv: 2604.11904 by Adriaan G. M. Janssen, Alessio Caratti o Garatti, Beno\^it Tabone, Brunella Nisini, Caroline Gieser, Daniel Harsono, Ewine F. van Dishoeck, Henrik Beuther, Jakobus M. Vorster, Katerina Slavicinska, Kay Justtanont, Korash Assani, Logan Francis, {\L}ukasz Tychoniec, Maria Gabriela Navarro, Martijn L. van Gelder, Nashanty G. C. Brunken, Pamela Klaassen, Pooneh Nazari, R. Devaraj, Simon Reyes, Teresa Giannini, Thomas P. Ray, Tyler Bourke, Valentin J. M. Le Gouellec, Yao-Lun Yang, Yuan Chen.

Figure 1
Figure 1. Figure 1: Top left: Large-scale view of the BHR71 globule in K s (red) and H (green) bands from the Persson’s Auxiliary Nasmyth Infrared Camera (PANIC; Martini et al. 2004) taken on 2009 January 17 and 18 and J (blue) band from Infrared Side-Port Imager (ISPI; van der Bliek et al. 2004) taken on 2009 June 11. See also (Tobin et al. 2010, 2019). Stars mark positions of the protostars from ALMA high-resolution images … view at source ↗
Figure 2
Figure 2. Figure 2: MIRI-MRS integrated Gaussian intensity maps of selected atomic and ionic emission lines. The coordinates are relative to the IR protostar source position indicated with the white star. The line quantum identifiers are listed in Table E.3. On the top left, regions identified as jet bullets are indicated. In the bottom-right corners, MIRI-MRS empirical FWHM of PSF (Law et al. 2023) is shown as a white circle… view at source ↗
Figure 3
Figure 3. Figure 3: MIRI-MRS integrated Gaussian intensity maps of representative H2 emission lines of the BHR71-IRS1 outflow. From left to right lines of increasing Eup: H2 v = 1 – 1 S(5), 10340 K; H2 v = 0 – 0 S(7), 7197 K; H2 v = 0 – 0 S(3), 2503 K; H2 v = 0 – 0 S(1), 1015 K. In the bottom-right corners, MIRI-MRS empirical FWHM of PSF (Law et al. 2023) is shown as a white circle. On the rightmost image, ALMA CO (2-1) integ… view at source ↗
Figure 4
Figure 4. Figure 4: MIRI-MRS integrated Gaussian intensity maps of molecular lines. For molecules with multiple lines detected, the lines are stacked. Left: rovibrational CO (v = 1 – 0) map is integrated from all unblended lines from Channel 1 in 4.9 to 5.1 µm range, in total 14 lines are stacked, with Eup = 4725 − 7338 K; Left middle: CO2, only the Q-branch at 14.98 µm (Eup = 964K) is imaged; Right middle: OH rotational line… view at source ↗
Figure 5
Figure 5. Figure 5: MIRI-MRS images of continuum emission from cubes centered on 6.3, 13.8, 19.6, 22.4, and 26.4 µm, from left to right. The images were generated by integrating cubes over 0.6 µm around the central wavelength. The cubes were PSF-subtracted before imaging. White star marks the position of the protostar in mid-IR. middle left). Both CO2 and the corresponding continuum com￾ponent match the width of the low-excit… view at source ↗
Figure 6
Figure 6. Figure 6: Line ratios of selected species at B1-B4 positions plotted as a function of the deprojected distance from the source. Left: Ratios of the same species sensitive to gas temperature: [Co ii] 14.5 µm to [Co II] 10.5µm (orange); [Fe ii] 17.9 µm to [Fe ii] 25.9 µm (blue); [Ni ii] 10.6 µm to [Ni ii] 6.6 µm (green); Right: Ratios sensitive to ionization fraction: [Fe ii] 25.9µm to [Fe i] 24.0 µm (grey); [Ne ii 12… view at source ↗
Figure 8
Figure 8. Figure 8: MIRI-MRS centroid velocity with respect to the vLS R of the pro￾tostar and corrected for source inclination. Top: Centroid velocity as a function of the upper energy level of the transition line. Color shows different positions along the jet (see [PITH_FULL_IMAGE:figures/full_fig_p007_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Fits to SED at the position B1. Extracted SED is smoothed with a Gaussian kernel. Dashed lines indicate the scattered, warm, and cold dust components, shown in blue, orange, and green, respectively. The solid red line shows the sum of the three blackbodies. The purple dashed line shows the sum of blackbody temperatures extincted with the KP5 extinction curve. of oscillation opposite to that seen in H2. Thi… view at source ↗
Figure 10
Figure 10. Figure 10 [PITH_FULL_IMAGE:figures/full_fig_p008_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Top: Properties of the shock along the jet. B1-B4 positions are indicated in [PITH_FULL_IMAGE:figures/full_fig_p009_11.png] view at source ↗
Figure 12
Figure 12. Figure 12 [PITH_FULL_IMAGE:figures/full_fig_p011_12.png] view at source ↗
read the original abstract

Protostellar winds can theoretically lift solids from the planet-forming disks, but direct evidence for launched dust has been scarce so far. Numerous atomic lines that are unique to mid-infrared (IR) wavelengths reveal refractories eroded from dust grains and provide information on wind properties in the earliest stages of the star formation process. We present JWST/MIRI-MRS spectral imaging of the inner 2000 au of the BHR71-IRS1 blueshifted side of the outflow. Atomic line intensities are compared to shock models to constrain the physical conditions and elemental abundances of the outflowing gas. Dust continuum maps are constructed from PSF-subtracted cubes, and the dust spectral energy distribution is analyzed. The ionized central jet of BHR71-IRS1 is spatially resolved and imaged for the first time, revealing a unique inventory of refractory, volatile, and noble-gas fine-structure lines (Fe, Ni, Co, Cl, S, Ne, Ar). The emission is concentrated along four bright knots that wiggle along the jet axis. PSF-subtracted continuum maps reveal extended mid-IR continuum emission co-spatial with the jet bullets and within the H$_2$-traced outflow cone. Spectral energy distributions along the jet are fit together with the extinction, revealing a warm (200-400 K) and a cold (70-90 K) dust component. Shock modeling constrained by the mid-IR lines indicates a decline in shock velocity from 70 to 35 km s$^{-1}$ and pre-shock density from $>$10$^5$ to $ 4\times 10^4$ cm$^{-3}$ with distance from the protostar. Gas-phase Fe and Ni are measurably depleted relative to Solar abundances, consistent with a substantial fraction of refractories remaining locked in grains in spite of the shocks. These JWST observations provide direct evidence that dust is launched in a Class 0 jet and at least partly survives shock processing.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper presents JWST/MIRI-MRS spectral imaging of the inner ~2000 au of the BHR71-IRS1 blueshifted outflow, identifying a rich set of atomic fine-structure lines (Fe, Ni, Co, Cl, S, Ne, Ar) concentrated in four wiggling knots. Line intensities are compared to shock models to derive declining shock velocities (70 to 35 km s^{-1}) and pre-shock densities (>10^5 to 4x10^4 cm^{-3}), along with measurable depletions of gas-phase Fe and Ni. PSF-subtracted continuum maps reveal extended mid-IR emission co-spatial with the jet bullets, whose SEDs are fit after extinction correction to yield warm (200-400 K) and cold (70-90 K) dust components. The authors conclude that these data furnish direct evidence that dust is launched in a Class 0 jet and at least partly survives shock processing.

Significance. If the continuum interpretation holds, the work supplies one of the first spatially resolved mid-IR constraints on dust entrainment and refractory survival in an embedded Class 0 jet, with direct relevance to the delivery of solids to planet-forming disks. The combination of multi-species line diagnostics with continuum mapping is a methodological strength and provides falsifiable predictions for depletion patterns that can be tested with future observations.

major comments (1)
  1. [Abstract and continuum-mapping section] Abstract and continuum-mapping section: The central claim of 'direct evidence' for dust launched in the jet rests on the extended mid-IR continuum being thermal emission from grains co-spatial with the jet bullets. The PSF subtraction of the bright central source can leave residuals or permit scattered-light contributions along the outflow cone; the manuscript must supply a quantitative description of the subtraction algorithm, residual maps, uncertainty propagation, and any tests (e.g., spatial correlation with H2 or multi-epoch checks) that exclude non-thermal origins. Without this, the 'launched dust' component of the conclusion lacks a secure observational basis even if the line-based depletion arguments remain intact.
minor comments (2)
  1. [Abstract] Abstract: the statement that 'Gas-phase Fe and Ni are measurably depleted' should be accompanied by the actual depletion factors and their uncertainties so readers can judge consistency with the shock-model grid.
  2. [Shock-modeling paragraph] Shock-modeling paragraph: the two free parameters (shock velocity and pre-shock density) are stated to decline with distance; the text should explicitly show how the full set of observed lines (including noble-gas species) jointly constrains these parameters and whether geometry or atomic-data uncertainties affect the derived depletions.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The single major comment concerns the robustness of the continuum interpretation and the need for more quantitative documentation of the PSF subtraction. We address this point directly below and will revise the manuscript to incorporate additional details and tests as outlined.

read point-by-point responses

  1. Referee: [Abstract and continuum-mapping section] Abstract and continuum-mapping section: The central claim of 'direct evidence' for dust launched in the jet rests on the extended mid-IR continuum being thermal emission from grains co-spatial with the jet bullets. The PSF subtraction of the bright central source can leave residuals or permit scattered-light contributions along the outflow cone; the manuscript must supply a quantitative description of the subtraction algorithm, residual maps, uncertainty propagation, and any tests (e.g., spatial correlation with H2 or multi-epoch checks) that exclude non-thermal origins. Without this, the 'launched dust' component of the conclusion lacks a secure observational basis even if the line-based depletion arguments remain intact.

    Authors: We agree that a more quantitative presentation of the continuum extraction is required to support the claim of launched dust. The current manuscript (Section 3.2) describes the use of the JWST pipeline with a custom PSF subtraction scaled to the central source flux, and notes the resulting extended emission is co-spatial with the H2-traced outflow and atomic-line knots. However, we will expand this section (and add an appendix) to include: (i) the precise subtraction algorithm, including the reference PSF selection, scaling factor determination, and any iterative masking; (ii) residual maps after subtraction, with quantitative metrics on residual levels relative to the extended emission; (iii) explicit propagation of subtraction uncertainties into the extracted SED fluxes and derived dust temperatures; and (iv) additional tests of thermal origin, such as the tight spatial correlation with both H2 and refractory lines (already shown in Figures 3 and 5) and the SED shape being inconsistent with scattered light or non-thermal processes. Multi-epoch checks are not possible with the existing single-epoch JWST data, but the combination of temperature fits (200-400 K warm component) and co-spatiality with shock tracers provides strong support for thermal dust emission. These revisions will be made in the next version, thereby securing the observational basis for the dust-launching conclusion while leaving the line-based depletion results unchanged. revision: yes

Circularity Check

0 steps flagged

No circularity: claims rest on new JWST data vs external shock models

full rationale

The paper's core claims derive from direct observational products (PSF-subtracted MIRI-MRS cubes, line intensities, and SED fits) compared against standard external shock models. No equations or steps reduce a 'prediction' or 'first-principles result' to the paper's own fitted inputs by construction. No self-citation chains, uniqueness theorems, or ansatzes are invoked to justify load-bearing choices. The derivation is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

Analysis depends on standard shock models and extinction corrections; two key parameters are fitted to match observed lines, and the dust interpretation rests on domain assumptions about continuum origin.

free parameters (2)
  • shock velocity = 70 to 35 km/s
    Decline from 70 to 35 km/s fitted to match atomic line intensities in shock models.
  • pre-shock density = >10^5 to 4e4 cm^-3
    Decline from >10^5 to 4e4 cm^-3 derived from the same line-fitting procedure.
axioms (2)
  • domain assumption Mid-IR fine-structure lines reliably trace gas-phase abundances released by dust erosion in shocks
    Invoked to convert line intensities into depletion factors and dust survival fraction.
  • domain assumption Standard shock models apply without major modifications to this protostellar jet environment
    Basis for deriving velocity, density, and abundance constraints.

pith-pipeline@v0.9.0 · 5822 in / 1514 out tokens · 81725 ms · 2026-05-10T15:09:34.392513+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. JOYS$+$: A JWST/MIRI survey of the evolution of H$_2$ winds and jets from low-mass protostars

    astro-ph.SR 2026-04 accept novelty 7.0

    JWST observations of 33 protostars show H2 winds broadening from ~20° to ~90° and warm mass-loss rates declining by 100x from Class 0 to Class I, matching MHD disk wind model predictions.

Reference graph

Works this paper leans on

6 extracted references · 6 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    JOYS$+$: A JWST/MIRI survey of the evolution of H$_2$ winds and jets from low-mass protostars

    Agra-Amboage, V ., Dougados, C., Cabrit, S., & Reunanen, J. 2011, A&A, 532, A59 Allen, C. W. 1973, Astrophysical quantities Anglada, G., López, R., Estalella, R., et al. 2007, AJ, 133, 2799 Argyriou, I., Glasse, A., Law, D. R., et al. 2023, A&A, 675, A111 Assani, K. D., Harsono, D., Ramsey, J. P., et al. 2024, A&A, 688, A26 Astropy Collaboration, Price-Wh...

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  2. [2]

    and would cause noticeable precession of the disk (Vioque et al. 2026). If instead the wiggle were due to orbital motion of the jet source itself, the observed single-jet geometry implies the launching star would be more massive than any companion. Ap- plying Eq. 7 of Lee et al. (2010) to the observed properties yields no physically plausible solution for...

    work page 2026

  3. [3]

    offer a powerful probe of inner disk structure and dynamical state. Disentangling precession, orbital Table C.1.Results of the H 2 line intensity analysis Reg, log(N warm) log(N hot) Nhot Nwarm Twarm Thot OPR cm−2 cm−2 % K K B1 21.2±0.2 19.2±0.3 0.9 674±75 2350±390 1.9±0.5 B2 20.2±0.1 19.0±0.2 6.0 881±80 2270±220 2.6±0.3 B3 19.8±0.1 18.4±0.1 3.8 813±52 22...

    work page 2080

  4. [4]

    We choose this extinction curve as it is based on mid-IR measurements fromSpitzerspanning a range of extinc- tions more relevant for the most embedded sources (Chapman et al

    with opacity magnitude as a free pa- rameter (τS(3)). We choose this extinction curve as it is based on mid-IR measurements fromSpitzerspanning a range of extinc- tions more relevant for the most embedded sources (Chapman et al. 2008). H2-based extinctions are reported in Table E.2. We denote this asτ S(3) rather thanτ 9.7 to avoid confusion with directly...

    work page 2008

  5. [5]

    Results are reported in Fig

    for undepleted species, xS =10 −5. Results are reported in Fig. D.2. The abundances found for the BHR71-71 shocks are predominantly lower than the solar values. For iron, the range is 0.1-0.4 of solar abundance; for nickel, 0.02-0.41; and for chlorine, 0.3-1.7. These results indi- cate a clear depletion of refractory species in the jet’s gas phase, sugges...

    work page 2013

  6. [6]

    E.1, dust properties as measured in 4.1

    In Tab. E.1, dust properties as measured in 4.1. The range of values is provided since different values are obtained based on density estimates from S and Cl. Tab. E.2 provides a summary of the properties of regions used for spectral analysis. Article number, page 18 of 21 Tychoniec et al.: JOYS: BHR71-IRS1 Fig. D.1.Extinction-corrected intensities of [Ne...

    work page 2080