We explore the production and escape of ionizing photons in young galaxies by investigating the ultraviolet and optical emission-line properties of models of ionization-bounded and density-bounded HII regions, active-galactic-nucleus (AGN) narrow-line regions and radiative shocks computed all using the same physically-consistent description of element abundances and depletion on to dust grains down to very low metallicities. We compare these models with a reference sample of metal-poor star-forming galaxies and Lyman-continuum (LyC) leakers at various redshifts, which allows the simultaneous exploration of more spectral diagnostics than typically available at once for individual subsamples. We confirm that current single- and binary-star population synthesis models do not produce hard-enough radiation to account for the high-ionization emission of the most metal-poor galaxies. Introducing either an AGN or radiative-shock component brings models into agreement with observations. A published model including X-ray binaries is an attractive alternative to reproduce the observed rise in HeII4686/Hbeta ratio with decreasing oxygen abundance in metal-poor star-forming galaxies, but not the high observed HeII4686/Hbeta ratios of galaxies with large EW(Hbeta). A source of harder ionizing radiation appears to be required in these extreme objects, such as an AGN or radiative-shock component, perhaps linked to an initial-mass-function bias toward massive stars at low metallicity. This would also account for the surprisingly high [OI]/[OIII] ratios of confirmed LyC leakers relative to ionization-bounded models. We find no simple by-eye diagnostic of the nature of ionizing sources and the escape of LyC photon, which require proper simultaneous fits of several lines to be discriminated against.

(Abridged) Long gamma-ray bursts (LGRB) have been suggested as promising tracers of star formation owing to their association with the core-collapse of massive stars. The goal of this work is to characterise the population of host galaxies of LGRBs at 1 < z < 2, investigate the conditions in which LGRBs form at these redshifts and assess their use as tracers of star formation. We perform a spectro-photometric analysis to determine the stellar mass, star formation rate, specific star formation rate and metallicity of the complete, unbiased host galaxy sample of the Swift/BAT6 LGRB sample at 1 < z < 2. We compare the distribution of these properties to the ones of typical star-forming galaxies from the MOSDEF and COSMOS2015 Ultra Deep surveys, within the same redshift range. We find that, similarly to z < 1, LGRBs do not directly trace star formation at 1 < z < 2, and they tend to avoid high-mass, high-metallicity host galaxies. We also find evidence for an enhanced fraction of starbursts among the LGRB host sample with respect to the star-forming population of galaxies. Nonetheless we demonstrate that the driving factor ruling the LGRB efficiency is metallicity. The LGRB host distributions can be reconciled with the ones expected from galaxy surveys by imposing a metallicity upper limit of 12+logOH ~ 8.55. Metallicity rules the LGRB production efficiency, which is stifled at Z > 0.7 Zsun. Under this hypothesis we can expect LGRBs to trace star formation at z > 3, once the bulk of the star forming galaxy population are characterised by metallicities below this limit. The moderately high metallicity threshold found is in agreement with the conditions necessary to rapidly produce a fast-rotating Wolf-Rayet star a in close binary system, and could be accommodated by single star models under chemically homogeneous mixing with very rapid rotation and weak magnetic coupling.

Measurements of the galaxy UV luminosity function at z>6 suggest that young stars hosted in low-mass star-forming galaxies produced the bulk of hydrogen-ionizing photons necessary to reionize the intergalactic medium (IGM) by redshift z~6. Whether star-forming galaxies dominated cosmic reionization, however, also depends on their stellar populations and interstellar medium properties, which set, among other things, the production rate of H-ionizing photons, $\xi_\text{ion}^\star$, and the fraction of these escaping into the IGM. Given the difficulty of constraining with existing observatories the physical properties of z>6 galaxies, in this work we focus on a sample of ten nearby objects showing UV spectral features comparable to those observed at z>6. We use the new-generation Beagle tool to model the UV-to-optical photometry and UV/optical emission lines of these Local 'analogues' of high-redshift galaxies, finding that our relatively simple, yet fully self-consistent, physical model can successfully reproduce the different observables considered. Our galaxies span a broad range of metallicities and are characterised by high ionization parameters, low dust attenuation, and very young stellar populations. Through our analysis, we derive a novel diagnostic of the production rate of H-ionizing photons per unit UV luminosity, $\xi_\text{ion}^\star$, based on the equivalent width of the bright $[\text{OIII}] \lambda 4959,5007$ doublet, which does not require measurements of H-recombination lines. This new diagnostic can be used to estimate $\xi_\text{ion}^\star$ from future direct measurements of the $[\text{OIII}] \lambda 4959,5007$ line using JWST/NIRSpec (out to z~9.5), and by exploiting the contamination by $\text{H}\beta + [\text{OIII}] \lambda 4959,5007$ of photometric observations of distant galaxies, for instance from existing Spitzer/IRAC data and from future ones with JWST/NIRCam.

We combine state-of-the-art models for the production of stellar radiation and its transfer through the interstellar medium (ISM) to investigate ultraviolet-line diagnostics of stars, the ionized and the neutral ISM in star-forming galaxies. We start by assessing the reliability of our stellar population synthesis modelling by fitting absorption-line indices in the ISM-free ultraviolet spectra of 10 Large-Magellanic-Cloud clusters. In doing so, we find that neglecting stochastic sampling of the stellar initial mass function in these young ($\sim10$-100 Myr), low-mass clusters affects negligibly ultraviolet-based age and metallicity estimates but can lead to significant overestimates of stellar mass. Then, we proceed and develop a simple approach, based on an idealized description of the main features of the ISM, to compute in a physically consistent way the combined influence of nebular emission and interstellar absorption on ultraviolet spectra of star-forming galaxies. Our model accounts for the transfer of radiation through the ionized interiors and outer neutral envelopes of short-lived stellar birth clouds, as well as for radiative transfer through a diffuse intercloud medium. We use this approach to explore the entangled signatures of stars, the ionized and the neutral ISM in ultraviolet spectra of star-forming galaxies. We find that, aside from a few notable exceptions, most standard ultraviolet indices defined in the spectra of ISM-free stellar populations are prone to significant contamination by the ISM, which increases with metallicity. We also identify several nebular-emission and interstellar-absorption features, which stand out as particularly clean tracers of the different phases of the ISM.

Nearby dwarf galaxies provide a unique laboratory in which to test stellar population models below $Z_\odot/2$. Such tests are particularly important for interpreting the surprising high-ionization UV line emission detected at $z>6$ in recent years. We present HST/COS ultraviolet spectra of ten nearby metal-poor star-forming galaxies selected to show He II emission in SDSS optical spectra. The targets span nearly a dex in gas-phase oxygen abundance ($7.8<12+\log\mathrm{O/H}<8.5$) and present uniformly large specific star formation rates (sSFR $\sim 10^2$ $\mathrm{Gyr}^{-1}$). The UV spectra confirm that metal-poor stellar populations can power extreme nebular emission in high-ionization UV lines, reaching C III] equivalent widths comparable to those seen in systems at $z\sim 6-7$. Our data reveal a marked transition in UV spectral properties with decreasing metallicity, with systems below $12+\log\mathrm{O/H}\lesssim 8.0$ ($Z/Z_\odot \lesssim 1/5$) presenting minimal stellar wind features and prominent nebular emission in He II and C IV. This is consistent with nearly an order of magnitude increase in ionizing photon production beyond the $\mathrm{He^+}$-ionizing edge relative to H-ionizing flux as metallicity decreases below a fifth solar, well in excess of standard stellar population synthesis predictions. Our results suggest that often neglected sources of energetic radiation such as stripped binary products and very massive O-stars produce a sharper change in the ionizing spectrum with decreasing metallicity than expected. Consequently, nebular emission in C IV and He II powered by these stars may provide useful metallicity constraints in the reionization era.

In this paper, we present a comprehensive analysis of star-forming galaxies (SFGs) at intermediate redshifts (z~1). We combine the ultra-deep optical spectro-photometric data from the Survey for High-z Absorption Red and Dead Sources (SHARDS) with deep UV-to-FIR observations in the GOODS-N field. Exploiting two of the 25 SHARDS medium-band filters, F687W17 and F823W17, we select [OII] emission line galaxies at z~0.84 and z~1.23 and characterize their physical properties. Their rest-frame equivalent widths (EW$_{\mathrm{rf}}$([OII])), line fluxes, luminosities, star formation rates (SFRs) and dust attenuation properties are investigated. The evolution of the EW$_{\mathrm{rf}}$([OII]) closely follows the SFR density evolution of the universe, with a trend of EW$_{\mathrm{rf}}$([OII])$\propto$(1+z)$^3$ up to redshift z~1, followed by a possible flattening. The SF properties of the galaxies selected on the basis of their [OII] emission are compared with complementary samples of SFGs selected by their MIR and FIR emission, and also with a general mass-selected sample of galaxies at the same redshifts. We demonstrate observationally that the UVJ diagram (or, similarly, a cut in the specific SFR) is only partially able to distinguish the quiescent galaxies from the SFGs. The SFR-M$_*$ relation is investigated for the different samples, yelding a logarithmic slope ~1, in good agreement with previous results. The dust attenuations derived from different SFR indicators (UV(1600), UV(2800), [OII], IR) are compared and show clear trends with respect to both the stellar mass and total SFR, with more massive and highly star-forming galaxies being affected by stronger dust attenuation.