Search

Context. Recent studies of the optical depth comparing [12CII] and [13CII] line profiles in Galactic star-forming regions revealed strong self-absorption in [12CII] by low excitation foreground material, implying a large column density of C+ corresponding to an equivalent AV of a few, up to about 10 mag. Aims. As the nature and origin of such a large column of cold C+ foreground gas are difficult to explain, it is essential to constrain the physical conditions of this material. Methods. We conducted high-resolution observations of [OI] 63 um and [OI] 145 um lines in M17 SW and Mon R2. The [OI] 145 um transition traces warm PDR-material, while the [OI] 63 um line traces foreground material as manifested by absorption dips. Results. Comparison of both [OI] line profiles with [CII] isotopic lines confirms warm PDR-origin background emission and a significant column of cold foreground material causing self-absorption visible in [12CII] and [OI] 63 um profiles. In M17 SW, the C+ and O column densities are comparable for both layers. Mon R2 exhibits larger O columns compared to C+, indicating additional material where the carbon is neutral or in molecular form. Small-scale spatial variation of the foreground absorption profiles and the large column density (around 1E18 cm-2 ) of the foreground material suggest emission from high-density regions associated with the cloud complex, not a uniform diffuse foreground cloud. Conclusions. The analysis confirms that the previously detected intense [CII] foreground absorption is attributable to a large column of low excitation dense atomic material, where carbon is ionized, and oxygen is in neutral atomic form., Comment: 30 pages, 24 figures

Using the upGREAT instrument on SOFIA, we have imaged [C II] 157.74 and [O I] 63.18 micron line emission from a bright photodissociation region (PDR) associated with an ionized ``bubble'' located in the Nessie Nebula, a filamentary infrared dark cloud. A comparison with ATCA data reveals a classic PDR structure, with a uniform progression from ionized gas, to photodissociated gas, and on to molecular gas from the bubble's interior to its exterior. [O I] line emission from the bubble's PDR reveals self-absorption features. Toward a FIR-bright protostar, both [O I] and [C II] show an absorption feature at a velocity of $-18$ km/s, the same velocity as an unrelated foreground molecular cloud. Since the gas density in typical molecular clouds is well below the [O I] and [C II] critical densities, the excitation temperatures for both lines are low (~20 K). The Meudon models demonstrate that the surface of a molecular cloud, externally illuminated by a standard G_0 = 1 interstellar radiation field, can produce absorption features in both transitions. Thus, the commonly observed [O I] and [C II] self-absorption and absorption features plausibly arise from the subthermally excited, externally illuminated, photodissociated envelopes of molecular clouds. The luminous young stellar object AGAL337.916-00.477, located precisely where the expanding bubble strikes the Nessie filament, is associated with two shock tracers: NH3 (3,3) maser emission and SiO 2-1 emission, indicating interaction between the bubble and the filament. The interaction of the expanding bubble with its parental dense filament has triggered star formation.

We explore the potential of the Gaussian Mixture Model (GMM), an unsupervised machine learning method, to identify coherent physical structures in the ISM. The implementation we present can be used on any kind of spatially and spectrally resolved data set. We provide a step-by-step guide to use these models on different sources and data sets. Following the guide, we run the models on NGC 1977, RCW 120 and RCW 49 using the [CII] 158 $\mu$m mapping observations from the SOFIA telescope. We find that the models identified 6, 4 and 5 velocity coherent physical structures in NGC 1977, RCW 120 and RCW 49, respectively, which are validated by analysing the observed spectra towards these structures and by comparison to earlier findings. In this work we demonstrate that GMM is a powerful tool that can better automate the process of spatial and spectral analysis to interpret mapping observations., Comment: 19 pages, 14 figures

We quantified the effects of stellar feedback in RCW 49 by determining the physical conditions in different regions using the [CII] 158 $\mu$m and [OI] 63 $\mu$m observations from SOFIA, the $^{12}$CO (3-2) observations from APEX and the H$_2$ line observations from Spitzer telescopes. Large maps of RCW 49 were observed with the SOFIA and APEX telescopes, while the Spitzer observations were only available towards three small areas. From our qualitative analysis, we found that the H$_2$ 0-0 S(2) emission line probes denser gas compared to the H$_2$ 0-0 S(1) line. In four regions ("northern cloud", "pillar", "ridge", and "shell"), we compared our observations with the updated PDR Toolbox models and derived the integrated far-ultraviolet flux between 6-13.6 eV ($G_{\rm 0}$), H nucleus density ($n$), temperatures and pressures. We found the ridge to have the highest $G_{\rm 0}$ (2.4 $\times$ 10$^3$ Habing units), while the northern cloud has the lowest $G_{\rm 0}$ (5 $\times$ 10$^2$ Habing units). This is a direct consequence of the location of these regions with respect to the Wd2 cluster. The ridge also has a high density (6.4 $\times$ 10$^3$ cm$^{-3}$), which is consistent with its ongoing star formation. Among the Spitzer positions, we found the one closest to the Wd2 cluster to be the densest, suggesting an early phase of star formation. Furthermore, the Spitzer position that overlaps with the shell was found to have the highest $G_{\rm 0}$ and we expect this to be a result of its proximity to an O9V star., Comment: 18 pages, 12 figures

We present [CII] 158 $\mu$m and [OI] 63 $\mu$m observations of the bipolar HII region RCW 36 in the Vela C molecular cloud, obtained within the SOFIA legacy project FEEDBACK, which is complemented with APEX $^{12/13}$CO(3-2) and Chandra X-ray (0.5-7 keV) data. This shows that the molecular ring, forming the waist of the bipolar nebula, expands with a velocity of 1 - 1.9 km s$^{-1}$. We also observe an increased linewidth in the ring indicating that turbulence is driven by energy injection from the stellar feedback. The bipolar cavity hosts blue-shifted expanding [CII] shells at 5.2$\pm$0.5$\pm$0.5 km s$^{-1}$ (statistical and systematic uncertainty) which indicates that expansion out of the dense gas happens non-uniformly and that the observed bipolar phase might be relatively short ($\sim$0.2 Myr). The X-ray observations show diffuse emission that traces a hot plasma, created by stellar winds, in and around RCW 36. At least 50 \% of the stellar wind energy is missing in RCW 36. This is likely due to leakage which is clearing even larger cavities around the bipolar RCW 36 region. Lastly, the cavities host high-velocity wings in [CII] which indicates relatively high mass ejection rates ($\sim$5$\times$10$^{-4}$ M$_{\odot}$ yr$^{-1}$). This could be driven by stellar winds and/or radiation pressure, but remains difficult to constrain. This local mass ejection, which can remove all mass within 1 pc of RCW 36 in 1-2 Myr, and the large-scale clearing of ambient gas in the Vela C cloud indicates that stellar feedback plays a significant role in suppressing the star formation efficiency (SFE)., Comment: 38 pages, 27 figures, 8 tables, accepted in ApJ

Crawford's et al. (2021) article on estimation of discrete choice models with unobserved or latent consideration sets, presents a unified framework to address the problem in practice by using "sufficient sets", defined as a combination of past observed choices. The proposed approach is sustained in a re-interpretation of a consistency result by McFadden (1978) for the problem of sampling of alternatives, but the usage of that result in Crawford et al. (2021) is imprecise in an important matter. It is stated that consistency would be attained if any subset of the true consideration set is used for estimation, but McFadden (1978) shows that, in general, one needs to do a sampling correction that depends on the protocol used to draw the choice set. This note derives the sampling correction that is required when the choice set for estimation is built from past choices. Then, it formalizes the conditions under which such correction would fulfill the uniform condition property and can therefore be ignored when building practical estimators, such as the ones analyzed by Crawford et al. (2021).

Context: How do expanding HII regions interact with their environmental cloud? This is one of the central questions driving the SOFIA legacy program FEEDBACK. Here, we present a case study toward the prototypical H{\sc ii} region NGC7538. Methods: With SOFIA we mapped an area of ~210'^2 around NGC7538 in the [CII] line at 1.9THz. Complementary observed atomic carbon [CI] and high-J CO(8-7) data as well as archival NIR/FIR, cm continuum, CO(3-2) and HI data are folded into the analysis. Results: While the overall [CII] morphology follows the general ionized gas, the channel maps show multiple bubble-like structures with sizes on the order of ~80-100" (~1.0-1.28pc). While at least one of them may be an individual feedback bubble driven by the main exciting sources of the region, the other bubble-morphologies may also be due to the intrinsically porous structure of the HII region. An analysis of the expansion velocities around 10km s^{-1} indicates that thermal expansion is not sufficient but that wind-driving from the central O-stars is required. The most blue-shifted [CII] component has barely any molecular or atomic counterparts. At the interface to the molecular cloud, we find a typical photon-dominated region (PDR) with a bar-shape. Ionized, atomic and molecular carbon show a layered structure in this PDR. The carbon in the PDR is dominated by its ionized form with atomic and molecular masses of ~0.45+-0.1M_{\odot} and ~1.2+-0.1M_{\odot}, respectively, compared to the ionized carbon in the range of 3.6-9.7M_{\odot}. Conclusions: The NGC7538 HII region exhibits a diverse set of sub-structures that interact with each other as well as with the adjacent cloud. Compared to other recent [CII] observations of HII regions (e.g., Orion Veil, RCW120, RCW49), bubble-shape morphologies revealed in [CII] emission, indicative of expanding shells, are recurring structures of PDRs., Comment: 21 pages, 20 figures, accepted for Astronomy & Astrophysics, a high-resolution version can be found at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

Revealing the 3D dynamics of HII regions and their associated molecular clouds is important for understanding the longstanding problem as to how stellar feedback affects the density structure and kinematics of the interstellar medium. We employed observations of the HII region RCW 120 in [CII], observed within the SOFIA legacy program FEEDBACK, and the $^{12}$CO and $^{13}$CO (3$\to$2) lines, obtained with APEX. In addition we used HI data from the Southern Galactic Plane Survey. Two radiative transfer models were used to fit the observed data. A line profile analysis with the 1D non-LTE radiative transfer code SimLine proves that the CO emission cannot stem from a spherically symmetric molecular cloud configuration. With a two-layer multicomponent model, we then quantified the amount of warm background and cold foreground gas. There is a deficit of CO emission along the line-of-sight toward the center of the HII region which indicates that the HII region is associated with a flattened molecular cloud. Self-absorption in the CO line may hide signatures of infalling and expanding molecular gas. The [CII] emission arises from an expanding [CII] bubble and from the PDRs. A significant part of [CII] emission is absorbed in a cool (~60-100 K), low-density (<500 cm$^{-3}$) atomic foreground layer with a thickness of a few parsec. We propose that the RCW 120 HII region formed in a flattened molecular cloud and is now bursting out of its parental cloud. The compressed surrounding molecular layer formed a torus around the spherically expanding HII bubble. This scenario can possibly be generalized for other HII bubbles and would explain the observed "flat" structure of molecular clouds associated with HII bubbles. We suggest that the [CII] absorption observed in many star-forming regions is at least partly caused by low-density, cool, HI-envelopes surrounding the molecular clouds.

We unveil the stellar wind driven shell of the luminous massive star-forming region of RCW 49 using SOFIA FEEDBACK observations of the [CII] 158 $\mu$m line. The complementary dataset of the $^{12}$CO and $^{13}$CO J = 3 - 2 transitions is observed by the APEX telescope and probes the dense gas toward RCW 49. Using the spatial and spectral resolution provided by the SOFIA and APEX telescopes, we disentangle the shell from a complex set of individual components of gas centered around RCW 49. We find that the shell of radius ~ 6 pc is expanding at a velocity of 13 km s$^{-1}$ toward the observer. Comparing our observed data with the ancillary data at X-Ray, infrared, sub-millimeter and radio wavelengths, we investigate the morphology of the region. The shell has a well defined eastern arc, while the western side is blown open and is venting plasma further into the west. Though the stellar cluster, which is ~ 2 Myr old gave rise to the shell, it only gained momentum relatively recently as we calculate the shell's expansion lifetime ~ 0.27 Myr, making the Wolf-Rayet star WR20a a likely candidate responsible for the shell's re-acceleration., Comment: 31 pages, 17 figures

FEEDBACK is a SOFIA legacy program dedicated to study the interaction of massive stars with their environment. It performs a survey of 11 galactic high mass star forming regions in the 158 $\mu$m (1.9 THz) line of CII and the 63 $\mu$m (4.7 THz) line of OI. We employ the 14 pixel LFA and 7 pixel HFA upGREAT instrument to spectrally resolve (0.24 MHz) these FIR structure lines. With an observing time of 96h, we will cover $\sim$6700 arcmin$^2$ at 14.1$''$ angular resolution for the CII line and 6.3$''$ for the OI line. The observations started in spring 2019 (Cycle 7). Our aim is to understand the dynamics in regions dominated by different feedback processes from massive stars such as stellar winds, thermal expansion, and radiation pressure, and to quantify the mechanical energy injection and radiative heating efficiency. The CII line provides the kinematics of the gas and is one of the dominant cooling lines of gas for low to moderate densities and UV fields. The OI line traces warm and high-density gas, excited in photodissociations regions with a strong UV field or by shocks. The source sample spans a broad range in stellar characteristics from single OB stars, to small groups of O stars, to rich young stellar clusters, to ministarburst complexes. It contains well-known targets such as Aquila, the Cygnus X region, M16, M17, NGC7538, NGC6334, Vela, and W43 as well as a selection of HII region bubbles, namely RCW49, RCW79, and RCW120. These CII maps, together with the less explored OI 63 $\mu$m line, provide an outstanding database for the community. They will be made publically available and will trigger further studies and follow-up observations.

Context. The [CII] 158 {\mu}m far-infrared (FIR) fine-structure line is one of the most important cooling lines of the star-forming interstellar medium (ISM). High spectral resolution observations have shown complex structures in the line profiles of the [CII] emission. Aims. Our aim is to determine whether the complex profiles observed in [^{12}CII] are due to individual velocity components along the line-of-sight or to self-absorption based on a comparison of the [^{12}CII] and isotopic [^{13}CII] line profiles. Methods. Deep integrations with the SOFIA/upGREAT 7-pixel array receiver in M43, Horsehead~PDR, Monoceros~R2, and M17~SW allow for the detection of optically thin [^{13}CII] emission lines, along with the [^{12}CII] emission lines, with a high signal-to-noise ratio. We first derived the [^{12}CII] optical depth and the [CII] column density from a single component model. However, the complex line profiles observed require a double layer model with an emitting background and an absorbing foreground. A multi-component velocity fit allows us to derive the physical conditions of the [CII] gas: column density and excitation temperature. Results. We find moderate to high [^{12}CII] optical depths in all four sources and self-absorption of [^{12}CII] in Mon R2 and M17 SW. The high column density of the warm background emission corresponds to an equivalent Av of up to 41 mag. The foreground absorption requires substantial column densities of cold and dense [CII] gas, with an equivalent Av ranging up to about 13 mag. Conclusions. The column density of the warm background material requires multiple photon-dominated region (PDR) surfaces stacked along the line of sight and in velocity. The substantial column density of dense and cold foreground [CII] gas detected in absorption cannot be explained with any known scenario and we can only speculate on its origins, Comment: 40 pages, 49 figures

We probe the chemical and energetic conditions in dense gas created by radiative feedback through observations of multiple CO, HCN and HCO$^+$ transitions toward the dense core of M17 SW. We used the dual band receiver GREAT on board the SOFIA airborne telescope to obtain maps of the $J=16-15$, $J=12-11$, and $J=11-10$ transitions of $^{12}$CO. We compare these maps with corresponding APEX and IRAM 30m telescope data for low- and mid-$J$ CO, HCN and HCO$^+$ emission lines, including maps of the HCN $J=8-7$ and HCO$^+$ $J=9-8$ transitions. The excitation conditions of $^{12}$CO, HCO$^+$ and HCN are estimated with a two-phase non-LTE radiative transfer model of the line spectral energy distributions (LSEDs) at four selected positions. The energy balance at these positions is also studied. We obtained extensive LSEDs for the CO, HCN and HCO$^+$ molecules toward M17 SW. The LSED shape, particularly the high-$J$ tail of the CO lines observed with SOFIA/GREAT, is distinctive for the underlying excitation conditions. The critical magnetic field criterion implies that the cold cloudlets at two positions are partially controlled by processes that create and dissipate internal motions. Supersonic but sub-Alfv\'enic velocities in the cold component at most selected positions indicates that internal motions are likely MHD waves. Magnetic pressure dominates thermal pressure in both gas components at all selected positions, assuming random orientation of the magnetic field. The magnetic pressure of a constant magnetic field throughout all the gas phases can support the total internal pressure of the cold components, but it cannot support the internal pressure of the warm components. If the magnetic field scales as $B \propto n^{2/3}$, then the evolution of the cold cloudlets at two selected positions, and the warm cloudlets at all selected positions, will be determined by ambipolar diffusion., Comment: 26 pages, 13 figures, A&A accepted