Search

The H$\alpha$ and H$\beta$ emission line luminosities measured in a single integrated spectrum are affected in non-trivial ways by point-to-point variations in dust attenuation in a galaxy. This work investigates the impact of this variation when estimating global H$\alpha$ luminosities corrected for the presence of dust by a global Balmer decrement. Analytical arguments show that the dust-corrected H$\alpha$ luminosity is always underestimated when using the global H$\alpha$/H$\beta$ flux ratio to correct for dust attenuation. We measure this effect on 156 face-on star-forming galaxies from the Mapping Nearby Galaxies at APO (MaNGA) survey. At 1-2 kpc spatial resolution, the effect is small but systematic, with the integrated dust-corrected H$\alpha$ luminosity underestimated by $2$-$4$ per cent (and typically not more than by $10$ per cent), and depends on the specific star formation rate of the galaxy. Given the spatial resolution of MaNGA, these are lower limits for the effect. From Multi Unit Spectroscopic Explorer (MUSE) observations of NGC 628 with a resolution of 36 pc we find the discrepancy between the globally and the point-by-point dust-corrected H$\alpha$ luminosity to be $14 \pm 1$ per cent, which may still underestimate the true effect. We use toy models and simulations to show that the true difference depends strongly on the spatial variance of the H$\alpha$/H$\beta$ flux ratio, and on the slope of the relation between H$\alpha$ luminosity and dust attenuation within a galaxy. Larger samples of higher spatial resolution observations are required to quantify the dependence of this effect as a function of galaxy properties., Comment: MNRAS accepted

The H$\alpha$ and H$\beta$ emission line luminosities measured in a single integrated spectrum are affected in non-trivial ways by point-to-point variations in dust attenuation in a galaxy. This work investigates the impact of this variation when estimating global H$\alpha$ luminosities corrected for the presence of dust by a global Balmer decrement. Analytical arguments show that the dust-corrected H$\alpha$ luminosity is always underestimated when using the global H$\alpha$/H$\beta$ flux ratio to correct for dust attenuation. We measure this effect on 156 face-on star-forming galaxies from the Mapping Nearby Galaxies at APO (MaNGA) survey. At 1-2 kpc spatial resolution, the effect is small but systematic, with the integrated dust-corrected H$\alpha$ luminosity underestimated by $2$-$4$ per cent (and typically not more than by $10$ per cent), and depends on the specific star formation rate of the galaxy. Given the spatial resolution of MaNGA, these are lower limits for the effect. From Multi Unit Spectroscopic Explorer (MUSE) observations of NGC 628 with a resolution of 36 pc we find the discrepancy between the globally and the point-by-point dust-corrected H$\alpha$ luminosity to be $14 \pm 1$ per cent, which may still underestimate the true effect. We use toy models and simulations to show that the true difference depends strongly on the spatial variance of the H$\alpha$/H$\beta$ flux ratio, and on the slope of the relation between H$\alpha$ luminosity and dust attenuation within a galaxy. Larger samples of higher spatial resolution observations are required to quantify the dependence of this effect as a function of galaxy properties., Comment: MNRAS accepted

We present a novel method to model galactic scale star formation and the resulting emission from star clusters and the multi-phase interstellar medium. We combine global parameters, such as SFR and CMF, with {\sc warpfield} which provides a description of the feedback-driven evolution of individual star-forming regions. Our approach includes stellar evolution, stellar winds, radiation pressure, supernovae, all of which couple to the dynamical evolution of the parental cloud in a highly non-linear fashion. The heating of diffuse gas and dust is calculated self-consistently with the age, mass and density dependent escape fractions of the local star-forming regions. From this we construct the interstellar radiation field at any point in the galaxy, and we employ the multi-frequency Monte Carlo radiative transfer code {\sc polaris} to produce synthetic emission maps for the one-to-one comparison with observational data. We demonstrate the capabilities of our approach by applying the method to a Milky Way like galaxy built-up in a high-resolution cosmological MHD simulation. We give three examples. First, we compute the multi-scale distribution of electron $n_{e^-}$ and $T_{e^{-}}$ and synthesize the MW all-sky H$\alpha$ emission. We use a multipole expansion method to show that the resulting maps are consistent with observations. Second, we predict the expected \SIII 9530~\AA\ emission. This line is a key target of several planned large survey programs. It suffers less extinction than other diagnostic lines and provides information about star formation in very dense environments that are otherwise observationally not readily accessible. Third, we explore the effects of differential extinction as seen by an extra-galactic observer, and discuss the consequences for the correct interpretation of \Ha emission as a star-formation rate tracer at different viewing angles.(abridged), Comment: 23 pages, 26 figures, submitted

Crucial progress in our understanding of star formation and feedback will depend on the ability to obtain spatially resolved spectroscopic observations of \ion{H}{ii} regions, from which reliable instantaneous measurements of their physical conditions can be obtained. Acquiring these datasets across full galactic systems will prove crucial for obtaining population samples that enable us to understand the time evolution of similar regions, and the variability of conditions among coeval regions. Separating the spatial and temporal dependencies in such way for different physical processes involved in star formation and the injection of feedback is crucial to overcome the inherit degeneracies associated with observing instantaneous snapshots of a dynamic ISM at any given time. Emission line diagnostics are at the core of measuring the physical condition in \ion{H}{ii} regions (e.g. dynamics, SFR, chemical abundances, dust extinction, ionization and excitation, etc.). These measurements require high spatial resolution, contiguous coverage across full galactic systems, and sensitivities significantly deeper than past efforts. The spatial scale required to resolve the \ion{H}{ii} regions of a few pc is only attainable in the Local Group where very large sky coverage is necessary., Comment: 7 pages, 3 figures

We present a novel method to model galactic scale star formation and the resulting emission from star clusters and the multi-phase interstellar medium. We combine global parameters, such as SFR and CMF, with {\sc warpfield} which provides a description of the feedback-driven evolution of individual star-forming regions. Our approach includes stellar evolution, stellar winds, radiation pressure, supernovae, all of which couple to the dynamical evolution of the parental cloud in a highly non-linear fashion. The heating of diffuse gas and dust is calculated self-consistently with the age, mass and density dependent escape fractions of the local star-forming regions. From this we construct the interstellar radiation field at any point in the galaxy, and we employ the multi-frequency Monte Carlo radiative transfer code {\sc polaris} to produce synthetic emission maps for the one-to-one comparison with observational data. We demonstrate the capabilities of our approach by applying the method to a Milky Way like galaxy built-up in a high-resolution cosmological MHD simulation. We give three examples. First, we compute the multi-scale distribution of electron $n_{e^-}$ and $T_{e^{-}}$ and synthesize the MW all-sky H$\alpha$ emission. We use a multipole expansion method to show that the resulting maps are consistent with observations. Second, we predict the expected \SIII 9530~\AA\ emission. This line is a key target of several planned large survey programs. It suffers less extinction than other diagnostic lines and provides information about star formation in very dense environments that are otherwise observationally not readily accessible. Third, we explore the effects of differential extinction as seen by an extra-galactic observer, and discuss the consequences for the correct interpretation of \Ha emission as a star-formation rate tracer at different viewing angles.(abridged), Comment: 23 pages, 26 figures, submitted

Crucial progress in our understanding of star formation and feedback will depend on the ability to obtain spatially resolved spectroscopic observations of \ion{H}{ii} regions, from which reliable instantaneous measurements of their physical conditions can be obtained. Acquiring these datasets across full galactic systems will prove crucial for obtaining population samples that enable us to understand the time evolution of similar regions, and the variability of conditions among coeval regions. Separating the spatial and temporal dependencies in such way for different physical processes involved in star formation and the injection of feedback is crucial to overcome the inherit degeneracies associated with observing instantaneous snapshots of a dynamic ISM at any given time. Emission line diagnostics are at the core of measuring the physical condition in \ion{H}{ii} regions (e.g. dynamics, SFR, chemical abundances, dust extinction, ionization and excitation, etc.). These measurements require high spatial resolution, contiguous coverage across full galactic systems, and sensitivities significantly deeper than past efforts. The spatial scale required to resolve the \ion{H}{ii} regions of a few pc is only attainable in the Local Group where very large sky coverage is necessary., Comment: 7 pages, 3 figures

We present images of radio recombination lines (RRLs) at wavelengths around 17 cm from the star-forming region W49A to determine the kinematics of ionized gas in the THOR survey (The HI/OH/Recombination line survey of the inner Milky Way) at an angular resolution of 16.8"x13.8". The distribution of ionized gas appears to be affected by feedback processes from the star clusters in W49A. The velocity structure of the RRLs shows a complex behavior with respect to the molecular gas. We find a shell-like distribution of ionized gas as traced by RRL emission surrounding the central cluster of OB stars in W49A. We describe the evolution of the shell with the recent feedback model code WARPFIELD that includes the important physical processes and has previously been applied to the 30 Doradus region in the Large Magellanic Cloud. The cloud structure and dynamics of W49A are in agreement with a feedback-driven shell that is re-collapsing. The shell may have triggered star formation in other parts of W49A. We suggest that W49A is a potential candidate for star formation regulated by feedback-driven and re-collapsing shells., Comment: 17 pages, 9 figures. Published in A&A. Updated to match published version

We present MUSE integral field data of two HII region complexes in the Large Magellanic Cloud (LMC), N44 and N180. Both regions consist of a main superbubble and a number of smaller, more compact HII regions that formed on the edge of the superbubble. For a total of 11 HII regions, we systematically analyse the radiative and mechanical feedback from the massive O-type stars on the surrounding gas. We exploit the integral field property of the data and the coverage of the HeII$\lambda$5412 line to identify and classify the feedback-driving massive stars, and from the estimated spectral types and luminosity classes we determine the stellar radiative output in terms of the ionising photon flux $Q_{0}$. We characterise the HII regions in terms of their sizes, morphologies, ionisation structure, luminosity and kinematics, and derive oxygen abundances via emission line ratios. We analyse the role of different stellar feedback mechanisms for each region by measuring the direct radiation pressure, the pressure of the ionised gas, and the pressure of the shock-heated winds. We find that stellar winds and ionised gas are the main drivers of HII region expansion in our sample, while the direct radiation pressure is up to three orders of magnitude lower than the other terms. We relate the total pressure to the star formation rate per unit area, $\Sigma_{SFR}$, for each region and find that stellar feedback has a negative effect on star formation, and sets an upper limit to $\Sigma_{SFR}$ as a function of increasing pressure., Comment: Accepted for publication in MNRAS, 27 pages, 21 figures

We present images of radio recombination lines (RRLs) at wavelengths around 17 cm from the star-forming region W49A to determine the kinematics of ionized gas in the THOR survey (The HI/OH/Recombination line survey of the inner Milky Way) at an angular resolution of 16.8"x13.8". The distribution of ionized gas appears to be affected by feedback processes from the star clusters in W49A. The velocity structure of the RRLs shows a complex behavior with respect to the molecular gas. We find a shell-like distribution of ionized gas as traced by RRL emission surrounding the central cluster of OB stars in W49A. We describe the evolution of the shell with the recent feedback model code WARPFIELD that includes the important physical processes and has previously been applied to the 30 Doradus region in the Large Magellanic Cloud. The cloud structure and dynamics of W49A are in agreement with a feedback-driven shell that is re-collapsing. The shell may have triggered star formation in other parts of W49A. We suggest that W49A is a potential candidate for star formation regulated by feedback-driven and re-collapsing shells., Comment: 17 pages, 9 figures. Published in A&A. Updated to match published version

We present MUSE integral field data of two HII region complexes in the Large Magellanic Cloud (LMC), N44 and N180. Both regions consist of a main superbubble and a number of smaller, more compact HII regions that formed on the edge of the superbubble. For a total of 11 HII regions, we systematically analyse the radiative and mechanical feedback from the massive O-type stars on the surrounding gas. We exploit the integral field property of the data and the coverage of the HeII$\lambda$5412 line to identify and classify the feedback-driving massive stars, and from the estimated spectral types and luminosity classes we determine the stellar radiative output in terms of the ionising photon flux $Q_{0}$. We characterise the HII regions in terms of their sizes, morphologies, ionisation structure, luminosity and kinematics, and derive oxygen abundances via emission line ratios. We analyse the role of different stellar feedback mechanisms for each region by measuring the direct radiation pressure, the pressure of the ionised gas, and the pressure of the shock-heated winds. We find that stellar winds and ionised gas are the main drivers of HII region expansion in our sample, while the direct radiation pressure is up to three orders of magnitude lower than the other terms. We relate the total pressure to the star formation rate per unit area, $\Sigma_{SFR}$, for each region and find that stellar feedback has a negative effect on star formation, and sets an upper limit to $\Sigma_{SFR}$ as a function of increasing pressure., Comment: Accepted for publication in MNRAS, 27 pages, 21 figures

We use the Spitzer SAGE survey of the Magellanic Clouds to evaluate the relationship between the 8-mic PAH emission, 24-mic hot dust emission, and HII region radiative transfer. We confirm that in the higher-metallicity Large Magellanic Cloud, PAH destruction is sensitive to optically thin conditions in the nebular Lyman continuum: objects identified as optically thin candidates based on nebular ionization structure show 6 times lower median 8-mic surface brightness (0.18 mJy arcsec^-2) than their optically thick counterparts (1.2 mJy arcsec^-2). The 24-mic surface brightness also shows a factor of 3 offset between the two classes of objects (0.13 vs 0.44 mJy arcsec^-2, respectively), which is driven by the association between the very small dust grains and higher density gas found at higher nebular optical depths. In contrast, PAH and dust formation in the low-metallicity Small Magellanic Cloud is strongly inhibited such that we find no variation in either 8-mic or 24-mic emission between our optically thick and thin samples. This is attributable to extremely low PAH and dust production together with high, corrosive UV photon fluxes in this low-metallicity environment. The dust mass surface densities and gas-to-dust ratios determined from dust maps using Herschel HERITAGE survey data support this interpretation., Comment: Accepted to ApJ, May 15, 2017. 10 pages, 9 figures

We use the Spitzer SAGE survey of the Magellanic Clouds to evaluate the relationship between the 8-mic PAH emission, 24-mic hot dust emission, and HII region radiative transfer. We confirm that in the higher-metallicity Large Magellanic Cloud, PAH destruction is sensitive to optically thin conditions in the nebular Lyman continuum: objects identified as optically thin candidates based on nebular ionization structure show 6 times lower median 8-mic surface brightness (0.18 mJy arcsec^-2) than their optically thick counterparts (1.2 mJy arcsec^-2). The 24-mic surface brightness also shows a factor of 3 offset between the two classes of objects (0.13 vs 0.44 mJy arcsec^-2, respectively), which is driven by the association between the very small dust grains and higher density gas found at higher nebular optical depths. In contrast, PAH and dust formation in the low-metallicity Small Magellanic Cloud is strongly inhibited such that we find no variation in either 8-mic or 24-mic emission between our optically thick and thin samples. This is attributable to extremely low PAH and dust production together with high, corrosive UV photon fluxes in this low-metallicity environment. The dust mass surface densities and gas-to-dust ratios determined from dust maps using Herschel HERITAGE survey data support this interpretation., Comment: Accepted to ApJ, May 15, 2017. 10 pages, 9 figures

We study the spectral energy distribution (SED) of the radio continuum emission from the KINGFISH sample of nearby galaxies to understand the energetics and origin of this emission. Effelsberg multi-wavelength observations at 1.4GHz, 4.8GHz, 8.5GHz, and 10.5GHz combined with archive data allow us, for the first time, to determine the mid-radio continuum (1-10 GHz, MRC) bolometric luminosities and further present calibration relations vs. the monochromatic radio luminosities. The 1-10 GHz radio SED is fitted using a Bayesian Markov Chain Monte Carlo (MCMC) technique leading to measurements for the nonthermal spectral index and the thermal fraction f_th with mean values of alpha_nt=0.97+-0.16 (0.79+-0.15 for the total spectral index) and f_th= 10% +- 9% at 1.4 GHz. The MRC luminosity changes over ~3 orders of magnitude in the sample. The thermal emission is responsible for ~23% of the MRC on average. We also compare the extinction-corrected diagnostics of star formation rate with the thermal and nonthermal radio tracers and derive the first star formation calibration relations using the MRC radio luminosity. The nonthermal spectral index flattens with increasing star formation rate surface density, indicating the effect of the star formation feedback on the cosmic ray electron population in galaxies. Comparing the radio and IR SEDs, we find that the FIR-to-MRC ratio could decrease with star formation rate, due to the amplification of the magnetic fields in star forming regions. This particularly implies a decrease in the ratio at high redshifts, where mostly luminous/star forming galaxies are detected., Comment: 29 pages, 15 figures, accepted for publication in Astrophysical Journal

We study the spectral energy distribution (SED) of the radio continuum emission from the KINGFISH sample of nearby galaxies to understand the energetics and origin of this emission. Effelsberg multi-wavelength observations at 1.4GHz, 4.8GHz, 8.5GHz, and 10.5GHz combined with archive data allow us, for the first time, to determine the mid-radio continuum (1-10 GHz, MRC) bolometric luminosities and further present calibration relations vs. the monochromatic radio luminosities. The 1-10 GHz radio SED is fitted using a Bayesian Markov Chain Monte Carlo (MCMC) technique leading to measurements for the nonthermal spectral index and the thermal fraction f_th with mean values of alpha_nt=0.97+-0.16 (0.79+-0.15 for the total spectral index) and f_th= 10% +- 9% at 1.4 GHz. The MRC luminosity changes over ~3 orders of magnitude in the sample. The thermal emission is responsible for ~23% of the MRC on average. We also compare the extinction-corrected diagnostics of star formation rate with the thermal and nonthermal radio tracers and derive the first star formation calibration relations using the MRC radio luminosity. The nonthermal spectral index flattens with increasing star formation rate surface density, indicating the effect of the star formation feedback on the cosmic ray electron population in galaxies. Comparing the radio and IR SEDs, we find that the FIR-to-MRC ratio could decrease with star formation rate, due to the amplification of the magnetic fields in star forming regions. This particularly implies a decrease in the ratio at high redshifts, where mostly luminous/star forming galaxies are detected., Comment: 29 pages, 15 figures, accepted for publication in Astrophysical Journal

We present new observations of the (6,6) and (9,9) inversion transitions of the hydronium ion toward Sagittarius B2(N) and W31C. Sensitive observations toward Sagittarius B2(N) show that the high, ~ 500 K, rotational temperatures characterizing the population of the highly-excited metastable H3O+ rotational levels are present over a wide range of velocities corresponding to the Sagittarius B2 envelope, as well as the foreground gas clouds between the Sun and the source. Observations of the same lines toward W31C, a line of sight that does not intersect the Central Molecular Zone, but instead traces quiescent gas in the Galactic disk, also imply a high rotational temperature of ~ 380 K, well in excess of the kinetic temperature of the diffuse Galactic interstellar medium. While it is plausible that some fraction of the molecular gas may be heated to such high temperatures in the active environment of the Galactic center, characterized by high X-ray and cosmic ray fluxes, shocks and high degree of turbulence, this is unlikely in the largely quiescent environment of the Galactic disk clouds. We suggest instead that the highly-excited states of the hydronium ion are populated mainly by exoergic chemical formation processes and temperature describing the rotational level population does not represent the physical temperature of the medium. The same arguments may be applicable to other symmetric top rotors, such as ammonia. This offers a simple explanation to the long-standing puzzle of the presence of a pervasive, hot molecular gas component in the central region of the Milky Way. Moreover, our observations suggest that this is a universal process, not limited to the active environments associated with galactic nuclei., Comment: Accepted for publication in The Astrophysical Journal, 10 pages, 4 figures, 3 tables

Ionization-parameter mapping (IPM) is a powerful technique for tracing the optical depth of Lyman continuum radiation from massive stars. Using narrow-band line-ratio maps, we examine trends in radiative feedback from ordinary HII regions of the Magellanic Clouds and nearby starburst galaxies. We find that the aggregate escape fraction for the Lyman continuum is sufficient to ionize the diffuse, warm ionized medium in the Magellanic Clouds, and that more luminous nebulae are more likely to be optically thin. We apply ionization-parameter mapping to entire starburst galaxies, revealing ionization cones in two nearby starbursts. Within the limits of our small sample, we examine the conditions for the propagation of ionizing radiation beyond the host galaxies., Comment: Invited paper to appear in "Massive Young Star Clusters Near and Far: From the Milky Way to Reionization", 2013 Guillermo Haro Conference, Eds. Y. D. Mayya, D. Rosa-Gonzalez & E. Terlevich, INAOE and AMC. 6 pages, 4 figures

We present new Herschel observations of the (6,6) and (9,9) inversion transitions of the hydronium ion toward Sagittarius B2(N) and W31C. Sensitive observations toward Sagittarius B2(N) show that the high, similar to 500 K, rotational temperatures characterizing the population of the highly excited metastable H3O+ rotational levels are present over a wide range of velocities corresponding to the Sagittarius B2 envelope, as well as the foreground gas clouds between the Sun and the source. Observations of the same lines toward W31C, a line of sight that does not intersect the Central Molecular Zone but instead traces quiescent gas in the Galactic disk, also imply a high rotational temperature of similar to 380 K, well in excess of the kinetic temperature of the diffuse Galactic interstellar medium. While it is plausible that some fraction of the molecular gas may be heated to such high temperatures in the active environment of the Galactic center, characterized by high X-ray and cosmic-ray fluxes, shocks, and high degree of turbulence, this is unlikely in the largely quiescent environment of the Galactic disk clouds. We suggest instead that the highly excited states of the hydronium ion are populated mainly by exoergic chemical formation processes and the temperature describing the rotational level population does not represent the physical temperature of the medium. The same arguments may be applicable to other symmetric top rotors, such as ammonia. This offers a simple explanation of the long-standing puzzle of the presence of a pervasive, hot molecular gas component in the central region of the Milky Way. Moreover, our observations suggest that this is a universal process not limited to the active environments associated with galactic nuclei.

We present new Herschel observations of the (6,6) and (9,9) inversion transitions of the hydronium ion toward Sagittarius B2(N) and W31C. Sensitive observations toward Sagittarius B2(N) show that the high, similar to 500 K, rotational temperatures characterizing the population of the highly excited metastable H3O+ rotational levels are present over a wide range of velocities corresponding to the Sagittarius B2 envelope, as well as the foreground gas clouds between the Sun and the source. Observations of the same lines toward W31C, a line of sight that does not intersect the Central Molecular Zone but instead traces quiescent gas in the Galactic disk, also imply a high rotational temperature of similar to 380 K, well in excess of the kinetic temperature of the diffuse Galactic interstellar medium. While it is plausible that some fraction of the molecular gas may be heated to such high temperatures in the active environment of the Galactic center, characterized by high X-ray and cosmic-ray fluxes, shocks, and high degree of turbulence, this is unlikely in the largely quiescent environment of the Galactic disk clouds. We suggest instead that the highly excited states of the hydronium ion are populated mainly by exoergic chemical formation processes and the temperature describing the rotational level population does not represent the physical temperature of the medium. The same arguments may be applicable to other symmetric top rotors, such as ammonia. This offers a simple explanation of the long-standing puzzle of the presence of a pervasive, hot molecular gas component in the central region of the Milky Way. Moreover, our observations suggest that this is a universal process not limited to the active environments associated with galactic nuclei.

We present new observations of the (6,6) and (9,9) inversion transitions of the hydronium ion toward Sagittarius B2(N) and W31C. Sensitive observations toward Sagittarius B2(N) show that the high, ~ 500 K, rotational temperatures characterizing the population of the highly-excited metastable H3O+ rotational levels are present over a wide range of velocities corresponding to the Sagittarius B2 envelope, as well as the foreground gas clouds between the Sun and the source. Observations of the same lines toward W31C, a line of sight that does not intersect the Central Molecular Zone, but instead traces quiescent gas in the Galactic disk, also imply a high rotational temperature of ~ 380 K, well in excess of the kinetic temperature of the diffuse Galactic interstellar medium. While it is plausible that some fraction of the molecular gas may be heated to such high temperatures in the active environment of the Galactic center, characterized by high X-ray and cosmic ray fluxes, shocks and high degree of turbulence, this is unlikely in the largely quiescent environment of the Galactic disk clouds. We suggest instead that the highly-excited states of the hydronium ion are populated mainly by exoergic chemical formation processes and temperature describing the rotational level population does not represent the physical temperature of the medium. The same arguments may be applicable to other symmetric top rotors, such as ammonia. This offers a simple explanation to the long-standing puzzle of the presence of a pervasive, hot molecular gas component in the central region of the Milky Way. Moreover, our observations suggest that this is a universal process, not limited to the active environments associated with galactic nuclei., Comment: Accepted for publication in The Astrophysical Journal, 10 pages, 4 figures, 3 tables

Ionization-parameter mapping (IPM) is a powerful technique for tracing the optical depth of Lyman continuum radiation from massive stars. Using narrow-band line-ratio maps, we examine trends in radiative feedback from ordinary HII regions of the Magellanic Clouds and nearby starburst galaxies. We find that the aggregate escape fraction for the Lyman continuum is sufficient to ionize the diffuse, warm ionized medium in the Magellanic Clouds, and that more luminous nebulae are more likely to be optically thin. We apply ionization-parameter mapping to entire starburst galaxies, revealing ionization cones in two nearby starbursts. Within the limits of our small sample, we examine the conditions for the propagation of ionizing radiation beyond the host galaxies., Comment: Invited paper to appear in "Massive Young Star Clusters Near and Far: From the Milky Way to Reionization", 2013 Guillermo Haro Conference, Eds. Y. D. Mayya, D. Rosa-Gonzalez & E. Terlevich, INAOE and AMC. 6 pages, 4 figures

We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows., Comment: Accepted for publication in ApJL 11/15/2013

We present a sample of 14 OB stars in the Small Magellanic Cloud that meet strong criteria for having formed under extremely sparse star-forming conditions in the field. These stars are a minimum of 28 pc in projection from other OB stars, and they are centered within symmetric, round HII regions. They show no evidence of bow shocks, implying that the targets are not transverse runaway stars. Their radial velocities relative to local HI also indicate that they are not line-of-sight runaway stars. A friends-of-friends analysis shows that 9 of the objects present a few low-mass companion stars, with typical mass ratios for the two highest-mass stars of around 0.1. This further substantiates that these OB stars formed in place, and that they can and do form in extremely sparse conditions. This poses strong constraints on theories of star formation and challenges proposed relations between cluster mass and maximum stellar mass., Comment: Accepted to ApJ, 12 pages

We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows., Comment: Accepted for publication in ApJL 11/15/2013

We present a sample of 14 OB stars in the Small Magellanic Cloud that meet strong criteria for having formed under extremely sparse star-forming conditions in the field. These stars are a minimum of 28 pc in projection from other OB stars, and they are centered within symmetric, round HII regions. They show no evidence of bow shocks, implying that the targets are not transverse runaway stars. Their radial velocities relative to local HI also indicate that they are not line-of-sight runaway stars. A friends-of-friends analysis shows that 9 of the objects present a few low-mass companion stars, with typical mass ratios for the two highest-mass stars of around 0.1. This further substantiates that these OB stars formed in place, and that they can and do form in extremely sparse conditions. This poses strong constraints on theories of star formation and challenges proposed relations between cluster mass and maximum stellar mass., Comment: Accepted to ApJ, 12 pages

The shape of the OB-star spectral energy distribution is a critical component in many diagnostics of the ISM and galaxy properties. We use single-star HII regions from the LMC to quantitatively examine the ionizing SEDs from widely available CoStar, TLUSTY, and WM-basic atmosphere grids. We evaluate the stellar atmosphere models by matching the emission-line spectra that they predict from CLOUDY photoionization simulations with those observed from the nebulae. The atmosphere models are able to reproduce the observed optical nebular line ratios, except at the highest energy transitions > 40 eV, assuming that the gas distribution is non-uniform. Overall we find that simulations using WM-basic produce the best agreement with the observed line ratios. The rate of ionizing photons produced by the model SEDs is consistent with the rate derived from the \Halpha\ luminosity for standard, log(g) = 4.0 models adopted from the atmosphere grids. However, there is a systematic offset between the rate of ionizing photons from different atmosphere models that is correlated with the relative hardness of the SEDs. In general WM-basic and TLUSTY atmosphere models predict similar effective temperatures, while CoStar predicts effective temperatures that are cooler by a few thousand degrees. We compare our effective temperatures, which depend on the nebular ionization balance, to conventional photospheric-based calibrations from the literature. We suggest that in the future, spectral type to effective temperature calibrations can be constructed from nebular data., Comment: Submitted to ApJ

We exploit ionization-parameter mapping as a powerful tool to measure the optical depth of star-forming HII regions. Our simulations using the photoionization code CLOUDY and our new, SURFBRIGHT surface brightness simulator demonstrate that this technique can directly diagnose most density-bounded, optically thin nebulae using spatially resolved emission line data. We apply this method to the Large and Small Magellanic Clouds, using the data from the Magellanic Clouds Emission Line Survey. We generate new HII region catalogs based on photoionization criteria set by the observed ionization structure in the [SII]/[OIII] ratio and Ha surface brightness. The luminosity functions from these catalogs generally agree with those from Ha-only surveys. We then use ionization-parameter mapping to crudely classify all the nebulae into optically thick vs optically thin categories, yielding fundamental new insights into Lyman continuum radiation transfer. We find that in both galaxies, the frequency of optically thin objects correlates with Ha luminosity, and that the numbers of these objects dominate above log L\geq37.0. The frequencies of optically thin objects are 40% and 33% in the LMC and SMC, respectively. Similarly, the frequency of optically thick regions correlates with H I column density, with optically thin objects dominating at the lowest N(HI). The integrated escape luminosity of ionizing radiation is dominated by the largest regions, and corresponds to luminosity-weighted, ionizing escape fractions from the H II region population of \geq0.42 and \geq0.40 in the LMC and SMC, respectively. These values correspond to global galactic escape fractions of 4% and 11%, respectively. This is sufficient to power the ionization rate of the observed diffuse ionized gas in both galaxies. Our results suggest the possibility of significant galactic escape fractions of Lyman continuum radiation., Comment: 31 pages, 24 figures, 6 tables, accepted to ApJ Errors in Tables B1 and B2 object types corrected. Erratum bundled with arxiv submission as "erratum_ms.pdf"

We exploit ionization-parameter mapping as a powerful tool to measure the optical depth of star-forming HII regions. Our simulations using the photoionization code CLOUDY and our new, SURFBRIGHT surface brightness simulator demonstrate that this technique can directly diagnose most density-bounded, optically thin nebulae using spatially resolved emission line data. We apply this method to the Large and Small Magellanic Clouds, using the data from the Magellanic Clouds Emission Line Survey. We generate new HII region catalogs based on photoionization criteria set by the observed ionization structure in the [SII]/[OIII] ratio and Ha surface brightness. The luminosity functions from these catalogs generally agree with those from Ha-only surveys. We then use ionization-parameter mapping to crudely classify all the nebulae into optically thick vs optically thin categories, yielding fundamental new insights into Lyman continuum radiation transfer. We find that in both galaxies, the frequency of optically thin objects correlates with Ha luminosity, and that the numbers of these objects dominate above log L\geq37.0. The frequencies of optically thin objects are 40% and 33% in the LMC and SMC, respectively. Similarly, the frequency of optically thick regions correlates with H I column density, with optically thin objects dominating at the lowest N(HI). The integrated escape luminosity of ionizing radiation is dominated by the largest regions, and corresponds to luminosity-weighted, ionizing escape fractions from the H II region population of \geq0.42 and \geq0.40 in the LMC and SMC, respectively. These values correspond to global galactic escape fractions of 4% and 11%, respectively. This is sufficient to power the ionization rate of the observed diffuse ionized gas in both galaxies. Our results suggest the possibility of significant galactic escape fractions of Lyman continuum radiation., Comment: 31 pages, 24 figures, 6 tables, accepted to ApJ Errors in Tables B1 and B2 object types corrected. Erratum bundled with arxiv submission as "erratum_ms.pdf"

The shape of the OB-star spectral energy distribution is a critical component in many diagnostics of the ISM and galaxy properties. We use single-star HII regions from the LMC to quantitatively examine the ionizing SEDs from widely available CoStar, TLUSTY, and WM-basic atmosphere grids. We evaluate the stellar atmosphere models by matching the emission-line spectra that they predict from CLOUDY photoionization simulations with those observed from the nebulae. The atmosphere models are able to reproduce the observed optical nebular line ratios, except at the highest energy transitions > 40 eV, assuming that the gas distribution is non-uniform. Overall we find that simulations using WM-basic produce the best agreement with the observed line ratios. The rate of ionizing photons produced by the model SEDs is consistent with the rate derived from the \Halpha\ luminosity for standard, log(g) = 4.0 models adopted from the atmosphere grids. However, there is a systematic offset between the rate of ionizing photons from different atmosphere models that is correlated with the relative hardness of the SEDs. In general WM-basic and TLUSTY atmosphere models predict similar effective temperatures, while CoStar predicts effective temperatures that are cooler by a few thousand degrees. We compare our effective temperatures, which depend on the nebular ionization balance, to conventional photospheric-based calibrations from the literature. We suggest that in the future, spectral type to effective temperature calibrations can be constructed from nebular data., Comment: Submitted to ApJ

The Orion Bar is one of the nearest and best-studied photodissociation or photon-dominated regions (PDRs). Observations reveal the presence of H2 lines from vibrationally or rotationally excited upper levels that suggest warm gas temperatures (400 to 700 K). However, standard models of PDRs are unable to reproduce such warm rotational temperatures. In this paper we attempt to explain these observations with new comprehensive models which extend from the H+ region through the Bar and include the magnetic field in the equation of state. We adopt the model parameters from our previous paper which successfully reproduced a wide variety of spectral observations across the Bar. In this model the local cosmic-ray density is enhanced above the galactic background, as is the magnetic field, and which increases the cosmic-ray heating elevating the temperature in the molecular region. The pressure is further enhanced above the gas pressure in the H+ region by the momentum transferred from the absorbed starlight. Here we investigate whether the observed H2 lines can be reproduced with standard assumptions concerning the grain photoelectric emission. We also explore the effects due to the inclusion of recently computed H2 + H2, H2 + H and H2 + He collisional rate coefficients., Comment: Accepted for publication in ApJ (34 pages, including 16 figures)

The Orion Bar is one of the nearest and best-studied photodissociation or photon-dominated regions (PDRs). Observations reveal the presence of H2 lines from vibrationally or rotationally excited upper levels that suggest warm gas temperatures (400 to 700 K). However, standard models of PDRs are unable to reproduce such warm rotational temperatures. In this paper we attempt to explain these observations with new comprehensive models which extend from the H+ region through the Bar and include the magnetic field in the equation of state. We adopt the model parameters from our previous paper which successfully reproduced a wide variety of spectral observations across the Bar. In this model the local cosmic-ray density is enhanced above the galactic background, as is the magnetic field, and which increases the cosmic-ray heating elevating the temperature in the molecular region. The pressure is further enhanced above the gas pressure in the H+ region by the momentum transferred from the absorbed starlight. Here we investigate whether the observed H2 lines can be reproduced with standard assumptions concerning the grain photoelectric emission. We also explore the effects due to the inclusion of recently computed H2 + H2, H2 + H and H2 + He collisional rate coefficients., Comment: Accepted for publication in ApJ (34 pages, including 16 figures)

Previous work has shown the Orion Bar to be an interface between ionized and molecular gas, viewed roughly edge on, which is excited by the light from the Trapezium cluster. Much of the emission from any star-forming region will originate from such interfaces, so the Bar serves as a foundation test of any emission model. Here we combine X-ray, optical, IR and radio data sets to derive emission spectra along the transition from H+ to H0 to H2 regions. We then reproduce the spectra of these layers with a simulation that simultaneously accounts for the detailed microphysics of the gas, the grains, and molecules, especially H2 and CO. The magnetic field, observed to be the dominant pressure in another region of the Orion Nebula, is treated as a free parameter, along with the density of cosmic rays. Our model successfully accounts for the optical, IR and radio observations across the Bar by including a significant magnetic pressure and also heating by an excess density of cosmic rays, which we suggest is due to cosmic rays being trapped in the compressed magnetic field. In the Orion Bar, as we had previously found in M17, momentum carried by radiation and winds from the newly formed stars pushes back and compresses the surrounding gas. There is a rough balance between outward momentum in starlight and the total pressure in atomic and molecular gas surrounding the H+ region. If the gas starts out with a weak magnetic field, the starlight from a newly formed cluster will push back the gas and compress the gas, magnetic field, and cosmic rays until magnetic pressure becomes an important factor., Comment: 32 pages, including 10 figures. ApJ (in press). Revised to correct 2 spelling errors

Previous work has shown the Orion Bar to be an interface between ionized and molecular gas, viewed roughly edge on, which is excited by the light from the Trapezium cluster. Much of the emission from any star-forming region will originate from such interfaces, so the Bar serves as a foundation test of any emission model. Here we combine X-ray, optical, IR and radio data sets to derive emission spectra along the transition from H+ to H0 to H2 regions. We then reproduce the spectra of these layers with a simulation that simultaneously accounts for the detailed microphysics of the gas, the grains, and molecules, especially H2 and CO. The magnetic field, observed to be the dominant pressure in another region of the Orion Nebula, is treated as a free parameter, along with the density of cosmic rays. Our model successfully accounts for the optical, IR and radio observations across the Bar by including a significant magnetic pressure and also heating by an excess density of cosmic rays, which we suggest is due to cosmic rays being trapped in the compressed magnetic field. In the Orion Bar, as we had previously found in M17, momentum carried by radiation and winds from the newly formed stars pushes back and compresses the surrounding gas. There is a rough balance between outward momentum in starlight and the total pressure in atomic and molecular gas surrounding the H+ region. If the gas starts out with a weak magnetic field, the starlight from a newly formed cluster will push back the gas and compress the gas, magnetic field, and cosmic rays until magnetic pressure becomes an important factor., Comment: 32 pages, including 10 figures. ApJ (in press). Revised to correct 2 spelling errors

The southwestern (SW) part of the Galactic H II region M17 contains an obscured ionization front that is most easily seen at infrared and radio wavelengths. It is nearly edge-on, thus offering an excellent opportunity to study the way in which the gas changes from fully ionized to molecular as radiation from the ionizing stars penetrates into the gas. M17 is also one of the very few H II regions for which the magnetic field strength can be measured in the photodissociation region ( PDR) that forms the interface between the ionized and molecular gas. Here we model an observed line of sight through the gas cloud, including the H+, H0 (PDR), and molecular layers, in a fully self-consistent single calculation. An interesting aspect of the M17 SW bar is that the PDR is very extended. We show that the strong magnetic field that is observed to be present inevitably leads to a very deep PDR, because the structure of the neutral and molecular gas is dominated by magnetic pressure, rather than by gas pressure, as previously had been supposed.We also show that a wide variety of observed facts can be explained if a hydrostatic geometry prevails, in which the gas pressure from an inner X-ray hot bubble and the outward momentum of the stellar radiation field compress the gas and its associated magnetic field in the PDR, as has already been shown to occur in the Orion Nebula. The magnetic field compression may also amplify the local cosmic-ray density. The pressure in the observed magnetic field balances the outward forces, suggesting that the observed geometry is a natural consequence of the formation of a star cluster within a molecular cloud., Comment: Published as 2007, ApJ,658,1119

The southwestern (SW) part of the Galactic H II region M17 contains an obscured ionization front that is most easily seen at infrared and radio wavelengths. It is nearly edge-on, thus offering an excellent opportunity to study the way in which the gas changes from fully ionized to molecular as radiation from the ionizing stars penetrates into the gas. M17 is also one of the very few H II regions for which the magnetic field strength can be measured in the photodissociation region ( PDR) that forms the interface between the ionized and molecular gas. Here we model an observed line of sight through the gas cloud, including the H+, H0 (PDR), and molecular layers, in a fully self-consistent single calculation. An interesting aspect of the M17 SW bar is that the PDR is very extended. We show that the strong magnetic field that is observed to be present inevitably leads to a very deep PDR, because the structure of the neutral and molecular gas is dominated by magnetic pressure, rather than by gas pressure, as previously had been supposed.We also show that a wide variety of observed facts can be explained if a hydrostatic geometry prevails, in which the gas pressure from an inner X-ray hot bubble and the outward momentum of the stellar radiation field compress the gas and its associated magnetic field in the PDR, as has already been shown to occur in the Orion Nebula. The magnetic field compression may also amplify the local cosmic-ray density. The pressure in the observed magnetic field balances the outward forces, suggesting that the observed geometry is a natural consequence of the formation of a star cluster within a molecular cloud., Comment: Published as 2007, ApJ,658,1119