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)

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.

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.

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.

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.

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.

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<z \lesssim 5$, properties that make them useful for a wide range of chemical evolution studies.

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.

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.

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.

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.

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.

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.

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.

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant LABs.

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.