Search

We introduce a probabilistic approach to select 6 ≤ $z$ ≤ 8 quasar candidates for spectroscopic follow-up, which is based on density estimation in the high-dimensional space inhabited by the optical and near-infrared photometry. Densities are modelled as Gaussian mixtures with principled accounting of errors using the extreme deconvolution (XD) technique, generalizing an approach successfully used to select lower redshift ($z$ ≤ 3) quasars. We train the probability density of contaminants on 1902 071 7-d flux measurements from the 1076 deg2 overlapping area from the Dark Energy Camera Legacy Survey (DECaLS) ($z$), VIKING (YJHKs), and unWISE (W1W2) imaging surveys, after requiring they dropout of DECaLS g and r, whereas the distribution of high-$z$ quasars are trained on synthetic model photometry. Extensive simulations based on these density distributions and current estimates of the quasar luminosity function indicate that this method achieves a completeness of $\ge 56{{\ \rm per\ cent}}$ and an efficiency of $\ge 5{{\ \rm per\ cent}}$ for selecting quasars at 6 < $z$ < 8 with JAB < 21.5. Among the classified sources are 8 known 6 < $z$ < 7 quasars, of which 2/8 are selected suggesting a completeness $\simeq 25{{\ \rm per\ cent}}$, whereas classifying the 6 known (JAB < 21.5) quasars at $z$ > 7 from the entire sky, we select 5/6 or a completeness of $\simeq 80{{\ \rm per\ cent}}$. The failure to select the majority of 6 < $z$ < 7 quasars arises because our quasar density model is based on an empirical quasar spectral energy distribution model that underestimates the scatter in the distribution of fluxes. This new approach to quasar selection paves the way for efficient spectroscopic follow-up of Euclid quasar candidates with ground-based telescopes and James Webb Space Telescope.

We build a multi-output generative model for quasar spectra and the properties of their black hole engines, based on a Gaussian process latent-variable model. This model treats every quasar as a vector of latent properties such that the spectrum and all physical properties of the quasar are associated with non-linear functions of those latent parameters; the Gaussian process kernel functions define priors on the function space. Our generative model is trained with a justifiable likelihood function that allows us to treat heteroscedastic noise and missing data correctly, which is crucial for all astrophysical applications. It can predict simultaneously unobserved spectral regions, as well as the physical properties of quasars in held-out test data. We apply the model to rest-frame ultraviolet and optical quasar spectra for which precise black hole masses (based on reverberation mapping measurements) are available. Unlike reverberation-mapping studies, which require multi-epoch data, our model predicts black hole masses from single-epoch spectra, even with limited spectral coverage. We demonstrate the capabilities of the model by predicting black hole masses and unobserved spectral regions. We find that we predict black hole masses at close to the best possible accuracy., accepted for publication in ApJ

The elemental abundances in the broad-line regions of high-redshift quasars trace the chemical evolution in the nuclear regions of massive galaxies in the early universe. In this work, we study metallicity-sensitive broad emission-line flux ratios in rest-frame UV spectra of 25 high-redshift (5.8 < z < 7.5) quasars observed with the VLT/X-shooter and Gemini/GNIRS instruments, ranging over $\log(M_{\rm{BH}}/M_{\odot})= 8.4-9.8$ in black hole mass and $\log(L_{\rm{bol}}/\rm{erg\, s}^{-1})= 46.7-47.7$ in bolometric luminosity. We fit individual spectra and composites generated by binning across quasar properties: bolometric luminosity, black hole mass, and blueshift of the \civ\, line, finding no redshift evolution in the emission-line ratios by comparing our high-redshift quasars to lower-redshift (2.0 < z < 5.0) results presented in the literature. Using Cloudy-based locally optimally-emitting cloud photoionisation model relations between metallicity and emission-line flux ratios, we find the observable properties of the broad emission lines to be consistent with emission from gas clouds with metallicity that are at least 2-4 times solar. Our high-redshift measurements also confirm that the blueshift of the CIV emission line is correlated with its equivalent width, which influences line ratios normalised against CIV. When accounting for the CIV blueshift, we find that the rest-frame UV emission-line flux ratios do not correlate appreciably with the black hole mass or bolometric luminosity., MNRAS accepted, 24 pages, 17 figures, 4 tables

The presence of excess scatter in the Ly-$\alpha$ forest at $z\sim 5.5$, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionisation. However, low data quality and systematic uncertainties complicate the use of Ly-$\alpha$ transmission as a precision probe of reionisation's end stages. In this paper, we assemble a sample of 67 quasar sightlines at $z>5.5$ with high signal-to-noise ratios of $>10$ per $\leq 15$ km s$^{-1}$ spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR $>20$ per pixel) spectra of 30 quasars at $z>5.7$. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-$\alpha$ systems. We present improved measurements of the mean Ly-$\alpha$ transmission over $4.93.5 \sigma$). Our results indicate that reionisation-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least $z=5.3$ ($t=1.1$ Gyr after the Big Bang). This is further evidence for a late end to reionisation., Comment: Accepted for publication in MNRAS. Full measurement datasets are available on the first author's website

We explore the kinematics of 27 z~6 quasar host galaxies observed in [CII]-158 micron ([CII]) emission with the Atacama Large Millimeter/sub-millimeter Array at a resolution of ~0.25''. We find that nine of the galaxies show disturbed [CII] emission, either due to a close companion galaxy or recent merger. Ten galaxies have smooth velocity gradients consistent with the emission arising from a gaseous disk. The remaining eight quasar host galaxies show no velocity gradient, suggesting that the gas in these systems is dispersion-dominated. All galaxies show high velocity dispersions with a mean of 129+-10 km/s. To provide an estimate of the dynamical mass within twice the half-light radius of the quasar host galaxy, we model the kinematics of the [CII] emission line using our publicly available kinematic fitting code, qubefit. This results in a mean dynamical mass of 5.0+-0.8(+-3.5) x 10^10 Msun. Comparison between the dynamical mass and the mass of the supermassive black hole reveals that the sample falls above the locally derived bulge mass--black hole mass relation at 2.4sigma significance. This result is robust even if we account for the large systematic uncertainties. Using several different estimators for the molecular mass, we estimate a gas mass fraction of >10%, indicating gas makes up a large fraction of the baryonic mass of z~6 quasar host galaxies. Finally, we speculate that the large variety in [CII] kinematics is an indication that gas accretion onto z~6 super massive black holes is not caused by a single precipitating factor., 33 pages. Kinematic fitting code is available at https://github.com/mneeleman/qubefit

X-ray emission from quasars has been detected up to redshift $z=7.5$, although only limited to a few objects at $z>6.5$. In this work, we present new Chandra observations of five $z>6.5$ quasars. By combining with archival Chandra observations of six additional $z>6.5$ quasars, we perform a systematic analysis on the X-ray properties of these earliest accreting supermassive black holes (SMBHs). We measure the black hole masses, bolometric luminosities ($L_{\rm bol}$), Eddington ratios ($\lambda_{\rm Edd}$), emission line properties, and infrared luminosities ($L_{\rm IR}$) of these quasars using infrared and sub-millimeter observations. Correlation analysis indicates that the X-ray bolometric correction (the factor that converts from X-ray luminosity to bolometric luminosity) decreases with increasing $L_{\rm bol}$, and that the UV/optical-to-X-ray ratio, $\alpha_{\rm ox}$, strongly correlates with $L_{\rm 2500}$, and moderately correlates with $\lambda_{\rm Edd}$ and blueshift of CIV emission lines. These correlations are consistent with those found in lower-$z$ quasars, indicating quasar accretion physics does not evolve with redshift. We also find that $L_{\rm IR}$ does not correlate with $L_{\rm 2-10 keV}$ in these luminous distant quasars, suggesting that the ratio of the SMBH growth rate and their host galaxy growth rate in these early luminous quasars are different from those of local galaxies. A joint spectral analysis of the X-ray detected $z>6.5$ quasars yields an average X-ray photon index of $ \Gamma=2.32^{+0.31}_{-0.30}$, steeper than that of low-$z$ quasars. By comparing it with the $\Gamma-\lambda_{\rm Edd}$ relation, we conclude that the steepening of $\Gamma$ for quasars at $z>6.5$ is mainly driven by their higher Eddington ratios., Comment: 16 pages, 10 figures, Accepted for publication in ApJ

We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range $6.3< z\le7.64$, including 32 quasars at $z>6.5$, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hole masses derived from the MgII emission lines are in the range $(0.3-3.6)\times10^{9}\,M_{\odot}$, which requires massive seed black holes with masses $\gtrsim10^{3-4}\,M_{\odot}$, assuming Eddington accretion since $z=30$. The Eddington ratio distribution peaks at $\lambda_{\rm Edd}\sim0.8$ and has a mean of 1.08, suggesting high accretion rates for these quasars. The CIV - MgII emission line velocity differences in our sample show an increase of CIV blueshift towards higher redshift, but the evolutionary trend observed from this sample is weaker than the previous results from smaller samples at similar redshift. The FeII/MgII flux ratios derived for these quasars up to $z=7.6$, compared with previous measurements at different redshifts, do not show any evidence of strong redshift evolution, suggesting metal-enriched environments in these quasars. Using this quasar sample, we create a quasar composite spectrum for $z>6.5$ quasars and find no significant redshift evolution of quasar broad emission lines and continuum slope, except for a blueshift of the CIV line. Our sample yields a strong broad absorption line quasar fraction of $\sim$24%, higher than the fractions in lower redshift quasar samples, although this could be affected by small sample statistics and selection effects., Comment: Accepted for publication in ApJ. 26 pages, 9 figures

Radio sources at the highest redshifts can provide unique information on the first massive galaxies and black holes, the densest primordial environments, and the epoch of reionization. The number of astronomical objects identified at z>6 has increased dramatically over the last few years, but previously only three radio-loud (R2500>10) sources had been reported at z>6, with the most distant being a quasar at z=6.18. Here we present the discovery and characterization of P172+18, a radio-loud quasar at z=6.823. This source has an MgII-based black hole mass of ~3x10^8 Msun and is one of the fastest accreting quasars, consistent with super-Eddington accretion. The ionized region around the quasar is among the largest measured at these redshifts, implying an active phase longer than the average lifetime of the z>6 quasar population. From archival data, there is evidence that its 1.4 GHz emission has decreased by a factor of two over the last two decades. The quasar's radio spectrum between 1.4 and 3.0 GHz is steep (alpha=-1.31) and has a radio-loudness parameter R2500~90. A second steep radio source (alpha=-0.83) of comparable brightness to the quasar is only 23.1" away (~120 kpc at z=6.82; projection probability, submitted to ApJ on Nov 29, 2020; accepted on Jan 31, 2021. See the companion paper by Momjian et al

Due to an error in the analysis of the ALMA data, the [C II] line fluxes Fline in 6 out of the 13 analyzed quasars were underestimated, and thus also based the inferred star formation rate (SFR) on the [C II] emission. We corrected these values and show the updated measurements for the six affected quasars in the corrected Table 2. The dust continuum estimates and their derived quantities were not affected by the error. Additionally, the observed peak frequency of the [C II] emission line nobs and thus the inferred systemic redshifts of the quasars were mildly offset, causing a shift of Dnobs ? 0.02 GHz and - Dzsub mm ? 0.0005. The derived proximity zone sizes were hence affected byDRp ? 0.04 pMpc. All updated values are shown in the corrected Tables 2 and 4. These changes did not affect any of the conclusions of the paper, nor did the figures in the paper change in any noticeable way. Furthermore, since the publication of our paper the offset between the peak of the Mg II emission line and the systemic redshift derived from submillimeter emission lines estimated by Schindler et al. (2020) changed slightly in their published version compared to the preliminary value we had used. Thus we use the now published velocity offset of Dv(Mg II-[C II]) = - - 391+ km s- 455 256 1 to estimate the systemic redshifts of the three quasars in our sample, for which no submillimeter data are available and therefore the redshift estimate is based on the Mg II emission. This resulted in minor changes of Dz ? 0.001, which are updated in Table 4.

We present the discovery of PSO J083.8371+11.8482, a weak emission line quasar with extreme star formation rate at $z=6.3401$. This quasar was selected from Pan-STARRS1, UHS, and unWISE photometric data. Gemini/GNIRS spectroscopy follow-up indicates a MgII-based black hole mass of $M_\mathrm{BH}=\left(2.0^{+0.7}_{-0.4}\right)\times10^9~M_\odot$ and an Eddington ratio of $L_\mathrm{bol}/L_\mathrm{Edd}=0.5^{+0.1}_{-0.2}$, in line with actively accreting supermassive black hole (SMBH) at $z\gtrsim6$. HST imaging sets strong constraint on lens-boosting, showing no relevant effect on the apparent emission. The quasar is also observed as a pure point-source with no additional emission component. The broad line region (BLR) emission is intrinsically weak and not likely caused by an intervening absorber. We found rest-frame equivalent widths of EW(Ly$\alpha$+NV) $=5.7\pm0.7$ Angstrom, EW(CIV) $\leq5.8$ Angstrom (3-sigma upper limit), and EW(MgII) $=8.7\pm0.7$ Angstrom. A small proximity zone size ($R_\mathrm{p}=1.2\pm0.4$ pMpc) indicates a lifetime of only $t_\mathrm{Q}=10^{3.4\pm0.7}$ years from the last quasar phase ignition. ALMA shows extended [CII] emission with a mild velocity gradient. The inferred far-infrared luminosity ($L_\mathrm{FIR}=(1.2\pm0.1)\times10^{13}\,L_\odot$) is one of the highest among all known quasar hosts at $z\gtrsim6$. Dust and [CII] emissions put a constraint on the star formation rate of SFR $=900-4900~M_\odot\,\mathrm{yr^{-1}}$, similar to that of hyper-luminous infrared galaxy. Considering the observed quasar lifetime and BLR formation timescale, the weak-line profile in the quasar spectrum is most likely caused by a BLR which is not yet fully formed rather than continuum boosting by gravitational lensing or a soft continuum due to super-Eddington accretion., Comment: 28 pages, 16 figures, 3 tables, accepted for publication in ApJ

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}$., Accepted for publication in ApJ

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, 21 pages, 8 figures, 3 tables, accepted for publication in ApJ

Observational studies are increasingly finding evidence against major mergers being the dominant mechanism responsible for triggering AGN. After studying the connection between major mergers and AGN with the highest Eddington ratios at z=2, we here expand our analysis to z0.3, and 25 mass- and redshift-matched control galaxies. To match the appearance of the two samples, we add synthetic point sources to the inactive comparison galaxies. The combined sample of AGN and inactive galaxies was independently ranked by 19 experts with respect to the degree of morphological distortion. We combine the resulting individual rankings into multiple overall rankings, from which we derive the respective major merger fractions of the two samples. With a best estimate of 0.41 $\pm$ 0.12 for the AGN host galaxies and 0.08 $\pm$ 0.06 for the inactive galaxies our results imply that our AGN host galaxies have a significantly higher merger rate, regardless of the observed wavelength or applied methodology. We conclude that although major mergers are an essential mechanism to trigger local high Eddington ratio AGNs at z=50% of this specific AGN subpopulation still remains unclear., Comment: 17 pages (plus 23 pages appendix), 8 figures, 3 tables; Accepted for publication in ApJ

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., Comment: 8 pages, 4 figures, published in the Astrophysical Journal. Revised to reflect accepted version

We investigate the distribution of companion galaxies around quasars using {\em Hubble Space Telescope} ({\em HST}) Advanced Camera for Surveys Wide Field Camera (ACS/WFC) archival images. Our master sample contains 532 quasars which have been observed by {\em HST} ACS/WFC, spanning a wide range of luminosity $(-3190\%)$ even for the faintest companions of interest. The number of physical companions is estimated by subtracting a background density from the number density of projected companions. We divide all the companions into three groups (faint, intermediate and bright) according to their fluxes. A control sample of galaxies is constructed to have similar redshift distribution and stellar mass range as the quasar sample using the data from {\em HST} deep fields. We find that quasars and control sample galaxies have similar numbers of faint and bright companions, while quasars show a $3.7\sigma$ deficit of intermediate companions compared to galaxies. The numbers of companions in all three groups do not show strong evolution with redshift, and the number of intermediate companions around quasars decreases with quasar luminosity. Assuming that merger-triggered quasars have entered the final coalescence stage during which individual companions are no longer detectable at large separations, our result is consistent with a picture in which a significant fraction of quasars is triggered by mergers., Comment: 18 pages, 10 figures, accepted by ApJ

We have designed the Extremely Luminous Quasar Survey (ELQS) to provide a highly complete census of unobscured UV-bright quasars during the cosmic noon, $z=2.8-5.0$. Here we report the discovery of 70 new quasars in the ELQS South Galactic Cap (ELQS-S) quasar sample, doubling the number of known extremely luminous quasars in $4,237.3\,\rm{deg}^2$ of the SDSS footprint. These observations conclude the ELQS and we present the properties of the full ELQS quasar catalog, containing 407 quasars over $11,838.5\,\rm{deg}^2$. Our novel ELQS quasar selection strategy resulted in unprecedented completeness at the bright end and allowed us to discover 109 new quasars in total. This marks an increase of $\sim36\%$ (109/298) to the known population at these redshifts and magnitudes, while we further are able to retain a selection efficiency of $\sim80\%$. On the basis of 166 quasars from the full ELQS quasar catalog, who adhere to the uniform criteria of the 2MASS point source catalog, we measure the bright-end quasar luminosity function (QLF) and extend it one magnitude brighter than previous studies. Assuming a single power law with exponential density evolution for the functional form of the QLF, we retrieve the best fit parameters from a maximum likelihood analysis. We find a steep bright-end slope of $\beta\approx-4.1$ and we can constrain the bright-end slope to $\beta\leq-3.4$ with $99\%$ confidence. The density is well modeled by the exponential redshift evolution, resulting in a moderate decrease with redshift ($\gamma\approx-0.4$)., Comment: 30 pages, 11 figures, ApJ accepted

This is the third paper in a series aimed at finding reionization-era quasars with the combination of DESI Legacy imaging Surveys (DELS), the Pan-STARRS1 (PS1) Survey, and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explorer (WISE) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at 6.4 ≲ z ≲ 6.9 from an area of ∼13,020 deg2, and present the quasar luminosity function (QLF) at z ∼ 6.7. The measured QLF follows {{Φ }}({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 < z> = 6.7 and M 1450 < -26.0 to be 0.39 ± 0.11 Gpc-3 and find the exponential density evolution parameter to be k = -0.78 ± 0.18 from z ∼ 6 to z ∼ 6.7, corresponding to a rapid decline by a factor of ∼6 per unit redshift toward earlier epochs. This indicates that the rapid decline of quasar spatial density at z > 5 that was found by previous works continues to z > 6, at a rate significantly faster than the average decline rate between z ∼ 3 and 5. We measured quasar comoving emissivity at z ∼ 6.7, which indicates that high-redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line quasar fraction at z ≳ 6.5 is measured to be ≳22%. In addition, we also report the discovery of six additional quasars at z ∼ 6 in the Appendix.

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., Comment: 53 pages, too many figures. Accepted for publication in ApJ. For access to the data and codes used in this work, please contact the authors or visit https://emastro.github.io/requiem/index.html

Distant quasars are unique tracers to study the formation of the earliest supermassive black holes (SMBHs) and the history of cosmic reionization. Despite extensive efforts, only two quasars have been found at $z\ge7.5$, due to a combination of their low spatial density and the high contamination rate in quasar selection. We report the discovery of a luminous quasar at $z=7.642$, J0313$-$1806, the most distant quasar yet known. This quasar has a bolometric luminosity of $3.6\times10^{13} L_\odot$. Deep spectroscopic observations reveal a SMBH with a mass of $(1.6\pm0.4) \times10^9M_\odot$ in this quasar. The existence of such a massive SMBH just $\sim$670 million years after the Big Bang challenges significantly theoretical models of SMBH growth. In addition, the quasar spectrum exhibits strong broad absorption line (BAL) features in CIV and SiIV, with a maximum velocity close to 20% of the speed of light. The relativistic BAL features, combined with a strongly blueshifted CIV emission line, indicate that there is a strong active galactic nucleus (AGN) driven outflow in this system. ALMA observations detect the dust continuum and [CII] emission from the quasar host galaxy, yielding an accurate redshift of $7.6423 \pm 0.0013$ and suggesting that the quasar is hosted by an intensely star-forming galaxy, with a star formation rate of $\rm\sim 200 ~M_\odot ~yr^{-1}$ and a dust mass of $\sim7\times10^7~M_\odot$. Followup observations of this reionization-era BAL quasar will provide a powerful probe of the effects of AGN feedback on the growth of the earliest massive galaxies., 4 figures, 1 table; accepted for publication in ApJL

We present X-SHOOTER near-infrared spectroscopy of a large sample of 38 luminous ($M_{1450}=-29.0$ to $-24.4$) quasars at $5.78, Comment: ApJ accepted

We present the North Galactic Cap sample of the Extremely Luminous Quasar Survey (ELQS-N), which targets quasars with $M_{1450}70\%$ over $z\sim3.0-5.0$ at $m_{i}\lesssim17.5$, limited toward fainter magnitudes by the depth of the Two Micron All Sky Survey (2MASS). The presented quasar catalog consists of a total of 270 objects, of which 39 are newly identified in this work with spectroscopy obtained at the Vatican Advanced Technology Telescope and the MMT $6.5\,\rm{m}$ telescope. In addition to the high completeness, which allowed us to discover new quasars in the already well-surveyed SDSS North Galactic Cap, the efficiency of our selection is relatively high at $\sim79\%$. Using 120 objects of this quasar sample we are able to extend the previously measured optical quasar luminosity function (QLF) by one magnitude toward the bright end at $2.8 \leq z \leq 4.5$. A first analysis of the QLF suggests a relatively steep bright-end slope of $\beta\approx-4$ for this sample. This result contrasts with previous results in the same redshift range, which find a much flatter slope around $\beta\sim-2.5$, but agrees with recent measurements of the bright-end slope at lower and higher redshifts. Our results constrain the bright-end slope at $z=2.8-4.5$ to $\beta, Comment: 23 pages, 13 figures, submitted to ApJ

Studies of the most luminous quasars at high redshift directly probe the evolution of the most massive black holes in the early Universe and their connection to massive galaxy formation. However, extremely luminous quasars at high redshift are very rare objects. Only wide area surveys have a chance to constrain their population. The Sloan Digital Sky Survey (SDSS) has so far provided the most widely adopted measurements of the quasar luminosity function (QLF) at $z>3$. However, a careful re-examination of the SDSS quasar sample revealed that the SDSS quasar selection is in fact missing a significant fraction of $z\gtrsim3$ quasars at the brightest end. We have identified the purely optical color selection of SDSS, where quasars at these redshifts are strongly contaminated by late-type dwarfs, and the spectroscopic incompleteness of the SDSS footprint as the main reasons. Therefore we have designed the Extremely Luminous Quasar Survey (ELQS), based on a novel near-infrared JKW2 color cut using WISE AllWISE and 2MASS all-sky photometry, to yield high completeness for very bright ($m_{\rm{i}} < 18.0$) quasars in the redshift range of $3.0\leq z\leq5.0$. It effectively uses random forest machine-learning algorithms on SDSS and WISE photometry for quasar-star classification and photometric redshift estimation. The ELQS will spectroscopically follow-up $\sim 230$ new quasar candidates in an area of $\sim12000\,\rm{deg}^2$ in the SDSS footprint, to obtain a well-defined and complete quasars sample for an accurate measurement of the bright-end quasar luminosity function at $3.0\leq z\leq5.0$. In this paper we present the quasar selection algorithm and the quasar candidate catalog., 16 pages, 8 figures, 9 tables; ApJ in press

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.65.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., Comment: Accepted by ApJS. Machine readable tables and an up-to-date census of z>5.6 quasars are available at https://users.obs.carnegiescience.edu/~ebanados/high-z-qsos.html

We present the results of the Extremely Luminous Quasar Survey in the $3\pi$ survey of the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS; PS1). This effort applies the successful quasar selection strategy of the Extremely Luminous Survey in the Sloan Digital Sky Survey footprint ($\sim12,000\,\rm{deg}^2$) to a much larger area ($\sim\rm{21486}\,\rm{deg}^2$). This spectroscopic survey targets the most luminous quasars ($M_{1450}\le-26.5$; $m_{i}\le18.5$) at intermediate redshifts ($z\ge2.8$). Candidates are selected based on a near-infrared JKW2 color cut using WISE AllWISE and 2MASS photometry to mainly reject stellar contaminants. Photometric redshifts ($z_{\rm{reg}}$) and star-quasar classifications for each candidate are calculated from near-infrared and optical photometry using the supervised machine learning technique random forests. We select 806 quasar candidates at $z_{\rm{reg}}\ge2.8$ from a parent sample of 74318 sources. After exclusion of known sources and rejection of candidates with unreliable photometry, we have taken optical identification spectra for 290 of our 334 good PS-ELQS candidates. We report the discovery of 190 new $z\ge2.8$ quasars and an additional 28 quasars at lower redshifts. A total of 44 good PS-ELQS candidates remain unobserved. Including all known quasars at $z\ge2.8$, our quasar selection method has a selection efficiency of at least $77\%$. At lower declinations $-30\le\rm{Decl.}\le0$ we approximately triple the known population of extremely luminous quasars. We provide the PS-ELQS quasar catalog with a total of 592 luminous quasars ($m_{i}\le18.5$, $z\ge2.8$). This unique sample will not only be able to provide constraints on the volume density and quasar clustering of extremely luminous quasars, but also offers valuable targets for studies of the intergalactic medium., Comment: 34 pages, 10 figures, accepted to ApJS

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., ApJL in press

At redshifts beyond $z{\sim}1$ measuring the black hole galaxy relations proves to be a difficult task. The bright light of the AGN aggravates deconvolution of black hole and galaxy properties. On the other hand high redshift data on these relations is vital to understand in what ways galaxies and black holes co-evolve and in what ways they don't. In this work we use black hole (BHMDs) and stellar mass densities (SMDs) to constrain the possible co-evolution of black holes with their host galaxies since $z{\sim}5$. The BHMDs are calculated from quasar luminosity functions (QLF) using the Soltan argument, while we use integrals over stellar mass functions (SMFs) or the star formation rate density to obtain values for the stellar mass density. We find that both quantities grow in lock-step below redshifts of $z{\sim}3$ with a non-evolving BHMD to SMD ratio. A fit to the data assuming a power law relation between the BHMD and the SMD yields exponents around unity ($1.0{-}1.5$). Up to $z{\sim}5$ the BHMD to SMD ratio doesn't show a strong evolution given the larger uncertainty in the completeness of high-redshift datasets. Our results, always applying the same analysis technique, seem to be consistent across all adopted data sets., 12 pages, 4 figures, submitted to the Astrophysical Journal

Asteroseismology of non-radial pulsations in Hot B Subdwarfs (sdB stars) offers a unique view into the interior of core-helium-burning stars. Ground-based and space-borne high precision light curves allow for the analysis of pressure and gravity mode pulsations to probe the structure of sdB stars deep into the convective core. As such asteroseismological analysis provides an excellent opportunity to test our understanding of stellar evolution. In light of the newest constraints from asteroseismology of sdB and red clump stars, standard approaches of convective mixing in 1D stellar evolution models are called into question. The problem lies in the current treatment of overshooting and the entrainment at the convective boundary. Unfortunately no consistent algorithm of convective mixing exists to solve the problem, introducing uncertainties to the estimates of stellar ages. Three dimensional simulations of stellar convection show the natural development of an overshooting region and a boundary layer. In search for a consistent prescription of convection in one dimensional stellar evolution models, guidance from three dimensional simulations and asteroseismological results is indispensable., Comment: 10 pages, 3 figures, proceedings of the Joint TASC2 - KASC9 Workshop - SPACEINN - HELAS8 Conference "Seismology of the Sun and the Distant Stars 2016". To be published by the EPJ Web of Conferences

We present a new algorithm to estimate quasar photometric redshifts (photo-$z$s), by considering the asymmetries in the relative flux distributions of quasars. The relative flux models are built with multivariate Skew-t distributions in the multi-dimensional space of relative fluxes as a function of redshift and magnitude. For 151,392 quasars in the SDSS, we achieve a photo-$z$ accuracy, defined as the fraction of quasars with the difference between the photo-$z$ $z_p$ and the spectroscopic redshift $z_s$, $|\Delta z| = |z_s-z_p|/(1+z_s)$ within 0.1, of 74%. Combining the WISE W1 and W2 infrared data with the SDSS data, the photo-$z$ accuracy is enhanced to 87%. Using the Pan-STARRS1 or DECaLS photometry with WISE W1 and W2 data, the photo-$z$ accuracies are 79% and 72%, respectively. The prior probabilities as a function of magnitude for quasars, stars and galaxies are calculated respectively based on (1) the quasar luminosity function; (2) the Milky Way synthetic simulation with the Besan\c{c}on model; (3) the Bayesian Galaxy Photometric Redshift estimation. The relative fluxes of stars are obtained with the Padova isochrones, and the relative fluxes of galaxies are modeled through galaxy templates. We test our classification method to select quasars using the DECaLS $g$, $r$, $z$, and WISE W1 and W2 photometry. The quasar selection completeness is higher than 70% for a wide redshift range $0.5, Comment: 22 pages, 17 figure, accepted by AJ. The code is available at https://doi.org/10.5281/zenodo.1014407

Stellar evolution calculations have had great success reproducing the observed atmospheric properties of different classes of stars. Recent detections of g-mode pulsations in evolved He burning stars allow a rare comparison of their internal structure with stellar models. Asteroseismology of subdwarf B stars suggests convective cores of $0.22-0.28\,M_\odot$, $\gtrsim 45\,\%$ of the total stellar mass. Previous studies found significantly smaller convective core masses ($\lesssim 0.19\,M_\odot$) at a comparable evolutionary stage. We evolved stellar models with MESA (Modules for Experiments in Stellar Astrophysics) to explore how well the interior structure inferred from asteroseismology can be reproduced by standard algorithms. Our qualitative evolutionary paths, position in the $\log g-T_{\rm{eff}}$ diagram and model timescales are consistent with previous results. SdB masses from our full evolutionary sequences fall within the range of the empirical sdB mass distribution, but are nearly always lower than the median. Using standard MLT with atomic diffusion we find convective core masses of $\sim 0.17-0.18 M_\odot$, averaged over the entire sdB lifetime. We can increase the convective core sizes to be as large as those inferred from asteroseismology, but only for extreme values of the overshoot parameter (overshoot gives numerically unstable and physically unrealistic behavior at the boundary). High resolution three-dimensional (3D) simulations of turbulent convection in stars suggest that the Schwarzschild criterion for convective mixing sytematically {\it underestimates the actual extent of mixing because a boundary layer forms.} Accounting for this would decrease the errors in both sdB total and convective core masses., Comment: This paper has been updated to a standard case with OPALII opacities. 13 pages, 6 figures. Submitted to ApJ, Accepted May 2015