pith. sign in

arxiv: 2605.18438 · v2 · pith:YWLWHIRJnew · submitted 2026-05-18 · 🌌 astro-ph.GA

Velocity-resolved [O I] 63,145 um, [C II] 158 um, and OH mapping along the Orion BN/KL explosive outflow and irradiated shocks

J. R. Goicoechea , R. G\"usten , B. Godard , H. Wiesemeyer , R. Higgins , A. Gusdorf , M. G. Santa-Maria , M. Zannese
show 1 more author
K. M. Menten
This is my paper

Pith reviewed 2026-05-20 09:14 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords Orion BN/KLexplosive outflow[O I] fine-structure linesdissociative J-type shocksUV irradiationmass-loss rateshock excitationOrion molecular cloud
0
0 comments X

The pith

Velocity-resolved [O I] and [C II] maps show Orion BN/KL outflow emission arises from 30-40 km/s dissociative J-type shocks under external UV irradiation.

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

The paper delivers the first velocity-resolved maps of the [O I] 63 and 145 micron lines across the wide-angle explosive outflow in Orion BN/KL, along with [C II] and OH data. It separates a broad outflow component with 20-30 km/s FWHM from the narrow quiescent cloud gas and finds that the line-wing intensity ratios reach unusually high values. These ratios are compared to grids of magnetized shock models that include UV irradiation. The comparison shows the data are consistent with dissociative J-type shocks at 30-40 km/s and preshock densities of a few 10^4 cm^{-3} illuminated by external UV from surrounding fast shocks or nearby massive protostars. The resulting [O I] luminosity implies a large outflow mass-loss rate and total mass.

Core claim

The [O I] and [C II] intensities are consistent with emission from dissociative J-type shocks with velocities of 30-40 km/s and preshock gas densities of a few 10^4 cm^{-3}, illuminated by external UV radiation generated by surrounding fast shocks and possibly by massive (proto)stars in the region.

What carries the argument

Grids of magnetized dissociative J-type shock models that include external UV irradiation, against which the observed [O I] 63/145 and [O I] 63/[C II] line-wing intensity ratios are compared.

If this is right

  • The total [O I] 63 and 145 um luminosity of 86.5 solar luminosities implies an outflow mass-loss rate of (9.1 +/- 2.6) x 10^{-3} solar masses per year.
  • The outflow mass is estimated at 3.3 to 5.9 solar masses.
  • Postshock gas reaches densities of 10^5 to 10^6 cm^{-3} and temperatures near 500 K.
  • The OH 119 um line shows a P-Cygni profile spanning 160 km/s, matching the velocity width of the broadest CO lines.

Where Pith is reading between the lines

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

  • Similar explosive outflows in other young clusters may be identified by searching for comparably high-velocity [O I] components with elevated [O I]/[C II] ratios.
  • External UV from overlapping shocks could set the overall cooling and ionization balance across larger portions of the Orion molecular cloud.
  • Higher-resolution mapping could isolate emission from individual shock surfaces embedded in the wide-angle outflow.

Load-bearing premise

The high observed line-wing ratios are produced by post-shock gas in magnetized dissociative J-type shocks rather than by UV-pumped PDRs or unresolved multiple velocity components.

What would settle it

Future spectra showing [O I] 63/145 or [O I] 63/[C II] ratios in the line wings that lie well outside the ranges predicted by the 30-40 km/s shock models at a few 10^4 cm^{-3} preshock density with added UV would rule out the shock interpretation.

Figures

Figures reproduced from arXiv: 2605.18438 by A. Gusdorf, B. Godard, H. Wiesemeyer, J. R. Goicoechea, K. M. Menten, M. G. Santa-Maria, M. Zannese, R. G\"usten, R. Higgins.

Figure 1
Figure 1. Figure 1: BN/KL outflow and the Trapezium cluster observed with JWST/NIRCam H2 F212N (McCaughrean & Pearson 2023). The dashed square marks the field of view mapped with SOFIA/GREAT. Cyan contours show the [O i] 63 µm redshifted line-wing emission over vLSR = 25–30 km s−1 , from 15 to 90 K km s−1 in steps of 15 K km s−1 . Snell et al. (1984), the outflow contains ∼ 8 M⊙ of molecular gas. About half of this mass is at… view at source ↗
Figure 2
Figure 2. Figure 2: Total line intensity maps (in erg s−1 cm−2 sr−1 = “cgs”) integrated over the complete line profile. The beam size is indicated in the bottom-right corner of each panel except for H2, where it is too small to display (0.1′′). The cyan triangle shows a position near the Trapezium, where the [O i] 63µm line intensity peaks (Table D.1). The H2 image refers to the JWST/NIRCam F212N image (McCaughrean & Pearson … view at source ↗
Figure 3
Figure 3. Figure 3: Total line intensity ratio maps (derived from line intensities in erg s−1 cm−2 sr−1 ). (a) [O i] 63/145 at a common angular resolution of 13′′ , (b) [O i] 63/[C ii] 158 at 15′′, and (c) [O i] 145/ [C ii]158 at 15′′. The dashed box shows the (100′′×80′′) area used to extract both the line luminosities in the outflow region ( [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Velocity-resolved spectra at representative positions, all from maps convolved to 15′′ (except CO 10–9, with a beam of 20′′). Left: Complete spectra, with offsets in arcseconds given in parenthesis; continuum levels of [Cii] 158 µm and [O i] 145 µm shifted for clarity. The OH 163 µm emission lines (yellow) are scaled by a factor of three or five. Right: Zoom on faint line wings and foreground line features… view at source ↗
Figure 5
Figure 5. Figure 5: P-Cygni profile of the OH 119 µm line toward the hot core region (the FIR continuum peak), compared to other line profiles. See Fig. B.3 for a zoom on the low-intensity features. “Northern Dark Lane” and the “Dark Bay”, which are apparent in optical images (O’Dell 2001). Because of the lower density in these components (nH ≈ 103 cm−3 ; Pabst et al. 2020), the fainter 158 µm continuum, and the much lower cr… view at source ↗
Figure 6
Figure 6. Figure 6: Possible compact outflow or wind bow shock around BN star. Left: Spatial distribution of the redshifted [O i] 63 µm emission in the vLSR range from +20 to +25 km s−1 . This map (in R.A. and Dec.) shows extended emission following the wide-angle H2 outflow and a possible spatially unresolved outflow around BN, revealed as a distinct broad spectral component. Right: [O i] 63 µm spectra toward Source I (red),… view at source ↗
Figure 7
Figure 7. Figure 7: FIR line intensity ratios as a function of velocity in a 15′′ beam toward Peak 1 (filled squares) and averaged over the 100′′ × 80′′ region of the BN/KL outflow (empty squares), binned in intervals of 5 km s−1 from maps convolved to a common angular resolution of 15′′. The dashed magenta line marks the LSR velocity of the quiescent gas in OMC-1. [O i] 63/145 line intensity ratio, 11.4 ± 0.1, occurs at vLSR… view at source ↗
Figure 8
Figure 8. Figure 8: [O i] 63 / 145 µm line intensity ratio versus [O i] 145 µm line inten￾sity. The light grey to black curves show a grid of radiative transfer mod￾els for four values of N(O): 2×1016 , 2×1017 , 2×1018, and 2×1019 cm−2 , from left to right, at 200 and 500 K. Tick marks denote log nH (in cm−3 ) and are separated by one decade along curves of constant N(O). The col￾ored squares show the observed values in LSR v… view at source ↗
Figure 9
Figure 9. Figure 9: compares the intensities of the [O i]63 µm, [O i]145 µm, and [C ii]158 µm lines observed in the red and blue wings toward Peak 1 ( [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: [O i] 63 / [C ii] 158 µm line intensity ratio as a function of shock velocity and preshock gas density in an externally irradiated shock with G ext 0 = 103 (from Lehmann et al. 2022). The horizontal blue and red lines indicate the observed line-wing intensities toward Peak 1. J-type shocks with velocities of 30–40 km s−1 , preshock densi￾ties of a few 104 cm−3 , and compression factors consistent with pos… view at source ↗
read the original abstract

Stellar mergers produce explosive outflows that serve as transient sources of IR line luminosity and inject mechanical energy early into the natal molecular cloud. We present the first velocity-resolved maps of the [O I] 63 and 145 um fine-structure line emission from the wide-angle outflow in Orion BN/KL, the nearest explosive outflow. The data were obtained with SOFIA and include sensitive [C II] 158 um and OH maps. They allowed us to disentangle the quiescent cloud gas from the outflow, traced by a broader [O I] component with a line FWHM of about 20-30 km/s and exhibiting a spatial distribution similar to that of the shock-excited H2 emission seen with JWST. The OH 119 um line shows a prominent P-Cygni profile covering 160 km/s, similar to the very broad CO lines. The total [O I] 63 and 145 line luminosity is remarkably high, 86.5 L_sun, comparable to the H2 and CO line luminosities, implying an outflow mass-loss rate of (9.1+/-2.6)x10^-3 M_sun/yr and a mass of 3.3-5.9 M_sun. The [O I] 63 / 145 and [O I] 63 / [C II] 158 intensity ratios reach very high values in the line wings (20-30 and 40-60, respectively), exceeding those found in PDRs. These ratios are consistent with the presence of dense (10^5 to 10^6 cm^-3 ) and warm (~500 K) postshock gas. We analyzed the fine-structure line-wing intensities using magnetized shock models that include UV irradiation, to which the [C II] 158 line intensity is particularly sensitive. We find that the [O I] and [C II] intensities are consistent with emission from dissociative J-type shocks with velocities of 30-40 km/s and preshock gas densities of a few 10^4 cm^-3, illuminated by external UV radiation generated by surrounding fast shocks and possibly by massive (proto)stars in the region. We also report a broad [O I] emission feature around the BN star, which we attribute to an unresolved outflow or wind bow shock.

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 the first velocity-resolved SOFIA maps of [O I] 63,145 μm, [C II] 158 μm, and OH lines toward the Orion BN/KL explosive outflow. It isolates a broad (FWHM ~20-30 km/s) outflow component whose spatial distribution matches JWST H2 shocks, reports a total [O I] luminosity of 86.5 L⊙, derives an outflow mass-loss rate of (9.1±2.6)×10^{-3} M⊙ yr^{-1}, and shows that the high line-wing ratios ([O I] 63/145 = 20-30; [O I] 63/[C II] = 40-60) are consistent with magnetized dissociative J-type shocks (v_shock = 30-40 km/s, n_pre ~ few×10^4 cm^{-3}) illuminated by external UV from surrounding shocks and/or massive protostars. A broad [O I] feature near BN is also noted.

Significance. If the shock interpretation holds, the work supplies direct observational constraints on the physical conditions and energy injection in the nearest explosive outflow, including quantitative mass-loss and luminosity comparisons to H2 and CO. The velocity-resolved data and explicit model-grid comparison constitute a clear advance over prior unresolved observations. The reported luminosities and ratios are reproducible from the new SOFIA spectra and can be tested against future higher-resolution observations or refined shock grids.

major comments (1)
  1. [model comparison and line-ratio discussion] The central claim that the observed wing ratios require (or are uniquely matched by) 30-40 km/s dissociative J-shocks rests on the statement that they exceed PDR values. An explicit side-by-side comparison to a PDR grid at the relevant high G0 (generated by surrounding fast shocks) and n = 10^4-10^6 cm^{-3} is not shown; without it, the exclusion of UV-pumped PDR or blended-velocity-component alternatives remains incompletely demonstrated (see the model-analysis paragraph following the ratio measurements).
minor comments (2)
  1. [Results] The OH 119 μm P-Cygni profile is described as covering 160 km/s; a quantitative comparison of its velocity extent to the CO line wings would strengthen the kinematic linkage.
  2. [Discussion] Notation for the preshock density range (“a few 10^4 cm^{-3}”) could be made more precise by quoting the exact grid points used in the shock models.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the positive overall assessment. We address the single major comment below and will incorporate the requested clarification in the revised version.

read point-by-point responses

  1. Referee: [model comparison and line-ratio discussion] The central claim that the observed wing ratios require (or are uniquely matched by) 30-40 km/s dissociative J-shocks rests on the statement that they exceed PDR values. An explicit side-by-side comparison to a PDR grid at the relevant high G0 (generated by surrounding fast shocks) and n = 10^4-10^6 cm^{-3} is not shown; without it, the exclusion of UV-pumped PDR or blended-velocity-component alternatives remains incompletely demonstrated (see the model-analysis paragraph following the ratio measurements).

    Authors: We agree that an explicit side-by-side comparison to PDR models at high G0 and the relevant densities would make the exclusion of PDR or blended-component interpretations more complete. In the revised manuscript we will add a new panel or table that directly overlays the observed wing ratios ([O I] 63/145 = 20-30 and [O I] 63/[C II] = 40-60) on predictions from standard PDR codes (Meudon PDR and Cloudy) for G0 = 10^4–10^6 and n = 10^4–10^6 cm^{-3}. This comparison will show that even under the strong external UV fields expected from surrounding shocks, PDR models yield [O I] 63/145 ratios below ~10. We will also note that the spatial coincidence with the JWST H2 shocks and the ~20–30 km s^{-1} line widths already argue against significant blending with quiescent PDR gas. These additions will be placed immediately after the ratio measurements. revision: yes

Circularity Check

0 steps flagged

No significant circularity; central inferences drawn from new observations matched to independent external model grids

full rationale

The paper reports new velocity-resolved SOFIA maps of [O I], [C II] and OH lines and directly compares the observed line-wing intensity ratios (e.g., [O I] 63/145 ≈ 20-30 and [O I] 63/[C II] ≈ 40-60) to published magnetized dissociative J-type shock model grids that incorporate external UV irradiation. These models pre-exist the present observations and are not constructed or fitted from the current dataset; the velocity and density inferences (30-40 km/s, few × 10^4 cm^{-3}) are therefore external benchmarks rather than self-derived quantities. No equations, parameters, or uniqueness theorems are defined in terms of the target results, and no self-citation chain is invoked to close the argument. The derivation chain remains self-contained against external, falsifiable shock-model predictions.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The analysis rests on standard astrophysical assumptions about shock excitation and external UV fields, plus two fitted parameters (shock velocity and preshock density) chosen to match the observed line ratios.

free parameters (2)
  • shock velocity = 30-40 km/s
    30-40 km/s range selected to reproduce the broad line wings and high [O I]/[C II] ratios in the magnetized shock models.
  • preshock density = few 10^4 cm^-3
    A few 10^4 cm^-3 chosen to place the postshock gas in the density regime that produces the observed [O I] 63/145 ratios of 20-30.
axioms (2)
  • domain assumption Line-wing emission arises from postshock gas in dissociative J-type shocks rather than PDR or other mechanisms
    Invoked when the observed ratios exceed typical PDR values and are compared to the shock-model grid.
  • domain assumption External UV radiation field is generated by surrounding fast shocks and massive (proto)stars
    Used to explain the sensitivity of the [C II] 158 μm line and to justify inclusion of UV irradiation in the models.

pith-pipeline@v0.9.0 · 6035 in / 1595 out tokens · 41384 ms · 2026-05-20T09:14:28.213767+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.