In a multi-wavelength survey of $13$ quasars at $5.8\lesssim z\lesssim6.5$, that were pre-selected to be potentially young, we find five objects with extremely small proximity zone sizes that may imply UV-luminous quasar lifetimes of $\lesssim 100,000$ years. Proximity zones are regions of enhanced transmitted flux in the vicinity of the quasars that are sensitive to the quasars' lifetimes because the intergalactic gas has a finite response time to their radiation. We combine sub-mm observations from the Atacama Large Millimetre Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA), as well as deep optical and near-infrared spectra from medium-resolution spectrograph on the Very Large Telescope (VLT) and on the Keck telescopes, in order to identify and characterize these new young quasars, which provide valuable clues about the accretion behavior of supermassive black holes (SMBHs) in the early universe, and pose challenges on current black hole formation models to explain the rapid formation of billion solar mass black holes. We measure the quasars' systemic redshifts, black hole masses, Eddington ratios, emission line luminosities, and star formation rates of their host galaxies. Combined with previous results we estimate the fraction of young objects within the high-redshift quasar population at large to be $5\%\lesssim f_{\rm young}\lesssim 10\%$. One of the young objects, PSO J158-14, shows a very bright dust continuum flux ($F_{\rm cont}=3.46\pm 0.02\,\rm mJy$), indicating a highly star-bursting host galaxy with a star formation rate of approximately $1420\,M_{\odot}\,\rm yr^{-1}$.

We present deep (9 hours) Gemini-N/GNIRS near-infrared spectroscopic observations of ULAS J1342+0928, a luminous quasar at z=7.54. Various broad emission lines were detected, as well as the underlying continuum and iron forests over the rest-frame wavelength 970-2930A. There is a clear trend that higher-ionization emission lines show larger blueshifts with CIV1549 exhibiting 5510^{+240}_{-110} km s-1 blueshift with respect to the systematic redshift from the far-infrared [CII] 158um emission line. Those high ionization lines have wide profiles with FWHM more than 10000 km s-1. A modest blueshift of 340^{+110}_{-80} km s-1 is also seen in MgII, the lowest ionization line identified in the spectrum. The updated MgII-based black hole mass of M_BH=9.1_{-1.3}^{+1.4} x 10^8 M_sun and the Eddington ratio of L_bol/L_Edd=1.1_{-0.2}^{+0.2} confirm that ULAS J1342+0928 is powered by a massive and actively accreting black hole. There is no significant difference in the emission line ratios such as SiIV/CIV and AlIII/CIV when compared to lower-redshift quasars in a similar luminosity range, which suggests early metal pollution of the broad-line-region clouds. This trend also holds for the FeII/MgII line ratio, known as a cosmic clock that traces the iron enrichment in the early universe. Different iron templates and continuum fitting ranges were used to explore how the FeII/MgII measurement changes as a function of spectral modeling. Quasars at even higher redshift or at fainter luminosity range (L_bol<10^46 erg s-1) are needed to probe the sites of early metal enrichment and a corresponding change in the FeII/MgII ratio.

We compare a sample of five high-resolution, high S/N Ly$\alpha$ forest spectra of bright $6<z \lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile fitting procedure for the inverted transmitted flux, $1-F$, similar to the widely used Voigt profile fitting of the transmitted flux $F$ at lower redshifts, to characterise the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominant in low temeperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution ($\leq $ 8 km/s) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of $T_0= (11000 \pm 1600,10500\pm 2100,12000 \pm 2200)$ K at $z= (5.4,5.6,5.8)$. This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to $z\sim5.3$, the width distribution of the observed transmission spikes is consistent with the range of $T_0$ caused by spatial fluctuations in the temperature-density relation.

Quasars at early redshifts ($z > 6$) with companion galaxies offer unique insights into the growth and evolution of the first supermassive black holes. Here, we report on a 150 ks Chandra observation of PSO J308.0416$-$21.2339, a $z=6.23$ quasar with a merging companion galaxy identified in [C II] and rest-frame UV emission. With $72.3^{+9.6}_{-8.6}$ net counts, we find that PSO J308.0416$-$21.2339 is powerful ($L_X = 2.31^{+1.14}_{-0.76} \times 10^{45}\ \textrm{erg}\,\textrm{s}^{-1}\,\textrm{cm}^{-2}$ in rest-frame $2.0-10.0$ keV) yet soft (spectral power-law index $\Gamma=2.39^{+0.37}_{-0.36}$ and optical-to-X-ray slope $\alpha_{\rm OX} = -1.41 \pm 0.11$). In addition, we detect three hard-energy photons $2.''0$ to the west of the main quasar, cospatial with the brightest UV emission of the merging companion. As no soft-energy photons are detected in the same area, this is potentially indicative of a highly-obscured source. With conservative assumptions, and accounting for both background fluctuations and the extended wings of the quasar's emission, these photons only have a probability $P=0.021$ of happening by chance. If confirmed by deeper observations, this system is the first high redshift quasar and companion individually detected in X-rays and is likely a dual AGN.

The discovery of quasars few hundred megayears after the Big Bang represents a major challenge to our understanding of black holes and galaxy formation and evolution. Their luminosity is produced by extreme gas accretion onto black holes, which already reached masses of 10$^9$ M$_\odot$ by z ~ 6. Simultaneously, their host galaxies form hundreds of stars per year, using up gas in the process. To understand which environments are able to sustain the rapid formation of these extreme sources we started a VLT/MUSE effort aimed at characterizing the surroundings of a sample of 5.7 < z < 6.6 quasars dubbed: the Reionization Epoch QUasar InvEstigation with MUSE (REQUIEM) survey. We here present results of our searches for extended Ly-Alpha halos around the first 31 targets observed as part of this program. Reaching 5-sigma surface brightness limits of 0.1-1.1 x 10$^{-17}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, we were able to unveil the presence of 12 Ly-Alpha nebulae, 8 of which are newly discovered. The detected nebulae show a variety of emission properties and morphologies with luminosities ranging from 8 x 10$^{42}$ to 2 x 10$^{44}$ erg/s, FWHMs between 300 and 1700 km/s, sizes < 30 pkpc, and redshifts consistent with those of the quasar host galaxies. As the first statistical and homogeneous investigation of the circum-galactic medium of massive galaxies at the end of the reionization epoch, the REQUIEM survey enables the study of the evolution of the cool gas surrounding quasars in the first 3 Gyr of the Universe. A comparison with the extended Ly-Alpha emission observed around bright (M$_{1450}$ < -25 mag) quasars at intermediate redshift indicates little variations on the properties of the cool gas from z ~ 6 to z ~ 3 followed by a decline in the average surface brightness down to z ~ 2.

Ratios of different ions of the same element encode ionization information independently from relative abundances in quasar absorption line systems, crucial for understanding the multiphase nature and origin of absorbing gas, particularly at $z>6$ where H I cannot be observed. Observational considerations have limited such studies to a small number of sightlines, with most surveys at $z>6$ focused upon the statistical properties of individual ions such as Mg II or C IV. Here we compare high- and low-ionization absorption within 69 intervening systems at $z>5$, including 16 systems at $z>6$, from Magellan/FIRE spectra of 47 quasars together with a Keck/HIRES spectrum of the `ultraluminous' $z=6.3$ quasar SDSSJ010013.02+280225.8. The highest redshift absorbers increasingly exhibit low-ionization species alone, consistent with previous single-ion surveys that show the frequency of Mg II is unchanging with redshift while C IV absorption drops markedly toward $z=6$. We detect no C IV or Si IV in half of all metal-line absorbers at $z>5.7$, with stacks not revealing any slightly weaker C IV just below our detection threshold, and most of the other half have $N_\mathrm{CII}>N_\mathrm{CIV}$. In contrast, only 20\% of absorbers at 5.0--5.7 lack high-ionization gas, and a search of 25 HIRES sightlines at $z\sim3$ yielded zero such examples. We infer these low-ionization high-redshift absorption systems may be analogous to metal-poor Damped Lyman-$\alpha$ systems ($\sim1\%$ of the absorber population at $z\sim3$), based on incidence rates and absolute and relative column densities. Simple photoionization models suggest that circumgalactic matter at redshift six has systematically lower chemical abundances and experiences a softer ionizing background relative to redshift three.

We use MUSE/VLT to conduct a survey of $z\sim3$ physical quasar pairs at close separation with a fast observation strategy. Our aim is twofold: (i) explore the Ly$\alpha$ glow around the faint-end of the quasar population; (ii) take advantage of the combined illumination of a quasar pair to unveil large-scale intergalactic structures extending between the two quasars. Here, we report the results for a quasar pair ($z=3.020,3.008$; $i=21.84,22.15$), separated by 11.6 arcsec (or 89 projected kpc). MUSE reveals filamentary Ly$\alpha$ structures extending between the two quasars with an average surface brightness of SB$_{\rm Ly\alpha}=1.8\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. Photoionization models of the constraints in the Ly$\alpha$, HeII, and CIV line emissions show that the emitting structures are intergalactic bridges with an extent between $\sim89$ and up to $\sim600$ kpc. Our models rule out the possibility that the structure extends for $\sim 2.9$ Mpc, i.e., the separation inferred from the uncertain systemic redshift difference of the quasars if the difference was only due to the Hubble flow. At the current spatial resolution and surface brightness limit, the average projected width of an individual bridge is about 35 kpc. We also detect a strong absorption in HI, NV, and CIV along the background sight-line at higher $z$, which we interpret as due to at least two components of cool, metal enriched, and relatively ionized CGM or IGM surrounding the quasar pair. Two additional HI absorbers are detected along both quasar sight-lines at $\sim -900$ and $-2800$ km s$^{-1}$ from the system, with the latter having associated CIV absorption only along the foreground quasar sight-line. The absence of galaxies in the MUSE field of view at the redshifts of these two absorbers suggests that they trace large-scale structures or expanding shells in front of the quasar pair.

Massive, quiescent galaxies are already observed at redshift $z\sim4$, i.e. $\sim$1.5 Gyr after the Big Bang. Current models predict them to be formed via massive, gas-rich mergers at $z>6$. Recent ALMA observations of the cool gas and dust in $z\gtrsim$6 quasars have discovered [CII]- and far infrared-bright galaxies adjacent to several quasars. In this work, we present sensitive imaging and spectroscopic follow-up observations, with HST/WFC3, Spitzer/IRAC, VLT/MUSE, Magellan/FIRE and LBT/LUCI-MODS, of ALMA-detected, dust-rich companion galaxies of four quasars at $z\gtrsim 6$, specifically acquired to probe their stellar content and unobscured star formation rate. Three companion galaxies do not show significant emission in the observed optical/IR wavelength range. The photometric limits suggest that these galaxies are highly dust-enshrouded, with unobscured star formation rates SFR$_{\rm UV}<$few M$_{\odot}$ yr$^{-1}$, and a stellar content of $M_{*}<$10$^{10}$ M$_{\odot}$ yr$^{-1}$. However, the companion to PJ167-13 shows bright rest-frame UV emission (F140W AB = 25.48). Its SED resembles that of a star-forming galaxy with a total SFR$\sim$50 M$_{\odot}$ yr$^{-1}$ and $M_{*}\sim 9 \times 10^{9}$ M$_{\odot}$. All the companion sources are consistent with residing on the galaxy main sequence at $z\sim$6. Additional, deeper data from future facilities, such as JWST, are needed in order to characterize these gas-rich sources in the first Gyr of cosmic history.

We report on ~0.35"(~2 kpc) resolution observations of the [CII] and dust continuum emission from five z>6 quasar host-companion galaxy pairs obtained with the Atacama Large Millimeter/submillimeter Array. The [CII] emission is resolved in all galaxies, with physical extents of 3.2-5.4 kpc. The dust continuum is on-average 40% more compact, which results in larger [CII] deficits in the center of the galaxies. However, the measured [CII] deficits are fully consistent with those found at lower redshifts. Four of the galaxies show [CII] velocity fields that are consistent with ordered rotation, while the remaining six galaxies show no clear velocity gradient. All galaxies have high (~80-200 km/s) velocity dispersions, consistent with the interpretation that the interstellar medium (ISM) of these high redshift galaxies is turbulent. By fitting the galaxies with kinematic models, we estimate the dynamical mass of these systems, which range between (0.3 -> 5.4) x 1E10 Msun. For the three closest separation galaxy pairs, we observe dust and [CII] emission from gas in between and surrounding the galaxies, which is an indication that tidal interactions are disturbing the gas in these systems. Although gas exchange in these tidal interactions could power luminous quasars, the existence of quasars in host galaxies without nearby companions suggests that tidal interactions are not the only viable method for fueling their active centers. These observations corroborate the assertion that accreting supermassive black holes do not substantially contribute to the [CII] and dust continuum emission of the quasar host galaxies, and showcase the diverse ISM properties of galaxies when the universe was less than one billion years old.

We present kpc-scale ALMA and HST imaging of the quasar PJ308-21 at $z$=$6.2342$, tracing dust, gas (via the [CII] 158 $\mu$m line) and young stars. At a resolution of $\sim0.3''$ ($\approx1.7$ kpc), the system is resolved over $>4''$ ($>$20 kpc). In particular, it features a main component, identified to be the quasar host galaxy, centered on the accreting supermassive black hole; and two other extended components on the West and East side, one redshifted and the other blueshifted relative to the quasar. The [CII] emission of the entire system stretches over $>$1500 km/s along the line of sight. All the components of the system are observed in dust, [CII], and rest-frame UV emission. The inferred [CII] luminosities [(0.9-4.6)$\times 10^9$ L$_\odot$], dust luminosities [(0.15-2.6)$\times10^{12}$ L$_\odot$], and rest-frame UV luminosities [(6.6-15)$\times10^{10}$ L$_\odot$], their ratios, and the implied gas/dust masses and star formation rates [11-290 M$_\odot$ yr$^{-1}$] are typical of high-redshift star-forming galaxies. A toy model of a single satellite galaxy that is tidally stripped by the interaction with the quasar host galaxy can account for the observed velocity and spatial extent of the two extended components. An outflow interpretation of the unique features in PJ308-21 is not supported by the data. PJ308-21 is thus one of the earliest galaxy mergers imaged at cosmic dawn.

The discoveries made over the past 20 years by Chandra and XMM-Newton surveys in conjunction with multiwavelength imaging and spectroscopic data available in the same fields have significantly changed the view of the supermassive black hole (SMBH) and galaxy connection. These discoveries have opened up several exciting questions that are beyond the capabilities of current X-ray telescopes and will need to be addressed by observatories in the next two decades. As new observatories peer into the early Universe, we will begin to understand the physics and demographics of SMBH infancy (at $z>6$) and investigate the influence of their accretion on the formation of the first galaxies ($\S$ 2.1). We will also be able to understand the accretion and evolution over the cosmic history (at $z\sim$1-6) of the full population of black holes in galaxies, including low accretion rate, heavily obscured AGNs at luminosities beyond the reach of current X-ray surveys ($\S$2.2 and $\S$2.3), enabling us to resolve the connection between SMBH growth and their environment.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

Despite extensive efforts, only two quasars have been found at $z>7$ to date due to a combination of low spatial density and high contamination from more ubiquitous Galactic cool dwarfs in quasar selection. This limits our current knowledge of the super-massive black hole (SMBH) growth mechanism and reionization history. In this letter, we report the discovery of a luminous quasar at $z=7.021$, DELS J003836.10$-$152723.6 (hereafter J0038$-$1527), selected using photometric data from DESI Legacy imaging Survey (DELS), Pan-STARRS1 (PS1) imaging Survey, as well as Wide-field Infrared Survey Explore ($WISE$) mid-infrared all-sky survey. With an absolute magnitude of $M_{1450}$=$-$27.1 and bolometric luminosity of $L_{\rm Bol}$=5.6$\times$10$^{13}$ $L_\odot$, J0038$-$1527 is the most luminous quasar known at $z>7$. Deep optical to near infrared spectroscopic observations suggest that J0038-1527 hosts a 1.3 billion solar mass BH accreting at the Eddington limit, with an Eddington ratio of 1.25$\pm$0.19. The CIV broad emission line of J0038$-$1527 is blue-shifted by more than 3000 km s$^{-1}$ to the systemic redshift. More detailed investigations of the high quality spectra reveal three extremely high velocity CIV broad absorption lines (BALs) with velocity from 0.08 to 0.14 times the speed of light and total balnicity index of more than 5000 km s$^{-1}$, suggesting the presence of relativistic outflows. J0038$-$1527 is the first quasar found at the epoch of reionization (EoR) with such strong outflows and provides a unique laboratory to investigate AGN feedback on the formation and growth of the most massive galaxies in the early universe.

This is the fourth paper in a series of publications aiming at discovering quasars at the epoch of reionization. In this paper, we expand our search for $z\sim 7$ quasars to the footprint of the Dark Energy Survey (DES) Data Release One (DR1), covering $\sim 5000$ deg$^2$ of new area. We select $z\sim 7$ quasar candidates using deep optical, near-infrared (near-IR) and mid-IR photometric data from the DES DR1, the VISTA Hemisphere Survey (VHS), the VISTA Kilo-degree Infrared Galaxy (VIKING) survey, the UKIRT InfraRed Deep Sky Surveys -- Large Area Survey (ULAS) and the unblurred coadds from the Wide-field Infrared Survey Explore ($WISE$) images (unWISE). The inclusion of DES and unWISE photometry allows the search to reach $\sim$ 1 magnitude fainter, comparing to our $z \gtrsim 6.5$ quasar survey in the northern sky (Wang et al. 2018). We report the initial discovery and spectroscopic confirmation of six new luminous quasars at $z>6.4$, including an object at $z=7.02$, the fourth quasar yet known at $z>7$, from a small fraction of candidates observed thus far. Based on the recent measurement of $z \sim 6.7 $ quasar luminosity function using the quasar sample from our survey in the northern sky, we estimate that there will be $\gtrsim$ 55 quasars at $z > 6.5$ at $M_{1450} < -24.5$ in the full DES footprint.

We present Atacama Large Millimeter Array 1mm observations of the rest-frame far-infrared (FIR) dust continuum in 27 quasars at redshifts 6.0 < z < 6.7. We detect FIR emission at >3sigma in all quasar host galaxies with flux densities at ~1900GHz in the rest-frame of 0.12 < S_rest,1900GHz < 5.9mJy, with a median (mean) flux density of 0.88mJy (1.59mJy). The implied FIR luminosities range from L_FIR = (0.27-13)x10^12 L_sun, with 74% of our quasar hosts having L_FIR > 10^12 L_sun. The estimated dust masses are M_dust = 10^7-10^9 M_sun. If the dust is heated only by star formation, then the star formation rates in the quasar host galaxies are between 50 and 2700 M_sun/yr. In the framework of the host galaxy-black hole coevolution model a correlation between ongoing black hole growth and star formation in the quasar host galaxy would be expected. However, combined with results from the literature to create a luminosity-limited quasar sample, we do not find a strong correlation between quasar UV luminosity (a proxy for ongoing black hole growth) and FIR luminosity (star formation in the host galaxy). The absence of such a correlation in our data does not necessarily rule out the coevolution model, and could be due to a variety of effects (including different timescales for black hole accretion and FIR emission).

This is the third paper in a series aims at finding reionzation-era quasars with the combination of DESI Legacy imaging Surveys (DELS) and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explore ($WISE$) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at $6.4\lesssim z \lesssim6.9$ from area of $\sim$13,020 deg$^2$, and present the quasar luminosity function (QLF) at $z\sim6.7$. The measured QLF follows $\Phi(L_{1450})\propto L_{1450}^{-2.35}$ in the magnitude range $27.6<M_{1450}<-25.5$. We determine the quasar comoving spatial density at $\langle z \rangle$=6.7 and $M_{1450}<-26.0$ to be $\rm 0.39\pm0.11 Gpc^{-3}$ and find that the exponential density evolution parameter to be $k=-0.78\pm0.18$ from $z\sim6$ to $z\sim6.7$, corresponding to a rapid decline by a factor of $\sim 6$ per unit redshift towards earlier epoch, a rate significantly faster than that at $z\sim 3- 5$. The cosmic time between $z\sim6$ and $z\sim6.7$ is only 121 Myrs. The quasar density declined by a factor of more than three within such short time requires that SMBHs must grow rapidly or they are less radiatively efficient at higher redshifts. We measured quasar comoving emissivity at $z\sim6.7$ which indicate that high redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line (BAL) quasar fraction at $z\gtrsim6.5$ is measured to be $\gtrsim$22%. In addition, we also report the discovery of additional five quasars at $z\sim6$ in the appendix.

We present the final results from our survey of luminous $z \sim $ 5.5 quasars. This is the first systematic quasar survey focusing on quasars at $z \sim$ 5.5, during the post-reionization epoch. It has been challenging to select quasars at $5.3 < z < 5.7$ using conventional color selections, due to their similar optical colors to those of late-type stars, especially M dwarfs. We developed a new selection technique for $z \sim$ 5.5 quasars based on optical, near-IR, and mid-IR photometry, using data from the Sloan Digital Sky Survey (SDSS), PanSTARR1 (PS1), the UKIRT Infrared Deep Sky Surveys - Large Area Survey, the UKIRT Hemisphere Survey, the VISTA Hemisphere Survey, and the Wide Field Infrared Survey Explorer ($WISE$), covering $\sim$ 11000 deg$^2$ of high galactic latitude sky. In this paper, we present the discovery of 15 new quasars at $z\sim 5.5$. Together with results from Yang et al. (2017), our survey provides a complete, flux-limited sample of 31 quasars at $5.3 \le z \le 5.7$. We measure the quasar spatial density at $z \sim 5.5$ and $M_{1450} < -26.2$. Our result is consistent with the rapid decline of quasar spatial density from $z = 5$ to 6, with $k=-0.66 \pm 0.05$ ($\rho(z) \propto 10^{kz}$). In addition, we present a new survey using optical colors only from the full PS1 area for luminous quasars at $z = 5.0 - 5.5$, beyond the SDSS footprint, and report the preliminary results from this survey, including 51 new quasars discovered at $4.61\le z \le5.71$.

We have used deep Optical Gravitational Lensing Experiment (OGLE-IV) images ($V \lesssim 23$ mag, $I \lesssim 23$ mag at $3\sigma$) of the Magellanic System, encompassing an area of $\sim$670 deg$^2$, to perform a search for high-$z$ quasar candidates. We combined the optical OGLE data with the mid-IR Wide-field Infrared Survey Explorer (WISE) 3.4/4.6/12 $\mu$m data, and devised a multi-color selection procedure. We have identified 33 promising sources and then spectroscopically observed the two most variable ones. We report the discovery of two high-$z$ quasars, OGLE J015531-752807 at a redshift $z=5.09$ and OGLE J005907-645016 at a redshift of $z=4.98$. The variability amplitude of both quasars at the rest-frame wavelength $\sim$1300\AA\ is much larger ($\sim$0.4 mag) than other quasars ($<0.15$ mag) at the same rest-frame wavelength but lower redshifts ($2<z<4$). To verify if there exist an increased variability amplitude in high-$z$ population of quasars, simply a larger sample of such sources with a decade long (or longer) light curves is necessary, which will be enabled by the Large Synoptic Survey Telescope (LSST) providing light curves for sources 3-4 mag fainter than OGLE.

Bright high-redshift quasars ($z>6$), hosting supermassive black holes ($\rm M_{BH} > 10^8 M_{\odot}$), are expected to reside in massive host galaxies embedded within some of the earliest and most massive galaxy overdensities. We analyze 1.2\,mm ALMA dust continuum maps of 35 bright quasars at $6<z<7$ and search the primary beam for excess dust continuum emission from sources with L$_{\rm IR} \gtrsim{ 10^{12}}$\,L$_{\odot}$ as evidence for early protoclusters. We compare the detection rates of continuum sources at $\geq5\sigma$ significance in the fields surrounding the quasars (A$_{\rm eff}$ = 4.3\,arcmin$^2$) with millimeter number counts in blank field surveys. We discover 15 millimeter sources in the fields excluding the quasars themselves, corresponding to an overdensity $\delta_{\rm gal} \equiv (N_{\rm gal} - N_{\rm exp})/N_{\rm exp} = -0.07\pm0.56$, consistent with no detected overdensity of dusty galaxies within 140 physical kpc of the quasars. However, the apparent lack of continuum overdensity does not negate the hypothesis that quasars live in overdense environments, as evidenced by strong [CII] overdensities found on the same scales to similarly-selected quasars. The small field of view of ALMA could miss a true overdensity if it exists on scales larger than 1\,cMpc, if the quasar is not centered in the overdensity, or if quasar feedback plays a role close to the quasar, but it is most likely that the large line of sight volume probed by a continuum survey will wash out a true overdensity signal. We discuss the necessary factors in determining the bias with which dusty star-forming galaxies trace true dark matter overdensities in order to improve upon overdensity searches in dust continuum.

Quasars are the most luminous non-transient objects known and as a result they enable studies of the Universe at the earliest cosmic epochs. Despite extensive efforts, however, the quasar ULAS J1120+0641 at z=7.09 has remained the only one known at z>7 for more than half a decade. Here we report observations of the quasar ULAS J134208.10+092838.61 (hereafter J1342+0928) at redshift z=7.54. This quasar has a bolometric luminosity of 4e13 times the luminosity of the Sun and a black hole mass of 8e8 solar masses. The existence of this supermassive black hole when the Universe was only 690 million years old---just five percent of its current age---reinforces models of early black-hole growth that allow black holes with initial masses of more than about 1e4 solar masses or episodic hyper-Eddington accretion. We see strong evidence of absorption of the spectrum of the quasar redwards of the Lyman alpha emission line (the Gunn-Peterson damping wing), as would be expected if a significant amount (more than 10 per cent) of the hydrogen in the intergalactic medium surrounding J1342+0928 is neutral. We derive a significant fraction of neutral hydrogen, although the exact fraction depends on the modelling. However, even in our most conservative analysis we find a fraction of more than 0.33 (0.11) at 68 per cent (95 per cent) probability, indicating that we are probing well within the reionization epoch of the Universe.

The recent detection of ULASJ1342+0928, a bright QSO at $z=7.54$, provides a powerful probe of the ionisation state of the intervening intergalactic medium, potentially allowing us to set strong constraints on the epoch of reionisation (EoR). Here we quantify the presence of Ly$\alpha$ damping wing absorption from the EoR in the spectrum of ULASJ1342+0928. Our Bayesian framework simultaneously accounts for uncertainties on: (i) the intrinsic QSO emission (obtained from reconstructing the Ly$\alpha$ profile from a covariance matrix of emission lines) and (ii) the distribution of HII regions during reionisation (obtained from three different 1.6$^3$ Gpc$^3$ simulations spanning the range of plausible EoR morphologies). Our analysis is complementary to that in the discovery paper (Ba\~nados et al.) and the accompanying method paper (Davies et al.) as it focuses solely on the damping wing imprint redward of Ly$\alpha$ ($1218 < \lambda < 1230$\AA), and uses a different methodology for (i) and (ii). We recover weak evidence for damping wing absorption. Our intermediate EoR model yields a volume-weighted neutral hydrogen fraction at $z=7.5$ of $\bar{x}_{\rm HI} = 0.21\substack{+0.17 \\ -0.19}$ (68 per cent). The constraints depend weakly on the EoR morphology. Our limits are lower than those presented previously, though they are consistent at ~1-1.5$\sigma$. We attribute this difference to: (i) a lower amplitude intrinsic Ly$\alpha$ profile obtained from our reconstruction pipeline, driven by correlations with other high-ionisation lines in the spectrum which are relatively weak; and (ii) only considering transmission redward of Ly$\alpha$ when computing the likelihood, which reduces the available constraining power but makes the results less model-dependent. Our results are consistent with previous estimates of the EoR history, and support the picture of a moderately extended EoR.

We present a 45 ks Chandra observation of the quasar ULAS J1342+0928 at z=7.54. We detect 14.0^{+4.8}_{-3.7} counts from the quasar in the observed-frame energy range 0.5-7.0 keV (6-sigma detection), representing the most distant non-transient astronomical source identified in X-rays to date. The present data are sufficient only to infer rough constraints on the spectral parameters. We find an X-ray hardness ratio of HR = -0.51^{+0.26}_{-0.28} between the 0.5-2.0 keV and 2.0-7.0 keV ranges and derive a power-law photon index of Gamma = 1.95^{+0.55}_{-0.53}. Assuming a typical value for high-redshift quasars of Gamma = 1.9, ULAS J1342+0928 has a 2-10 keV rest-frame X-ray luminosity of L_{2-10} = 11.6^{+4.3}_{-3.5} x 10^{44} erg/s. Its X-ray-to-optical power-law slope is alpha_{OX}=-1.67^{+0.16}_{-0.10}, consistent with the general trend indicating that the X-ray emission in the most bolometrically powerful quasars is weaker relative to their optical emission.

Quasars (QSOs) hosting supermassive black holes are believed to reside in massive halos harboring galaxy overdensities. However, many observations revealed average or low galaxy densities around $z\gtrsim6$ QSOs. This could be partly because they measured galaxy densities in only tens of arcmin$^2$ around QSOs and might have overlooked potential larger scale galaxy overdensities. Some previous studies also observed only Lyman break galaxies (LBGs, massive older galaxies) and missed low mass young galaxies like Ly$\alpha$ emitters (LAEs) around QSOs. Here we present observations of LAE and LBG candidates in $\sim700$ arcmin$^2$ around a $z=6.61$ luminous QSO using Subaru Telescope Suprime-Cam with narrow/broadbands. We compare their sky distributions, number densities and angular correlation functions with those of LAEs/LBGs detected in the same manner and comparable data quality in our control blank field. In the QSO field, LAEs and LBGs are clustering in 4-20 comoving Mpc angular scales, but LAEs show mostly underdensity over the field while LBGs are forming $30\times60$ comoving Mpc$^2$ large scale structure containing 3-$7\sigma$ high density clumps. The highest density clump includes a bright (23.78 mag in the narrowband) extended ($\gtrsim 16$ kpc) Ly$\alpha$ blob candidate, indicative of a dense environment. The QSO could be part of the structure but is not located exactly at any of the high density peaks. Near the QSO, LAEs show underdensity while LBGs average to $4\sigma$ excess densities compared to the control field. If these environments reflect halo mass, the QSO may not be in the most massive halo, but still in a moderately massive one.

During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral Universe, and enables measurements of the hydrogen neutral fraction $x_{\rm HI}(z)$ at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Ly$\alpha$. We combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift Universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its own ionizing radiation. Combining this detailed model with our new technique for predicting the quasar continuum and its associated uncertainty, we introduce a Bayesian statistical method to jointly constrain the neutral fraction of the Universe and the quasar lifetime from individual quasar spectra. We apply this methodology to the spectra of the two highest redshift quasars known, ULAS J1120+0641 and ULAS J1342+0928, and measured volume-averaged neutral fractions $\langle x_{\rm HI} \rangle(z=7.09)=0.48^{+0.26}_{-0.26}$ and $\langle x_{\rm HI} \rangle(z=7.54)=0.60^{+0.20}_{-0.23}$ (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of $10^3 \leq t_{\rm q} \leq 10^8$ years.

Measuring the proximity effect and the damping wing of intergalactic neutral hydrogen in quasar spectra during the epoch of reionization requires an estimate of the intrinsic continuum at rest-frame wavelengths $\lambda_{\rm rest}\sim1200$-$1260$ {\AA}. In contrast to previous works which used composite spectra with matched spectral properties or explored correlations between parameters of broad emission lines, we opted for a non-parametric predictive approach based on principal component analysis (PCA) to predict the intrinsic spectrum from the spectral properties at redder (i.e. unabsorbed) wavelengths. We decomposed a sample of $12764$ spectra of $z\sim2$-$2.5$ quasars from SDSS/BOSS into 10 red-side ($1280$ {\AA} $<\lambda_{\rm rest}<2900$ {\AA}) and 6 blue-side ($1180$ {\AA} $<\lambda_{\rm rest}<1280$ {\AA}) PCA basis spectra, and constructed a projection matrix to predict the blue-side coefficients from a fit to the red-side spectrum. We found that our method predicts the blue-side continuum with $\sim6$-$12\%$ precision and $\lesssim1\%$ bias by testing on the full training set sample. We then computed predictions for the blue-side continua of the two quasars currently known at $z>7$: ULAS J1120+0641 ($z=7.09$) and ULAS J1342+0928 ($z=7.54$). Both of these quasars are known to exhibit extreme emission line properties, so we individually calibrated the precision of the continuum predictions from similar quasars in the training set. We find that both $z>7$ quasars, and in particular ULAS J1342+0928, show signs of damping wing-like absorption at wavelengths redward of Ly$\alpha$.

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 [CII] 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 [CII] 158micron 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~6), making this source by far the FIR-brightest galaxy known at z>7.5. The [CII] emission is spatially unresolved, with an upper limit on the diameter of 7 kpc. Together with the measured FWHM of the [CII] line, this yields a dynamical mass of the host of <1.5x10^11 M_sun. Using standard assumptions about the dust temperature and emissivity, the NOEMA measurements give a dust mass of (0.6-4.3)x10^8 M_sun. The brightness of the [CII] luminosity, together with the high dust mass, imply active ongoing star formation in the quasar host. Using [CII]-SFR scaling relations, we derive star formation rates of 85-545 M_sun/yr 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 ~10^8 M_sun of dust implied by the observations.

We present final statistics from a survey for intervening MgII absorption towards 100 quasars with emission redshifts between $z=3.55$ and $z=7.08$. Using infrared spectra from Magellan/FIRE, we detect 279 cosmological MgII absorbers, and confirm that the incidence rate of $W_r>0.3 \AA$ MgII absorption per comoving path length does not evolve measurably between $z=0.25$ and $z=7$. This is consistent with our detection of seven new MgII systems at $z>6$, a redshift range that was not covered in prior searches. Restricting to relatively strong MgII systems ($W_r>1$\AA), there is significant evidence for redshift evolution. These systems roughly double in number density between $z=0$ and $z=2$-$3$, but decline by an order of magnitude from this peak by $z\sim 6$. This evolution mirrors that of the global star formation rate density, which could reflect a connection between star formation feedback and strong MgII absorbers. We compared our results to the Illustris cosmological simulation at $z=2$-$4$ by assigning absorption to catalogued dark-matter halos and by direct extraction of spectra from the simulation volume. To reproduce our results using the halo catalogs, we require circumgalactic (CGM) MgII envelopes within halos of progressively smaller mass at earlier times. This occurs naturally if we define the lower integration cutoff using SFR rather than mass. MgII profiles calculated directly from the Illustris volume yield far too few strong absorbers. We argue that this arises from unresolved phase space structure of CGM gas, particularly from turbulent velocities on sub-mesh scales. The presence of CGM MgII at $z>6$-- just $\sim 250$ Myr after the reionization redshift implied by Planck--suggests that enrichment of intra-halo gas may have begun before the presumed host galaxies' stellar populations were mature and dynamically relaxed. [abridged]

On 2017 August 17, gravitational waves were detected from a binary neutron star merger, GW170817, along with a coincident short gamma-ray burst, GRB170817A. An optical transient source, Swope Supernova Survey 17a (SSS17a), was subsequently identified as the counterpart of this event. We present ultraviolet, optical and infrared light curves of SSS17a extending from 10.9 hours to 18 days post-merger. We constrain the radioactively-powered transient resulting from the ejection of neutron-rich material. The fast rise of the light curves, subsequent decay, and rapid color evolution are consistent with multiple ejecta components of differing lanthanide abundance. The late-time light curve indicates that SSS17a produced at least ~0.05 solar masses of heavy elements, demonstrating that neutron star mergers play a role in r-process nucleosynthesis in the Universe.

Quasars are galaxies hosting accreting supermassive black holes; due to their brightness, they are unique probes of the early universe. To date, only few quasars have been reported at $z > 6.5$ ($<$800 Myr after the Big Bang). In this work, we present six additional $z \gtrsim 6.5$ quasars discovered using the Pan-STARRS1 survey. We use a sample of 15 $z \gtrsim 6.5$ quasars to perform a homogeneous and comprehensive analysis of this highest-redshift quasar population. We report four main results: (1) the majority of $z\gtrsim$6.5 quasars show large blueshifts of the broad CIV 1549\AA$\,$emission line compared to the systemic redshift of the quasars, with a median value $\sim$3$\times$ higher than a quasar sample at $z\sim$1; (2) we estimate the quasars' black hole masses (M$\rm_{BH}\sim$0.3$-$5 $\times$ 10$^{9}$ M$_{\odot}$) via modeling of the MgII 2798\AA$\,$emission line and rest-frame UV continuum; we find that quasars at high redshift accrete their material (with $\langle (L_{\mathrm{bol}}/L_{\mathrm{Edd}}) \rangle = 0.39$) at a rate comparable to a luminosity-matched sample at lower$-$redshift, albeit with significant scatter ($0.4$ dex); (3) we recover no evolution of the FeII/MgII abundance ratio with cosmic time; (4) we derive near zone sizes; together with measurements for $z\sim6$ quasars from recent work, we confirm a shallow evolution of the decreasing quasar near zone sizes with redshift. Finally, we present new millimeter observations of the [CII] 158 $\mu$m emission line and underlying dust continuum from NOEMA for four quasars, and provide new accurate redshifts and [CII]/infrared luminosities estimates. The analysis presented here shows the large range of properties of the most distant quasars.

We utilize the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) to search for extended Lyman-Alpha emission around the z~6.6 QSO J0305-3150. After carefully subtracting the point-spread-function, we reach a nominal 5-sigma surface brightness limit of SB = 1.9x10$^{-18}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, collapsing 5 wavelength slices centered at the expected location of the redshifted Lyman-Alpha emission (i.e. at 9256 Ang.). Current data suggest the presence (5-sigma, accounting for systematics) of a Lyman-Alpha nebula that extends for 9 kpc around the QSO. This emission is displaced and redshifted by 155 km/s with respect to the location of the QSO host galaxy traced by the [CII] emission line. The total luminosity is L = 3.0x10$^{42}$ erg/s. Our analysis suggests that this emission is unlikely to rise from optically thick clouds illuminated by the ionizing radiation of the QSO. It is more plausible that the Lyman-Alpha emission is due to fluorescence of the highly ionized optically thin gas. This scenario implies a high hydrogen volume density of n$_H$ ~ 6 cm$^{-3}$. In addition, we detect a Lyman-Alpha emitter (LAE) in the immediate vicinity of the QSO: i.e., with a projected separation of 12.5 kpc and a line-of-sight velocity difference of 560 km/s. The luminosity of the LAE is L = 2.1x10$^{42}$ erg/s and its inferred star-formation-rate is SFR ~ 1.3 M$_\odot$/yr. The probability of finding such a close LAE is one order of magnitude above the expectations based on the QSO-galaxy cross-correlation function. This discovery is in agreement with a scenario where dissipative interactions favour the rapid build-up of super-massive black holes at early Cosmic times.

The existence of massive ($10^{11}$ solar masses) elliptical galaxies by redshift z~4 (when the Universe was 1.5 billion years old) necessitates the presence of galaxies with star-formation rates exceeding 100 solar masses per year at z>6 (corresponding to an age of the Universe of less than 1 billion years). Surveys have discovered hundreds of galaxies at these early cosmic epochs, but their star-formation rates are more than an order of magnitude lower. The only known galaxies with very high star-formation rates at z>6 are, with only one exception, the host galaxies of quasars, but these galaxies also host accreting supermassive (more than $10^9$ solar masses) black holes, which probably affect the properties of the galaxies. Here we report observations of an emission line of singly ionized carbon ([CII] at a wavelength of 158 micrometres) in four galaxies at z>6 that are companions of quasars, with velocity offsets of less than 600 kilometers per second and linear offsets of less than 600 kiloparsecs. The discovery of these four galaxies was serendipitous; they are close to their companion quasars and appear bright in the far-infrared. On the basis of the [CII] measurements, we estimate star-formation rates in the companions of more than 100 solar masses per year. These sources are similar to the host galaxies of the quasars in [CII] brightness, linewidth and implied dynamical masses, but do not show evidence for accreting supermassive black holes. Similar systems have previously been found at lower redshift. We find such close companions in four out of twenty-five z>6 quasars surveyed, a fraction that needs to be accounted for in simulations. If they are representative of the bright end of the [CII] luminosity function, then they can account for the population of massive elliptical galaxies at z~4 in terms of cosmic space density.

Bright quasars, observed when the Universe was less than one billion years old (z>5.5), are known to host massive black holes (~10$^{9}$ M$_{\odot}$), and are thought to reside in the center of massive dark matter overdensities. In this picture, overdensities of galaxies are expected around high redshift quasars. However, observations based on the detection of Lyman Break Galaxies (LBGs) around these quasars do not offer a clear picture: this may be due to the uncertain redshift constraints of LBGs, which are selected through broad-band filters only. To circumvent such uncertainties, we here perform a search for Lyman Alpha Emitting galaxies (LAEs) in the field of the quasar PSO J215.1512-16.0417 at z~5.73, through narrow band, deep imaging with FORS2 at the VLT. We study an area of 37 arcmin$^{2}$, i.e. ~206 comoving Mpc$^{2}$ at the redshift of the quasar. We find no evidence for an overdensity of LAEs in the quasar field with respect to blank field studies. Possible explanations for these findings include that our survey volume is too small, or that the strong ionizing radiation from the quasar hinders galaxy formation in its immediate proximity. Another possibility is that these quasars are not situated in the dense environments predicted by some simulations.

Luminous quasars at z>5.6 can be studied in detail with the current generation of telescopes and provide us with unique information on the first gigayear of the universe. Thus far these studies have been statistically limited by the number of quasars known at these redshifts. Such quasars are rare and therefore wide-field surveys are required to identify them and multiwavelength data are needed to separate them efficiently from their main contaminants, the far more numerous cool dwarfs. In this paper, we update and extend the selection for z~6 quasars presented in Banados et al. (2014) using the Pan-STARRS1 (PS1) survey. We present the PS1 distant quasar sample, which currently consists of 124 quasars in the redshift range 5.6<z<6.7 that satisfy our selection criteria. Seventy-seven of these quasars have been discovered with PS1, and 63 of them are newly identified in this paper. We present composite spectra of the PS1 distant quasar sample. This sample spans a factor of ~20 in luminosity and shows a variety of emission line properties. The number of quasars at z>5.6 presented in this work almost double the quasars previously known at these redshifts, marking a transition phase from studies of individual sources to statistical studies of the high-redshift quasar population, which was impossible with earlier, smaller samples.

We present a survey of $z\sim3$ Ly$\alpha$ emitters (LAEs) within the fields of the VLT LBG Redshift Survey. The data encompasses 5 independent survey fields co-spatial with spectroscopic LBG data and covering a larger total area than previously analysed for LAE number counts and clustering. This affords an improved analysis over previous work by minimising the effects of cosmic variance and allowing the cross-clustering analysis of LAEs and LBGs. Our photometric sample consists of $\approx600$ LAE candidates, over an area of 1.07~deg$^2$, with equivalent widths of $\gtrsim65$~\AA\ and a flux limit of $\approx2\times10^{-17}$~erg~cm$^{-2}$~s$^{-1}$. From spectroscopic follow-up, we measured a success rate of $78\pm18\%$. We find the $R$-band continuum luminosity function to be $\sim10\times$ lower than the luminosity function of LBGs at this redshift, consistent with previous studies. Exploiting the large area of the survey, we estimate the LAE auto-correlation function and find a clustering length of $r_0=2.86\pm0.33~h^{-1}$~Mpc, low compared to the $z\sim3$ LBG population, but somewhat higher than previous LAE measurements. This corresponds to a median halo mass of $M_{\rm DM}=10^{11.0\pm0.3}~h^{-1}~$M$_{\odot}$. We present an analysis of clustering length versus continuum magnitude and find that the measurements for LAEs and LBGs are consistent at faint magnitudes. Our combined dataset of LAEs and LBGs allows us to measure, for the first time, the LBG-LAE cross-correlation, finding a clustering length of $r_0=3.29\pm0.57~h^{-1}$~Mpc and a LAE halo mass of $10^{11.1\pm0.4}~h^{-1}$~M$_{\odot}$. Overall, we conclude that LAEs inhabit primarily low mass halos, but form a relatively small proportion of the galaxy population found in such halos.

We present the results of our first year of quasar search in the on-going ESO public Kilo Degree Survey (KiDS) and VISTA Kilo-Degree Infrared Galaxy (VIKING) surveys. These surveys are among the deeper wide-field surveys that can be used to uncovered large numbers of z~6 quasars. This allows us to probe a more common population of z~6 quasars that is fainter than the well-studied quasars from the main Sloan Digital Sky Survey. From this first set of combined survey catalogues covering ~250 deg^2 we selected point sources down to Z_AB=22 that had a very red i-Z (i-Z>2.2) colour. After follow-up imaging and spectroscopy, we discovered four new quasars in the redshift range 5.8<z<6.0. The absolute magnitudes at a rest-frame wavelength of 1450 A are between -26.6 < M_1450 < -24.4, confirming that we can find quasars fainter than M^*, which at z=6 has been estimated to be between M^*=-25.1 and M^*=-27.6. The discovery of 4 quasars in 250 deg^2 of survey data is consistent with predictions based on the z~6 quasar luminosity function. We discuss various ways to push the candidate selection to fainter magnitudes and we expect to find about 30 new quasars down to an absolute magnitude of M_1450=-24. Studying this homogeneously selected faint quasar population will be important to gain insight into the onset of the co-evolution of the black holes and their stellar hosts.

The [CII] 158$\mu$m fine-structure line is known to trace regions of active star formation and is the main coolant of the cold, neutral atomic medium. In this \textit{Letter}, we report a strong detection of the [CII] line in the host galaxy of the brightest quasar known at $z>6.5$, the Pan-STARRS1 selected quasar PSO J036.5078+03.0498 (hereafter P036+03), using the IRAM NOEMA millimeter interferometer. Its [CII] and total far-infrared luminosities are $(5.8 \pm 0.7) \times 10^9 \,L_\odot$ and $(7.6\pm1.5) \times 10^{12}\,L_\odot$, respectively. This results in a $L_{[CII]} /L_{TIR}$ ratio of $\sim 0.8\times 10^{-3}$, which is at the high end for those found for active galaxies, though it is lower than the average found in typical main sequence galaxies at $z\sim 0$. We also report a tentative additional line which we identify as a blended emission from the $3_{22} - 3_{13}$ and $5_{23} - 4_{32}$ H$_2$O transitions. If confirmed, this would be the most distant detection of water emission to date. P036+03 rivals the current prototypical luminous J1148+5251 quasar at $z=6.42$, in both rest-frame UV and [CII] luminosities. Given its brightness and because it is visible from both hemispheres (unlike J1148+5251), P036+03 has the potential of becoming an important laboratory for the study of star formation and of the interstellar medium only $\sim 800\,$Myr after the Big Bang.

Radio-loud Active Galactic Nuclei at z~2-4 are typically located in dense environments and their host galaxies are among the most massive systems at those redshifts, providing key insights for galaxy evolution. Finding radio-loud quasars at the highest accessible redshifts (z~6) is important to study their properties and environments at even earlier cosmic time. They would also serve as background sources for radio surveys intended to study the intergalactic medium beyond the epoch of reionization in HI 21 cm absorption. Currently, only five radio-loud ($R=f_{\nu,5{\rm GHz}}/f_{\nu,4400\AA}>10$) quasars are known at z~6. In this paper we search for 5.5 < z < 7.2 quasars by cross-matching the optical Pan-STARRS1 and radio FIRST surveys. The radio information allows identification of quasars missed by typical color-based selections. While we find no good 6.4 < z <7.2 quasar candidates at the sensitivities of these surveys, we discover two new radio-loud quasars at z~6. Furthermore, we identify two additional z~6 radio-loud quasars which were not previously known to be radio-loud, nearly doubling the current z~6 sample. We show the importance of having infrared photometry for z>5.5 quasars to robustly classify them as radio-quiet or radio-loud. Based on this, we reclassify the quasar J0203+0012 (z=5.72), previously considered radio-loud, to be radio-quiet. Using the available data in the literature, we constrain the radio-loud fraction of quasars at z~6, using the Kaplan--Meier estimator, to be $8.1^{+5.0}_{-3.2}\%$. This result is consistent with there being no evolution of the radio-loud fraction with redshift, in contrast to what has been suggested by some studies at lower redshifts.

Luminous distant quasars are unique probes of the high redshift intergalactic medium (IGM) and of the growth of massive galaxies and black holes in the early universe. Absorption due to neutral Hydrogen in the IGM makes quasars beyond a redshift of z~6.5 very faint in the optical $z$-band, thus locating quasars at higher redshifts require large surveys that are sensitive above 1 micron. We report the discovery of three new z>6.5 quasars, corresponding to an age of the universe of <850 Myr, selected as z-band dropouts in the Pan-STARRS1 survey. This increases the number of known z>6.5 quasars from 4 to 7. The quasars have redshifts of z=6.50, 6.52, and 6.66, and include the brightest z-dropout quasar reported to date, PSO J036.5078+03.0498 with M_1450=-27.4. We obtained near-infrared spectroscopy for the quasars and from the MgII line we estimate that the central black holes have masses between 5x10^8 and 4x10^9 M_sun, and are accreting close to the Eddington limit (L_Bol/L_Edd=0.13-1.2). We investigate the ionized regions around the quasars and find near zone radii of R_NZ=1.5-5.2 proper Mpc, confirming the trend of decreasing near zone sizes with increasing redshift found for quasars at 5.7<z<6.4. By combining R_NZ of the PS1 quasars with those of 5.7<z<7.1 quasars in the literature, we derive a luminosity corrected redshift evolution of R_NZ,corrected=(7.2+/-0.2)-(6.1+/-0.7)x(z-6) Mpc. However, the large spread in R_NZ in the new quasars implies a wide range in quasar ages and/or a large variation in the neutral Hydrogen fraction along different lines of sight.

High-redshift quasars are currently the only probes of the growth of supermassive black holes and potential tracers of structure evolution at early cosmic time. Here we present our candidate selection criteria from the Panoramic Survey Telescope & Rapid Response System 1 and follow-up strategy to discover quasars in the redshift range 5.7<z<6.2. With this strategy we discovered eight new 5.7<z<6.0 quasars, increasing the number of known quasars at z>5.7 by more than 10%. We additionally recovered 18 previously known quasars. The eight quasars presented here span a large range of luminosities (-27.3 < M_{1450} < -25.4; 19.6 < z_ps1 < 21.2) and are remarkably heterogeneous in their spectral features: half of them show bright emission lines whereas the other half show a weak or no Ly$\alpha$ emission line (25% with rest-frame equivalent width of the Ly$\alpha$ + Nv line lower than 15{\AA}). We find a larger fraction of weak-line emission quasars than in lower redshift studies. This may imply that the weak-line quasar population at the highest redshifts could be more abundant than previously thought. However, larger samples of quasars are needed to increase the statistical significance of this finding.

High-redshift quasars are believed to reside in massive halos in the early universe and should therefore be located in fields with overdensities of galaxies, which are thought to evolve into galaxy clusters seen in the local universe. However, despite many efforts, the relationship between galaxy overdensities and z~6 quasars is ambiguous. This can possibly be attributed to the difficulty of finding galaxies with accurate redshifts in the vicinity of z~6 quasars. So far, overdensity searches around z~6 quasars have been based on studies of Lyman break galaxies (LBGs), which probe a redshift range of Delta z ~ 1. This range is large enough to select galaxies that may not be physically related to the quasar. We use deep narrow- and broadband imaging to study the environment of the z=5.72 quasar ULAS J0203+0012. The redshift range probed by our narrow-band selection of Lyman alpha emitters (LAEs) is Delta z ~ 0.1, which is significantly narrower than the LBG searches. This is the first time that LAEs were searched for near a z~6 quasar, in an effort to provide clues about the environments of quasars at the end of the epoch of reionization. We find no enhancement of LAEs in the surroundings of ULAS J0203+0012 in comparison with blank fields. We explore different explanations and interpretations for this non-detection of a galaxy overdensity, including that (1) the strong ionization from the quasar may prevent galaxy formation in its immediate vicinity and (2) high-redshift quasars may not reside in the center of the most massive dark matter halos.