# The SOFIA FEEDBACK [CII] Legacy Survey: Rapid molecular cloud dispersal in RCW 79

L Bonne, S Kabanovic, N Schneider, A. Zavagno, E Keilmann, R Simon, C Buchbender, R Güsten, A M Jacob, K Jacobs, et al.

► To cite this version:

L Bonne, S Kabanovic, N Schneider, A. Zavagno, E Keilmann, et al. The SOFIA FEEDBACK [CII] Legacy Survey: Rapid molecular cloud dispersal in RCW 79. Astronomy & Astrophysics - A&A, 2023, 679, 10.1051/0004-6361/202347721. [hal-04457387]

HAL Id: hal-04457387

https://hal.science/hal-04457387v1

Submitted on 14 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
LETTER TO THE EDITOR

# The SOFIA FEEDBACK [CII] Legacy Survey: Rapid molecular cloud dispersal in RCW 79*

L. Bonne<fnref n="1" />, S. Kabanovic<fnref n="2" />, N. Schneider<fnref n="2" />, A. Zavagno<fnref n="3" />,<fnref n="4" />, E. Keilmann<fnref n="2" />, R. Simon<fnref n="2" />, C. Buchbender<fnref n="2" />, R. Güsten<fnref n="5" />,
A. M. Jacob<fnref n="5" />,<fnref n="6" />, K. Jacobs<fnref n="2" />, U. Kavak<fnref n="1" />, F. L. Polles<fnref n="1" />, M. Tiwari<fnref n="5" />, F. Wyrowski<fnref n="5" />, and A. G. G. M. Tielens<fnref n="7" />,<fnref n="8" />

Received 14 August 2023 / Accepted 2 October 2023

ABSTRACT

It has long been discussed whether stellar feedback in the form of winds and/or radiation can shred the nascent molecular cloud, thereby controlling the star formation rate. However, directly probing and quantifying the impact of stellar feedback on the neutral gas of the nascent clouds is challenging. We present an investigation of this impact toward the RCW 79 H II region using the ionized carbon line at 158 µm ([C II]) from the FEEDBACK Legacy Survey. We combine this data with information on the dozen ionizing O stars responsible for the evolution of the region, and observe in [C II] for the first time both blue- and redshifted high-velocity gas that reaches velocities of up to 25 km s⁻<fnref n="1" /> relative to the bulk emission of the molecular cloud. This high-velocity gas mostly contains neutral gas, and partly forms a fragmented shell, similar to recently found shells in a few Galactic H II regions. However, this shell does not account for all of the observed neutral high-velocity gas. We also find high-velocity gas streaming out of the nascent cloud through holes, and obtain a range of dynamical timescales below 1.0 Myr for the high-velocity gas that is well below the 2.3 ± 0.5 Myr age of the OB cluster. This suggests a different scenario for the evolution of RCW 79, where the high-velocity gas does not solely stem from a spherical expanding bubble, but also from gas recently ablated at the edge of the turbulent molecular cloud into the surrounding interstellar medium through low-pressure holes or chimneys. The resulting mass ejection rate estimate for the cloud is 0.9–3.5 × 10⁻<fnref n="2" /> M☉ yr⁻<fnref n="1" />, which leads to short erosion timescales (<5 Myr) for the nascent molecular cloud. This finding provides direct observational evidence of rapid molecular cloud dispersal.

**Key words.** HII regions – ISM: clouds – ISM: bubbles

**1. Introduction**

Feedback from OB stars, in the form of radiation and stellar winds, injects copious amounts of thermal and mechanical energy in the interstellar medium (ISM). As a result, feedback plays a central role in the evolution of the ISM, and might even maintain the low star formation rate observed in molecular clouds (Gao & Solomon 2004; Evans et al. 2009; Lada et al. 2010).

The [C II] fine-structure line at 158 µm is ideally suited to study the impact of stellar feedback on the embedding molecular cloud as it mostly probes photodissociation regions (PDRs) where far-ultraviolet (FUV) photons between 6 and 13.6 eV dominate the physical and chemical state of the gas (Hollenbach & Tielens 1999; Wolfire et al. 2022). This has made the [C II] spectral line a prominent target for observing pro-

grams on far-infrared observatories (e.g., Stacey et al. 1993; Pineda et al. 2013, 2014; Goicoechea et al. 2015). Recently, the Stratospheric Observatory for Infrared Astronomy (SOFIA) facilitated efficient velocity-resolved [C II] mapping around H II regions. This particularly revealed the presence of expanding [C II] shell morphologies at high velocity (up to ~15 km s⁻<fnref n="1" />). However, these shells were typically found to expand in one direction, mostly toward the observer (Pabst et al. 2019, 2020; Luisi et al. 2021; Tiwari et al. 2021; Kabanovic et al. 2022; Beuther et al. 2022; Bonne et al. 2022; Tram et al. 2023). Early work by Pabst et al. (2019) and Luisi et al. (2021) attributed these observed dynamic features in [C II] to a single roughly coherent expanding shell or bubble that is typically considered in analytical models (e.g., Spitzer 1978; Weaver et al. 1977). However, observations in Beuther et al. (2022) and Bonne et al. (2022) showed that [C II] emission at high velocities can also have a more complex morphology. Constraining the nature of high-velocity [C II] gas is thus important in order to understand the impact of stellar feedback on molecular cloud

<footnote n="1">¹ SOFIA Science Center, USRA, NASA Ames Research Center, Moffett Field, CA 94 045, USA
e-mail: lbonne@usra.edu</footnote>

<footnote n="2">² I. Physik. Institut, University of Cologne, Zülpicher Str. 77, 50937 Cologne, Germany</footnote>

<footnote n="3">³ Aix Marseille Univ, CNRS, CNES, LAM, Marseille, France</footnote>

<footnote n="4">⁴ Institut Universitaire de France, 1 rue Descartes, Paris, France</footnote>

<footnote n="5">⁵ Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany</footnote>

<footnote n="6">⁶ Department of Physics & Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA</footnote>

<footnote n="7">⁷ University of Maryland, Department of Astronomy, College Park, MD 20742-2421, USA</footnote>

<footnote n="8">⁸ Leiden Observatory, PO Box 9513, 2300 RA Leiden, The Netherlands</footnote>

<footnote>* Movies and interactive figures associated to Fig. 1 are available at
https://www.aanda.org</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Integrated intensity map and kinematic structure of RCW 79 in [C II]. Top left: [C II] integrated intensity map of RCW 79 between −70 and −20 km s⁻¹ with the integrated ¹²CO(3–2) emission at 25, 50, 75, and 100 K km s⁻¹ overlaid by white contours. The white plus signs and black crosses respectively give the locations of the previously identified O stars and YSOs in the region. The red arrows outline the cuts used for the PV diagrams below, with the arrow indicating a positive radius. The cyan circle delimits the region used for the average spectrum presented in Fig. 2 and the magenta circle indicates the compact H II region. Movies of PV cuts over the map are found online and at https://hera.phl.uni-koeln.de/~nschneid/rcw79_animations.html. Top right: 2D projection of an interactive 3D isocontour plot, which is found at the same url as above. Bottom left and right: [C II] PV diagrams with the ¹²CO(3–2) emission indicated by white contours. Both blue- and redshifted high-velocity [C II] emission is detected. The high-velocity [C II] gas is typically found toward the edges of the H II region, and flowing out through the opening in the northwest (see PV cut 2).</figcaption></figure>

evolution, which will have implications for molecular cloud lifetimes and the long-standing debate of whether clouds evolve dynamically or in a quasi-static fashion (e.g., Shu et al. 1987; Elmegreen 2000; Hartmann et al. 2001; Krumholz & McKee 2005; Schneider et al. 2023).

Here we present a currently unique combination of [C II] observations from the SOFIA legacy survey FEEDBACK (Schneider et al. 2020) and the results from a detailed study of the ionizing O stars (Martins et al. 2010) for the H II region RCW 79 (see Fig. 1). RCW 79 has been located at distances between 4.0 and 4.3 kpc (Russeil et al. 1998; Mège et al. 2021); here we update the distance to 3.9 ± 0.4 kpc using the Gaia DR3 release (see Appendix A). The region is heated by a cluster of 12 identified O stars, including O4-6V/III candidates, that have an estimated age of 2.3 ± 0.5 Myr (Martins et al. 2010). Gas compression by stellar feedback has created a fragmented quasi-circular ring of dense molecular gas (Zavagno et al. 2006;

Liu et al. 2017) with a radius of 6–7 pc. This ring hosts young stellar objects (YSOs) and one compact H II region in the south-east that were proposed by Zavagno et al. (2006) to be the result of the collect and collapse mechanism (Elmegreen & Lada 1977). Figure 1 shows that the O stars are located in the southwestern area of the H II region.

## 2. Observations

### 2.1. SOFIA

The [C II] ³P₃/₂ → ³P₁/₂ line at 158 µm was observed in parallel with the [O I] 63 µm line with the upgraded German REceiver for Astronomy at Terahertz frequencies (upGREAT; Risacher et al. 2018) on board SOFIA (Young et al. 2012). An area of ~470 arcmin² was mapped in the on-the-fly (OTF) mode and atmospheric calibration was done with the GREAT pipeline
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Average [C II], <sup>12</sup>CO(3–2), and <sup>13</sup>CO(3–2) spectrum toward RCW 79. The gray dashed vertical lines indicate the escape velocities of the molecular cloud relative to its central velocity of −47 km s<sup>−1</sup>. The [C II] spectrum has evident high-velocity wings outside the escape velocity range of the molecular cloud.</figcaption></figure>

(Guan et al. 2012). Including calibration, this required a total observing time of 10 h. The antenna temperature was converted to main beam temperature using a forward efficiency η<sub>for</sub> = 0.97 and a main beam efficiency η<sub>mb</sub> = 0.65. To improve the data quality and remove scanning effects in the form of stripes, we employed a method based on principal component analysis (PCA) that was also used for the data presented in other FEEDBACK papers (Tiwari et al. 2021; Kabanovic et al. 2022; Schneider et al. 2023; Bonne et al. 2023b). The nominal angular resolution of the [C II] and [O I] data is 14.1" and 6", respectively, but here we convolve the [C II] data cube to an angular resolution of 20" and a spectral binning of 0.5 km s<sup>−1</sup>. This slight smoothing and rebinning improves the quality of the data cube, for example by reducing the noise rms and striping effects. The noise rms in one channel is typically 0.5–0.6 K (see Schneider et al. 2020 for more observational details, such as the OFF position).

## 2.2. APEX

The <sup>12</sup>CO(3–2) and <sup>13</sup>CO(3–2) transitions at 345.796 GHz and 330.588 GHz, respectively, were mapped over the same region as the [C II] map with the LAsMA receiver on the APEX telescope (Güsten et al. 2006). The observations were done in the total power OTF mode, and the antenna temperature was converted to main beam temperature using η<sub>mb</sub> = 0.68. The final data cubes were produced with first-order baseline removal with an angular resolution of ~20" and a spectral binning of 0.3 km s<sup>−1</sup>. This results in a typical noise rms of 0.4–0.5 K in one channel.

# 3. Results

## 3.1. Integrated intensity map and spectra

The [C II] integrated intensity map of RCW 79 in Fig. 1 (top left) is dominated by a circular morphology very similar to that observed in the Spitzer infrared maps presented in Zavagno et al. (2006). Most of the [C II] emission is located at the interface of the H II region and the <sup>12</sup>CO emission from the nascent molecular cloud, and has a brightness peak at the location of a compact H II

region in the southeast. However, there is also emission outside of the ring, which is discussed below, and low-surface brightness emission inside the [C II] ring. Figure 2 presents the spatially averaged spectrum over the H II region which reveal prominent high-velocity wings in [C II] between −70 and −52.6 km s<sup>−1</sup> (blueshifted) and −41.4 and −25 km s<sup>−1</sup> (redshifted), that are mostly undetected in <sup>12</sup>CO. These velocities are outside the escape velocity<sup>1</sup> range of the molecular cloud. Such [C II] high-velocity wings were detected early on in spectrally resolved [C II] observations of H II regions (Simon et al. 2012; Schneider et al. 2018), but their nature has not been discussed in detail.

The different integrated emission features are also clearly identified when applying the Gaussian mixture model (GMM) methodology, introduced in Kabanovic et al. (2022), to the [C II] data cube (see Fig. 3). In short, GMM sorts spectra into clusters with similar line profiles by describing the spectral parameters<sup>2</sup> of all spectra in the data cube with a linear combination of Gaussian distributions (see Sect. 4.3.1 in Kabanovic et al. 2022 for a detailed description of the method). The brightness range of all spectra is normalized so that GMM clusters the spectra based on the line profile regardless of the absolute intensity. We then employ the Bayesian information criterion (BIC) to determine the best fitting number of clusters. Figure 3 displays the spatial distribution of the seven identified clusters in RCW 79, and shows the spatial association of kinematic substructure in the region. The bright [C II] ring is separated into three clusters: cluster 1 (blue), 2 (orange), and 3 (green) where the spectra from clusters 2 and 3, which trace the inner regions of the ring, have high-velocity wings. The wings in cluster 2 are predominantly blueshifted, while cluster 3 has equally bright blue- and redshifted wings. Clusters 4 (red) and 5 (purple) are located outside of the PDR ring and have slightly more blueshifted velocities. Clusters 6 (brown) and 7 (pink) probe the region inside the ring with brightness peaks at more extreme blue- and redshifted velocities (up to −64 km s<sup>−1</sup> and −34 km s<sup>−1</sup>, respectively).

In Fig. 4 we present the line emission integrated over the blue- and redshifted high-velocity wings. This emission is unevenly distributed for the blue- and redshifted gas along the ring that identifies the RCW 79 bubble, indicating that their structure is affected by the 3D morphology of the nascent cloud. We also note that the blue- and the redshifted high-velocity gas is observed outside of the H II region. The blueshifted gas is mostly observed south of the region, while the redshifted gas is breaking through the northwestern and northeastern holes in the ring.

## 3.2. Position-velocity diagrams

The position-velocity (PV) diagrams and movies in Fig. 1 present the [C II] and <sup>12</sup>CO(3–2) kinematics over the H II region. The PV cut 1 in Fig. 1 shows the high-velocity gas structure in [C II] that is expected for a spherical expanding shell. It is also the first time that we simultaneously clearly observe both a blue- and redshifted shell morphology. Cut 2, on the other hand, was chosen as the axis goes through a region in the northwest where the ring has clearly broken open. This shows that the [C II] high-velocity gas is not spherically expanding and that

<footnote n="1">1 $v_{\text{esc},0} = \sqrt{2G M_{\text{cloud},0}/r_{\text{cloud},0}}$ with G the gravitational constant, $M_{\text{cloud},0} (= 5.5 \times 10^4 M_{\odot})$ the original mass of the cloud (i.e., the current mass of the cloud corrected for the mass ejection that has been going on for 2.3 ± 0.5 Myr; see Sect. 4.2), and $r_{\text{cloud},0}$ (= 15 pc) the estimated original size of the cloud.</footnote>

<footnote n="2">2 The spectral parameters are the intensity of each spectral bin in the spectra.</footnote>
<figure><img src="image_4.png" /></figure>

**Fig. 3.** Output of the unsupervised GMM clustering algorithm on the [C II] data of RCW 79. Top: Cluster locations in RCW 79, representing zones with similar [C II] spectral line profiles, identified with the GMM. The contours indicate the total [C II] integrated intensity, starting at 25 K km s⁻¹ with increments of 25 K km s⁻¹, revealing that the regions identified based on the normalized [C II] spectral morphology with the GMM are clearly correlated to spatial features identified in the intensity map of the region. Bottom: Corresponding average spectra for each cluster identified by the GMM calculations on RCW 79.

high-velocity gas is flowing out of the region through the opening. The interactive 3D isocontour plot in Fig. 1 and the PV movies show the [C II] velocity structure in more detail, which further demonstrates the complex structure of the high-velocity gas. This complexity is also demonstrated by the GMM analysis that identifies four distinct clusters with high-velocity gas (see previous section). In addition, we also find that the spatial distribution of the high-velocity [C II] emission is inconsistent with a single spherical expanding bubble (see Appendix B).

# 4. Discussion

## 4.1. High-velocity mass ablation in RCW 79

To further discuss the observed [C II] high-velocity gas, we need to constrain the 3D distance of the [C II] high-velocity emission region from the ionizing O star cluster. For this we use the integrated [C II] brightness in the wings, which has val-

<figure><img src="image_5.png" /></figure>

**Fig. 4.** Spatial distribution of the blue- and redshifted [C II] high-velocity gas in RCW 79. Top: Intensity map, integrated over the blueshifted high-velocity wing (i.e., -70 to -52.6 km s⁻¹). The high-velocity gas is also observed to the south of RCW 79. The white contours indicate the total integrated [C II] emission (presented in Fig. 1) starting at 50 K km s⁻¹ with increments of 50 K km s⁻¹. Bottom: Same, but for the redshifted high-velocity wing (i.e., -41.4 to -25 km s⁻¹). The high-velocity gas is observed passing through the openings of the H II region in the northeast and northwest.

ues spread over the observed region from the 3σ detection limit (~8 K km s⁻¹) up to 50 K km s⁻¹ for each wing. Using the PDR Toolbox (Pound & Wolfire 2023), this points to a typical strength of G₀ ≈ 100–1000 for the FUV field that excites the [C II] emission in the wings (assuming a density n_H ~ 10²–10⁴ cm⁻<fnref n="3" />).

The stellar parameters from Martins et al. (2010) allow us to calculate the FUV field strength as a function of distance from the cluster<fnref n="3" />. This results in 3D distances from the O-star

<footnote n="3">³ This method follows the $r^{-2}$ decrease in FUV intensity, integrated between 6 and 13.6 eV without dust attenuation, from the surface of the stellar objects (see Bonne et al. 2020 for more details).</footnote>
<figure><img src="image_6.png" /></figure>

**Fig. 5.** Schematic illustration of the proposed scenario for RCW 79, seen edge-on (i.e., from in the plane of the sky). The green ring indicates the expanding fragmented shells with multiple holes. Through these holes, high-velocity neutral gas, indicated with the blue arrows, is escaping into the ambient ISM. Also shown (in orange) is the ring-like structure of molecular gas, wrapping the H II region, seen in CO from an observing position in the plane of the sky.

cluster of 5–12 pc for the [C II] emission in the wings<fnref n="4" />. Based on the observed plane-of-the-sky (POS) radius this gives a distance d in the line of sight (LOS) between 0 and 10 pc, which allows us to estimate the dynamical timescales ($t_{dyn} = d_{LOS}/v_{wing}$). Using a typical $d_{LOS} = 5$ pc and the velocity range of the [C II] wings relative to the cloud for $v_{wing}$ then results in a full range of dynamical timescales of ≈0.2–1.0 Myr. This is significantly shorter than the observationally calculated cluster age of 2.3 ± 0.5 Myr (Martins et al. 2010) and an age estimate of 2.2 ± 0.1 Myr by Tremblin et al. (2014). This suggests that the high-velocity gas is continuously replenished as it flows away, and that [C II] thus only traces the gas that is currently being ejected from the region (<1 Myr old). Combined with the observed complex morphology of the high-velocity gas, we argue that the high-velocity gas is not a single expanding shell in RCW 79, in contrast to what has been observed for several other H II regions (Luisi et al. 2021; Pabst et al. 2019). We propose that the high-velocity [C II] emission also originates from gas that is continuously entrained at the interface of the H II region and the dense ring and turbulent molecular cloud. This entrained gas escapes from the molecular clouds through low-density holes in the turbulent cloud (see Fig. 5). In this context it is interesting to note that observations in the Orion cloud found that protostellar outflows can also create holes in an expanding shell (Kavak et al. 2022b,a). In addition to the turbulent structure of the cloud, protostellar outflows could thus also contribute at the early feedback stages in the creation of these holes. However, there are currently no clear indications of an important role for protostellar outflows in RCW 79. We also note that Zavagno et al. (2006) suggested that ionized gas from

the H II region is flowing outward based on their $H\alpha$ observations. The scenario for RCW 79 is conceptually similar to that of mass erosion from a ring of dense gas that was recently described in Whitworth et al. (2022).

## 4.2. Mass ejection rate and cloud erosion timescale

With the estimate of the dynamic timescales for the high-velocity wings and the calculated $C^+$ column density in the velocity channel of the wings we derive a mass ejection rate of $0.9–3.5 \times 10^{-2} M_{\odot} \text{yr}^{-1}$ associated with the [C II] wings (see Appendix C for details). This range of values is in agreement with the mass ejection rate ($\sim 10^{-2} M_{\odot} \text{yr}^{-1}$) found in the simulations by Walch et al. (2012, 2013), with cluster and cloud properties that closely match RCW 79. Interestingly, these simulations also contain a variety of expanding structures that fit the observed range of dynamical timescales of the observed mass ejection. From the Herschel dust column density map<fnref n="5" /> in Liu et al. (2017), we know that most of the dense gas is located in the swept-up ring at a radius of 6–7 pc with ongoing star formation. However, we consider a radius of 15 pc to include most of the mass of the molecular cloud, which gives a current total cloud mass of $2.1 \times 10^4 M_{\odot}$. Combined with the mass ejection rates, this gives a current cloud erosion timescale of 0.6–2.3 Myr. With the estimated cluster age of 2.3 ± 0.5 Myr, we thus obtain an erosion timescale for the original molecular cloud that is <5 Myr. Focusing only on the swept-up ring, within a radius of 8 pc, gives a current erosion timescale of 0.2–0.9 Myr (which adds up to a total erosion time of 2.5–3.2 Myr with the cluster age). This assumes that the high-velocity mass ejection from the cloud is continuous, which is supported by the presence of high-velocity wings in all the targets of the FEEDBACK survey, and in other H II regions (Bonne et al., in prep.). It is interesting to note that this direct estimate of the molecular cloud erosion timescale fits with numerical simulations of stellar feedback in turbulent clouds (Kim et al. 2018, 2021b) and indirect observational estimates of the molecular cloud erosion timescales in the Milky Way and nearby galaxies based on extinction measurements and the molecular gas distribution near OB clusters (e.g., Leisawitz & Hauser 1988; Leisawitz et al. 1989; Hannon et al. 2019; Kim et al. 2021a; Chevance et al. 2022).

## 5. Conclusion

We presented high spectral resolution [C II] observations of the H II region RCW 79 from the SOFIA legacy project FEEDBACK. The [C II] emission traces the PDRs at the edge of the H II region and a compact H II region in the southeast. Most importantly, the map reveals a large quantity of high-velocity gas. This gas displays an expanding shell morphology over the region along a few axes, where for the first time we identify both a blue- and a redshifted expanding [C II] shell. However, the overall kinematic structure of this high-velocity gas varies over the H II region with clear evidence of high-velocity gas escaping from the bubble through holes in the nascent cloud that gives rise to a complex dynamic morphology. These proposed holes can naturally emerge in a turbulent cloud, but it cannot be excluded that protostellar outflows might also play a role at earlier stages in their creation. From the velocities of

<footnote n="4">⁴ This fits with the observation of high-velocity wings in RCW 36. These observations show that the 3D extent of the high-velocity wings was within a factor two of the radius of the swept-up ring (Bonne et al. 2022).</footnote>

<footnote n="5">⁵ The Herschel dust column density map at 18'' angular resolution of the RCW 79 region (Liu et al. 2017) was produced from the HOBYS survey (Motte et al. 2010) using the method described in Palmeirim et al. (2013).</footnote>
the high-velocity gas we derive a range of dynamical timescales below 1.0 Myr. As these timescales are significantly shorter than the estimated cluster age of $2.3 \pm 0.5$ Myr (Martins et al. 2010), we propose that the [C II] wings only trace the most recent gas that is ablated from the dense ring and ejected from the turbulent cloud. Quantifying the associated mass ejection rate ($0.9-3.5 \times 10^{-2} M_{\odot} \text{yr}^{-1}$) implies relatively short estimated erosion timescales for the initial molecular cloud that are <5 Myr, which provides direct observational evidence for rapid molecular cloud evolution in this region, consistent with extragalactic studies. However, the current data does not allow us to draw a conclusion on the physical process that drives the molecular cloud erosion (i.e., whether stellar wind or radiation dominates). Additional observations will be required to address this question.

**Acknowledgements.** We thank the anonymous referee for constructive and insightful comments that improved the clarity of this Letter. We thank F. Comeron, D. Russeil and F. Martins for fruitful discussions on the determination of the new distance of RCW 79 and how it affects the deduced properties of the ionizing O stars. We also thank L. Lancaster for interesting discussions, and J-G. Kim for insightful comments on a draft version. This study was based on observations made with the NASA/DLR Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is jointly operated by the Universities Space Research Association Inc. (USRA), under NASA contract NNA17BF53C, and the Deutsches SOFIA Institut (DSI), under DLR contract 50 OK 0901 to the University of Stuttgart. upGREAT is a development by the MPIfR and the KOSMA/University Cologne, in cooperation with the DLR Institut für Optische Sensorsysteme. Financial support for FEEDBACK at the University of Maryland was provided by NASA through award SOF070077 issued by USRA. The FEEDBACK project is supported by the BMWI via DLR, project number 50 OR 2217 (FEEDBACK-plus). S.K. acknowledges support by the BMWI via DLR, project number 50OR2311 (Orion-Legacy). A.Z. acknowledges the support of the IUF.

## References

Beuther, H., Schneider, N., Simon, R., et al. 2022, A&A, 659, A77

Bonne, L., Schneider, N., Bontemps, S., et al. 2020, A&A, 641, A17

Bonne, L., Schneider, N., García, P., et al. 2022, ApJ, 935, 171

Bonne, L., Andersson, B. G., Minchin, R., et al. 2023a, AJ, 165, 243

Bonne, L., Bontemps, S., Schneider, N., et al. 2023b, ApJ, 951, 39

Chevance, M., Kruijssen, J. M. D., Krumholz, M. R., et al. 2022, MNRAS, 509, 272

Elmegreen, B. G. 2000, ApJ, 530, 277

Elmegreen, B. G., & Lada, C. J. 1977, ApJ, 214, 725

Evans, N. J. I., Dunham, M. M., Jørgensen, J. K., et al. 2009, ApJS, 181, 321

Ferland, G. J., Chatzikos, M., Guzmán, F., et al. 2017, Rev. Mex. Astron. Astrofis., 53, 385

Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1

Gao, Y., & Solomon, P. M. 2004, ApJ, 606, 271

Goicoechea, J. R., Teyssier, D., Etxaluze, M., et al. 2015, ApJ, 812, 75

Goldsmith, P. F., Langer, W. D., Pineda, J. L., & Velusamy, T. 2012, ApJS, 203, 13

Guan, X., Stutzki, J., Graf, U. U., et al. 2012, A&A, 542, L4

Güsten, R., Nyman, L. Å., Schilke, P., et al. 2006, A&A, 454, L13

Hannon, S., Lee, J. C., Whitmore, B. C., et al. 2019, MNRAS, 490, 4648

Hartmann, L., Ballesteros-Paredes, J., & Bergin, E. A. 2001, ApJ, 562, 852

Hollenbach, D. J., & Tielens, A. G. G. M. 1999, Rev. Mod. Phys., 71, 173

Kabanovic, S., Schneider, N., Ossenkopf-Okada, V., et al. 2022, A&A, 659, A36

Kavak, Ü., Bally, J., Goicoechea, J. R., et al. 2022a, A&A, 663, A117

Kavak, Ü., Goicoechea, J. R., Pabst, C. H. M., et al. 2022b, A&A, 660, A109

Kim, J.-G., Kim, W.-T., & Ostriker, E. C. 2018, ApJ, 859, 68

Kim, J., Chevance, M., Kruijssen, J. M. D., et al. 2021a, MNRAS, 504, 487

Kim, J.-G., Ostriker, E. C., & Filippova, N. 2021b, ApJ, 911, 128

Krumholz, M. R., & McKee, C. F. 2005, ApJ, 630, 250

Lada, C. J., Lombardi, M., & Alves, J. F. 2010, ApJ, 724, 687

Leisawitz, D., & Hauser, M. G. 1988, ApJ, 332, 954

Leisawitz, D., Bash, F. N., & Thaddeus, P. 1989, ApJS, 70, 731

Liu, H.-L., Figueira, M., Zavagno, A., et al. 2017, A&A, 602, A95

Luisi, M., Anderson, L. D., Schneider, N., et al. 2021, Sci. Adv., 7, eabe9511

Martins, F., Pomarès, M., Deharveng, L., Zavagno, A., & Bouret, J. C. 2010, A&A, 510, A32

Mège, P., Russeil, D., Zavagno, A., et al. 2021, A&A, 646, A74

Motte, F., Zavagno, A., Bontemps, S., et al. 2010, A&A, 518, L77

Pabst, C., Higgins, R., Goicoechea, J. R., et al. 2019, Nature, 565, 618

Pabst, C. H. M., Goicoechea, J. R., Teyssier, D., et al. 2020, A&A, 639, A2

Palmeirim, P., André, P., Kirk, J., et al. 2013, A&A, 550, A38

Pineda, J. L., Langer, W. D., Velusamy, T., & Goldsmith, P. F. 2013, A&A, 554, A103

Pineda, J. L., Langer, W. D., & Goldsmith, P. F. 2014, A&A, 570, A121

Pound, M. W., & Wolfire, M. G. 2023, AJ, 165, 25

Risacher, C., Güsten, R., Stutzki, J., et al. 2018, J. Astron. Instrum., 7, 1840014

Russeil, D., Georgelin, Y. M., Amram, P., et al. 1998, A&AS, 130, 119

Schneider, N., Röllig, M., Simon, R., et al. 2018, A&A, 617, A45

Schneider, N., Simon, R., Guevara, C., et al. 2020, PASP, 132, 104301

Schneider, N., Bonne, L., Bontemps, S., et al. 2023, Nat. Astron., 7, 546

Shu, F. H., Adams, F. C., & Lizano, S. 1987, ARA&A, 25, 23

Simon, R., Schneider, N., Stutzki, J., et al. 2012, A&A, 542, L12

Sofia, U. J., Lauroesch, J. T., Meyer, D. M., & Cartledge, S. I. B. 2004, ApJ, 605, 272

Spitzer, L. 1978, Physical Processes in the Interstellar Medium (New York: Wiley)

Stacey, G. J., Jaffe, D. T., Geis, N., et al. 1993, ApJ, 404, 219

Tiwari, M., Karim, R., Pound, M. W., et al. 2021, ApJ, 914, 117

Tram, L. N., Bonne, L., Hu, Y., et al. 2023, ApJ, 946, 8

Tremblin, P., Anderson, L. D., Didelon, P., et al. 2014, A&A, 568, A4

Walch, S. K., Whitworth, A. P., Bisbas, T., Wünsch, R., & Hubber, D. 2012, MNRAS, 427, 625

Walch, S., Whitworth, A. P., Bisbas, T. G., Wünsch, R., & Hubber, D. A. 2013, MNRAS, 435, 917

Weaver, R., McCray, R., Castor, J., Shapiro, P., & Moore, R. 1977, ApJ, 218, 377

Whitworth, A. P., Priestley, F. D., & Geen, S. T. 2022, MNRAS, 517, 4940

Wolfire, M. G., Vallini, L., & Chevance, M. 2022, ARA&A, 60, 247

Young, E. T., Becklin, E. E., Marcum, P. M., et al. 2012, ApJ, 749, L17

Zavagno, A., Deharveng, L., Comerón, F., et al. 2006, A&A, 446, 171
# Appendix A: Distance update of RCW 79 with Gaia DR3

With Gaia data release 3 (DR3; Gaia Collaboration 2023), we updated the distance of the RCW 79 ionizing cluster. From a list of all stars in the region of the RCW 79 cluster with parallax measurements, we selected stars that lie within 1" of the O stars from Martins et al. (2010) and have a distance within 3.2 and 5.2 kpc (which is the largest error bar for the region). This provided us with two candidate stars that correspond to an O star identified in Martins et al. (2010). These two stars have a consistent distance of 3.87±0.07 kpc and 3.93±0.06 kpc, a good astrometric solution, and only a velocity difference of 6 km s⁻¹ that might be associated with internal motion in the cluster. Furthermore, the magnitudes of the stars are consistent with reddened hot massive stars. Based on these values, and taking into account a possible parallax zero-point offset, we propose a new distance of 3.9±0.4 kpc. As the spectral types in Martins et al. (2010) were derived based on spectral classification, the new distance does not affect these results nor does it affect the upper limit on the stellar wind rates. Only the ionizing photon flux requires a 14% correction.

# Appendix B: Limb brightening in a spherical shell

The sketch in Fig. B.1 illustrates the limb brightening effect in a spherical shell. For the calculations we assume that the observed line is optically thin such that the observed intensity is given by $S_v\tau_v = S_v\kappa_v\Delta y = j_v\Delta y$, where $S_v$ is the specific intensity, $\tau_v$ the optical depth, $\kappa_v$ the attenuation coefficient, and $\Delta y$ the distance. At the center of the shell (x=0) there is no limb brightening, which results in an observed intensity of $j_v r$. At the position with peak limb brightening, the observed intensity is $j_v y_{peak}$. At this position, $y_{peak}$ can be described by

$$y_{peak} = R \sqrt{1 - \left(\frac{R-r}{R}\right)^2} \qquad (B.1)$$

using the Pythagorean theorem, where R is the radius of the shell and r the 3D thickness of the shell. The observed limb brightening is thus given by

$$\frac{I_{peak}}{I_{cent}} = \frac{R}{r} \sqrt{1 - \left(\frac{R-r}{R}\right)^2} \qquad (B.2)$$

In an isolated region, the observed shell width would be simply given by r as there is no emission outside the shell. However, shells are embedded in interstellar gas which means that the above definition is not always applicable. Therefore, we decided to define the width at half of the maximum intensity which is defined by $\Delta x_{WHM}$ in Fig. B.1. We do not employ the full width at half maximum (FWHM) because the projected shell intensity profile is not symmetric around the peak emission (see Fig. 14 in Kabanovic et al. 2022). The associated length through the shell is then given by $y_{WHM} = y_{peak}/2$ so that $\theta_{WHM}$ is defined by

$$\theta_{WHM} = \sin^{-1}\left(\frac{1}{2}\sqrt{1 - \left(\frac{R-r}{R}\right)^2}\right) \qquad (B.3)$$

On the other hand, we also have

$$R \cos\theta_{WHM} = R - r + \Delta x_{WHM} \qquad (B.4)$$

Using the relation

$$\cos(\sin^{-1}(\phi)) = \sqrt{1 - \phi^2}, \qquad (B.5)$$

<table><caption>Table B.1. Inferred shell widths for RCW 79</caption><thead><tr><th></th><th>r<sub>I<sub>peak</sub>/I<sub>cent</sub></sub></th><th>r<sub>Δx<sub>WHM</sub></sub></th></tr></thead><tbody><tr><td>Total</td><td>0.05-0.13 pc</td><td>0.5-1.1 pc</td></tr><tr><td>Blue wing</td><td>0.08-0.14 pc</td><td>0.7-0.8 pc</td></tr><tr><td>Red wing</td><td>0.26-0.93 pc</td><td>0.7-1.1 pc</td></tr></tbody></table>

Note: The inferred width of the shell (r) in [C II] based on the limb brightening ($r_{I_{peak}/I_{cent}}$) and the observed half width of the shell ($r_{\Delta x_{WHM}}$). The results based on the two methods are clearly inconsistent.

<figure><img src="image_7.png" /></figure>

Fig. B.1. Spherical shell or bubble observed along the y-axis, as indicated by the telescope illustration. The shell (in blue) surrounding a bubble (white) is shown with a constant specific intensity (S), radius R, and thickness r. The maximum observed limb brightening is at x = R-r, and is indicated by $y_{peak}$. The observed width at half the intensity in the shell ($y_{WHM}$) is indicated with $\Delta x_{WHM}$.

we then obtain

$$R - r + \Delta x_{WHM} = R \sqrt{\frac{3}{4} + \frac{1}{4} \left(\frac{R-r}{R}\right)^2}, \qquad (B.6)$$

which results in

$$\frac{\Delta x_{WHM}}{R} = \frac{r}{R} - 1 + \sqrt{\frac{3}{4} + \frac{1}{4} \left(\frac{R-r}{R}\right)^2} \qquad (B.7)$$

These equations are also employed in the analysis of Bonne et al. (2023a). The observed $I_{peak}/I_{cent}$ and $\Delta x/R$ thus provide two independent methods for calculating r, using Brent's method,<fnref n="6" /> and for evaluating whether the results are consistent (which is expected for a spherical expanding shell). To this end, we created velocity-integrated intensity cuts over the center of the RCW 79 H II region for the blue- and redshifted wing and the total velocity range (-25 km s⁻¹ to -70 km s⁻¹); this was done in intervals of 45°. To make sure we were not working with upper limits, we only worked with the cases where $I_{cent} \ge 3\sigma$. To determine $\Delta x_{WHM}$ we focused on the pixels with a higher radius than the peak intensity to find the first pixel with an observed intensity <$I_{peak}/2$. With linear interpolation between this and the previous pixel, we determined $\Delta x_{WHM}$. Based on these observed values

<footnote n="6">⁶ https://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.brentq.html</footnote>
<figure><img src="image_8.png" /><figcaption>**Fig. B.2.** Limb brightening factor and projected shell FWHM assuming a bubble geometry with a radius of 6.5 pc for the total intensity of RCW 79, based on the formalism in Kabanovic et al. (2022). To be a self-consistent spherical bubble, the deduced geometric shell width from the limb brightening and the projected shell FWHM must be the same within the error, which is not the case. This shows that the emission distribution of RCW 79 is not consistent with a spherical shell morphology. In this plot the results for the highest projected shell FWHM (1.8 pc and 2.3 pc) were not included as the FWHM formalism does not allow a solution for these values.</figcaption></figure>

we then inferred the actual shell width (r). The results are summarized in Table B.1; they show that the resulting shell width (r) from the limb brightening and $\Delta x_{WHM}$ measurements are inconsistent. To explain the observed limb brightening, a significantly smaller width of the shell would be required than the observed width of the shell. We also note that this should not be a spatial resolution effect since our spatial resolution of ~0.4 pc is smaller than the observed half width of the shell ($\Delta x_{WHM}$).

The model we presented here is based on three assumptions: the emission is optically thin, the shell has a constant local emissivity ($j_v$), and the shell is spherical. Optical depth effects for the limb brightening would give results opposite to the observations. Therefore, the high-velocity emission has to be highly inhomogeneous, and thus is likely non-spherical. Based on the very bright [C II] emission at the edge of the H II region, an elegant solution would be that the high-velocity gas is predominantly gas ablated from the edge of the ring. This appears to fit with the observed kinematics, and would also provide a solution for the timescale issue.

Recently, the same model was also discussed in Kabanovic et al. (2022) in the context of RCW 120. We independently confirmed our result for RCW 79 with the formalism for the geometry of a spherically expanding bubble presented in Kabanovic et al. (2022). Figure B.2 shows the results from this independent analysis and demonstrates again that the observed limb brightening in RCW 79 is inconsistent with the observed shell width. Although in essence the model of Kabanovic et al. (2022) is the same as the one above, it builds on the observational definition of a FWHM (instead of WHM in the formalism above) for the observed shell width. Applying this formalism thus verifies that the results above are not the result of a bias in our observational definition.

## Appendix C: Column density in the [C II] high-velocity wings

The C$^+$ column density can be estimated in each velocity channel using the equation from Goldsmith et al. (2012)

$$ \Delta T_A = 3.43 \times 10^{-16} \times \left[ 1 + 0.5 \times e^{91.25/T_{kin}} \left( 1 + \frac{2.4 \times 10^{-6}}{C_{ul}} \right) \right]^{-1} \frac{N(C^+)}{\delta v} \quad (C.1) $$

Here $T_{kin}$ is the kinetic temperature, N(C$^+$) the C$^+$ column density, $\delta v$ the spectral bin width, and $C_{ul}$ the collisional de-excitation rate given by

$$ C_{ul} = n \times R_{ul}, \quad (C.2) $$

where n is the density and $R_{ul}$ is the de-excitation rate coefficient for atomic hydrogen given by

$$ R_{ul} = 7.6 \times 10^{-10} \text{ cm}^3 \text{ s}^{-1} \left( \frac{T_{kin}}{100 \text{ K}} \right)^{0.14}. \quad (C.3) $$

Based on the results of the PDR Toolbox for the FUV field strength, we use $T_{kin}$ values between 50 and 250 K. The quoted mass, momentum, and energy ejection intervals in the paper are based on this temperature range. As there is no direct estimate of the density, we assume $n_H = 10^3 \text{ cm}^{-3}$. This value is reasonable seeing $^{12}$CO is basically undetected in the high-velocity wings (but does show some hints of localized clumps), and it is close to the predicted density for the PDR assuming thermal equilibrium with the ionized gas phase. Using lower density values would only lead to higher mass, momentum, and energy ejection rates. To convert the C$^+$ column density to hydrogen column density we used [C$^+$/[H] = $1.6 \times 10^{-4}$ from Sofia et al. (2004) and a mean atomic mass of $2.2 \times 10^{-27}$ kg. Based on the column density maps, we can calculate the total mass associated with the blue- and redshifted high-velocity wings. For temperatures between 50 and 250 K, this results in masses between $5.9 \times 10^3 \text{ M}_\odot$ and $1.5 \times 10^3 \text{ M}_\odot$, respectively, for the blueshifted gas and $6.5 \times 10^3 \text{ M}_\odot$ and $1.6 \times 10^3 \text{ M}_\odot$, respectively, for the redshifted gas.

To estimate a timescale, which is necessary to calculate the mass, momentum, and energy ejection rates, we work with the observed LOS velocity relative to the bulk motion of the cloud and the estimated typical distance in the LOS (i.e., 5 pc). Lastly, it should be noted that the calculated mass, momentum, and energy ejection rates have several factors of uncertainty, for example projection uncertainties, density uncertainty, and the assumed [C$^+$/[H] value. This makes it difficult to estimate the error bars, but suggests that at least a factor two should be considered for the uncertainty. As the calculations are plagued by many uncertainties, future work will explore numerical simulations to estimate correction factors for potential biases.

Lastly, we note that we verified whether the high-velocity [C II] emission could originate from the H II region instead of the PDR. To do this, we ran several constant pressure CLOUDY models (Ferland et al. 2017) for a H II region resembling RCW 79, using the O-star luminosity from Martins et al. (2010), a blackbody radiation spectrum of $3.9 \times 10^4 \text{ K}$, radii for the H II gas between 1 and 5 pc, and Orion nebula abundances. This showed that <14% of the total [C II] emission originates from the H II region itself. This implies that <30% of the emission in the high-velocity [C II] wings can be attributed to the H II region. [C II] excitation in the ionized gas phase can give rise to slightly lower masses in the high-velocity wings per unit of emission (by a factor 2-3). However, as the [C II] emission coming from the ionized gas phase is limited, this would lead to a correction of <15%, which is well below the other uncertainties that go into the calculation with currently available data.