---
title: 'DR21 Outflow: Explosive Dynamics in Cygnus-X'
url: https://www.emergentmind.com/topics/dr21-outflow
type: topic
---

# DR21 Outflow: Explosive Dynamics in Cygnus-X

The DR21 outflow is the powerful molecular outflow associated with the high-mass star formation region DR21 in Cygnus-X. It has long been regarded as one of the most extreme outflows in the Milky Way in mass and size, but its dynamical classification has remained contested. In one line of work, the outflow is interpreted as an explosive dispersal outflow whose CO streamers radiate from a common center with Hubble–Lemaître-like expansion; in another, it is described as a strongly collimated bipolar protostellar outflow whose red- and blue-shifted overlap can be explained without invoking an explosion [2406.03564] [2309.09687].

## 1. Astrophysical setting and nomenclature

DR21 lies in the Cygnus-X star-forming complex at a distance of **1.5 kpc**, within a larger filamentary environment that includes DR21(OH), W75S, and additional dense clumps and H II regions [2509.02679] [2309.09687]. The broader DR21 filament is not dynamically simple. Large-scale molecular-line studies found widespread infall signatures in optically thick HCO\(^+\) and \(^{12}\)CO, a typical infall speed of **\(\sim 0.6\) km s\(^{-1}\)**, and mass accretion rates of the order of a few \(10^{-3}\,M_\odot\,{\rm yr}^{-1}\) for the two main clumps, supporting a state of global gravitational collapse rather than equilibrium turbulence [1003.4198]. In the neighboring DR21(OH) clump, high-resolution N\(_2\)H\(^+\) and CH\(_3\)CN observations were interpreted in terms of convergent flows and low-velocity shocks, not classical high-velocity protostellar outflows [1108.4451].

This nomenclature matters because the phrase “DR21 outflow” is used in two ways in the literature. In the strict sense, it denotes the large DR21 Main outflow with eastern and western lobes and very large energetics [2309.09687]. In a broader complex-wide sense, it can also refer to distinct outflows in DR21(OH), including the methanol maser outflow and several compact CO/SiO flows associated with MM1/MM2 [1109.3153]. A common misconception is that every kinematic feature in the DR21 complex is outflow-driven; the filament studies show that infall, convergent flows, and sub-filamentary accretion are also major contributors to the observed velocity field [1003.4198] [1108.4451].

## 2. Morphological phenomenology

Historically, DR21 was known as a very energetic bipolar flow with eastern and western extensions seen in molecular gas and shocked H\(_2\) emission [1005.2517] [1005.2874]. Herschel-SPIRE maps of the DR21 core showed a **prominent central peak** with **extensions to the east and west** aligned with the known outflow axis, and the high-\(J\) CO ladder was interpreted as tracing warm outflow-associated gas rather than only a compact core [1005.2874]. Near-infrared echelle spectroscopy of the DR21/W75N region detected **79** discrete H\(_2\) components; within that survey, object **A9-1** was explicitly identified as a **collimated flow** corresponding to the Main DR21 outflow, and several DR21-related components exhibited radial velocities in excess of **80 km s\(^{-1}\)** [1406.5771].

The morphology changed substantially when interferometric CO(\(2-1\)) imaging began to isolate high-velocity structure. SMA observations at \(\sim 2''\) resolution revealed **about 25 high-velocity molecular filaments** made up of roughly **300 compact emission features**, arranged in a “finger” pattern radiating from a common center [1302.0208]. The authors argued that the apparent east–west bipolarity seen in H\(_2\) was misleading and could be shaped by extinction from a dense north–south dusty lane [1302.0208].

ALMA Band 6 observations sharpened the picture. A **1.3 mm continuum mosaic** at about **0.74″ \(\times\) 0.29″** resolution and rms **0.06 mJy beam\(^{-1}\)** revealed the known compact H II regions plus **five previously uncatalogued compact millimeter sources (MM1–MM5)** [2406.03564]. The decisive dataset was the CO(\(2-1\)) cube at about **0.66″ \(\times\) 0.24″** resolution. By manual channel-by-channel inspection outside the heavily contaminated cloud-velocity range near the systemic velocity, the study identified **eighteen CO streamers** distributed into narrow, coherent linear chains that radiate **quasi-isotropically from one common origin** [2406.03564]. These streamers are interspersed on the sky, include both red- and blue-shifted structures, and are not organized as a standard bipolar pair [2406.03564].

The common origin was refined to about **RA \(20^{\rm h}39^{\rm m}00^{\rm s}.8\), Dec \(+42^\circ19'37.62''\) (J2000)**, within the cometary H II region containing sources **A, B, and C** [2406.03564]. The absence of an embedded protostar or compact millimeter core exactly at that origin is a central morphological argument against a normal steady protostellar jet [2406.03564].

## 3. Kinematics and energetics

The outflow kinematics are extreme. In the SMA high-velocity CO data, the radial velocities span approximately **\(+3\) to \(+50\) km s\(^{-1}\)** for redshifted gas and **\(+3\) down to \(-80\) km s\(^{-1}\)** for blueshifted gas, with most filaments extrapolating back to a common systemic velocity of about **\(-3\) km s\(^{-1}\)**, matching the ambient molecular velocity of DR21 [1302.0208]. The ALMA analysis extended this picture: the blue-shifted emission spans roughly **\(-25\) to \(-101\) km s\(^{-1}\)** and the red-shifted emission roughly **\(+20\) to \(+69\) km s\(^{-1}\)** relative to the parent cloud, giving a total observed velocity spread from about **\(-100\) to \(+70\) km s\(^{-1}\)** [2406.03564].

A hallmark of the explosive interpretation is that individual streamers show approximately linear velocity gradients along their length, so that radial velocity increases with projected distance from the center [2406.03564]. The ALMA analysis expressed the bulk physical estimates through the usual relations
\[
M = \mu m_{\rm H}\,\Omega\,N_{\rm tot}/X_{\rm CO},
\]
\[
P = \sum_i M_i\,v_i,\qquad E_k = \tfrac{1}{2}\sum_i M_i\,v_i^2,
\]
and
\[
t_{\rm dyn} \sim \frac{R}{v} \approx 8600\ {\rm yr}.
\]
Assuming LTE and optically thin CO emission, and using excitation temperatures of **70–140 K**, the total mass of the streamer ensemble was estimated as **120–210 \(M_\odot\)**, with total momentum of order **\(10^3\,M_\odot\) km s\(^{-1}\)** and kinetic energy of about **\(10^{48}\) erg** [2406.03564]. The adopted most distant extent was about **120″** and the characteristic velocity about **100 km s\(^{-1}\)**, yielding a dynamical age of **8600 yr** [2406.03564].

The classical DR21 outflow literature also contains larger-scale, tracer-dependent energetic estimates. The 2013 SMA paper reiterated earlier values of **mass \(>3000\,M_\odot\)** and **energy \(>10^{48}\) erg**, and framed DR21 as about **20 times more energetic** than Orion BN/KL [1302.0208]. By contrast, the CASCADE 3.6 mm study, using HCO\(^+\) as the main outflow tracer and the separation method, derived a total outflow **mass of \(124\,M_\odot\)**, **momentum of \(1928\,M_\odot\,\mathrm{km\,s^{-1}}\)**, **kinetic energy of \(5.0\times10^{47}\) erg**, and dynamical times of roughly **3100–4900 yr** [2309.09687]. In that same analysis, the opening angles of both lobes decreased systematically with velocity, from about **\(\sim 80^\circ\)** to **\(\sim 20^\circ\)** over the range **5 to 45 km s\(^{-1}\)** relative to the source velocity [2309.09687]. This suggests that different tracers and velocity selections isolate different components of the DR21 flow.

## 4. Explosive dispersal outflow versus bipolar protostellar outflow

The explosive interpretation entered the DR21 literature explicitly with the SMA study titled “A 10,000 Years Old Explosion in DR21” [1302.0208]. Its argument combined five elements: a radial filamentary morphology rather than broad bipolar lobes, Hubble-like position–velocity behavior, a common kinematic origin near \(-3\) km s\(^{-1}\), the lack of a clear central powering source, and close similarity to the Orion BN/KL explosive outflow [1302.0208]. The inferred explosion center was given as
\[
\alpha(\mathrm{J2000}) = 20^h\,39^m\,1.1^s \pm 0.3^s,\qquad
\delta(\mathrm{J2000}) = +42^\circ\,19'\,37.9'' \pm 5'',
\]
and the authors suggested that **DR21-D** could be a displaced remnant of the disrupted system [1302.0208]. The proposed physical origin was the disintegration of a massive non-hierarchical stellar system, analogous in kind to Orion BN/KL but larger in scale and older in dynamical age [1302.0208].

For several years DR21 remained a debated case. The G5.89 explosive outflow paper described DR21 as the **“marginal case”** that had stood beside Orion BN/KL as the only suggested but not fully secure example [2010.13835]. That historical status changed with the ALMA DR21 study, which stated that the morphology and kinematics of the CO streamers **confirm the presence of an explosive dispersal outflow at the heart of DR21** [2406.03564]. Its strongest evidence was the combination of many narrow CO streamers, a common center, radial rather than bipolar geometry, linear velocity–distance gradients, and very large mass, momentum, and energy concentrated into a compact region [2406.03564].

The competing view is represented most clearly by the CASCADE analysis of DR21 Main [2309.09687]. Using combined IRAM 30 m and NOEMA 3.6 mm data at \(\sim 3''\) resolution, that study found **integrated intensity maps of the HCO\(^+\)** emission revealing a **strongly collimated bipolar outflow** with substantial red/blue overlap, hollow cavities in both lobes, and no evidence for elongated, filament-like structures of the type expected in confirmed explosive outflows [2309.09687]. It emphasized that red/blue overlap is not diagnostic by itself, because a bipolar flow close to the plane of the sky can show the same behavior [2309.09687]. The paper also argued that several weak structures are better interpreted as cavity walls than explosion debris, and that north–south HCO\(^+\) emission is ridge material rather than a second orthogonal outflow [2309.09687].

The literature therefore contains two explicit and incompatible characterizations of DR21 Main. A plausible implication is that the disagreement is sensitive to tracer choice, angular filtering, and the morphological criteria used to distinguish explosive ejecta from cavity walls. That inference is consistent with the fact that the explosive case is built primarily from high-velocity CO streamer identification, whereas the bipolar case is built primarily from HCO\(^+\), HCN, and dense-gas/shock tracers over the full outflow extent [2406.03564] [2309.09687].

## 5. Thermal, chemical, and shock structure

The DR21 outflow is chemically and thermally stratified. Herschel-SPIRE spectroscopy over **196–671 \(\mu\)m** detected **19 molecular lines** plus atomic lines, including the CO ladder from **\(J=4\!-\!3\) through \(J=13\!-\!12\)** [1005.2874]. LVG/RADEX modeling of the SPIRE CO and \(^{13}\)CO data yielded a warm component with **\(T_{\rm kin} \sim 125\) K**, **\(n_{\mathrm{H}_2} \sim 7\times10^4\,\mathrm{cm^{-3}}\)**, **\(N(\mathrm{CO}) \sim 3.5\times10^{18}\,\mathrm{cm^{-2}}\)**, and filling factor **\(\sim 12\%\)**, interpreted as a warm, moderately dense, clumpy component associated with the outflow [1005.2874]. A lower-excitation component with **\(T \approx 78\) K** and **\(N(\mathrm{CO}) \approx 4.5\times10^{21}\,\mathrm{cm^{-2}}\)** emerged when the SPIRE data were combined with ground-based CO observations, indicating that the outflow environment is multiphase [1005.2874].

Herschel-HIFI observations refined the heating problem. One HIFI study concluded that, despite the known shock-excitation of H\(_2\) and the clearly visible strong outflow, the emission of all lines up to **\(> 2\) THz** can be explained by **purely radiative heating** in two clumpy PDR ensembles, with hot dense clumps close to the central cluster probably dynamically affected by the outflow and a more widespread distribution of cooler dense clumps [1005.2517]. This is an important corrective to the assumption that a powerful outflow must dominate the thermal budget of every warm molecular component. Another HIFI study focused on water abundance and found broad outflow emission in both \(^{13}\)CO \(10\!-\!9\) and H\(_2\)O near **1100 GHz**, with the H\(_2\)O broad emission having **FWHM \(23.8\) km s\(^{-1}\)** centered at **\(-1.0\) km s\(^{-1}\)** [1005.2903]. The outflow **para-H\(_2\)O abundance** was derived as
\[
x(\mathrm{p\text{-}H_2O}) \approx 7\times10^{-7},
\]
far above the abundances inferred for the dense core and the foreground cloud, and interpreted as the result of warm gas liberating water from icy grain mantles [1005.2903].

The ionized and irradiated interfaces are also well constrained. HIFI detection of ortho-H\(_2\)O\(^+\) absorption at **1115 GHz** toward DR21 yielded a lower limit
\[
N(\mathrm{H_2O^+}) \ge 2.3 \times 10^{13}\ \mathrm{cm^{-2}},
\]
and the velocity distribution closely matched OH 6 cm absorption and [C II] \(158\,\mu\)m, indicating a thin layer of hot gas on clump surfaces directly facing the blue-shifted blister outflow [1005.2521]. At \(\sim 0.1\) pc resolution, para-H\(_2\)CO mapping of the DR21 filament found **24–114 K** kinetic temperatures with a mean of **\(48.3\pm0.5\) K**; it found **no direct evidence** that dense gas in DR21 is heated by shocks from the putative \(\sim 10^4\)-yr-old explosion, but it did identify **DR21W1** as a shock-dominated region associated with the western DR21 flow, with temperatures in the range **33–88 K**, exceptionally high non-thermal linewidths, and high Mach numbers [2405.18767].

Spatially resolved far-infrared spectroscopy with SOFIA FIFI-LS added a lobe-by-lobe view [2503.15059]. The **western lobe** is dominated by **high-\(J\) CO, [O I] \(63.18\,\mu\)m, and OH**, while the **eastern lobe** is dominated by **[C II] \(157.74\,\mu\)m, [O I] \(145.53\,\mu\)m, and [O III]** [2503.15059]. Non-LTE modeling of [O I] line ratios suggested **\(n({\rm H}_2)\sim10^5\,\mathrm{cm^{-3}}\)** in the west and **\(\sim10^4\,\mathrm{cm^{-3}}\)** in the east, consistent with a denser shock-dominated western interaction region and a more PDR-like or ionized eastern cavity [2503.15059]. Assuming that the bulk of the far-infrared line emission arises in shocks, the study derived an outflow power of **\(4.3\!-\!4.8\times10^2\,L_\odot\)** and a mass-loss rate of **\(3.3\!-\!3.7\times10^{-3}\,M_\odot\,\mathrm{yr}^{-1}\)** [2503.15059].

## 6. Related outflows in DR21(OH) and the surrounding complex

A full account of “DR21 outflow” usage must distinguish DR21 Main from the neighboring DR21(OH) outflows. DR21(OH), about **3 arcmin north of DR21**, hosts a cluster of compact millimeter sources and multiple compact outflows rather than one dominant flow [1811.01880] [1109.3153]. SMA continuum observations at \(\sim1''\) resolution resolved **nine compact sources** in MM1/MM2 and found that the most distinctive outflow there is a **well-collimated east–west low-velocity methanol/formaldehyde outflow** with **blueshifted emission to the east** and **redshifted emission to the west**, centered on **SMA4** in MM2 [1109.3153]. Its characteristic velocities are near the systemic value, roughly **\(-8\) to \(+5\) km s\(^{-1}\)**, and the integrated CH\(_3\)OH and H\(_2\)CO fluxes of **80 to 150 Jy km s\(^{-1}\)** correspond to isotropic luminosities of about **\(10^{-4}\,L_\odot\)** [1109.3153]. The same SMA study showed that high-velocity CO and SiO trace **other compact outflows** from MM1/MM2 and are **not related** to the low-velocity methanol maser outflow [1109.3153].

The maser phenomenology in DR21(OH) is correspondingly rich. EVLA and SMA imaging showed that the distribution of **36 GHz** class I methanol masers in the DR21(OH) outflow is similar to that of the other class I transitions, with numerous multitransition overlaps; all **14** of the 36 GHz masers in the outflow region are coincident, within positional uncertainty, with 44 GHz masers, and class I masers at **36 and 229 GHz** occur in virtual overlap with class II **6.7 GHz** masers at the main continuum source [1012.5771]. A separate SMA line survey of the CH\(_3\)OH maser outflow derived a methanol rotational temperature of **\(28\pm2.5\) K** and column density of **\(6.0\pm0.9\times10^{15}\,\mathrm{cm^{-2}}\)** at the best-constrained outflow spot, while not detecting complex molecules such as CH\(_3\)OCHO, \((\mathrm{CH_3})_2\mathrm{CO}\), and CH\(_3\)CH\(_2\)CN [1811.01880]. More recently, wide-band LTE and RADEX modeling of DR21(OH)/N44 identified a distinct broad outflow component in CH\(_3\)OH and H\(_2\)CO at **\(v_{\rm LSR}=-8.0\) km s\(^{-1}\)** with **\(\Delta v = 15\) km s\(^{-1}\)**, reinforcing the view that DR21(OH) is a highly dynamical, actively star-forming environment with multiple flows and shocks [2411.12916].

These DR21(OH) results are not measurements of the DR21 Main outflow itself. Their importance is classificatory: they show that the DR21 complex contains several physically distinct outflow systems, so shorthand references to “the DR21 outflow” must be interpreted with care.

## 7. Galactic context and high-energy significance

DR21 now occupies a central place in the emerging Galactic class of explosive dispersal outflows. The 2024 ALMA study stated that **five** dispersal explosive outflows associated with massive star-forming regions have been confirmed in the Galaxy: **Orion BN/KL, G5.89-0.39, S106-IR, IRAS 16076-5134, and IRAS 12326-6245**, with DR21 added as another clear case [2406.03564]. It also stressed that their **frequency of occurrence in the Galaxy and the originating nature are still uncertain** [2406.03564]. In the G5.89 ALMA paper, DR21 served as the benchmark older system whose kinematic age helped motivate an occurrence estimate of about **one explosion every \(\sim 100\) yr**, close to the supernova rate; in that narrative DR21 moved from being an isolated marginal exception to a member of a small but growing empirical class [2010.13835].

The broader significance of DR21 expanded further with the report of gamma-ray emission spatially coincident with the outflow [2509.02679]. Using **15 years of Fermi-LAT data** in the **0.2–500 GeV** range, the study found a detection at **greater than \(35\sigma\)**, best described spatially by an **extended radial Gaussian** with **\(\sigma = 0.65^\circ \pm 0.03^\circ\)** and spectrally by a **power law with an exponential cutoff** with
\[
\Gamma_1 = 2.08 \pm 0.02,\qquad E_{\rm c} = 10089 \pm 2963\,\mathrm{MeV}.
\]
Assuming the **1.5 kpc** distance, it inferred a gamma-ray luminosity
\[
L_{\gamma} \simeq (2.17\pm 0.15)\times 10^{35}\,\mathrm{erg\,s^{-1}}
\]
in the **0.1–500 GeV** band, and estimated that **no more than about 15\%** of the outflow’s kinetic power is converted into cosmic-ray acceleration [2509.02679]. The gamma-ray morphology overlaps dense gas tracers and allWISE mid-infrared emission, and the paper interpreted the emission in terms of hadronic \(pp\) collisions in a dense star-forming environment [2509.02679]. If that interpretation holds, explosive dispersal outflows become not only a dynamical subclass of massive-star feedback, but also a newly recognized class of Galactic high-energy sources.

Source: https://www.emergentmind.com/topics/dr21-outflow