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.

We present resolved [CI] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from $J_\mathrm{up} = 1$ to 7 to interpret what phase of the interstellar medium the [CI] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4-3) and [CI](2-1) lines; we hypothesize this is due to the similar upper level temperature of these two lines. We modeled the [CI] and CO line emission using large velocity gradient models combined with an empirical template. According to this modeling, the [CI](1-0) line is clearly dominated by the low-excitation component. We determine [CI] to molecular mass conversion factors for both the [CI](1-0) and [CI](2-1) lines, with mean values of $\alpha_{\mathrm{[CI](1-0)}} = 7.3$ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ and $\alpha_{\mathrm{[CI](2-1)}} = 34 $ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ with logarithmic root-mean-square spreads of 0.20 and 0.32 dex, respectively. The similar spread of $\alpha_{\mathrm{[CI](1-0)}}$ to $\alpha_{\mathrm{CO}}$ (derived using the CO(2-1) line) suggests that [CI](1-0) may be just as good a tracer of cold molecular gas as CO(2-1) in galaxies of this type. On the other hand, the wider spread of $\alpha_{\mathrm{[CI](2-1)}}$ and the tight relation found between [CI](2-1) and CO(4-3) suggest that much of the [CI](2-1) emission may originate in warmer molecular gas.

We explore to what extent simple algebraic models can be used to describe H II regions when winds, radiation pressure, gravity and photon breakout are included. We a) develop algebraic models to describe the expansion of photoionised H II regions under the influence of gravity and accretion in power-law density fields with $\rho \propto r^{-w}$, b) determine when terms describing winds, radiation pressure, gravity and photon breakout become significant enough to affect the dynamics of the H II region where $w=2$, and c) solve these expressions for a set of physically-motivated conditions. We find that photoionisation feedback from massive stars is the principal mode of feedback on molecular cloud scales, driving accelerating outflows from molecular clouds in cases where the peaked density structure around young massive stars is considered at radii between $\sim$0.1 and 10-100 pc. Under a large range of conditions the effect of winds and radiation on the dynamics of H II regions is around 10% of the contribution from photoionisation. The effect of winds and radiation pressure are most important at high densities, either close to the star or in very dense clouds such as those in the Central Molecular Zone of the Milky Way. Out to $\sim$0.1 pc they are the principal drivers of the H II region. Lower metallicities make the relative effect of photoionisation even stronger as the ionised gas temperature is higher.

We present the {\sc warpfield} emission predictor, {\sc warpfield-emp}, which couples the 1D stellar feedback code {\sc warpfield} with the {\sc cloudy} \hii region/PDR code and the {\sc polaris} radiative transfer code, in order to make detailed predictions for the time-dependent line and continuum emission arising from the H{\sc ii} region and PDR surrounding an evolving star cluster. {\sc warpfield-emp} accounts for a wide range of physical processes (stellar winds, supernovae, radiation pressure, gravity, thermal conduction, radiative cooling, dust extinction etc.) and yet runs quickly enough to allow us to explore broad ranges of different cloud parameters. We compare the results of an extensive set of models with SITELLE observations of a large sample of \hii regions in NGC~628 and find very good agreement, particularly for the highest signal-to-noise observations. We show that our approach of modeling individual clouds from first principles (instead of in terms of dimensionless quantities such as the ionization parameter) allows us to avoid long-standing degeneracies in the interpretation of \hii region diagnostics and enables us to relate these diagnostics to important physical parameters such as cloud mass or cluster age. Finally, we explore the implications of our models regarding the reliability of simple metallicity diagnostics, the properties of long-lived embedded clusters, and the role played by winds and supernovae in regulating \hii region and PDR line emission.

SIGNALS, the Star formation, Ionized Gas, and Nebular Abundances Legacy Survey, is a large observing program designed to investigate massive star formation and HII regions in a sample of local extended galaxies. The program will use the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope. Over 355 hours (54.7 nights) have been allocated beginning in fall 2018 for eight consecutive semesters. Once completed, SIGNALS will provide a statistically reliable laboratory to investigate massive star formation, including over 50 000 resolved HII regions : the largest, most complete, and homogeneous database of spectroscopically and spatially resolved extragalactic HII regions ever assembled. For each field observed, three datacubes covering the spectral bands of the filters SN1 (363 -386 nm), SN2 (482 - 513 nm), and SN3 (647 - 685 nm) are gathered. The spectral resolution selected for each spectral band is 1000, 1000, and 5000, respectively. As defined, the project sample will facilitate the study of small-scale nebular physics and many other phenomena linked to star formation at a mean spatial resolution of 20 pc. This survey also has considerable legacy value for additional topics including planetary nebulae, diffuse ionized gas, andsupernova remnants. The purpose of this paper is to present a general outlook of the survey, notably the observing strategy, galaxy sample, and science requirements.

We present an updated version of POLARIS, a well established code designated for dust polarisation and line radiative transfer (RT) in arbitrary astrophysical environments. We extend the already available capabilities with a synchrotron feature for polarised emission. Here, we combine state-of-the-art solutions of the synchrotron RT coefficients with numerical methods for solving the complete system of equations of the RT problem, including Faraday rotation (FR) as well as Faraday conversion (FC). We validate the code against Galactic and extragalactic observations by performing a statistical analysis of synthetic all-sky synchrotron maps for positions within the galaxy and for extragalactic observations. For these test scenarios we apply a model of the Milky Way based on sophisticated magneto-hydrodynamic (MHD) simulations and population-synthesis post-processing techniques.We explore different parameters for modeling the distribution of free electrons and for a turbulent magnetic field component. We find that a strongly fluctuating field is necessary for simulating synthetic synchrotron observations on small scales, we argue that Faraday rotation alone can account for the depolarisation of the synchrotron signal, and we discuss the importance of the observer position within the Milky Way. Altogether, we conclude that POLARIS is a highly reliable tool for predicting synchrotron emission and polarisation, including Faraday rotation in a realistic galactic context. It can thus contribute to better understand the results from current and future observational missions.

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.

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)

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.

Star formation is an inefficient process and in general only a small fraction of the gas in a giant molecular cloud (GMC) is turned into stars. This is partly due to the negative effect of stellar feedback from young massive star clusters. In a recent paper, we introduced a novel 1D numerical treatment of the effects of stellar feedback from young massive clusters on their natal clouds, which we named WARPFIELD. Here, we present version 2 of the WARPFIELD code, containing improved treatments of the thermal evolution of the gas and the fragmentation of the feedback-driven shell. As part of this update, we have produced new cooling and heating tables that account for the combined effects of photoionization and collisional ionization on the cooling rate of the gas, which we now make publically available. We employ our updated version of WARPFIELD to investigate the impact of stellar feedback on GMCs with a broad range of masses and surface densities and a variety of density profiles. We show that the minimum star formation efficiency $\epsilon_{\mathrm{min}}$, i.e. the star formation efficiency above which the cloud is destroyed by feedback and further star formation is shut off, is mainly set by the average cloud surface density. A star formation efficiency of 1-6 % is generally sufficient to destroy a GMC. We also find star formation efficiencies per free-fall time $\epsilon_{\mathrm{ff}} \sim 0.3$ %, in good agreement with recent observations. Our results imply that stellar feedback alone is sufficient to explain the low observed star formation efficiencies of GMCs. Finally, we show that very massive clouds with steep density profiles - possible proxies of the giant clumps observed in galaxies at $z \approx 2$ - are more resilient to feedback than typical GMCs, with $\epsilon_{\mathrm{min}}$ between 1 and 12 %.

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.

The C$^{+}$ ion is an important coolant of interstellar gas, and so the [CII] fine structure line is frequently observed in the interstellar medium. However, the physical and chemical properties of the [CII]-emitting gas are still unclear. We carry out non-LTE radiative transfer simulations with RADMC-3D to study the [CII] line emission from a young, turbulent molecular cloud before the onset of star formation, using data from the SILCC-Zoom project. The [CII] emission is optically thick over 40% of the observable area with $I_{[\textrm{CII}]} > 0.5$ K km s$^{-1}$. To determine the physical properties of the [CII] emitting gas, we treat the [CII] emission as optically thin. We find that the [CII] emission originates primarily from cold, moderate density gas ($40 \lesssim T \lesssim 65$ K and $50 \lesssim n \lesssim 440$ cm$^{-3}$), composed mainly of atomic hydrogen and with an effective visual extinction between $\sim 0.50$ and $\sim 0.91$. Gas dominated by molecular hydrogen contributes only $\lesssim$20% of the total [CII] line emission. Thus, [CII] is not a good tracer for CO-dark H$_2$ at this early phase in the cloud's lifetime. We also find that the total gas, H and C$^+$ column densities are all correlated with the integrated [CII] line emission, with power law slopes ranging from 0.5 to 0.7. Further, the median ratio between the total column density and the [CII] line emission is $Y_{{\rm CII}}\approx 1.1 \times 10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, and $Y_{{\rm CII}}$ scales with $I_{[\textrm{CII}]}^{-0.3}$. We expect $Y_{{\rm CII}}$ to change in environments with a lower or higher radiation field than simulated here.

We use the POLARIS radiative transport code to generate predictions of the two main observables directly sensitive to the magnetic field morphology and strength in filaments: dust polarization and gas Zeeman line splitting. We simulate generic gas filaments with power-law density profiles assuming two density-field strength dependencies, six different filament inclinations, and nine distinct magnetic field morphologies, including helical, toroidal, and warped magnetic field geometries. We present idealized spatially resolved dust polarization and Zeeman-derived field strengths and directions maps. Under the assumption that dust grains are aligned by radiative torques (RATs), dust polarization traces the projected plane-of-the-sky magnetic field morphology. Zeeman line splitting delivers simultaneously the intensity-weighted line-of-sight field strength and direction. We show that linear dust polarization alone is unable to uniquely constrain the 3D field morphology. We demonstrate that these ambiguities are ameliorated or resolved with the addition of the Zeeman directional information. Thus, observations of both the dust polarization and Zeeman splitting together provide the most promising means for obtaining constraints of the 3D magnetic field configuration. We find that the Zeeman-derived field strengths are at least a factor of a few below the input field strengths due to line-of-sight averaging through the filament density gradient. Future observations of both dust polarization and Zeeman splitting are essential for gaining insights into the role of magnetic fields in star and cluster forming filaments.

The Small Magellanic Cloud (SMC) provides the only laboratory to study the structure of molecular gas at high resolution and low metallicity. We present results from the Herschel Spectroscopic Survey of the SMC (HS$^{3}$), which mapped the key far-IR cooling lines [CII], [OI], [NII], and [OIII] in five star-forming regions, and new ALMA 7m-array maps of $^{12}$CO and $^{13}$CO $(2-1)$ with coverage overlapping four of the five HS$^{3}$ regions. We detect [CII] and [OI] throughout all of the regions mapped. The data allow us to compare the structure of the molecular clouds and surrounding photodissociation regions using $^{13}$CO, CO, [CII], and [OI] emission at $<10$" ($<3$ pc) scales. We estimate Av using far-IR thermal continuum emission from dust and find the CO/[CII] ratios reach the Milky Way value at high A$_{V}$ in the centers of the clouds and fall to $\sim{1/5-1/10}\times$ the Milky Way value in the outskirts, indicating the presence of translucent molecular gas not traced by bright CO emission. We estimate the amount of molecular gas traced by bright [CII] emission at low A$_{V}$ and bright CO emission at high A$_{V}$. We find that most of the molecular gas is at low A$_{V}$ and traced by bright [CII] emission, but that faint CO emission appears to extend to where we estimate the H$_{2}$-to-HI transition occurs. By converting our H$_{2}$ gas estimates to a CO-to-H$_{2}$ conversion factor ($X_{CO}$), we show that $X_{CO}$ is primarily a function of A$_{V}$, consistent with simulations and models of low metallicity molecular clouds.

${\bf Aim:}$ To test the hypothesis that radiation pressure from star clusters acting on dust is the dominant feedback agent disrupting the largest star-forming molecular clouds and thus regulating the star-formation process. ${\bf Methods:}$ We perform multi-frequency, 3D, RT calculations including scattering, absorption, and re-emission to longer wavelengths for clouds with masses of $10^4$-$10^7\,$M$_{\odot}$, with embedded clusters and a star formation efficiencies of 0.009%-91%, and varying maximum grain sizes up to 200$\,\mu$m. We calculate the ratio between radiative force and gravity to determine whether radiation pressure can disrupt clouds. ${\bf Results:}$ We find that radiation acting on dust almost never disrupts star-forming clouds. UV and optical photons to which the cloud is optically thick do not scatter much. Instead, they quickly get absorbed and re-emitted by at thermal wavelengths. As the cloud is typically optically thin to far-IR radiation, it promptly escapes, depositing little momentum. The resulting spectrum is more narrowly peaked than the corresponding Planck function with an extended tail at longer wavelengths. As the opacity drops significantly across the sub-mm and mm, the resulting radiative force is even smaller than for the corresponding single-temperature black body. The force from radiation pressure falls below the strength of gravitational attraction by an order of magnitude or more for either Milky Way or starbust conditions. For unrealistically large maximum grain sizes, and star formation efficiencies far exceeding 50% do we find that the strength of radiation pressure can exceed gravity. ${\bf Conclusions:}$ We conclude that radiation pressure acting on dust does not disrupt star-forming molecular clouds in any Local Group galaxies. Radiation pressure thus appears unlikely to regulate the star-formation process on either local or global scales.

SDSS-V will be an all-sky, multi-epoch spectroscopic survey of over six million objects. It is designed to decode the history of the Milky Way, trace the emergence of the chemical elements, reveal the inner workings of stars, and investigate the origin of planets. It will also create an integral-field spectroscopic map of the gas in the Galaxy and the Local Group that is 1,000x larger than the current state of the art and at high enough spatial resolution to reveal the self-regulation mechanisms of galactic ecosystems. SDSS-V will pioneer systematic, spectroscopic monitoring across the whole sky, revealing changes on timescales from 20 minutes to 20 years. The survey will thus track the flickers, flares, and radical transformations of the most luminous persistent objects in the universe: massive black holes growing at the centers of galaxies. The scope and flexibility of SDSS-V will be unique among extant and future spectroscopic surveys: it is all-sky, with matched survey infrastructures in both hemispheres; it provides near-IR and optical multi-object fiber spectroscopy that is rapidly reconfigurable to serve high target densities, targets of opportunity, and time-domain monitoring; and it provides optical, ultra-wide-field integral field spectroscopy. SDSS-V, with its programs anticipated to start in 2020, will be well-timed to multiply the scientific output from major space missions (e.g., TESS, Gaia, eROSITA) and ground-based projects. SDSS-V builds on the 25-year heritage of SDSS's advances in data analysis, collaboration infrastructure, and product deliverables. The project is now refining its science scope, optimizing the survey strategies, and developing new hardware that builds on the SDSS-IV infrastructure. We present here an overview of the current state of these developments as we seek to build our worldwide consortium of institutional and individual members.

The 30 Doradus Nebula in the Large Magellanic Cloud (LMC) contains the massive starburst cluster NGC 2070 with a massive and probably younger stellar sub clump at its center: R136. It is not clear how such a massive inner cluster could form several million years after the older stars in NGC 2070, given that stellar feedback is usually thought to expel gas and inhibit further star formation. Using the recently developed 1D feedback scheme WARPFIELD to scan a large range of cloud and cluster properties, we show that an age offset of several million years between the stellar populations is in fact to be expected given the interplay between feedback and gravity in a giant molecular cloud (GMC) with a density $\gtrsim 500$ cm$^{-3}$ due to re-accretion of gas onto the older stellar population. Neither capture of field stars nor gas retention inside the cluster have to be invoked in order to explain the observed age offset in NGC 2070 as well as the structure of the interstellar medium (ISM) around it.

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.

Star clusters interact with the interstellar medium (ISM) in various ways, most importantly in the destruction of molecular star-forming clouds, resulting in inefficient star formation on galactic scales. On cloud scales, ionizing radiation creates \hii regions, while stellar winds and supernovae drive the ISM into thin shells. These shells are accelerated by the combined effect of winds, radiation pressure and supernova explosions, and slowed down by gravity. Since radiative and mechanical feedback is highly interconnected, they must be taken into account in a self-consistent and combined manner, including the coupling of radiation and matter. We present a new semi-analytic one-dimensional feedback model for isolated massive clouds ($\geq 10^5\,M_{\odot}$) to calculate shell dynamics and shell structure simultaneously. It allows us to scan a large range of physical parameters (gas density, star formation efficiency, metallicity) and to estimate escape fractions of ionizing radiation $f_{\rm{esc,i}}$, the minimum star formation efficiency $\epsilon_{\rm{min}}$ required to drive an outflow, and recollapse time scales for clouds that are not destroyed by feedback. Our results show that there is no simple answer to the question of what dominates cloud dynamics, and that each feedback process significantly influences the efficiency of the others. We find that variations in natal cloud density can very easily explain differences between dense-bound and diffuse-open star clusters. We also predict, as a consequence of feedback, a $4-6$ Myr age difference for massive clusters with multiple generations.

The [CII] 158um fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, [CII] balances the heating, including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of [CII] emission remains unclear, because C+ can be found in multiple phases of the interstellar medium. Here we measure the fractions of [CII] emission originating in the ionized and neutral gas phases of a sample of nearby galaxies. We use the [NII] 205um fine-structure line to trace the ionized medium, thereby eliminating the strong density dependence that exists in the ratio of [CII]/[NII] 122um. Using the FIR [CII] and [NII] emission detected by the KINGFISH and Beyond the Peak Herschel programs, we show that 60-80% of [CII] emission originates from neutral gas. We find that the fraction of [CII] originating in the neutral medium has a weak dependence on dust temperature and the surface density of star formation, and a stronger dependence on the gas-phase metallicity. In metal-rich environments, the relatively cooler ionized gas makes substantially larger contributions to total [CII] emission than at low abundance, contrary to prior expectations. Approximate calibrations of this metallicity trend are provided.

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.

In this study we focus on the $z\sim6.6$ Lyman-$\alpha$ CR7 consisting of clump A that is host to a potential direct collapse black hole (DCBH), and two metal enriched star forming clumps B and C. In contrast to claims that signatures of metals rule out the existence of DCBHs, we show that metal pollution of A from star forming clumps clumps B and C is inevitable, and that A can form a DCBH well before its metallicity exceeds the critical threshold of $10^{-5}-10^{-6}\ \rm Z_{\odot}$. Assuming metal mixing happens instantaneously, we derive the metallicity of A based on the star formation history of B and C. We find that treating a final accreting black hole of $10^6-10^7\ \rm M_{\odot}$ in A for nebular emission already pushes its $H_{160}$ - [3.6] and [3.6]-[4.5] colours into the 3$\sigma$ limit of observations. Hence, we show that the presence of metals in DCBH hosts is inevitable, and that it is the coevolution of the LW radiation field and metals originating from neighbouring galaxies that governs DCBH formation in a neighbouring {initially} pristine atomic cooling haloes.

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.

We present three-dimensional radiation-hydrodynamical simulations of the impact of stellar winds, photoelectric heating, photodissociating and photoionising radiation, and supernovae on the chemical composition and star formation in a stratified disc model. This is followed with a sink-based model for star clusters with populations of individual massive stars. Stellar winds and ionising radiation regulate the star formation rate at a factor of ~10 below the simulation with only supernova feedback due to their immediate impact on the ambient interstellar medium after star formation. Ionising radiation (with winds and supernovae) significantly reduces the ambient densities for most supernova explosions to rho < 10^-25 g cm^-3, compared to 10^-23 g cm^-3 for the model with only winds and supernovae. Radiation from massive stars reduces the amount of molecular hydrogen and increases the neutral hydrogen mass and volume filling fraction. Only this model results in a molecular gas depletion time scale of 2 Gyr and shows the best agreement with observations. In the radiative models, the Halpha emission is dominated by radiative recombination as opposed to collisional excitation (the dominant emission in non-radiative models), which only contributes ~1-10 % to the total Halpha emission. Individual massive stars (M >= 30 M_sun) with short lifetimes are responsible for significant fluctuations in the Halpha luminosities. The corresponding inferred star formation rates can underestimate the true instantaneous star formation rate by factors of ~10.

CR7 is the brightest Lyman-$\alpha$ emitter observed at $z>6$, which shows very strong Lyman-$\alpha$ and HeII 1640\AA\ line luminosities, but no metal line emission. Previous studies suggest that CR7 hosts either young primordial stars with a total stellar mass of $\sim 10^7\,\mathrm{M}_\odot$ or a black hole of $\gtrsim 10^6\,\mathrm{M}_\odot$. Here, we explore different formation scenarios for CR7 with a semianalytical model, based on the random sampling of dark matter merger trees. We are unable to reproduce the observational constraints with a primordial stellar source, given our model assumptions, due to the short stellar lifetimes and the early metal enrichment. Black holes that are the remnants of the first stars are either not massive enough, or reside in metal-polluted haloes, ruling out this possible explanation of CR7. Our models instead suggest that direct collapse black holes, which form in metal-free haloes exposed to large Lyman-Werner fluxes, are more likely the origin of CR7. However, this result is derived under optimistic assumptions and future observations are necessary to further constrain the nature of CR7.

The [NII] 122 and 205 \mu m transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range, which corresponds to electron gas densities $n_e$~1-300 cm$^{-3}$, with a median value of $n_e$=30 cm$^{-3}$. Variations in the electron density within individual galaxies can be as a high as a factor of ~50, frequently with strong radial gradients. We find that $n_e$ increases as a function of infrared color, dust-weighted mean starlight intensity, and star formation rate surface density ($\Sigma_{SFR}$). As the intensity of the [NII] transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the star formation rate (SFR). We derive relations between the [NII] emission and SFR in the low-density limit and in the case of a log-normal distribution of densities. The scatter in the correlation between [NII] surface brightness and $\Sigma_{SFR}$ can be understood as a property of the $n_e$ distribution. For regions with $n_e$ close to or higher than the [NII] line critical densities, the low-density limit [NII]-based SFR calibration systematically underestimates the SFR since [NII] emission is collisionally quenched. Finally, we investigate the relation between [NII] emission, SFR, and $n_e$ by comparing our observations to predictions from the MAPPINGS-III code.

The Magellanic Clouds provide the only laboratory to study the effect of metallicity and galaxy mass on molecular gas and star formation at high (~20 pc) resolution. We use the dust emission from HERITAGE Herschel data to map the molecular gas in the Magellanic Clouds, avoiding the known biases of CO emission as a tracer of H$_{2}$. Using our dust-based molecular gas estimates, we find molecular gas depletion times of ~0.4 Gyr in the LMC and ~0.6 SMC at 1 kpc scales. These depletion times fall within the range found for normal disk galaxies, but are shorter than the average value, which could be due to recent bursts in star formation. We find no evidence for a strong intrinsic dependence of the molecular gas depletion time on metallicity. We study the relationship between gas and star formation rate across a range in size scales from 20 pc to ~1 kpc, including how the scatter in molecular gas depletion time changes with size scale, and discuss the physical mechanisms driving the relationships. We compare the metallicity-dependent star formation models of Ostriker, McKee, and Leroy (2010) and Krumholz (2013) to our observations and find that they both predict the trend in the data, suggesting that the inclusion of a diffuse neutral medium is important at lower metallicity.

We present a method to simultaneously model the dust far-infrared spectral energy distribution (SED) and the total infrared $-$ carbon monoxide (CO) integrated intensity $(S_{\rm IR}-I_{\rm CO})$ relationship. The modelling employs a hierarchical Bayesian (HB) technique to estimate the dust surface density, temperature ($T_{\rm eff}$), and spectral index at each pixel from the observed far-infrared (FIR) maps. Additionally, given the corresponding CO map, the method simultaneously estimates the slope and intercept between the FIR and CO intensities, which are global properties of the observed source. The model accounts for correlated and uncorrelated uncertainties, such as those present in Herschel observations. Using synthetic datasets, we demonstrate the accuracy of the HB method, and contrast the results with common non-hierarchical fitting methods. As an initial application, we model the dust and gas on 100 pc scales in the Magellanic Clouds from Herschel FIR and NANTEN CO observations. The slopes of the $\log S_{\rm IR}-\log I_{\rm CO}$ relationship are similar in both galaxies, falling in the range 1.1$-$1.7. However, in the SMC the intercept is nearly 3 times higher, which can be explained by its lower metallicity than the LMC, resulting in a larger $S_{\rm IR}$ per unit $I_{\rm CO}$. The HB modelling evidences an increase in $T_{\rm eff}$ in regions with the highest $I_{\rm CO}$ in the LMC. This may be due to enhanced dust heating in the densest molecular regions from young stars. Such simultaneous dust and gas modelling may reveal variations in the properties of the ISM and its association with other galactic characteristics, such as star formation rates and/or metallicities.

To compute the SFR of galaxies from the rest-frame UV it is essential to take into account the obscuration by dust. To do so, one of the most popular methods consists in combining the UV with the emission from the dust itself in the IR. Yet, different studies have derived different estimators, showing that no such hybrid estimator is truly universal. In this paper we aim at understanding and quantifying what physical processes drive the variations between different hybrid estimators. Doing so, we aim at deriving new universal UV+IR hybrid estimators to correct the UV for dust attenuation, taking into account the intrinsic physical properties of galaxies. We use the CIGALE code to model the spatially-resolved FUV to FIR SED of eight nearby star-forming galaxies drawn from the KINGFISH sample. This allows us to determine their local physical properties, and in particular their UV attenuation, average SFR, average specific SFR (sSFR), and their stellar mass. We then examine how hybrid estimators depend on said properties. We find that hybrid UV+IR estimators strongly depend on the stellar mass surface density (in particular at 70 and 100 micron) and on the sSFR (in particular at 24 micron and the TIR). Consequently, the IR scaling coefficients for UV obscuration can vary by almost an order of magnitude. This result contrasts with other groups who found relatively constant coefficients with small deviations. We exploit these variations to construct a new class of hybrid estimators based on observed UV to near-IR colours and near-IR luminosity densities per unit area. We find that they can reliably be extended to entire galaxies. The new estimators provide better estimates of attenuation-corrected UV emission than classical hybrid estimators. Naturally taking into account the variable impact of dust heated by old stellar populations, they constitute a step towards universal estimators.

More complete knowledge of galaxy evolution requires understanding the process of star formation and interaction between the interstellar radiation field and the interstellar medium in galactic environments traversing a wide range of physical parameter space. Here we focus on the impact of massive star formation on the surrounding low metallicity ISM in 30 Doradus in the Large Magellanic Cloud. A low metal abundance, as is the case of some galaxies of the early universe, results in less ultra-violet shielding for the formation of the molecular gas necessary for star formation to proceed. The half-solar metallicity gas in this region is strongly irradiated by the super star cluster R136, making it an ideal laboratory to study the structure of the ISM in an extreme environment. Our spatially resolved study investigates the gas heating and cooling mechanisms, particularly in the photo-dissociation regions where the chemistry and thermal balance are regulated by far-ultraviolet photons (6 eV< h\nu <13.6 eV). We present Herschel observations of far-infrared fine-structure lines obtained with PACS and SPIRE/FTS. We have combined atomic fine-structure lines from Herschel and Spitzer observations with ground-based CO data to provide diagnostics on the properties and the structure of the gas by modeling it with the Meudon PDR code. We derive the spatial distribution of the radiation field, the pressure, the size, and the filling factor of the photodissociated gas and molecular clouds. We find a range of pressure of ~ 10^5 - 1.7x10^6 cm^{-3} K and a range of incident radiation field G_UV ~ 10^2 - 2.5x10^4 through PDR modeling. Assuming a plane-parallel geometry and a uniform medium, we find a total extinction of 1-3 mag , which correspond to a PDR cloud size of 0.2 to 3pc, with small CO depth scale of 0.06 to 0.5pc. We also determine the three dimensional structure of the gas. (Abridged)

We combine Spitzer and Herschel data of the star-forming region N11 in the Large Magellanic Cloud to produce detailed maps of the dust properties in the complex and study their variations with the ISM conditions. We also compare APEX/LABOCA 870um observations with our model predictions in order to decompose the 870um emission into dust and non-dust (free-free emission and CO(3-2) line) contributions. We find that in N11, the 870um can be fully accounted for by these 3 components. The dust surface density map of N11 is combined with HI and CO observations to study local variations in the gas-to-dust mass ratios. Our analysis leads to values lower than those expected from the LMC low-metallicity as well as to a decrease of the gas-to-dust mass ratio with the dust surface density. We explore potential hypotheses that could explain the low observed gas-to-dust mass ratios (variations in the XCO factor, presence of CO-dark gas or of optically thick HI or variations in the dust abundance in the dense regions). We finally decompose the local SEDs using a Principal Component Analysis (i.e. with no a priori assumption on the dust composition in the complex). Our results lead to a promising decomposition of the local SEDs in various dust components (hot, warm, cold) coherent with that expected for the region. Further analysis on a larger sample of galaxies will follow in order to understand how unique this decomposition is or how it evolves from one environment to another.

We present the first fully calibrated H$_2$, 1-0 S(1) image of the entire 30 Doradus nebula. The observations were conducted using the NOAO Extremely Wide-Field Infrared Imager on the CTIO 4-meter Blanco Telescope. Together with a NEWFIRM Br$\gamma$ image of 30 Doradus, our data reveal the morphologies of the warm molecular gas and ionized gas in 30 Doradus. The brightest H$_2$-emitting area, which extends from the northeast to the southwest of R136, is a photodissociation region viewed face-on, while many clumps and pillar features located at the outer shells of 30 Doradus are photodissociation regions viewed edge-on. Based on the morphologies of H$_2$, Br$\gamma$, $^{12}$CO, and 8$\mu$m emission, the H$_2$ to Br$\gamma$ line ratio and Cloudy models, we find that the H$_2$ emission is formed inside the photodissociation regions of 30 Doradus, 2 - 3 pc to the ionization front of the HII region, in a relatively low-density environment $<$ 10$^4$ cm$^{-3}$. Comparisons with Br$\gamma$, 8$\mu$m, and CO emission indicate that H$_2$ emission is due to fluorescence, and provide no evidence for shock excited emission of this line.

We analyze the applicability of far-infrared fine-structure lines [CII] 158 micron, [OI] 63 micron and [OIII] 88 micron to reliably trace the star formation rate (SFR) in a sample of low-metallicity dwarf galaxies from the Herschel Dwarf Galaxy Survey and compare with a broad sample of galaxies of various types and metallicities in the literature. We study the trends and scatter in the relation between the SFR (as traced by GALEX FUV and MIPS 24 micron) and far-infrared line emission, on spatially resolved and global galaxy scales, in dwarf galaxies. We assemble far-infrared line measurements from the literature and infer whether the far-infrared lines can probe the SFR (as traced by the total-infrared luminosity) in a variety of galaxy populations. In metal-poor dwarfs, the [OI] and [OIII] lines show the strongest correlation with the SFR with an uncertainty on the SFR estimates better than a factor of 2, while the link between [CII] emission and the SFR is more dispersed (uncertainty factor of 2.6). The increased scatter in the SFR-L([CII]) relation towards low metal abundances, warm dust temperatures, large filling factors of diffuse, highly ionized gas suggests that other cooling lines start to dominate depending on the density and ionization state of the gas. For the literature sample, we evaluate the correlations for a number of different galaxy populations. The [CII] and [OI] lines are considered to be reliable SFR tracers in starburst galaxies, recovering the star formation activity within an uncertainty of factor 2. [Abridged]

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.

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.

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.

The metal content of a galaxy is a key property for distinguishing between viable galaxy evolutionary scenarios, and it strongly influences many of the physical processes in the interstellar medium. An absolute and robust determination of extragalactic metallicities is essential in constraining models of chemical enrichment and chemical evolution, however, current gas phase abundance determinations from optical fine-structure lines are uncertain to 0.8 dex as conversion of these optical line fluxes to abundances is strongly dependent on the electron temperature of the ionized gas. In contrast, the far-IR emission lines can be used to derive an O++ abundance that is relatively insensitive to temperature, while the ratio of the optical to far-IR lines provides a consistent temperature to be used in the derivation of an O$^+$ abundance. We present observations of the [O III] 88 \micron\ fine-structure line in NGC 628 that were obtained as part of the KINGFISH program. These data are combined with optical IFU data to derive oxygen abundances for seven HII regions. We find the abundance of these regions to all lie between the high and low values of strong line calibrations and in agreement with estimates that assume temperature fluctuations are present in the HII regions.

We present kiloparsec (kpc) spatial resolution maps of the CO-to-H2 conversion factor (alpha_co) and dust-to-gas ratio (DGR) in 26 nearby, star-forming galaxies. We have simultaneously solved for alpha_co and DGR by assuming that the DGR is approximately constant on kpc scales. With this assumption, we can combine maps of dust mass surface density, CO integrated intensity and HI column density to solve for both alpha_co and DGR with no assumptions about their value or dependence on metallicity or other parameters. Such a study has just become possible with the availability of high resolution far-IR maps from the Herschel key program KINGFISH, 12CO J=(2-1) maps from the IRAM 30m large program HERACLES and HI 21-cm line maps from THINGS. We use a fixed ratio between the (2-1) and (1-0) lines to present our alpha_co results on the more typically used 12CO J=(1-0) scale and show using literature measurements that variations in the line ratio do not effect our results. In total, we derive 782 individual solutions for alpha_co and DGR. On average, alpha_co = 3.1 Msun pc^-2 (K km s^-1)^-1 for our sample with a standard deviation of 0.3 dex. Within galaxies we observe a generally flat profile of alpha_co as a function of galactocentric radius. However, most galaxies exhibit a lower alpha_co in the central kpc---a factor of ~2 below the galaxy mean, on average. In some cases, the central alpha_co value can be factors of 5 to 10 below the standard Milky Way (MW) value of alpha_co,MW =4.4 Msun pc^-2 (K km s^-1)^-1. While for alpha_co we find only weak correlations with metallicity, DGR is well-correlated with metallicity, with an approximately linear slope. Finally, we present several recommendations for choosing an appropriate alpha_co for studies of nearby galaxies.

We present observations of the nearby spiral galaxy IC342 with the Herschel Spectral and Photometric Imaging Receiver (SPIRE) Fourier Transform Spectrometer. The spectral range afforded by SPIRE, 196-671 microns, allows us to access a number of 12CO lines from J=4--3 to J=13--12 with the highest J transitions observed for the first time. In addition we present measurements of 13CO, [CI] and [NII]. We use a radiative transfer code coupled with Bayesian likelihood analysis to model and constrain the temperature, density and column density of the gas. We find two 12CO components, one at 35 K and one at 400 K with CO column densities of 6.3x10^{17} cm^{-2} and 0.4x10^{17} cm^{-2} and CO gas masses of 1.26x10^{7} Msolar and 0.15x10^{7} Msolar, for the cold and warm components, respectively. The inclusion of the high-J 12CO line observations, indicate the existence of a much warmer gas component (~400 K) confirming earlier findings from H_{2} rotational line analysis from ISO and Spitzer. The mass of the warm gas is 10% of the cold gas, but it likely dominates the CO luminosity. In addition, we detect strong emission from [NII] 205microns and the {3}P_{1}->{3}P_{0} and {3}P_{2} ->{3}P_{1} [CI] lines at 370 and 608 microns, respectively. The measured 12CO line ratios can be explained by Photon-dominated region (PDR) models although additional heating by e.g. cosmic rays cannot be excluded. The measured [CI] line ratio together with the derived [C] column density of 2.1x10^{17} cm^{-2} and the fact that [CI] is weaker than CO emission in IC342 suggests that [CI] likely arises in a thin layer on the outside of the CO emitting molecular clouds consistent with PDRs playing an important role.

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.

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.

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.

We report the discovery of a new candidate ultraluminous X-ray source (ULX) in the nearby edge-on spiral galaxy NGC 891. The source, which has an absorbed flux of F_X ~ 10^-12 erg/s/cm^2 (corresponding to L_X > 10^40 erg/s at 9 Mpc), must have begun its outburst in the past 5 years as it is not detected in prior X-ray observations between 1986 and 2006. We try empirical fits to the XMM-Newton spectrum, finding that the spectrum is fit very well as emission from a hot disk, a cool irradiated disk, or blurred reflection from the innermost region of the disk. The simplest physically motivated model with an excellent fit is a hot disk around a stellar-mass black hole (a super-Eddington outburst), but equally good fits are found for each model. We suggest several follow-up experiments that could falsify these models.

We use our new optical-imaging and spectrophotometric survey of key diagnostic emission lines in 30 Doradus, together with CLOUDY photoionization models, to study the physical conditions and ionization mechanisms along over 4000 individual lines of sight at points spread across the face of the extended nebula, out to a projected radius 75 pc from R136 at the center of the ionizing cluster NGC 2070. We focus on the physical conditions, geometry and importance of radiation pressure on a point-by-point basis, with the aim of setting observational constraints on important feedback processes. We find that the dynamics and large scale structure of 30 Dor are set by a confined system of X-ray bubbles in rough pressure equilibrium with each other and with the confining molecular gas. Although the warm (10,000K) gas is photoionized by the massive young stars in NGC 2070, the radiation pressure does not currently play a major role in shaping the overall structure. The completeness of our survey also allows us to create a composite spectrum of 30 Doradus, simulating the observable spectrum of a spatially-unresolved, distant giant extragalactic H II region. We find that the highly simplified models used in the "strong line" abundance technique do in fact reproduce our observed lines strengths and deduced chemical abundances, in spite of the more than one order of magnitude range in the ionization parameter and density of the actual gas in 30 Dor.

We have carried out a narrow-band survey of the Local Group galaxy, M33, in the HeII4686 emission line, to identify HeII nebulae in this galaxy. With spectroscopic follow-up observations, we confirm three of seven candidate objects, including identification of two new HeII nebulae, BCLMP651, HBW673. We also obtain spectra of associated ionizing stars for all the HII regions, identifying two new WN stars. We demonstrate that the ionizing source for the known HeII nebula, MA 1, is consistent with being the early-type WN star MC8 (M33-WR14), by carrying out a combined stellar and nebular analysis of MC8 and MA1. We were unable to identify the helium ionizing sources for HBW 673 and BCLMP 651, which do not appear to be Wolf-Rayet stars. According to the [OIII]5007/Hbeta vs [NII]6584/Halpha diagnostic diagram, excitation mechanisms apart from hot stellar continuum are needed to account for the nebular emission in HBW 673, which appears to have no stellar source at all.

We have completed a a new optical imaging and spectrophotometric survey of a 140 x 80 pc$^2$ region of 30 Doradus centered on R136, covering key optical diagnostic emission lines including \Ha, \Hb, \Hg, [O III] $\lambda\lambda$4363, 4959, 5007, [N II] $\lambda\lambda$6548, 6584, [S II] $\lambda\lambda$6717, 6731 [S III] $\lambda $6312 and in some locations [S III] $\lambda$9069. We present maps of fluxes and intensity ratios for these lines, and catalogs of isolated ionizing stars, elephant-trunk pillars, and edge-on ionization fronts. The final science-quality spectroscopic data products are available to the public. Our analysis of the new data finds that, while stellar winds and supernovae undoubtedly produce shocks and are responsible for shaping the nebula, there are no global spectral signatures to indicate that shocks are currently an important source of ionization. We conclude that the considerable region covered by our survey is well described by photoionization from the central cluster where the ionizing continuum is dominated by the most massive O stars. We show that if 30 Dor were viewed at a cosmological distance, its integrated light would be dominated by its extensive regions of lower surface-brightness rather than by the bright, eye-catching arcs.

Light curves in the B, V, and I_c passbands have been obtained for the type II Cepheids V154 in M3 and V42 and V84 in M5. Alternating cycle behavior, similar to that seen among RV Tauri variables, is confirmed for V84. Old and new observations, spanning more than a century, show that V154 has increased in period while V42 has decreased in period. V84, on the other hand, has shown large, erratic changes in period that do not appear to reflect the long term evolution of V84 through the HR diagram.

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.

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.

Emission-line abundances have been uncertain for more than a decade due to unexplained discrepancies in the relative intensities of the forbidden lines and weak permitted recombination lines in planetary nebulae (PNe) and H II regions. The observed intensities of forbidden and recombination lines originating from the same parent ion differ from their theoretical values by factors of more than an order of magnitude in some of these nebulae. In this study we observe UV resonance line absorption in the central stars of PNe produced by the nebular gas, and from the same ions that emit optical forbidden lines. We then compare the derived absorption column densities with the emission measures determined from ground-based observations of the nebular forbidden lines. We find for our sample of PNe that the collisionally excited forbidden lines yield column densities that are in basic agreement with the column densities derived for the same ions from the UV absorption lines. A similar comparison involving recombination line column densities produces poorer agreement, although near the limits of the formal uncertainties of the analyses. An additional sample of objects with larger abundance discrepancy factors will need to be studied before a stronger statement can be made that recombination line abundances are not correct.

The s-process should occur in all but the lower mass progenitor stars of planetary nebulae, and this should be reflected in the chemical composition of the gas which is expelled to create the current planetary nebula shell. Weak forbidden emission lines are expected from several s-process elements in these shells, and have been searched for and in some cases detected in previous investigations. Here we extend these studies by combining very high signal-to-noise echelle spectra of a sample of PNe with a critical analysis of the identification of the emission lines of Z>30 ions. Emission lines of Br, Kr, Xe, Rb, Ba, and Pb are detected with a reasonable degree of certainty in at least some of the objects studied here, and we also tentatively identify lines from Te and I, each in one object. The strengths of these lines indicate enhancement of s-process elements in the central star progenitors, and we determine the abundances of Br, Kr, and Xe, elements for which atomic data relevant for abundance determination have recently become available. As representative elements of the ``light'' and ``heavy'' s-process peaks Kr and Xe exhibit similar enhancements over solar values, suggesting that PNe progenitors experience substantial neutron exposure.

We identify emission lines of post-iron peak elements in very high signal-to-noise spectra of a sample of planetary nebulae. Analysis of lines from ions of Kr and Xe reveals enhancements in most of the PNe, in agreement with the theories of s-process in AGB star. Surprisingly, we did not detect lines from Br even though s-process calculations indicate that it should be produced with Kr at detectable levels.

Nuclei with magic numbers serve as important benchmarks in nuclear theory. In addition, neutron-rich nuclei play an important role in the astrophysical rapid neutron-capture process (r-process). 78Ni is the only doubly-magic nucleus that is also an important waiting point in the r-process, and serves as a major bottleneck in the synthesis of heavier elements. The half-life of 78Ni has been experimentally deduced for the first time at the Coupled Cyclotron Facility of the National Superconducting Cyclotron Laboratory at Michigan State University, and was found to be 110 (+100 -60) ms. In the same experiment, a first half-life was deduced for 77Ni of 128 (+27 -33) ms, and more precise half-lives were deduced for 75Ni and 76Ni of 344 (+20 -24) ms and 238 (+15 -18) ms respectively.