Search

We present the localization and host galaxy of FRB 20190208A, a repeating source of fast radio bursts (FRBs) discovered using CHIME/FRB. As part of the PRECISE repeater localization program on the EVN, we monitored FRB 20190208A for 65.6 hours at $\sim1.4$ GHz and detected a single burst, which led to its VLBI localization with 260 mas uncertainty (2$\sigma$). Follow-up optical observations with the MMT Observatory ($i\gtrsim 25.7$ mag (AB)) found no visible host at the FRB position. Subsequent deeper observations with the GTC, however, revealed an extremely faint galaxy ($r=27.32 \pm0.16$ mag), very likely ($99.95 \%$) associated with FRB 20190208A. Given the dispersion measure of the FRB ($\sim580$ pc cm$^{-3}$), even the most conservative redshift estimate ($z_{\mathrm{max}}\sim0.83$) implies that this is the lowest-luminosity FRB host to date ($\lesssim10^8L_{\odot}$), even less luminous than the dwarf host of FRB 20121102A. We investigate how localization precision and the depth of optical imaging affect host association, and discuss the implications of such a low-luminosity dwarf galaxy. Unlike the other repeaters with low-luminosity hosts, FRB 20190208A has a modest Faraday rotation measure of a few tens of rad m$^{-2}$, and EVN plus VLA observations reveal no associated compact persistent radio source. We also monitored FRB 20190208A for 40.4 hours over 2 years as part of the \'ECLAT repeating FRB monitoring campaign on the Nan\c{c}ay Radio Telescope, and detected one burst. Our results demonstrate that, in some cases, the robust association of an FRB with a host galaxy will require both high localization precision, as well as deep optical follow-up., Comment: submitted to ApJL

Fast radio bursts (FRBs) are unique probes of extragalactic ionized baryonic structure as each signal, through its burst properties, holds information about the ionized matter it encounters along its sightline. FRB 20200723B is a burst with a scattering timescale of $\tau_\mathrm{400\,MHz} >$1 second at 400 MHz and a dispersion measure of DM $\sim$ 244 pc cm$^{-3}$. Observed across the entire CHIME/FRB frequency band, it is the single-component burst with the largest scattering timescale yet observed by CHIME/FRB. The combination of its high scattering timescale and relatively low dispersion measure present an uncommon opportunity to use FRB 20200723B to explore the properties of the cosmic web it traversed. With an $\sim$arcminute-scale localization region, we find the most likely host galaxy is NGC 4602 (with PATH probability $P(O|x)=0.985$), which resides $\sim$30 Mpc away within a sheet filamentary structure on the outskirts of the Virgo Cluster. We place an upper limit on the average free electron density of this filamentary structure of $\langle n_e \rangle < 4.6^{+9.6}_{-2.0} \times 10^{-5}$ cm$^{-3}$, broadly consistent with expectations from cosmological simulations. We investigate whether the source of scattering lies within the same galaxy as the FRB, or at a farther distance from an intervening structure along the line of sight. Comparing with Milky Way pulsar observations, we suggest the scattering may originate from within the host galaxy of FRB 20200723B., Comment: 20 pages, 6 figures, submitted. Comments welcome!

A common issue for machine learning models applied to time-series forecasting is the temporal evolution of the data distributions (i.e., concept drift). Because most of the training data does not reflect such changes, the models present poor performance on the new out-of-distribution scenarios and, therefore, the impact of such events cannot be reliably anticipated ahead of time. We present and publicly release CounterfacTS, a tool to probe the robustness of deep learning models in time-series forecasting tasks via counterfactuals. CounterfacTS has a user-friendly interface that allows the user to visualize, compare and quantify time series data and their forecasts, for a number of datasets and deep learning models. Furthermore, the user can apply various transformations to the time series and explore the resulting changes in the forecasts in an interpretable manner. Through example cases, we illustrate how CounterfacTS can be used to i) identify the main features characterizing and differentiating sets of time series, ii) assess how the model performance depends on these characateristics, and iii) guide transformations of the original time series to create counterfactuals with desired properties for training and increasing the forecasting performance in new regions of the data distribution. We discuss the importance of visualizing and considering the location of the data in a projected feature space to transform time-series and create effective counterfactuals for training the models. Overall, CounterfacTS aids at creating counterfactuals to efficiently explore the impact of hypothetical scenarios not covered by the original data in time-series forecasting tasks., Comment: Submitted. Code publicly available

We present LIMFAST, a semi-numerical code for simulating high-redshift galaxy formation and cosmic reionization as revealed by multi-tracer line intensity mapping (LIM) signals. LIMFAST builds upon and extends the 21cmFAST code widely used for 21 cm cosmology by implementing state-of-the-art models of galaxy formation and evolution. The metagalactic radiation background, including the production of various star-formation lines, together with the 21 cm line signal tracing the neutral intergalactic medium (IGM), are self-consistently described by photoionization modeling and stellar population synthesis coupled to the galaxy formation model. We introduce basic structure and functionalities of the code, and demonstrate its validity and capabilities by showing broad agreements between the predicted and observed evolution of cosmic star formation, IGM neutral fraction, and metal enrichment. We also present the LIM signals of 21 cm, Ly$\alpha$, H$\alpha$, H$\beta$, [OII], and [OIII] lines simulated by LIMFAST, and compare them with results from the literature. We elaborate on how several major aspects of our modeling framework, including models of star formation, chemical enrichment, and photoionization, may impact different LIM observables and thus become testable once applied to observational data. LIMFAST aims at being an efficient and resourceful tool for intensity mapping studies in general, exploring a wide range of scenarios of galaxy evolution and reionization and frequencies over which useful cosmological signals can be measured., Comment: 24 pages, 12 figures; to appear in ApJ after being significantly updated from v1

The epoch of reionization (EoR) offers a unique window into the dawn of galaxy formation, through which high-redshift galaxies can be studied by observations of both themselves and their impact on the intergalactic medium. Line intensity mapping (LIM) promises to explore cosmic reionization and its driving sources by measuring intensity fluctuations of emission lines tracing the cosmic gas in varying phases. Using LIMFAST, a novel semi-numerical tool designed to self-consistently simulate LIM signals of multiple EoR probes, we investigate how building blocks of galaxy formation and evolution theory, such as feedback-regulated star formation and chemical enrichment, might be studied with multi-tracer LIM during the EoR. On galaxy scales, we show that the star formation law and the feedback associated with star formation can be indicated by both the shape and redshift evolution of LIM power spectra. For a baseline model of metal production that traces star formation, we find that lines highly sensitive to metallicity are generally better probes of galaxy formation models. On larger scales, we demonstrate that inferring ionized bubble sizes from cross-correlations between tracers of ionized and neutral gas requires a detailed understanding of the astrophysics that shape the line luminosity--halo mass relation. Despite various modeling and observational challenges, wide-area, multi-tracer LIM surveys will provide important high-redshift tests for the fundamentals of galaxy formation theory, especially the interplay between star formation and feedback by accessing statistically the entire low-mass population of galaxies as ideal laboratories, complementary to upcoming surveys of individual sources by new-generation telescopes., Comment: 24 pages, 15 figures; to appear in ApJ

We demonstrate the potential of line-intensity mapping to place constraints on the initial mass function (IMF) of Population III stars via measurements of the mean He II 1640A/H$\alpha$ line-intensity ratio. We extend the 21cmFAST code with modern high-redshift galaxy-formation and photoionization models, and estimate the line emission from Population II and Population III galaxies at redshifts $5 \le z \le 20$. In our models, mean ratio values of He II/H$\alpha \gtrsim 0.1$ indicate top-heavy Population III IMFs with stars of several hundred solar masses, reached at $z \gtrsim 10$ when Population III stars dominate star formation. A next-generation space mission with capabilities moderately superior to those of CDIM will be able to probe this scenario by measuring the He II and H$\alpha$ fluctuation power spectrum signals and their cross-correlation at high significance up to $z\sim 20$. Moreover, regardless of the IMF, a ratio value of He II/H$\alpha \lesssim 0.01$ indicates low Population III star formation and, therefore, it signals the end of the period dominated by this stellar population. However, a detection of the corresponding He II power spectrum may be only possible for top-heavy Population III IMFs or through cross-correlation with the stronger H$\alpha$ signal. Finally, ratio values of $0.01 \lesssim$ He II/H$\alpha$ $\lesssim 0.1$ are complex to interpret because they can be driven by several competing effects. We discuss how various measurements at different redshifts and the combination of the line-intensity ratio with other probes can assist in constraining the Population III IMF in this case., Comment: Updated to match the accepted ApJ version. Main results in Fig.3 and Fig.6

We forecast the reionization history constraints, inferred from Lyman-alpha damping wing absorption features, for a future sample of $\sim 20$ $z \geq 6$ gamma-ray burst (GRB) afterglows. We describe each afterglow spectrum by a three-parameter model. First, L characterizes the size of the ionized region (the "bubble size") around a GRB host halo. Second, $\langle{x_{\rm HI}\rangle}$ is the volume-averaged neutral fraction outside of the ionized bubble around the GRB, which is approximated as spatially uniform. Finally, $N_{\mathrm{HI}}$ denotes the column-density of a local damped Lyman-alpha absorber (DLA) associated with the GRB host galaxy. The size distribution of ionized regions is extracted from a numerical simulation of reionization, and evolves strongly across the Epoch of Reionization (EoR). The model DLA column densities follow the empirical distribution determined from current GRB afterglow spectra. We use a Fisher matrix formalism to forecast the $\langle{x_{\rm HI}(z)\rangle}$ constraints that can be obtained from follow-up spectroscopy of afterglows with SNR = 20 per R=3,000 resolution element at the continuum. We find that the neutral fraction may be determined to better than 10-15\% (1-$\sigma$) accuracy from this data across multiple independent redshift bins at $z \sim 6-10$, spanning much of the EoR, although the precision degrades somewhat near the end of reionization. A more futuristic survey with $80$ GRB afterglows at $z \geq 6$ can improve the precision here by a factor of $2$ and extend measurements out to $z \sim 14$. We further discuss how these constraints may be combined with estimates of the escape fraction of ionizing photons, derived from the DLA column density distribution towards GRBs extracted at slightly lower redshift. This combination will help in testing whether we have an accurate census of the sources that reionized the universe., Comment: To be submitted to ApJ

We present a formalism that incorporates hydrogen Lyman-alpha (Ly$\alpha$) polarization arising from the scattering of radiation in galaxy halos into the intensity mapping approach. Using the halo model, and emission profiles based on simulations and observations, we calculate auto and cross power spectra at redshifts $3\leq z \leq13$ for the Ly$\alpha$ total intensity, $I$, polarized intensity, $P$, degree of polarization, $\Pi=P/I$, and two new quantities, the astrophysical $E$ and $B$ modes of Ly$\alpha$ polarization. The one-halo terms of the $\Pi$ power spectra show a turnover that signals the average extent of the polarization signal, and thus the extent of the scattering medium. The position of this feature depends on redshift, and on the specific emission profile shape and extent, in our formalism. Therefore, the comparison of various Ly$\alpha$ polarization quantities and redshifts can break degeneracies between competing effects, and it can reveal the true shape of the emission profiles, and thus, the physical properties of the cool gas in halos. Furthermore, measurements of Ly$\alpha$ $E$ and $B$ modes may be used as probes of galaxy evolution, because they are related to the average degree of anisotropy in the emission and in the halo gas distribution across redshifts. The detection of the polarization signal at $z \sim 3-5$ requires improvements in the sensitivity of current ground-based experiments by a factor of $\sim 10$, and of $\sim 100$ for space-based instruments targeting the redshifts $z\sim 9-10$, the exact values depending on the specific redshift and experiment. Interloper contamination in polarization is expected to be small, because the interlopers need to also be polarized. Overall, Ly$\alpha$ polarization boosts the amount of physical information retrievable on galaxies and their surroundings, most of it not achievable with total emission alone. (abridged), Comment: Phys. Rev. D, accepted version. Main plots: figs. 3, 8, 9

The Cosmic Dawn and Reionization epochs remain a fundamental but challenging frontier of astrophysics and cosmology. We advocate a large-scale, multi-tracer approach to develop a comprehensive understanding of the physics that led to the formation and evolution of the first stars and galaxies. We highlight the line intensity mapping technique to trace the multi-phase reionization topology on large scales, and measure reionization history in detail. Besides 21cm, we advocate for Lya tomography mapping during the epoch of Wouthuysen-Field coupling as an additional probe of the cosmic dawn era., Comment: Astro2020 Science White Paper; 5 pages, 1 figure

We report on quasar outflow properties revealed by analyzing more than 60 composite outflow spectra built from $\sim 60\,000$ CIV absorption troughs in the SDSS-III/BOSS DR12QBAL catalog. We assess the dependences of the equivalent widths of many outflow metal absorption features on outflow velocity, trough width and position, and quasar magnitude and redshift. The evolution of the equivalent widths of the OVI and NV lines with outflow velocity correlates with that of the mean absorption-line width, the outflow electron density, and the strength of lines arising from collisionally-excited meta-stable states. None of these correlations is found for the other high- or low-ionization species, and different behaviors with trough width are also suggested. We find no dependence on quasar magnitude or redshift in any case. All the observed trends can be reconciled by considering a multiphase stratified outflow structure, where inner regions are colder, denser and host lower-ionization species. Given the prevalence of radiative acceleration in quasar outflows found by Mas-Ribas & Mauland (2019), we suggest that radiation pressure sweeps up and compresses the outflowing gas outwards, creating waves or filaments where the multiphase stratified structure could take form. This scenario is supported by the suggested correlation between electron density and outflow velocity, and the similar behavior observed for the line and line-locking components of the absorption features. We show that this outflow structure is also consistent with other X-ray, radiative transfer, and polarization results, and discuss the implications of our findings for future observational and numerical quasar outflow studies., Comment: Main results Figs. 3 and 7. ApJ accepted

Observational evidence revealing the main mechanisms that accelerate quasar outflows has proven difficult to obtain, due to the complexity of the absorption features that this gas produces in the spectra of the emission sources. We build 36 composite outflow spectra, covering a large range of outflow and quasar parameters, by stacking broad ($>450\,{\rm km\,s^{-1}}$) absorption line systems in the spectra of SDSS-III/BOSS DR12 quasars. The two lines of the atomic doublet of CIV, with a separation of $\approx 497\,{\rm km\,s^{-1}}$, as well as those of other species appear well resolved in most of our composites. This agrees with broad outflow troughs consisting in the superposition of narrow absorbers. We also report on the ubiquitous detection of the radiative-acceleration signature known as line locking in all our composite outflow spectra, including one spectrum strictly built from broad absorption line (BAL) systems. This is the first line-locking detection in BAL composite spectra. Line locking is driven by the CIV atomic doublet, and is visible on the blue side of most strong absorption transitions. Similar effects from the doublets of OVI, SiIV, or NV, however, seem to not be present. Our results confirm that radiation pressure is a prevalent mechanism for accelerating outflows in quasars., Comment: Accepted for publication at ApJ. Data at https://github.com/lluism/BALs

We measure the cosmological bias factor of DLAs from their cross-correlation with the Ly$\alpha$ forest absorption, as a function of the DLA metal strength, defined from an average of equivalent widths of the strongest detectable low-ionization metal lines. A clear increase of the bias factor with metal strength is detected, as expected from a relation of metallicity and velocity dispersion with host halo mass. The relation is stronger after the metal strength is corrected for the HI column density, to make it more related to metallicity instead of metal column density. After correcting for the effects of measurement errors of the metal strength parameter, we find that the bias factor of DLAs with the weakest metal lines is close to unity, consistent with an origin in dwarf galaxies with host halo masses $\sim10^{10} M_{\odot}$, whereas the most metal rich DLAs have a bias factor as large as $b_{\rm DLA} \sim 3$, indicative of massive galaxies or galaxy groups in host halos with masses $\sim 10^{12} M_{\odot}$. Our result confirms the physical origin of the relation of bias factors measured from cross-correlation studies to the host halos of the absorbers., Comment: 18 pages, 5 figures, 2 tables, accepted to MNRAS

We present the first paper of the series Origin of Metals around Galaxies (OMG) aimed to study the origin of the metals observed in the circumgalactic and intergalactic media. In this work we extract and build the catalogs of metal absorbers that will be used in future analyses, and make our results publicly available to the community. We design a fully automatic algorithm to search for absorption metal-line doublets of the species CIV, NV, SiIV and MgII in high-resolution ($R\gtrsim30\,000$) quasar spectra without human intervention, and apply it to the high-resolution and signal-to-noise ratio spectra of 690 quasars, observed with the UVES and HIRES instruments. We obtain $5\,656$ CIV doublets, $7\,919$ doublets of MgII, $2\,258$ of SiIV, and 239 of NV, constituting the largest high-resolution metal-doublet samples to date, and estimate the dependence of their completeness and purity on various doublet parameters such as equivalent width and redshift, using real and artificial quasar spectra. The catalogs include doublets with rest-frame line equivalent widths down to a few ${\rm m\AA}$, all detected at a significance above 3$\sigma$, and covering the redshifts between $1

The baryonic content and physical properties of the warm and hot ($10^5\lesssim T\lesssim 10^7$ K) phases of the circumgalactic medium (CGM) are poorly constrained, owing to the lack of observables probing the requisite range of temperature, spatial scale, halo mass, and redshift. The radiation from a luminous quasar produces a spatially extended emission halo resulting from Thomson scattering off of free electrons in the CGM, which can be used to measure the electron density profile, and therefore, the amount of warm and hot baryonic matter present. We predict the resulting surface brightness profiles and show that they are easily detectable in a three hour integration with the James Webb Space Telescope (JWST), out to $\sim 100$ physical kpc from the centers of individual hyper-luminous quasars. This electron scattering surface brightness is redshift independent, and the signal-to-noise ratio depends only very weakly on redshift, in principle allowing measurements of the warm and hot CGM into the Epoch of Reionization at $z\sim 6.5$. We consider a litany of potential contaminants, and find that for fainter quasars at $z\lesssim1$, extended stellar halos might be of comparable surface brightness. At $z>2$, JWST mid-IR observations start to probe rest-frame optical/UV wavelengths implying that scattering by dust grains in the CGM becomes significant, although multi-color observations should be able to distinguish these scenarios given that Thomson scattering is achromatic., Comment: 22 pages, 8 figures. Version accepted for publication at ApJ. Busy reader, see Fig. 6 (quasi-constant detectability of Thomson scattering halos up to Reionization)

The Baryon Oscillation Spectroscopic Survey of SDSS-III has provided an unprecedentedly large sample of Damped \lya systems (DLAs), the largest repositories of neutral hydrogen in the Universe. This DLA sample has been used to determine the DLA bias factor from their cross-correlation with the \lya forest absorption in \cite{FontRibera2012,Perez2018}, showing that DLAs are associated with relatively massive halos. However, the low resolution and signal-to-noise of BOSS spectra do not allow precise measurements of the DLA metal lines. We define a metal strength parameter, $S$, based on combining equivalent widths of 17 metal lines, that can be measured with an optimal signal-to-noise ratio for individual DLAs in BOSS spectra, allowing for the classification of the DLA population into subgroups of different $S$. We present the distribution of this DLA metal strength and the dependence of its mean value on $N_{\rm HI}$ and redshift. We search for systematic effects and variations in the catalogue purity by examining the dependence of the $S$ distribution on the spectral signal-to-noise and the estimated error on $S$. A catalogue of DLAs with measured equivalent widths for the selected 17 metal lines and the value of $S$ is made publicly available, which will be used to measure the dependence of the DLA bias factor on the $S$ parameter. The relation of the metal strength on the gas metal abundances and velocity dispersion can be constrained by studying the stacked metal absorption spectra of DLAs as a function of $S$, allowing for future determinations of the dependence of the bias factor on the metallicity and velocity dispersion of DLAs., Comment: 12 pages, 10 figures, link to online catalogue (https://github.com/andreuandreu/Catalogue_Metal_Strength_DLAs_SDSS-BOSS_DR12)

We present a measurement of the DLA mean bias from the cross-correlation of DLA and the Ly$\alpha$ forest, updating earlier results of Font-Ribera et al. 2012 with the final BOSS Data Release and an improved method to address continuum fitting corrections. Our cross-correlation is well fitted by linear theory with the standard $\Lambda CDM$ model, with a DLA bias of $b_{\rm DLA} = 1.99\pm 0.11$; a more conservative analysis, which removes DLA in the Ly$\beta$ forest and uses only the cross-correlation at $r> 10{\rm h^{-1}\,Mpc}$, yields $b_{\rm DLA} = 2.00\pm 0.19$. This assumes the cosmological model from \cite{Planck2015} and the Ly$\alpha$ forest bias factors of Bautista et al. 2017, and includes only statistical errors obtained from bootstrap analysis. The main systematic errors arise from possible impurities and selection effects in the DLA catalogue, and from uncertainties in the determination of the Ly$\alpha$ forest bias factors and a correction for effects of high column density absorbers. We find no dependence of the DLA bias on column density or redshift. The measured bias value corresponds to a host halo mass $\sim 4\cdot10^{11} {\rm M_{\odot}}$ if all DLA were hosted in halos of a similar mass. In a realistic model where host halos over a broad mass range have a DLA cross section $\Sigma(M_h) \propto M_h^{\alpha}$ down to $M_h > M_{\rm min} =10^{8.5} {\rm M_{\odot}}$, we find that $\alpha > 1$ is required to have $b_{\rm DLA}> 1.7$, implying a steeper relation or higher value of $M_{\rm min}$ than is generally predicted in numerical simulations of galaxy formation., Comment: accepted to MNRAS, 23 pages, 8 figures, 4 tables

We present a new method to quantify the value of the escape fraction of ionizing photons, and the existence of ultra-faint galaxies clustered around brighter objects during the epoch of cosmic reionization, using the diffuse Ly$\alpha$, continuum and H$\alpha$ emission observed around galaxies at $z\sim6$. We model the surface brightness profiles of the diffuse halos considering the fluorescent emission powered by ionizing photons escaping from the central galaxies, and the nebular emission from satellite star-forming sources, by extending the formalisms developed in Mas-Ribas & Dijkstra (2016) and Mas-Ribas et al. (2017). The comparison between our predicted profiles and Ly$\alpha$ observations at $z=5.7$ and $z=6.6$ favors a low ionizing escape fraction, $f_{\rm esc}^{\rm ion}\sim5\%$, for galaxies in the range $-19\gtrsim M_{\rm UV} \gtrsim -21.5$. However, uncertainties and possible systematics in the observations do not allow for firm conclusions. We predict H$\alpha$ and rest-frame visible continuum observations with JWST, and show that JWST will be able to detect extended (a few tens of kpc) fluorescent H$\alpha$ emission powered by ionizing photons escaping from a bright, $L\gtrsim 5L^*$, galaxy. Such observations can differentiate fluorescent emission from nebular emission by satellite sources. We discuss how observations and stacking of several objects may provide unique constraints on the escape fraction for faint galaxies and/or the abundance of ultra-faint radiation sources., Comment: 9 pages, 4 figures, re-submitted after referee report to ApJ

Lyman alpha halos are observed ubiquitously around star-forming galaxies at high redshift, but their origin is still a matter of debate. We demonstrate that the emission from faint unresolved satellite sources, $M_{\rm UV} \gtrsim -17$, clustered around the central galaxies may play a major role in generating spatially extended Ly$\alpha$, continuum (${\rm UV + VIS}$) and H$\alpha$ halos. We apply the analytic formalism developed in Mas-Ribas & Dijkstra (2016) to model the halos around Lyman Alpha Emitters (LAEs) at $z=3.1$, for several different satellite clustering prescriptions. In general, our UV and Ly$\alpha$ surface brightness profiles match the observations well at $20\lesssim r \lesssim 40$ physical kpc from the centers of LAEs. We discuss how our profiles depend on various model assumptions and how these can be tested and constrained with future H$\alpha$ observations by the James Webb Space Telescope (JWST). Our analysis shows how spatially extended halos constrain (i) the presence of otherwise undetectable satellite sources, (ii) the integrated, volumetric production rates of Ly$\alpha$ and LyC photons, and (iii) their population-averaged escape fractions. These quantities are all directly relevant for understanding galaxy formation and evolution and, for high enough redshifts, cosmic reionization., Comment: 13 pages, 6 figures, edited to match accepted ApJ version. Results unaffected. New descriptive flow-chart figure (Fig.6)

Direct collapse black holes forming in pristine, atomically-cooling haloes at $z \approx 10-20$ may act as the seeds of supermassive black holes (BH) at high redshifts. In order to create a massive BH seed, the host halo needs to be prevented from forming stars. H$_2$ therefore needs to be irradiated by a large flux of Lyman-Werner (LW) UV photons in order to suppress H$_2$ cooling. A key uncertainty in this scenario is the escape fraction of LW radiation from first galaxies, the dominant source of UV photons at this epoch. To better constrain this escape fraction, we have performed radiation-hydrodynamical simulations of the growth of HII regions and their associated photodissociation regions in the first galaxies using the ZEUS-MP code. We find that the LW escape fraction crucially depends on the propagation of the ionisation front (I-front). For an R-type I-front overrunning the halo, the LW escape fraction is always larger than 95%. If the halo recombines later from the outside--in, due to a softened and weakened spectrum, the LW escape fraction in the rest-frame of the halo (the near-field) drops to zero. A detailed and careful analysis is required to analyse slowly moving, D-type I-fronts, where the escape fraction depends on the microphysics and can be as small as 3% in the near-field and 61% in the far-field or as large as 100% in both the near-field and the far-field., Comment: 14 pages, 10 figures, 7 tables, accepted for publication in MNRAS

We study the mean absorption spectrum of the Damped Lyman alpha population at $z\sim 2.6$ by stacking normalized, rest-frame shifted spectra of $\sim 27\,000$ DLAs from the DR12 of BOSS/SDSS-III. We measure the equivalent widths of 50 individual metal absorption lines in 5 intervals of DLA hydrogen column density, 5 intervals of DLA redshift, and overall mean equivalent widths for an additional 13 absorption features from groups of strongly blended lines. The mean equivalent width of low-ionization lines increases with $N_{\rm HI}$, whereas for high-ionization lines the increase is much weaker. The mean metal line equivalent widths decrease by a factor $\sim 1.1-1.5$ from $z\sim2.1$ to $z \sim 3.5$, with small or no differences between low- and high-ionization species. We develop a theoretical model, inspired by the presence of multiple absorption components observed in high-resolution spectra, to infer mean metal column densities from the equivalent widths of partially saturated metal lines. We apply this model to 14 low-ionization species and to AlIII, SIII, SiIII, CIV, SiIV, NV and OVI. We use an approximate derivation for separating the equivalent width contributions of several lines to blended absorption features, and infer mean equivalent widths and column densities from lines of the additional species NI, ZnII, CII${}^{*}$, FeIII, and SIV. Several of these mean column densities of metal lines in DLAs are obtained for the first time; their values generally agree with measurements of individual DLAs from high-resolution, high signal-to-noise ratio spectra when they are available., Comment: Resubmitted after referee revision. Added evolution of metal-line equivalent widths with redshift (Section 5). Added assessment of result dependencies on sample and methodology. Comparison of relative abundances of DLAs vs Milky Way ISM and halo (Figure 16). Publicly available videos of composite quasar and DLA spectra realizations here: https://github.com/lluism/

We revisit calculations of nebular hydrogen Lya and HeII 1640A line strengths for population III galaxies, undergoing continuous and bursts of star formation. We focus on initial mass functions (IMFs) motivated by recent theoretical studies, which generally span a lower range of stellar masses than earlier works. We also account for case-B departures and the stochastic sampling of the IMF. In agreement with previous works, we find that departures from case-B can enhance the Lya flux by a factor of a few, but we argue that this enhancement is driven mainly by collisional excitation and ionization, and not due to photoionization from the n = 2 state of atomic hydrogen. The increased sensitivity of the Lya flux to the high-energy end of the galaxy spectrum makes it more subject to stochastic sampling of the IMF. The latter introduces a dispersion in the predicted nebular line fluxes around the deterministic value by as much as a factor of ~4. In contrast, the stochastic sampling of the IMF has less impact on the emerging Lyman Werner (LW) photon flux. When case-B departures and stochasticity effects are combined, nebular line emission from population III galaxies can be up to one order of magnitude brighter than predicted by 'standard' calculations that do not include these effects. This enhances the prospects for detection with future facilities such as JWST and large, groundbased telescopes., Comment: 12 pages, 6 figures, SEDs publicly available. Accepted for publication in ApJ

We quantify the contribution of Ly$\alpha$ fluorescence to observed spatially extended Ly$\alpha$ halos around Ly$\alpha$ emitters (LAE) at redshift ${\rm z=3.1}$. The key physical quantities that describe the fluorescent signal include (${\it i}$) the distribution of cold gas in the circum-galactic medium (CGM); we explore simple analytical models and fitting functions to recent hydrodynamical simulations; (${\it ii}$) local variations in the ionizing background due to ionizing sources that cluster around the central galaxy. We account for clustering by boosting the observationally inferred volumetric production rate of ionizing photons, $\epsilon_{\rm LyC}$, by a factor of $1+\xi_{\rm LyC}(r)$, in which $\xi_{\rm LyC}(r)$ quantifies the clustering of ionizing sources around the central galaxy. We compute $\xi_{\rm LyC}(r)$ by assigning an 'effective' bias parameter to the ionizing sources. This novel approach allows us to quantify our ignorance of the population of ionizing sources in a simple parametrized form. We find a maximum enhancement in the local ionizing background in the range $50-200$ at $r \sim 10$ physical kpc. For spatially uncorrelated ionizing sources and fluorescing clouds we find that fluorescence can contribute up to $\sim 50-60\%$ of the observed spatially extended Ly$\alpha$ emission. We briefly discuss how future observations can shed light on the nature of Ly$\alpha$ halos around star forming galaxies., Comment: 13 pages, 6 Figures, ApJ in Press. Improved Rahmati et al. model. Results unaffected