This paper provides an update of our previous scaling relations between galaxy-integrated molecular gas masses, stellar masses, and star formation rates (SFRs), in the framework of the star formation main sequence (MS), with the main goal of testing for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and similar to 1 mm dust photometry, in a large sample of 1444 star-forming galaxies between z = 0 and 4. The sample covers the stellar mass range log(M-*/M-circle dot) = 9.0-11.8, and SFRs relative to that on the MS, delta MS = SFR/SFR (MS), from 10(-1.3) to 10(2.2). Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero-point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time t(depl), defined as the ratio of molecular gas mass to SFR, scales as (1 + z)(-0.6) x (delta MS)(-0.44) and is only weakly dependent on stellar mass. The ratio of molecular to stellar mass mu(gas) depends on (1+ z)(2.5) x (delta MS)(0.52) x (M-*)(-0.36), which tracks the evolution of the specific SFR. The redshift dependence of mu(gas) requires a curvature term, as may the mass dependences of t(depl) and mu(gas). We find no or only weak correlations of t(depl) and mu(gas) with optical size R or surface density once one removes the above scalings, but we caution that optical sizes may not be appropriate for the high gas and dust columns at high z.

We investigate the relationship between the black hole accretion rate (BHAR) and star formation rate (SFR) for Milky Way (MW) and Andromeda (M31)-mass progenitors from z = 0.2 to 2.5. We source galaxies from the K-s-band-selected ZFOURGE survey, which includes multiwavelength data spanning 0.3-160 mu m. We use decomposition software to split the observed spectral energy distributions (SEDs) of our galaxies into their active galactic nuclei (AGNs) and star-forming components, which allows us to estimate BHARs and SFRs from the infrared (IR). We perform tests to check the robustness of these estimates, including a comparison with BHARs and SFRs derived from X-ray stacking and far-IR analysis, respectively. We find that, as the progenitors evolve their relative black hole-galaxy growth (i.e. their BHAR/SFR ratio) increases from low to high redshift. The MW-mass progenitors exhibit a log-log slope of 0.64 +/- 0.11, while the M31-mass progenitors are 0.39 +/- 0.08. This result contrasts with previous studies that find an almost flat slope when adopting X-ray-/AGN-selected or mass-limited samples and is likely due to their use of a broad mixture of galaxies with different evolutionary histories. Our use of progenitor-matched samples highlights the potential importance of carefully selecting progenitors when searching for evolutionary relationships between BHAR/SFRs. Additionally, our finding that BHAR/SFR ratios do not track the rate at which progenitors quench casts doubts over the idea that the suppression of star formation is predominantly driven by luminous AGN feedback (i.e. high BHARs).

We present a spectroscopic survey of high-redshift, luminous galaxies over four square degrees on the sky, aiming to build a large and homogeneous sample of Ly alpha emitters (LAEs) at z approximate to 5.7 and 6.5, and Lyman-break galaxies (LBGs) at 5.5 < z < 6.8. The fields that we choose to observe are well studied, such as by the Subaru XMM-Newton Deep Survey and COSMOS. They have deep optical imaging data in a series of broad and narrow bands, allowing for the efficient selection of galaxy candidates. Spectroscopic observations are being carried out using the multi-object spectrograph M2FS on the Magellan Clay telescope. M2FS is efficient enough to identify high-redshift galaxies, owing to its 256 optical fibers deployed over a circular field of view 30' in diameter. We have observed similar to 2.5 square degrees. When the program is completed, we expect to identify more than 400 bright LAEs at z approximate to 5.7 and 6.5, and a substantial number of LBGs at z >= 6. This unique sample will be used to study a variety of galaxy properties and to search for large protoclusters. Furthermore, the statistical properties of these galaxies will be used to probe cosmic reionization. We describe the motivation, program design, target selection, and M2FS observations. We also outline our science goals, and present a sample of the brightest LAEs at z approximate to 5.7 and 6.5. This sample contains 32 LAEs with Ly alpha luminosities higher than 10(43) erg s(-1). A few of them reach >= 3 x 10(43) erg s(-1), comparable to the two most luminous LAEs known at z >= 6, "CR7" and "COLA1." These LAEs provide ideal targets to study extreme galaxies in the distant universe.

Context. The reionization of the Universe is one of the most important topics of present-day astrophysical research. The most plausible candidates for the reionization process are star-forming galaxies, which according to the predictions of the majority of the theoretical and semi-analytical models should dominate the H I ionizing background at z greater than or similar to 3. Aims. We measure the Lyman continuum escape fraction, which is one of the key parameters used to compute the contribution of star-forming galaxies to the UV background. It provides the ratio between the photons produced at lambda <= 912 angstrom rest-frame and those that are able to reach the inter-galactic medium, i.e. that are not absorbed by the neutral hydrogen or by the dust of the galaxy's inter-stellar medium. Methods. We used ultra-deep U-band imaging (U = 30.2 mag at 1 sigma) from Large Binocular Camera at the Large Binocular Telescope (LBC/LBT) in the CANDELS/GOODS-North field and deep imaging in the COSMOS and EGS fields in order to estimate the Lyman continuum escape fraction of 69 star-forming galaxies with secure spectroscopic redshifts at 3.27 <= z <= 3.40 to faint magnitude limits (L = 0.2L*, or equivalently M-1500 similar to -19). The narrow redshift range implies that the LBC U-band filter exclusively samples the lambda <= 912 angstrom rest-frame wavelengths. Results. We measured through stacks a stringent upper limit (<1.7% at 1 sigma) for the relative escape fraction of H I ionizing photons from bright galaxies (L > L*), while for the faint population (L = 0.2L*) the limit to the escape fraction is less than or similar to 10%. We computed the contribution of star-forming galaxies to the observed UV background at z similar to 3 and find that it is not sufficient to keep the Universe ionized at these redshifts unless their escape fraction increases significantly (>= 10%) at low luminosities (M-1500 >= -19). Conclusions. We compare our results on the Lyman continuum escape fraction of high-z galaxies with recent estimates in the literature, and discuss future prospects to shed light on the end of the Dark Ages. In the future, strong gravitational lensing will be fundamental in order to measure the Lyman continuum escape fraction down to faint magnitudes (M-1500 similar to -16) that are inaccessible with the present instrumentation on blank fields. These results will be important in order to quantify the role of faint galaxies to the reionization budget.

A brief Chandra observation of the ultraluminous quasar SDSS J010013.02+280225.8 at redshift 6.326 showed it to be a relatively bright, soft X-ray source with a count rate of about 1 count ks(-1). In this article, we present results for the quasar from a 65-ks XMM-Newton observation, which constrains its spectral shape well. The quasar is clearly detected with a total of similar to 460 net counts in the 0.2-10 keV band. The spectrum is characterized by a simple power-law model with a photon index of Gamma = 2.30(-0.10)(+0.10) and the intrinsic 2-10 keV luminosity is 3.14 x 10(45) erg s(-1). The 1 sigma upper limit to any intrinsic absorption column density is N-H = 6.07 x 10(22) cm(-2). No significant iron emission lines were detected. We derive an X-ray-to- optical flux ratio alpha(ox) of -1.74 +/- 0.01, consistent with the values found in other quasars of comparable ultraviolet luminosity. We did not detect significant flux variations either in the XMM-Newton exposure or between XMM-Newton and Chandra observations, which are separated by similar to 8 months. The X-ray observation enables the bolometric luminosity to be calculated after modelling the spectral energy distribution: the accretion rate is found to be sub-Eddington.

Cosmological simulations suggest a strong correlation between high optical-depth Lya absorbers, which arise from the intergalactic medium, and 3D mass overdensities on scales of 10-30 h(-1) comoving Mpc. By examining the absorption spectra of similar to 80,000 QSO sight lines over a volume of 0.1 Gpc(3) in the Sloan Digital Sky Survey III (SDSS-III), we have identified an extreme overdensity, BOSS1441, which contains a rare group of strong Lya absorbers at z = 2.32 +/- 0.02. This absorber group is associated with six QSOs at the same redshift on a 30 comoving Mpc scale. Using Mayall/MOSAIC narrowband and broadband imaging, we detect Lya emitters (LAEs) down to 0.7 x L-Ly alpha ' * a and reveal a large-scale structure of LAEs in this field. Our follow-up LBT observations have spectroscopically confirmed 20 galaxies in the density peak. We show that BOSS1441 has an LAE overdensity of 10.8 +/- 2.6 on a 15 comoving Mpc scale, which could collapse to a massive cluster with M greater than or similar to 10(15) M-circle dot at z similar to 0. This overdensity is among the most massive large-scale structures at z similar to 2 discovered to date.

We present the first discoveries from a survey of z greater than or similar to 6 quasars using imaging data from the DECam Legacy Survey (DECaLS) in the optical, the UKIRT Deep Infrared Sky Survey (UKIDSS) and a preliminary version of the UKIRT Hemisphere Survey (UHS) in the near-IR, and ALLWISE in the mid-IR. DECaLS will image 9000 deg(2) of sky down to z(AB) similar to 23.0, and UKIDSS and UHS will map the northern sky at 0 < decl. < + 60 degrees, reaching J(VEGA) similar to 19.6 (5-sigma). The combination of these data sets allows us to discover quasars at redshift z greater than or similar to 7 and to conduct a complete census of the faint quasar population at z greater than or similar to 6. In this paper, we report on the selection method of our search, and on the initial discoveries of two new, faint z greater than or similar to 6 quasars and one new z = 6.63 quasar in our pilot spectroscopic observations. The two new z similar to 6 quasars are at z = 6.07 and z = 6.17 with absolute magnitudes at rest-frame wavelength 1450 angstrom being M-1450 = -25.83 and M-1450 = -25.76, respectively. These discoveries suggest that we can find quasars close to or fainter than the break magnitude of the Quasar Luminosity Function (QLF) at z greater than or similar to 6. The new z = 6.63 quasar has an absolute magnitude of M-1450 = -25.95. This demonstrates the potential of using the combined DECaLS and UKIDSS/UHS data sets to find z greater than or similar to 7 quasars. Extrapolating from previous QLF measurements, we predict that these combined data sets will yield similar to 200 z similar to 6 quasars to z(AB) < 21.5, similar to 1000 z similar to 6 quasars to z(AB) < 23, and similar to 30 quasars at z > 6.5 to J(VEGA) < 19.5.

Imaging with the James Webb Space Telescope (JWST) will allow observations of the bulk of distant galaxies at the epoch of reionization. The recovery of their properties, such as age, color excess E(B-V), specific star formation rate (sSFR), and stellar mass, will mostly rely on spectral energy distribution fitting, based on the data provided by JWST's two imager cameras, namely the Near Infrared Camera (NIRCam) and the Mid Infrared Imager (MIRI). In this work we analyze the effect of choosing different combinations of NIRCam and MIRI broadband filters, from 0.6 to 7.7 mu m, on the recovery of these galaxy properties. We performed our tests on a sample of 1542 simulated galaxies, with known input properties, at z = 7-10. We found that, with only eight NIRCam broadbands, we can recover the galaxy age within 0.1 Gyr and the color excess within 0.06 mag for 70% of the galaxies. Additionally, the stellar masses and sSFR are recovered within 0.2 and 0.3 dex, respectively, at z = 7-9. Instead, at z = 10, no NIRCam band traces purely the lambda > 4000 angstrom regime and the percentage of outliers in stellar mass (sSFR) increases by > 20% (> 90%), in comparison to z = 9. The MIRI F560W and F770W bands are crucial to improve the stellar mass and the sSFR estimation at z = 10. When nebular emission lines are present, deriving correct galaxy properties is challenging at any redshift and with any band combination. In particular, the stellar mass is systematically overestimated in up to 0.3 dex on average with NIRCam data alone and including MIRI observations only marginally improves the estimation.

We present a measurement of the abundance of carbon monoxide in the early universe, utilizing the final results from the CO Power Spectrum Survey (COPSS). Between 2013 and 2015, we performed observations with the Sunyaev-Zel'dovich Array to measure aggregate CO emission from z similar to 3 galaxies with the intensity mapping technique. Data were collected on 19 fields, covering an area of 0.7 square degrees, over the frequency range 27-35 GHz. With these data, along with data analyzed in COPSS I, we are able to observe the CO(1-0) transition within the redshift range z = 2.3-3.3 for spatial frequencies between k = 0.5-10 h Mpc(-1), spanning a comoving volume of 4.9 x 10(6) h(-3) Mpc(3). We present estimates of contributions from continuum sources and ground illumination within our measurement. We constrain the amplitude of the CO power spectrum to P-CO = 3.0(-1.3)(+1.3) x 10(3) mu K-2(h(-1) Mpc)(3), or Delta(2)(CO)(k=1 h Mpc(-1)) = 1.5(-0.7)(+0.7) x 10(3) mu K-2, at 68% confidence, and PCO > 0 at 98.9% confidence. These results are a factor of 10 improvement in sensitivity compared to those of COPSS I. With this measurement, we constrain on the CO(1-0) galaxy luminosity function at z similar to 3. Assuming that CO emission is proportional to halo mass and using theoretical estimates of the scatter in this relationship, we constrain the ratio of CO(1-0) luminosity to halo mass to A(CO) = 6.3(-2.1)(+1.4) x 10 (7) L circle dot M circle dot-1. Assuming a Milky Way-like linear relationship between CO luminosity and molecular gas mass, we estimate a mass fraction of molecular gas of f(H2) = 5.5(-2.2)(+3.4) x 10(-2) for halos with masses of similar to 10(12)M(circle dot). Using theoretical estimates for the scaling of molecular gas mass fraction and halo mass, we estimate the cosmic molecular gas density to be rho(z similar to 3) (H-2) = 1.1(-0.4)(+0.7) x 10(8) M(circle dot)Mpc(-3).

We use cosmological hydrodynamic simulations to study the impact of out-flows and radiative feedback on high-redshift galaxies. For outflows, we consider simulations that assume (i) no winds, (ii) a .constant-wind. model in which the mass-loading factor and outflow speed are constant, and (iii) "momentum driven" winds in which both parameters vary smoothly with mass. In order to treat radiative feedback, we develop a moment-based radiative transfer technique that operates in both post-processing and coupled radiative hydrodynamic modes. We first ask how outflows impact the broadband spectral energy distributions (SEDs) of six observed reionization-epoch galaxies. Simulations reproduce five regardless of the outflow prescription, while the sixth suggests an unusually bursty star formation history. We conclude that (i) simulations broadly account for available constraints on reionization-epoch galaxies, (ii) individual SEDs do not constrain outflows, and (iii) SED comparisons efficiently isolate objects that challenge simulations. We next study how outflows impact the galaxy mass metallicity relation (MZR). Momentum-driven outflows uniquely reproduce observations at z = 2. In this scenario, galaxies obey two equilibria: (i) The rate at which a galaxy processes gas into stars and outflows tracks its inflow rate; and (ii) The gas enrichment rate owing to star formation balances the dilution rate owing to inflows. Combining these conditions indicates that the MZR is dominated by the (instantaneous) variation of outflows with mass, with more-massive galaxies driving less gas into outflows per unit stellar mass formed. Turning to radiative feedback, we use post-processing simulations to study the topology of reionization. Reionization begins in overdensities and then .leaks. directly into voids, with filaments reionizing last owing to their high density and low emissivity. This result conflicts with previous findings that voids ionize last. We argue that it owes to the uniqely-biased emissivity field produced by our star formation prescriptions, which have previously been shown to reproduce numerous post-reionization constraints. Finally, preliminary results from coupled radiative hydrodynamic simulations indicate that reionization suppresses the star formation rate density by at most 10.20% by z = 5. This is much less than previous estimates, which we attribute to our unique reionization topology although confirmation will have to await more detailed modeling.