We used the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST) to observe the semi-forbidden CIII] emission in Green Pea galaxies at 0.13 < z < 0.3. We detect CIII] emission in 7/10 galaxies with CIII] equivalent widths that range from 2-10\AA~. The observed CIII] emission line strengths are consistent with the predictions from photoionization models which incorporate the effects of binary stellar evolution with young stellar ages < 3-5 Myrs, and high ionization parameters (logU > -2). The hard ionizing radiation from young massive stars, and high nebular temperatures at low-metallicities can account for the observed high equivalent widths of CIII] and [OIII] emission lines. The Green Pea galaxies do not show a significant correlation between the Ly$\alpha$ and CIII] equivalent widths, and the observed scatter is likely due to the variations in the optical depth of Ly$\alpha$ to the neutral gas. Green Pea galaxies are likely to be density-bounded, and we examined the dependence of CIII] emission on the Lyman continuum optical depth. The potential LyC leaker galaxies in our sample have high CIII] equivalent widths that can only be reproduced by starburst ages as young as < 3 Myrs and harder ionizing spectra than the non-leakers. Among the galaxies with similar metallicities and ionization parameters, the CIII] equivalent width appears to be stronger for those with higher optical depth to LyC, as expected from the photoionization models. Further investigation of a larger sample of CIII]-emitters is necessary to calibrate the dependence of CIII] emission on the escape of LyC radiation, and to enable application of the CIII] diagnostics to galaxies in the reionization epoch.

We investigate the ionised gas kinematics, physical properties and chemical abundances of SDSS$J$142947, a Green Pea galaxy at redshift z$\sim$0.17 with strong, double-peak Ly$\alpha$ emission and indirect evidence of Lyman continuum (LyC) leakage. Using high-dispersion spectroscopy, we perform a multi-component analysis of emission-line profiles. Our model consistently fits all lines as a narrow component with intrinsic velocity dispersion $\sigma \sim$40 km s$^{-1}$, and two broader blue-shifted components with $\sigma \sim$90 km s$^{-1}$ and $\sim$250 km s$^{-1}$. We find electron densities and temperatures, ionisation conditions, and direct O/H and N/O abundances for each component. A highly ionised, metal-poor, young and compact starburst dominates narrow emission, showing evidence of hard radiation fields and elevated N/O. The blue-shifted broader components are consistent with highly turbulent, possibly clumpy ionised gas at the base of a strong photoionised outflow, which accounts for $\gtrsim$50% of the integrated emission-line fluxes. The outflow is dense and metal-enriched compared to the HII regions, with expansion velocities larger than those obtained from UV interstellar absorption lines under standard assumptions. Some of these metals may be able to escape, with outflows loading factors comparable to those found in high-$z$ galaxies of similar SFR/Area. Our findings depict a two-stage starburst picture; hard radiation fields from young star clusters illuminate a turbulent and clumpy ISM that has been eroded by SNe feedback. While UV data suggest an extended Ly$\alpha$ halo with high average HI column density, LyC photons could only escape from SDSS$J$142947 through low HI density channels or filaments in the ISM approaching density-bounded conditions, traced by outflowing gas.

We report on the HST detection of the Lyman-continuum (LyC) radiation emitted by a galaxy at redshift z=3.794, dubbed Ion1 (Vanzella et al. 2012). The LyC from Ion1 is detected at rest-frame wavelength 820$\sim$890 \AA with HST WFC3/UVIS in the F410M band ($m_{410}=27.60\pm0.36$ magnitude (AB), peak SNR = 4.17 in a circular aperture with radius r = 0.12'') and at 700$\sim$830 \AA with the VLT/VIMOS in the U-band ($m_U = 27.84\pm0.19$ magnitude (AB), peak SNR = 6.7 with a r = 0.6'' aperture). A 20-hr VLT/VIMOS spectrum shows low- and high-ionization interstellar metal absorption lines, the P-Cygni profile of CIV and Ly$\alpha$ in absorption. The latter spectral feature differs from what observed in known LyC emitters, which show strong Ly$\alpha$ emission. An HST far-UV color map reveals that the LyC emission escapes from a region of the galaxy that is bluer than the rest, presumably because of lower dust obscuration. The F410M image shows that the centroid of the LyC emission is offset from the centroid of the non-ionizing UV emission by 0.12''$\pm$0.03'', corresponding to 0.85$\pm$0.21 kpc (physical), and that its morphology is likely moderately resolved. These morphological characteristics favor a scenario where the LyC photons produced by massive stars escape from low HI column-density "cavities" in the ISM, possibly carved by stellar winds and/or supernova. We also collect the VIMOS U-band images of a sample of 107 Lyman-break galaxies with spectroscopic redshifts at $3.40<z<3.95$, i.e. sampling the LyC, and stack them with inverse-variance weights. No LyC emission is detected in the stacked image, resulting in a 32.5 magnitude (AB) flux limit (1$\sigma$) and an upper limit of absolute LyC escape fraction $f_{esc}^{abs} < 0.63\%$. LyC emitters like Ion1 are very likely at the bright-end of the LyC luminosity function.

As some of the only Lyman continuum (LyC) emitters at z~0, Green Pea (GP) galaxies are possible analogs of the sources that reionized the universe. We present HST COS spectra of 13 of the most highly ionized GPs, with [O III]/[O II]=6-35, and investigate correlations between Ly-alpha, galaxy properties, and low-ionization UV lines. Galaxies with high [O III]/[O II] have higher H-alpha equivalent widths (EWs), and high intrinsic Ly-alpha production may explain the prevalence of high Ly-alpha EWs among GPs. While Ly-alpha escape fraction is closely linked to low gas covering fractions, implying a clumpy gas geometry, narrow Ly-alpha velocity peak separation (delta_v,LyA) correlates with the ionization state, suggesting a density-bounded geometry. We therefore suggest that delta_v,LyA may trace the residual transparency of low-column-density pathways. Metallicity is associated with both [O III]/[O II] and delta_v,LyA. This trend may result from catastrophic cooling around low-metallicity star clusters, which generates a compact geometry of dense clouds within a low-density inter-clump medium. We find that the relative strength of low-ionization UV emission to absorption correlates with Ly-alpha emission strength and is related to Ly-alpha profile shape. However, as expected for optically thin objects, the GPs with the lowest delta_v,LyA show both weak low-ionization emission and weak absorption. The strengths of the low-ionization absorption and emission lines in a stacked spectrum do not correspond to any individual spectrum. Galaxies with high [O III]/[O II] contain a high fraction of LyC emitter candidates, but [O III]/[O II] alone is an insufficient diagnostic of LyC escape.

Quantifying the physical conditions that allow radiation emitted shortward of the hydrogen ionization edge at 911.7 {\AA} to escape the first collapsed objects and ultimately reionize the universe is a compelling problem for astrophysics. The escape of LyC emission from star-forming galaxies and AGN is intimately tied to the emergence and sustenance of the metagalactic ionizing background that pervades the universe to the present day and in turn is tied to the emergence of structure at all epochs. JWST was built in part to search for the source(s) responsible for reionization, but it cannot observe LyC escape directly, because of the progressive increase in the mean transmission of the intergalactic medium towards the epoch of reionization. Remarkable progress has been made to date in directly detecting LyC leaking from star-forming galaxies using space-based and the ground-based observatories, but there remain significant gaps in our redshift coverage of the phenomenon. Ongoing projects to measure LyC escape at low- and intermediate-z will provide guidance to JWST investigations by analyzing the robustness of a set of proposed LyC escape proxies, and also provide a closeup examination of the physical conditions that favor LyC escape. However, currently available facilities are inadequate for deeply probing LyC escape at the faint end of the galaxy luminosity function. Doing so will require facilities that can detect LyC emission in the restframe to limiting magnitudes approaching 28 $< m^*_{(1+z)900} <$ 32 for $M^*_{(1+z)1500}$ galaxies. The goal of acquiring statistically robust samples for determining LyC luminosity functions across cosmic time will require multi-object spectroscopy from spacebased flagship class and groundbased ELT class telescopes along with ancillary panchromatic imaging and spectroscopy spanning the far-UV to the mid-IR.

Diagnosing the physical and chemical conditions within star-forming galaxies (SFGs) is of paramount importance to understanding key components of galaxy formation and evolution: star-formation, gas enrichment, outflows, and accretion. Well established optical emission-line diagnostics used to discern such properties (i.e., metal content, density, strength/shape of ionizing radiation) will be observationally inaccessible for the earliest galaxies, emphasizing the need for robust, reliable interstellar medium (ISM) diagnostics at ultraviolet (UV) wavelengths. Calibrating these UV diagnostics requires a comprehensive comparison of the UV and optical emission lines in nearby SFGs. Optical integral field unit (IFU) surveys have revealed the inhomogeneous nature of the ISM in SFGs, which leads to non-systematic biases in the interpretation of unresolved sources. Spatial variations are especially important to consider at UV wavelengths, where the strongest emission features originate from only the highest excitation regions of the nebula and are challenging to distinguish from competing high-ionization sources (e.g., shocks, AGN, etc.). Since surveys collecting large-scale optical integral field unit (IFU) spectroscopy are already underway, this white paper calls for an IFU or multi-object far-UV (FUV) spectroscopic instrument with high sensitivity, high spatial resolution, and large field of view (FoV). Given the impact of large-scale optical IFU surveys over the past decade, this white paper emphasizes the scientific need for a comparable foundation of spatially-resolved far-UV spectroscopy survey of nearby galaxies that will lay the foundation of diagnostics critical to the interpretation of the distant universe.

Mechanisms regulating the escape of Ly$\alpha$ photons and ionizing radiation remain poorly understood. To study these processes we analyze VLA 21cm observations of one Green Pea (GP), J160810+352809 (hereafter J1608), and HST COS spectra of 17 GP galaxies at $z<0.2$. All are highly ionized: J1608 has the highest [O III] $\lambda5007$/[O II] $\lambda3727$ for star-forming galaxies in SDSS, and the 17 GPs have [O III]/[O II] $\geq6.6$. We set an upper limit on J1608's HI mass of $\log M_{HI}/M_\odot=8.4$, near or below average compared to similar mass dwarf galaxies. In the COS sample, eight GPs show Ly$\alpha$ absorption components, six of which also have Ly$\alpha$ emission. The HI column densities derived from Ly$\alpha$ absorption are high, $\log N_{HI}/$cm$^{-2}=19-21$, well above the LyC optically thick limit. Using low-ionization absorption lines, we measure covering fractions ($f_{\mbox{cov}}$) of $0.1-1$, and find that $f_{\mbox{cov}}$ strongly anti-correlates with Ly$\alpha$ escape fraction. Low covering fractions may facilitate Ly$\alpha$ and LyC escape through dense neutral regions. GPs with $f_{\mbox{cov}}\sim1$ all have low neutral gas velocities, while GPs with lower $f_{\mbox{cov}}=0.2-0.6$ have a larger range of velocities. Conventional mechanical feedback may help establish low $f_{\mbox{cov}}$ in some cases, whereas other processes may be important for GPs with low velocities. Finally, we compare $f_{\mbox{cov}}$ with proposed indicators of LyC escape. Ionizing photon escape likely depends on a combination of neutral gas geometry and kinematics, complicating the use of emission-line diagnostics for identifying LyC emitters.

Low redshift, spatially resolved Lyman continuum (LyC) emitters allow us to clarify the processes for LyC escape from these starburst galaxies. We use Hubble Space Telescope (HST) WFC3 and ACS imaging of the confirmed low-redshift LyC emitter Tol 1247-232 to study the ionization structure of the gas and its relation to the ionizing star clusters. We perform ionization parameter mapping (IPM) using [O III]4959, 5007 and [O II]3727 imaging as the high- and low-ionization tracers, revealing broad, large-scale, optically thin regions originating from the center, and reaching the outskirts of the galaxy, consistent with LyC escape. We carry out stellar population synthesis modeling of the 26 brightest clusters using our HST photometry. Combining these data with the nebular photometry, we find a global LyC escape fraction of f_esc = 0.12, with uncertainties also consistent with zero escape and with all measured f_esc values for this galaxy. Our analysis suggests that, similar to other candidate LyC emitters, a two-stage starburst has taken place in this galaxy, with a 12 Myr old, massive, central cluster likely having pre-cleared regions in and around the center, and the second generation of 2 - 4 Myr old clusters dominating the current ionization, including some escape from the galaxy.

Recent observations have discovered an escape of ionizing Lyman-continuum radiation from a population of compact, actively star-forming galaxies at redshift z~0.3, also known as "green peas". We here study the peculiar green pea Lyman-alpha (Lya) spectra, which are mostly double-peaked, unlike in any other galaxy sample. We select a sample of twelve archival green peas and we apply numerical radiative transfer models to reproduce the observed Lya spectral profiles, using the geometry of expanding, homogeneous spherical shells. We use ancillary optical and ultraviolet data to constrain the model parameters, and we evaluate the match between the models and the observed Lya spectra. As a second step, we allow all the fitting parameters to be free, and examine the agreement between the interstellar medium parameters derived from the models and those from ancillary data. The green pea double-peaked Lya line profiles are not correctly reproduced by the constrained shell models. Conversely, unconstrained models fit the spectra, but parameters derived from the best-fitting models are not in agreement with the ancillary data. In particular: 1) the best-fit systemic redshifts are larger by 10 - 250 km/s than those derived from optical emission lines, 2) the double-peaked Lya profiles are best reproduced with low-velocity (<150 km/s) outflows that contradict the observed ultraviolet absorption lines of low-ionization-state elements with characteristic velocities as large as 300 km/s, and 3) the models need to consider intrinsic Lya profiles that are on average three times broader than the observed Balmer lines. Differences between the modelled and observed velocities are larger for targets with prominent Lya blue peaks. The blue peak position and flux appear to be connected to low column densities of neutral hydrogen, leading to Lya and Lyman-continuum escape.

By clearing neutral gas away from a young starburst, superwinds may regulate the escape of Lyman continuum (LyC) photons from star-forming galaxies. However, models predict that superwinds may not launch in the most extreme, compact starbursts. We explore the role of outflows in generating low optical depths in the Green Peas (GPs), the only known star-forming population with several confirmed and candidate LyC-leaking galaxies. With Hubble Space Telescope UV spectra of 25 low-redshift GPs, including new observations of 13 of the most highly ionized GPs, we compare the kinematics of UV absorption lines with indirect HI optical depth diagnostics: Ly-alpha escape fraction, Ly-alpha peak separation, or low-ionization absorption line equivalent width. The data suggest that high ionization kinematics tracing superwind activity may correlate with low optical depth in some objects. However, the most extreme GPs, including many of the best candidate LyC emitters with weak low-ionization absorption and strong, narrow Ly-alpha profiles, show the lowest velocities. These results are consistent with models for suppressed superwinds, which suggests that outflows may not be the only cause of LyC escape from galaxies.

We report the detection of CO(J=2-1) coincident with the super star cluster (SSC) Mrk 71-A in the nearby Green Pea analog galaxy, NGC 2366. Our NOEMA observations reveal a compact, ~7 pc, molecular cloud whose mass (10^5 M_sun) is similar to that of the SSC, consistent with a high star-formation efficiency, on the order of 0.5. There are two, spatially distinct components separated by 11 km/s. If expanding, these could be due to momentum-driven, stellar wind feedback. Alternatively, we may be seeing the remnant infalling, colliding clouds responsible for triggering the SSC formation. The kinematics are also consistent with a virialized system. These extreme, high-density, star-forming conditions inhibit energy-driven feedback; the co-spatial existence of a massive, molecular cloud with the SSC supports this scenario, and we quantitatively confirm that any wind-driven feedback in Mrk 71-A is momentum-driven, rather than energy-driven. Since Mrk 71-A is a candidate Lyman continuum emitter, this implies that energy-driven superwinds may not be a necessary condition for the escape of ionizing radiation. In addition, the detection of the nebular continuum emission yields an accurate astrometric position for the Mrk 71-A. We also detect four other massive, molecular clouds in this giant star-forming complex.

Identifying the mechanism by which high energy Lyman continuum (LyC) photons escaped from early galaxies is one of the most pressing questions in cosmic evolution. Haro 11 is the best known local LyC leaking galaxy, providing an important opportunity to test our understanding of LyC escape. The observed LyC emission in this galaxy presumably originates from one of the three bright, photoionizing knots known as A, B, and C. It is known that Knot C has strong Ly$\alpha$ emission, and Knot B hosts an unusually bright ultraluminous X-ray source, which may be a low-luminosity AGN. To clarify the LyC source, we carry out ionization-parameter mapping (IPM) by obtaining narrow-band imaging from the Hubble Space Telescope WFC3 and ACS cameras to construct spatially resolved ratio maps of [OIII]/[OII] emission from the galaxy. IPM traces the ionization structure of the interstellar medium and allows us to identify optically thin regions. To optimize the continuum subtraction, we introduce a new method for determining the best continuum scale factor derived from the mode of the continuum-subtracted, image flux distribution. We find no conclusive evidence of LyC escape from Knots B or C, but instead, we identify a high-ionization region extending over at least 1 kpc from Knot A. Knot A shows evidence of an extremely young age ($\lesssim 1$ Myr), perhaps containing very massive stars ($>100$ M$_\odot$). It is weak in Ly$\alpha$, so if it is confirmed as the LyC source, our results imply that LyC emission may be independent of Ly$\alpha$ emission.

We present the remarkable discovery that the dwarf irregular galaxy NGC 2366 is an excellent analog of the Green Pea (GP) galaxies, which are characterized by extremely high ionization parameters. The similarities are driven predominantly by the giant H II region Markarian 71 (Mrk 71). We compare the system with GPs in terms of morphology, excitation properties, specific star-formation rate, kinematics, absorption of low-ionization species, reddening, and chemical abundance, and find consistencies throughout. Since extreme GPs are associated with both candidate and confirmed Lyman continuum (LyC) emitters, Mrk 71/NGC 2366 is thus also a good candidate for LyC escape. The spatially resolved data for this object show a superbubble blowout generated by mechanical feedback from one of its two super star clusters (SSCs), Knot B, while the extreme ionization properties are driven by the <1 Myr-old, enshrouded SSC Knot A, which has ~ 10 times higher ionizing luminosity. Very massive stars (> 100 Msun) may be present in this remarkable object. Ionization-parameter mapping indicates the blowout region is optically thin in the LyC, and the general properties also suggest LyC escape in the line of sight. Mrk 71/NGC 2366 does differ from GPs in that it is 1 - 2 orders of magnitude less luminous. The presence of this faint GP analog and candidate LyC emitter (LCE) so close to us suggests that LCEs may be numerous and commonplace, and therefore could significantly contribute to the cosmic ionizing budget. Mrk 71/NGC 2366 offers an unprecedentedly detailed look at the viscera of a candidate LCE, and could clarify the mechanisms of LyC escape.

The increasing neutrality of the intergalactic medium at z>6 suppresses Ly-alpha emission, and spectroscopic confirmation of galaxy redshifts requires detecting alternative UV lines. The strong [C III] 1907 + C III] 1909 doublet frequently observed in low-metallicity, actively star-forming galaxies is a promising emission feature. We present CLOUDY photoionization model predictions for C III] equivalent widths (EWs) and line ratios as a function of starburst age, metallicity, and ionization parameter. Our models include a range of C/O abundances, dust content, and gas density. We also examine the effects of varying the nebular geometry and optical depth. Only the stellar models that incorporate binary interaction effects reproduce the highest observed C III] EWs. The spectral energy distributions from the binary stellar population models also generate observable C III] over a longer timescale relative to single-star models. We show that diagnostics using C III] and nebular He II 1640 can separate star-forming regions from shock-ionized gas. We also find that density-bounded systems should exhibit weaker C III] EWs at a given ionization parameter, and C III] EWs could therefore select candidate Lyman continuum-leaking systems. In almost all models, C III] is the next strongest line at < 2700 Angstroms after Ly-alpha, and C III] reaches detectable levels for a wide range of conditions at low metallicity. C III] may therefore serve as an important diagnostic for characterizing galaxies at z>6.

We analyze archival Ly$\alpha$ spectra of 12 "Green Pea" galaxies observed with the Hubble Space Telescope, model their Ly$\alpha$ profiles with radiative transfer models, and explore the dependence of Ly$\alpha$ escape fraction on various properties. Green Pea galaxies are nearby compact starburst galaxies with [OIII]$\lambda$5007 equivalent widths of hundreds of \AA. All 12 Green Pea galaxies in our sample show Ly$\alpha$ lines in emission, with a Ly$\alpha$ equivalent width distribution similar to high redshift Ly$\alpha$ emitters. Combining the optical and UV spectra of Green Pea galaxies, we estimate their Ly$\alpha$ escape fractions and find correlations between Ly$\alpha$ escape fraction and kinematic features of Ly$\alpha$ profiles. The escape fraction of Ly$\alpha$ in these galaxies ranges from 1.4% to 67%. We also find that the Ly$\alpha$ escape fraction depends strongly on metallicity and moderately on dust extinction. We compare their high-quality Ly$\alpha$ profiles with single HI shell radiative transfer models and find that the Ly$\alpha$ escape fraction anti-correlates with the derived HI column densities. Single shell models fit most Ly$\alpha$ profiles well, but not the ones with highest escape fractions of Ly$\alpha$. Our results suggest that low HI column density and low metallicity are essential for Ly$\alpha$ escape, and make a galaxy a Ly$\alpha$ emitter.

HI in galaxies traces the fuel for future star formation and reveals the effects of feedback on neutral gas. Using a statistically uniform, HI-selected sample of 565 galaxies from the ALFALFA H-alpha survey, we explore HI properties as a function of star formation activity. ALFALFA H-alpha provides R-band and H-alpha imaging for a volume-limited subset of the 21-cm ALFALFA survey. We identify eight starbursts based on H-alpha equivalent width and six with enhanced star formation relative to the main sequence. Both starbursts and non-starbursts have similar HI to stellar mass ratios (MHI/M*), which suggests that feedback is not depleting the starbursts' HI. Consequently, the starbursts do have shorter HI depletion times (t_dep), implying more efficient HI-to-H2 conversion. While major mergers likely drive this enhanced efficiency in some starbursts, the lowest mass starbursts may experience periodic bursts, consistent with enhanced scatter in t_dep at low M*. Two starbursts appear to be pre-coalescence mergers; their elevated MHI/M* suggest that HI-to-H2 conversion is still ongoing at this stage. By comparing with the GASS sample, we find that t_dep anti-correlates with stellar surface density for disks, while spheroids show no such trend. Among early-type galaxies, t_dep does not correlate with bulge-to-disk ratio; instead, the gas distribution may determine the star formation efficiency. Finally, the weak connection between galaxies' specific star formation rates and MHI/M* contrasts with the well-known correlation between MHI/M* and color. We show that dust extinction can explain the HI-color trend, which may arise from the relationship between M*, MHI, and metallicity.

We suggest that low optical depth in the Lyman continuum (LyC) may relate the Ly-alpha emission, C II and Si II absorption, and C II* and Si II* emission seen in high-redshift galaxies. We base this analysis on Hubble Space Telescope COS spectra of four Green Pea (GP) galaxies, which may be analogs of z>2 Ly-alpha emitters (LAEs). In the two GPs with the strongest Ly-alpha emission, the Ly-alpha line profiles show reduced signs of resonant scattering. Instead, the Ly-alpha profiles resemble the H-alpha line profiles of evolved star ejecta, suggesting that the Ly-alpha emission originates from a low column density and similar outflow geometry. The weak C II absorption and presence of non-resonant C II* emission in these GPs support this interpretation and imply a low LyC optical depth along the line of sight. In two additional GPs, weak Ly-alpha emission and strong C II absorption suggest a higher optical depth. These two GPs differ in their Ly-alpha profile shapes and C II* emission strengths, however, indicating different inclinations of the outflows to our line of sight. With these four GPs as examples, we explain the observed trends linking Ly-alpha, C II, and C II* in stacked LAE spectra, in the context of optical depth and geometric effects. Specifically, in some galaxies with strong Ly-alpha emission, a low LyC optical depth may allow Ly-alpha to escape with reduced scattering. Furthermore, C II absorption, C II* emission, and Ly-alpha profile shape can reveal the optical depth, constrain the orientation of neutral outflows in LAEs, and identify candidate LyC emitters.

Our understanding of radiative feedback and star formation in galaxies at high redshift is hindered by the rarity of similar systems at low redshift. However, the recently identified Green Pea (GP) galaxies are similar to high-redshift galaxies in their morphologies and star formation rates and are vital tools for probing the generation and transmission of ionizing photons. The GPs contain massive star clusters that emit copious amounts of high-energy radiation, as indicated by intense [OIII] 5007 emission and HeII 4686 emission. We focus on six GP galaxies with high ratios of [O III] 5007,4959/[O II] 3727 ~10 or more. Such high ratios indicate gas with a high ionization parameter or a low optical depth. The GP line ratios and ages point to chemically homogeneous massive stars, Wolf-Rayet stars, or shock ionization as the most likely sources of the He II emission. Models including shock ionization suggest that the GPs may have low optical depths, consistent with a scenario in which ionizing photons escape along passageways created by recent supernovae. The GPs and similar galaxies can shed new light on cosmic reionization by revealing how ionizing photons propagate from massive star clusters to the intergalactic medium.

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 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.

Although Lyman continuum (LyC) radiation from star-forming galaxies likely drove the reionization of the Universe, observations of star-forming galaxies at low redshift generally indicate low LyC escape fractions. However, the extreme [O III]/[O II] ratios of the z=0.1-0.3 Green Pea galaxies may be due to high escape fractions of ionizing radiation. To analyze the LyC optical depths and ionizing sources of these rare, compact starbursts, we compare nebular photoionization and stellar population models with observed emission lines in the Peas' SDSS spectra. We focus on the six most extreme Green Peas, the galaxies with the highest [O III]/[O II] ratios and the best candidates for escaping ionizing radiation. The Balmer line equivalent widths and He I {\lambda}3819 emission in the extreme Peas support young ages of 3-5 Myr, and He II {\lambda}4686 emission in five extreme Peas signals the presence of hard ionizing sources. Ionization by active galactic nuclei or high-mass X-ray binaries is inconsistent with the Peas' line ratios and ages. Although stacked spectra reveal no Wolf-Rayet (WR) features, we tentatively detect WR features in the SDSS spectra of three extreme Peas. Based on the Peas' ages and line ratios, we find that WR stars, chemically homogeneous O stars, or shocks could produce the observed He II emission. If hot stars are responsible, the Peas' optical depths are ambiguous. However, accounting for emission from shocks lowers the inferred optical depth and suggests that the Peas may be optically thin. The Peas' ages likely optimize the escape of Lyman-continuum radiation; they are old enough for supernovae and stellar winds to reshape the interstellar medium, but young enough to possess large numbers of UV-luminous O or WR stars.

The hot, X-ray-emitting gas in superbubbles imparts energy and enriched material to the interstellar medium (ISM) and generates the hot ionized medium, the ISM's high-temperature component. The evolution of superbubble energy budgets is not well understood, however, and the processes responsible for enhanced X-ray emission in superbubbles remain a matter of debate. We present Chandra ACIS-S observations of two X-ray-bright superbubbles in the Large Magellanic Cloud (LMC), DEM L50 (N186) and DEM L152 (N44), with an emphasis on disentangling the true superbubble X-ray emission from non-related diffuse emission and determining the spatial origin and spectral variation of the X-ray emission. An examination of the superbubble energy budgets shows that on the order of 50% of the X-ray emission comes from regions associated with supernova remnant (SNR) impacts. We find some evidence of mass-loading due to swept-up clouds and metallicity enrichment, but neither mechanism provides a significant contribution to the X-ray luminosities. We also find that one of the superbubbles, DEM L50, is likely not in collisional ionization equilibrium. We compare our observations to the predictions of the standard Weaver et al. model and to 1-D hydrodynamic simulations including cavity supernova impacts on the shell walls. Our observations show that mass-loading due to thermal evaporation from the shell walls and SNR impacts are the dominant source of enhanced X-ray luminosities in superbubbles. These two processes should affect most superbubbles, and their contribution to the X-ray luminosity must be considered when determining the energy available for transport to the ISM.

We have obtained spectrophotometric observations of 41 anticenter planetary nebulae (PNe) located in the disk of the Milky Way. Electron temperatures and densities, as well as chemical abundances for He, N, O, Ne, S, Cl, and Ar were determined. Incorporating these results into our existing database of PN abundances yielded a sample of 124 well-observed objects with homogeneously-determined abundances extending from 0.9-21 kpc in galactocentric distance. We performed a detailed regression analysis which accounted for uncertainties in both oxygen abundances and radial distances in order to establish the metallicity gradient across the disk to be: 12+log(O/H)=(9.09+/-.05) - (0.058+/-.006) x Rg, with Rg in kpc. While we see some evidence that the gradient steepens at large galactocentric distances, more objects toward the anticenter need to be observed in order to confidently establish the true form of the metallicity gradient. We find no compelling evidence that the gradient differs between Peimbert Types I and II, nor is oxygen abundance related to the vertical distance from the galactic plane. Our gradient agrees well with analogous results for H II regions but is steeper than the one recently published by Stanghellini & Haywood (2010) over a similar range in galactocentric distance. A second analysis using PN distances from a different source implied a flatter gradient, and we suggest that we have reached a confusion limit which can only be resolved with greatly improved distance measurements and an understanding of the natural scatter in oxygen abundances. Finally, a consideration of recently published chemical evolution models of the Galactic disk suggests that reconciling the current range in published oxygen gradients is necessary for adequately constraining parameters such as the surface density threshold for star formation and the characteristic timescale for disk formation.