Search

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 doing exactly that toward the RCW 79 HII region using the ionized carbon line at 158 $\mu$m ([CII]) 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 [CII] for the first time both blue- and red-shifted mostly neutral high-velocity gas which reaches velocities up to 25 km s$^{-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 HII 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 which is well below the 2.3$\pm$0.5 Myr age of the OB cluster. This suggests a different scenario for the evolution of RCW 79, where the high-velocity gas is not solely stemming 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$\times$10$^{-2}$ M$_{\odot}$~yr$^{-1}$, which leads to short erosion timescales, i.e. $<$5 Myr, for the nascent molecular cloud. This finding provides direct observational evidence of rapid molecular cloud dispersal., Comment: 9 pages, 7 figures, accepted in A&A Letters

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

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.

IRAS 20319+3958 in Cygnus X South is a rare example of a free-floating globule (mass ~240 Msun, length ~1.5 pc) with an internal HII region created by the stellar feedback of embedded intermediate-mass stars, in particular, one Herbig Be star. Here, we present a Herschel/HIFI CII 158 mu map of the whole globule and a large set of other FIR lines (mid-to high-J CO lines observed with Herschel/PACS and SPIRE, the OI 63 mu line and the CO 16-15 line observed with upGREAT on SOFIA), covering the globule head and partly a position in the tail. The CII map revealed that the whole globule is probably rotating. Highly collimated, high-velocity CII emission is detected close to the Herbig Be star. We performed a PDR analysis using the KOSMA-tau PDR code for one position in the head and one in the tail. The observed FIR lines in the head can be reproduced with a two-component model: an extended, non-clumpy outer PDR shell and a clumpy, dense, and thin inner PDR layer, representing the interface between the HII region cavity and the external PDR. The modelled internal UV field of ~2500 Go is similar to what we obtained from the Herschel FIR fluxes, but lower than what we estimated from the census of the embedded stars. External illumination from the ~30 pc distant Cyg OB2 cluster, producing an UV field of ~150-600 G0 as an upper limit, is responsible for most of the CII emission. For the tail, we modelled the emission with a non-clumpy component, exposed to a UV-field of around 140 Go., Comment: accepted by Astronomy & Astrophysics

Context. The [CII] 158micron far-infrared fine-structure line is one of the dominant cooling lines of the star-forming interstellar medium (ISM). Hence [CII] emission originates in and thus can be used to trace a range of ISM processes. Velocity-resolved large-scale mapping of [CII] in star-forming regions provides a unique perspective of the kinematics of these regions and their interactions with the exciting source of radiation. Aims. We explore the scientific applications of large-scale mapping of velocity-resolved [CII] observations. With the [CII] observations, we investigate the effect of stellar feedback on the ISM. We present the details of observation, calibration, and data reduction using a heterodyne array receiver mounted on an airborne observatory. Results. A square-degree [CII] map with a spectral resolution of 0.3 km/s is presented. The scientific potential of this data is summarized with discussion of mechanical and radiative stellar feedback, filament tracing using [CII], [CII] opacity effects, [CII] and carbon recombination lines, and [CII] interaction with the large molecular cloud. The data quality and calibration is discussed in detail, and new techniques are presented to mitigate the effects of unavoidable instrument deficiencies (e.g. baseline stability) and thus to improve the data quality. A comparison with a smaller [CII] map taken with the Herschel/Heterodyne Instrument for the Far-Infrared (HIFI) spectrometer is presented.

Ionized interstellar gas is an important component of the interstellar medium and its lifecycle. The recent evidence for a widely distributed highly ionized warm interstellar gas with a density intermediate between the warm ionized medium (WIM) and compact HII regions suggests that there is a major gap in our understanding of the interstellar gas. Here we investigate the properties of the dense warm ionized medium (D-WIM) in the Milky Way using spectrally resolved SOFIA GREAT [NII] 205 micron line emission and Green Bank Telescope hydrogen radio recombination lines (RRL) data, supplemented by Herschel PACS [NII] 122 micron data, and spectrally resolved 12CO. We observed eight lines of sight in the 20deg 50% of the observed [CII] intensity along these LOS. The kinetic temperatures we derive are too low to explain the presence of N+ resulting from electron collisional ionization and/or proton charge transfer of atomic nitrogen. Rather, these regions most likely are ionized by extreme ultraviolet radiation., Comment: 16 pages, 10 figures

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.

Observations with the Herschel Space Telescope have established that most of the star forming gas is organised in interstellar filaments, a finding that is supported by numerical simulations of the supersonic interstellar medium (ISM) where dense filamentary structures are ubiquitous. We aim to understand the formation of these dense structures by performing observations covering the $^{12}$CO(4-3), $^{12}$CO(3-2), and various CO(2-1) isotopologue lines of the Musca filament, using the APEX telescope. The observed CO intensities and line ratios cannot be explained by PDR (photodissociation region) emission because of the low ambient far-UV field that is strongly constrained by the non-detections of the [C II] line at 158 $\mu$m and the [O I] line at 63 $\mu$m, observed with the upGREAT receiver on SOFIA, as well as a weak [C I] 609 $\mu$m line detected with APEX. We propose that the observations are consistent with a scenario in which shock excitation gives rise to warm and dense gas close to the highest column density regions in the Musca filament. Using shock models, we find that the CO observations can be consistent with excitation by J-type low-velocity shocks. A qualitative comparison of the observed CO spectra with synthetic observations of dynamic filament formation simulations shows a good agreement with the signature of a filament accretion shock that forms a cold and dense filament from a converging flow. The Musca filament is thus found to be dense molecular post-shock gas. Filament accretion shocks that dissipate the supersonic kinetic energy of converging flows in the ISM may thus play a prominent role in the evolution of cold and dense filamentary structures., Comment: 16 pages, 17 figures

The aim of the present work is to study structure and gas kinematics in the photodissociation regions (PDRs) around the compact HII regions S235A and S235C. We observe the [CII], [13CII] and [OI] line emission, using SOFIA/upGREAT and complement them by data of HCO+ and CO. We use the [13CII] line to measure the optical depth of the [CII] emission, and find that the [CII] line profiles are influenced by self-absorption, while the [13CII] line remains unaffected by these effects. Hence, for dense PDRs, [13CII] emission is a better tracer of gas kinematics. The optical depth of the [CII] line is up to 10 in S235A. We find an expanding motion of the [CII]-emitting layer of the PDRs into the front molecular layer in both regions. Comparison of the gas and dust columns shows that gas components visible neither in the [CII] nor in low-J CO lines may contribute to the total column across S235A. We test whether the observed properties of the PDRs match the predictions of spherical models of expanding HII region + PDR + molecular cloud. Integrated intensities of the [13CII], [CII] and [OI] lines are well-represented by the model, but the models do not reproduce the double-peaked [CII] line profiles due to an insufficient column density of C+. The model predicts that the [OI] line could be a more reliable tracer of gas kinematics, but the foreground self-absorbing material does not allow using it in the considered regions., Comment: accepted to MNRAS

Strong winds and ultraviolet (UV) radiation from O-type stars disrupt and ionize their molecular core birthplaces, sweeping up material into parsec-size shells. Owing to dissociation by starlight, the thinnest shells are expected to host low molecular abundances and therefore little star formation. Here, we expand previous maps taken with the IRAM 30m telescope and present square-degree 12CO and 13CO (J=2-1) maps of the wind-driven "Veil bubble'' that surrounds the Trapezium cluster and its natal Orion molecular core (OMC). Although widespread and extended CO emission is largely absent from the Veil, we show that several CO "globules'' exist and are embedded in the [CII]158um-bright shell that confines the bubble. This includes the first detection of quiescent CO at negative LSR velocities in Orion. Given the harsh UV irradiation conditions in this translucent material, the detection of CO globules is surprising. These globules are small (R=7,100 AU), not massive (M=0.3M_Sun), and are moderately dense: n_ H=4x10^4 cm^-3 (median values). They are confined by the external pressure of the shell, P_ext/k~10^7 cm^-3 K, and are likely magnetically supported. They are either transient objects formed by instabilities or have detached from pre-existing molecular structures, sculpted by the passing shock associated with the expanding shell and by UV radiation from the Trapezium. Some represent the first stages in the formation of small pillars, others of isolated small globules. Although their masses do not suggest they will form stars, one globule matches the position of a known YSO. The lack of extended CO in the "Veil shell'' demonstrates that feedback from massive stars expels, agitates, and reprocesses most of the disrupted molecular cloud gas, thereby limiting the star-formation rate in the region. The presence of globules is a result of this feedback., Comment: Accepted for publication in A&A. (Abridged abstract. English edited in this v2. Some figures have been bitmapped to lower resolution)

We present the performance of the upGREAT heterodyne array receivers on the SOFIA telescope after several years of operations. This instrument is a multi-pixel high resolution (R > 10^7) spectrometer for the Stratospheric Observatory for Far-Infrared Astronomy (SOFIA). The receivers use 7-pixel subarrays configured in a hexagonal layout around a central pixel. The low frequency array receiver (LFA) has 2x7 pixels (dual polarization), and presently covers the 1.83-2.06 THz frequency range, which allows to observe the [CII] and [OI] lines at 158 um and 145 um wavelengths. The high frequency array (HFA) covers the [OI] line at 63 um and is equipped with one polarization at the moment (7 pixels, which can be upgraded in the near future with a second polarization array). The 4.7 THz array has successfully flown using two separate quantum-cascade laser local oscillators from two different groups. NASA completed the development, integration and testing of a dual-channel closed-cycle cryocooler system, with two independently operable He compressors, aboard SOFIA in early 2017 and since then, both arrays can be operated in parallel using a frequency separating dichroic mirror. This configuration is now the prime GREAT configuration and has been added to SOFIA's instrument suite since observing cycle 6., Comment: Accepted to the Journal of Astronomical Instrumentation (SOFIA Special Edition) on 12th November 2018

We present a new multi-pixel high resolution (R >10^7) spectrometer for the Stratospheric Observatory for Far-Infrared Astronomy (SOFIA). The receiver uses 2 x 7-pixel subarrays in orthogonal polarization, each in an hexagonal array around a central pixel. We present the first results for this new instrument after commissioning campaigns in May and December 2015 and after science observations performed in May 2016 . The receiver is designed to ultimately cover the full 1.8-2.5 THz frequency range but in its first implementation, the observing range was limited to observations of the [CII] line at 1.9 THz in 2015 and extended to 1.83-2.07 THz in 2016. The instrument sensitivities are state-of-the-art and the first scientific observations performed shortly after the commissioning confirm that the time efficiency for large scale imaging is improved by more than an order of magnitude as compared to single pixel receivers. An example of large scale mapping around the Horsehead Nebula is presented here illustrating this improvement. The array has been added to SOFIA's instrument suite already for ongoing observing cycle 4., Comment: 7 pages, 11 figures. Accepted for publication in Astronomy & Astrophysics

We aim to better understand the heating of the gas by observing the prominent gas cooling line [CII] at 158um in the low-metallicity environment of the Local Group spiral galaxy M33 at scales of 280pc. In particular, we aim at describing the variation of the photoelectric heating efficiency with galactic environment. In this unbiased study, we used ISO/LWS [CII] observations along the major axis of M33, in combination with Herschel PACS and SPIRE continuum maps, IRAM 30m CO 2-1 and VLA HI data to study the variation of velocity integrated intensities. The ratio of [CII] emission over the far-infrared continuum is used as a proxy for the heating efficiency, and models of photon-dominated regions are used to study the local physical densities, FUV radiation fields, and average column densities of the molecular clouds. The heating efficiency stays constant at 0.8% in the inner 4.5kpc radius of the galaxy where it starts to increase to reach values of ~3% in the outskirts at about 6kpc radial distance. The rise of efficiency is explained in the framework of PDR models by lowered volume densities and FUV fields, for optical extinctions of only a few magnitudes at constant metallicity. In view of the significant fraction of HI emission stemming from PDRs, and for typical pressures found in the Galactic cold neutral medium (CNM) traced by HI emission, the CNM contributes ~15% to the observed [CII] emission in the inner 2kpc radius of M33. The CNM contribution remains largely undetermined in the south, while positions between 2 and 7.3kpc radial distance in the north of M33 show a contribution of ~40%+-20%., Comment: (16 pages, a version with higher resolution figures may be downloaded from http://www.iram.es/IRAMES/hermesWiki/HermesPublications)

Within the framework of the Herschel M 33 extended survey HerM33es we study the Spectral Energy Distribution (SED) of a set of HII regions in M 33 as a function of the morphology. We present a catalogue of 119 HII regions morphologically classified: 9 filled, 47 mixed, 36 shell, and 27 clear shell HII regions. For each object we extract the photometry at twelve available wavelength bands (from FUV-1516A to IR-250mi) and obtain the SED. We also obtain emission line profiles across the regions to study the location of the stellar, ionised gas, and dust components. We find trends for the SEDs related to the morphology, showing that the star and gas-dust configuration affects the ratios of the emission in different bands. The mixed and filled regions show higher emission at 24mi than the shells and clear shells, which could be due to the proximity of the dust to the stellar clusters in the case of mixed and filled regions. The FIR peak for shells and clear shells seems to be located towards longer wavelengths, indicating that the dust is colder for this type of objects.The logarithmic 100/70mi ratio for filled and mixed regions remains constant over one order of magnitude in Halpha and FUV surface brightness, while the shells and clear shells exhibit a wider range of values of almost two orders of magnitude. We derive dust masses and temperatures fitting the individual SEDs with dust models proposed in the literature. The derived dust mass range is between 10^2-10^4 Msun and the cold dust temperature spans T(cold)~12-27 K. The spherical geometrical model proposed for the Halpha clear shells is confirmed by the emission profile obtained from the observations and is used to infer the electron density within the envelope: the typical electron density is 0.7+-0.3 cm^-3, while filled regions can reach values two to five times higher., Comment: 21 pages, 19 figures. Accepted to A&A

We aim to better understand the emission of molecular tracers of the diffuse and dense gas in giant molecular clouds and the influence that metallicity, optical extinction, density, far-UV field, and star formation rate have on these tracers. Using the IRAM 30m telescope, we detected HCN, HCO+, 12CO, and 13CO in six GMCs along the major axis of M33 at a resolution of ~ 114pc and out to a radial distance of 3.4kpc. Optical, far-infrared, and submillimeter data from Herschel and other observatories complement these observations. To interpret the observed molecular line emission, we created two grids of models of photon-dominated regions, one for solar and one for M33-type subsolar metallicity. The observed HCO+/HCN line ratios range between 1.1 and 2.5. Similarly high ratios have been observed in the Large Magellanic Cloud. The HCN/CO ratio varies between 0.4% and 2.9% in the disk of M33. The 12CO/13CO line ratio varies between 9 and 15 similar to variations found in the diffuse gas and the centers of GMCs of the Milky Way. Stacking of all spectra allowed HNC and C2H to be detected. The resulting HCO+/HNC and HCN/HNC ratios of ~ 8 and 6, respectively, lie at the high end of ratios observed in a large set of (ultra-)luminous infrared galaxies. HCN abundances are lower in the subsolar metallicity PDR models, while HCO+ abundances are enhanced. For HCN this effect is more pronounced at low optical extinctions. The observed HCO+/HCN and HCN/CO line ratios are naturally explained by subsolar PDR models of low optical extinctions between 4 and 10 mag and of moderate densities of n = 3x10^3 - 3x10^4 cm^-3, while the FUV field strength only has a small effect on the modeled line ratios. The line ratios are almost equally well reproduced by the solar-metallicity models, indicating that variations in metallicity only play a minor role in influencing these line ratios., Comment: 21 pages, 9 figures, published in A&A

We study the far-infrared emission from the nearby spiral galaxy M33 in order to investigate the dust physical properties such as the temperature and the luminosity density across the galaxy. Taking advantage of the unique wavelength coverage (100, 160, 250, 350 and 500 micron) of the Herschel Space Observatory and complementing our dataset with Spitzer-IRAC 5.8 and 8 micron and Spitzer-MIPS 24 and 70 micron data, we construct temperature and luminosity density maps by fitting two modified blackbodies of a fixed emissivity index of 1.5. We find that the 'cool' dust grains are heated at temperatures between 11 and 28 K with the lowest temperatures found in the outskirts of the galaxy and the highest ones in the center and in the bright HII regions. The infrared/submillimeter total luminosity (5 - 1000 micron) is estimated to be 1.9x10^9 Lsun. 59% of the total luminosity of the galaxy is produced by the 'cool' dust grains (~15 K) while the rest 41% is produced by 'warm' dust grains (~55 K). The ratio of the cool-to-warm dust luminosity is close to unity (within the computed uncertainties), throughout the galaxy, with the luminosity of the cool dust being slightly enhanced in the center of the galaxy. Decomposing the emission of the dust into two components (one emitted by the diffuse disk of the galaxy and one emitted by the spiral arms) we find that the fraction of the emission in the disk in the mid-infrared (24 micron) is 21%, while it gradually rises up to 57% in the submillimeter (500 micron). We find that the bulk of the luminosity comes from the spiral arm network that produces 70% of the total luminosity of the galaxy with the rest coming from the diffuse dust disk. The 'cool' dust inside the disk is heated at a narrow range of temperatures between 18 and 15 K (going from the center to the outer parts of the galaxy)., Comment: 12 pages, 14 figures, accepted for publication in A&A

Power spectra of de-projected images of late-type galaxies in gas and/or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power-laws, a shallow one at large scales (larger than 500 pc) and a steeper one at small scales, with the break between the two corresponding to the line-of-sight thickness of the galaxy disk. We present a thorough analysis of the power spectra of the dust and gas emission at several wavelengths in the nearby galaxy M33. In particular, we use the recently obtained images at five wavelengths by PACS and SPIRE onboard Herschel. The large dynamical range (2-3 dex in scale) of most images allow us to determine clearly the change in slopes from -1.5 to -4, with some variations with wavelength. The break scale is increasing with wavelength, from 100 pc at 24 and 100micron to 350 pc at 500micron, suggesting that the cool dust lies in a thicker disk than the warm dust, may be due to star formation more confined to the plane. The slope at small scale tends to be steeper at longer wavelength, meaning that the warmer dust is more concentrated in clumps. Numerical simulations of an isolated late-type galaxy, rich in gas and with no bulge, like M33, are carried out, in order to better interpret these observed results. Varying the star formation and feedback parameters, it is possible to obtain a range of power-spectra, with two power-law slopes and breaks, which nicely bracket the data. The small-scale power-law is indeed reflecting the 3D behaviour of the gas layer, steepening strongly while the feedback smoothes the structures, by increasing the gas turbulence. M33 appears to correspond to a fiducial model with an SFR of $\sim$ 0.7 Mo/yr, with 10% supernovae energy coupled to the gas kinematics., Comment: 11 pages, 24 figures, accepted in Astronomy & Astrophysics

Context: The emission line of [CII] at 158 micron is one of the strongest cooling lines of the interstellar medium (ISM) in galaxies. Aims: Disentangling the relative contributions of the different ISM phases to [CII] emission, is a major topic of the HerM33es program, a Herschel key project to study the ISM in the nearby spiral galaxy M33. Methods: Using PACS, we have mapped the emission of [CII] 158 micron, [OI] 63 micron, and other FIR lines in a 2'x2' region of the northern spiral arm of M33, centered on the HII region BCLMP302. At the peak of H-alpha emission, we have observed in addition a velocity resolved [CII] spectrum using HIFI. We use scatterplots to compare these data with PACS 160 micron continuum maps, and with maps of CO and HI data, at a common resolution of 12 arcsec or 50 pc. Maps of H-alpha and 24 micron emission observed with Spitzer are used to estimate the SFR. We have created maps of the [CII] and [OI] 63 micron emission and detected [NII] 122 micron and NIII 57 micron at individual positions. Results: The [CII] line observed with HIFI is significantly broader than that of CO, and slightly blue-shifted. In addition, there is little spatial correlation between [CII] observed with PACS and CO over the mapped region. There is even less spatial correlation between [CII] and the atomic gas traced by HI. Detailed comparison of the observed intensities towards the HII region with models of photo ionization and photon dominated regions, confirms that a significant fraction, 20--30%, of the observed [CII] emission stems from the ionized gas and not from the molecular cloud. The gas heating efficiency, using the ratio between [CII] and the TIR as a proxy, varies between 0.07 and 1.5%, with the largest variations found outside the HII region., Comment: Accepted for publication in A&A. 12 pages, 12 figures

Within the key project "Herschel M33 extended survey" (HerM33es), we are studying the physical and chemical processes driving star formation and galactic evolution in the nearby galaxy M33, combining the study of local conditions affecting individual star formation with properties only becoming apparent on global scales. Here, we present recent results obtained by the HerM33es team. Combining Spitzer and Herschel data ranging from 3.6um to 500um, along with HI, Halpha, and GALEX UV data, we have studied the dust at high spatial resolutions of 150pc, providing estimators of the total infrared (TIR) brightness and of the star formation rate. While the temperature of the warm dust at high brightness is driven by young massive stars, evolved stellar populations appear to drive the temperature of the cold dust. Plane-parallel models of photon dominated regions (PDRs) fail to reproduce fully the [CII], [OI], and CO maps obtained in a first spectroscopic study of one 2'x2' subregion of M33, located on the inner, northern spiral arm and encompassing the HII region BCLMP302., Comment: 6 pages, to appear in the proceedings of the 5th Zermatt ISM Symposium "Conditions and impact of star formation: New results with Herschel and beyond"

We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $\sigma$ using SPIRE and recent CO and \HI\ observations; a variation in $\sigma$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33., Comment: 5 pages, 5 figures Accepted for publication in Astronomy and Astrophysics

Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33., Comment: 5 pages, 3 figures, accepted for publication in the A&A Herschel Special Issue

Within the framework of the HERM33ES Key Project, using the high resolution and sensitivity of the Herschel photometric data, we study the compact emission in the Local Group spiral galaxy M33 to investigate the nature of the compact SPIRE emission sources. We extracted a catalogue of sources at 250um in order to investigate the nature of this compact emission. Taking advantage of the unprecedented Herschel resolution at these wavelengths, we also focus on a more precise study of some striking Halpha shells in the northern part of the galaxy. We present a catalogue of 159 compact emission sources in M33 identified by SExtractor in the 250um SPIRE band that is the one that provides the best spatial resolution. We also measured fluxes at 24um and Halpha for those 159 extracted sources. The morphological study of the shells also benefits from a multiwavelength approach including Halpha, far-UV from GALEX, and infrared from both Spitzer IRAC 8um and MIPS 24um in order to make comparisons. For the 159 compact sources selected at 250um, we find a very strong Pearson correlation coefficient with the MIPS 24um emission (r24 = 0.94) and a rather strong correlation with the Halpha emission, although with more scatter (rHa = 0.83). The morphological study of the Halpha shells shows a displacement between far-ultraviolet, Halpha, and the SPIRE bands. The cool dust emission from SPIRE clearly delineates the Halpha shell structures. The very strong link between the 250um compact emission and the 24um and Halpha emissions, by recovering the star formation rate from standard recipes for HII regions, allows us to provide star formation rate calibrations based on the 250um compact emission alone. The different locations of the Halpha and far-ultraviolet emissions with respect to the SPIRE cool dust emission leads to a dynamical age of a few Myr for the Halpha shells and the associated cool dust., Comment: 4 pages, 3 figures, Accpeted for publication in the A&A Herschel Special Issue

Over the past few years several studies have provided estimates of the SFR (star-formation rate) or the total infrared luminosity from just one infrared band. However these relations are generally derived for entire galaxies, which are known to contain a large scale diffuse emission that is not necessarily related to the latest star-formation episode. We provide new relations to estimate the SFR from resolved star-forming regions at 100 mum and 160 mum. We select individual star-forming regions in the nearby (840 kpc) galaxy M33. We estimate the SFR combining the emission in Halpha and at 24 mum to calibrate the emission at 100 mum and 160 mum as SFR estimators, as mapped with PACS/Herschel. The data are obtained in the framework of the HERM33ES open time key project. There is less emission in the HII regions at 160 mum than at 100 mum. Over a dynamic range of almost 2 dex in Sigma(SFR) we find that the 100 mum emission is a nearly linear estimator of the SFR, whereas that at 160 mum is slightly superlinear. The behaviour of individual star-forming regions is surprisingly similar to that of entire galaxies. At high Sigma(SFR), star formation drives the dust temperature, whereas uncertainties and variations in radiation-transfer and dust-heated processes dominate at low Sigma(SFR). Detailed modelling of both galaxies and individual star forming regions will be needed to interpret similarities and differences between the two and assess the fraction of diffuse emission in galaxies., Comment: 5 pages, 3 figures, accepted for publication in the A&A Herschel special issue