We present IRAM/NOEMA and JVLA observations of the quasar J1342+0928 at z = 7.54 and report detections of copious amounts of dust and [C Pi] emission in the interstellar medium (ISM) of its host galaxy. At this redshift, the age of the universe is 690 Myr, about 10% younger than the redshift of the previous quasar record holder. Yet, the ISM of this new quasar host galaxy is significantly enriched by metals, as evidenced by the detection of the [C 158 mu m cooling line and the underlying far-infrared (FIR) dust continuum emission. To the first order, the FIR properties of this quasar host are similar to those found at a slightly lower redshift (z similar to 6), making this source by far the FIR-brightest galaxy known at z greater than or similar to 7.5. The [C Pi]emission is spatially unresolved, with an upper limit on the diameter of 7 kpc. Together with the measured FWHM of the [C Pi]line, this yields a dynamical mass of the host of <1.5 x 10(11) M-circle dot Using standard assumptions about the dust temperature and emissivity, the NOEMA measurements give a dust mass of (0.6-4.3) x 10(8) M-circle dot The brightness of the [C Pi] luminosity, together with the high dust mass, imply active ongoing star formation in the quasar host. Using [C Pi]-SFR scaling relations, we derive star formation rates of 85-545 M-circle dot yr(-1) in the host, consistent with the values derived from the dust continuum. Indeed, an episode of such past high star formation is needed to explain the presence of similar to 10(8) M-circle dot of dust implied by the observations.

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.

Understanding the time evolution of the star formation in the Universe is one of the main aims of observational astronomy. Since a significant portion of the UV starlight is being absorbed by dust and re-emitted in the IR, we need to understand both of those regimes to properly describe the cosmic star formation history. In UV, the depth and the resolution of the data permits calculations of the star formation rate densities out to very high redshifts (z ∼ 8 − 9). In IR however, the large beam sizes and the relatively shallow data limits these calculations to z ∼ 2. In this thesis, I explore the SMA and PdBI high-resolution follow-up of 30 bright sources originally selected by AzTEC and LABOCA instruments at 1.1 mm and 870 μm respectively in conjunction with the SCUBA-2 Cosmology Legacy Survey (S2CLS) deep COSMOS and wide UDS maps, where 106 and 283 sources were detected, with the signal-to-noise ratio of > 5 and > 3.5 at 850 μm respectively. I find that the (sub)mm-selected galaxies reside and the mean redshifts of ¯z ≃ 2.5±0.05 with the exception of the brightest sources which seem to lie at higher redshifts (¯z ≃ 3.5 ± 0.2), most likely due to the apparent correlation of the (sub)mm flux with redshift, where brighter sources tend to lie at higher redshifts. Stellar masses, M⋆, and star formation rates, SFRs, were found (M⋆ & 1010M⊙ and SFR & 100M⊙ yr−1) and used to calculate the specific SFRs. I determine that the (sub)mm-selected sources mostly lie on the high-mass end of the star formation ‘main-sequence’ which makes them a high-mass extension of normal star forming galaxies. I also find that the specific SFR slightly evolves at redshifts 2−4, suggesting that the efficiency of the star formation seems to be increasing at these redshifts. Using the S2CLS data, the bolometric IR luminosity functions (IR LFs) were found for a range of redshifts z = 1.2 − 4.2 and the contribution of the SMGs to the total star formation rate density (SFRD) was calculated. The IR LFs were found to evolve out to redshift ∼ 2.5. The star formation activity in the Universe was found to peak at z ≃ 2 followed by a slight decline. Assuming the IR to total SFRD correction found in the literature the SFRD found in this work closely follows the best-fitting function of Madau & Dickinson (2014).

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.