Search

High redshift quasars ($z>5$) that also shine brightly at radio wavelengths are unique signposts of supermassive black hole activity in the early universe. However, bright radio sources at $z\ge5$ are extremely rare and therefore we have started a campaign to search for new high-$z$ quasars by combining an optical dropout selection driven by the $g$, $r$, and $z$ bands from the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys with low-frequency radio observations from the LOFAR Two-metre Sky Survey (LoTSS). Currently, LoTSS covers a large fraction of the northern sky (5720 deg$^2$) to such a depth (median noise level of 83 $\mu$Jy beam$^{-1}$) that about 30% of the general quasar population is detected $-$ which is a factor of 5-10 more than previous large sky radio surveys such as NVSS and FIRST, respectively. In this paper, we present the discovery of 20 new quasars (and the independent confirmation of 4) between $4.9\leq z\leq 6.6$. Out of the 24 quasars, 21 satisfy the traditional radio-loudness criterion of $R=f_{5\text{GHz}}/f_{4400A} > 10$, with the full sample spanning $R\sim$6-1000, thereby more than doubling the sample of known radio-loud quasars at $z \ge 5$. Our radio detection requirement strongly decreases the contamination of stellar sources and allows one to select these quasars in a broad redshift range. Despite selecting our quasar candidates using fewer and less conservative colour restrictions, both the optical and near-infrared colours, Ly$\alpha$ emission line properties, and dust reddening, $E(B-V)$, measurements of our quasar sample do not deviate from the known radio-quiet quasar population, suggesting similar optical quasar properties of the radio-loud and radio-quiet quasar population at high-$z$. Our campaign demonstrates the potential for discovering new high-$z$ quasar populations through next generation radio continuum surveys., Comment: 19 pages, 11 figures, Tab A1 corrected, accepted for publication in A&A

In this work we present a new hybrid method to simulate the thermal effects of the reionization in cosmological hydrodynamical simulations. The method improves upon the standard approach used in simulations of the intergalactic medium (IGM) and galaxy formation without a significant increase of the computational cost allowing for efficient exploration of the parameter space. The method uses a small set of phenomenological input parameters and combines a semi-numerical reionization model to solve for the topology of reionization and an approximate model of how reionization heats the IGM, with the massively parallel \texttt{Nyx} hydrodynamics code, specifically designed to solve for the structure of diffuse IGM gas. We have produced several large-scale high resolution cosmological hydrodynamical simulations ($2048^3$, $L_{\rm box} = 40$ Mpc/h) with different instantaneous and inhomogeneous HI reionization models that use this new methodology. We study the IGM thermal properties of these models and find that large scale temperature fluctuations extend well beyond the end of reionization. Analyzing the 1D flux power spectrum of these models, we find up to $\sim 50\%$ differences in the large scale properties (low modes, $k\lesssim0.01$ s/km) of the post-reionization power spectrum due to the thermal fluctuations. We show that these differences could allow one to distinguish between different reionization scenarios already with existing Ly$\alpha$ forest measurements. Finally, we explore the differences in the small-scale cutoff of the power spectrum and we find that, for the same heat input, models show very good agreement provided that the reionization redshift of the instantaneous reionization model happens at the midpoint of the inhomogeneous model., Comment: 24 pages, 16 figures. Accepted by MNRAS. Minor changes to match published version

Motivated by the recent discovery of rare Enormous Lyman-Alpha Nebulae (ELAN) around z~2 quasars, we have initiated a long-term observational campaign with the MUSE instrument to directly uncover the astrophysics of the gas around quasars. We present here the first 61 targets of our effort under the acronym QSO MUSEUM (Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMission). These quasars have a median redshfit of z=3.17, absolute $i$ magnitude in the range $-29.67\leq M_i(z=2)\leq-27.03$, and different levels of radio-loudness. This sample unveils diverse specimens of Ly$\alpha$ nebulosities extending for tens of kiloparsecs around these quasars (on average out to a maximum projected distance of 80 kpc) above a surface brightness SB$>8.8\times10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ ($2\sigma$). The bulk of the extended Ly$\alpha$ emission is within R<50 kpc, and is characterized by relatively quiescent kinematics, with average velocity dispersions of $\langle \sigma_{\rm Ly\alpha}\rangle < 400$~km~s$^{-1}$. Therefore, the motions within all these Ly$\alpha$ nebulosities have amplitudes consistent with gravitational motions expected in dark matter halos hosting quasars at these redshifts, possibly reflecting the complexity in propagating a fast wind on large scales. Our current data suggest a combination of photoionization and resonant scattering as powering mechanisms of the Ly$\alpha$ emission. We discover the first $z\sim3$ ELAN, which confirms a very low probability ($\sim1\%$) of occurrence of such extreme systems at these cosmic epochs. Finally, we discuss the redshift evolution currently seen in extended Ly$\alpha$ emission around radio-quiet quasars from $z\sim3$ to $z\sim2$, concluding that it is possibly linked to a decrease of cool gas mass within the quasars' CGM from $z\sim3$ to $z\sim2$, and thus to the balance of cool vs hot media., Comment: 48 pages, 3 page-long tables, 28 figures, 3 appendices. Submitted to MNRAS. Comments are welcome

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., Comment: Updated to match the final journal version

The latest measurements of CMB electron scattering optical depth reported by Planck significantly reduces the allowed space of HI reionization models, pointing toward a later ending and/or less extended phase transition than previously believed. Reionization impulsively heats the intergalactic medium (IGM) to $\sim10^4$ K, and owing to long cooling and dynamical times in the diffuse gas, comparable to the Hubble time, memory of reionization heating is retained. Therefore, a late ending reionization has significant implications for the structure of the $z\sim5-6$ Lyman-$\alpha$ (ly$\alpha$) forest. Using state-of-the-art hydrodynamical simulations that allow us to vary the timing of reionization and its associated heat injection, we argue that extant thermal signatures from reionization can be detected via the Ly$\alpha$ forest power spectrum at $5< z<6$. This arises because the small-scale cutoff in the power depends not only the the IGM temperature at these epochs, but is also particularly sensitive to the pressure smoothing scale set by the IGM full thermal history. Comparing our different reionization models with existing measurements of the Ly$\alpha$ forest flux power spectrum at $z=5.0-5.4$, we find that models satisfying Planck's $\tau_e$, constraint favor a moderate amount of heat injection consistent with galaxies driving reionization, but disfavoring quasar-driven scenarios. We explore the impact of different reionization histories and heating models on the shape of the power spectrum, and find that they can produce similar effects, but argue that this degeneracy can be broken with high enough quality data. We study the feasibility of measuring the flux power spectrum at $z\simeq 6$ using mock quasar spectra and conclude that a sample of $\sim10$ high-resolution spectra with an attainable signal-to-noise ratio will allow distinguishing between different reionization scenarios., Comment: 14 pages, 7 Figures. Matches version accepted by ApJ. UVB models used in this work can be downloaded from http://thermal.joseonorbe.com/uvbtables.html

We present the most precise estimate to date of the clustering of quasars on very small scales, based on a sample of 47 binary quasars with magnitudes of $g<20.85$ and proper transverse separations of $\sim 25\,h^{-1}$\,kpc. Our sample of binary quasars, which is about 6 times larger than any previous spectroscopically confirmed sample on these scales, is targeted using a Kernel Density Estimation technique (KDE) applied to Sloan Digital Sky Survey (SDSS) imaging over most of the SDSS area. Our sample is "complete" in that all of the KDE target pairs with $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc in our area of interest have been spectroscopically confirmed from a combination of previous surveys and our own long-slit observational campaign. We catalogue 230 candidate quasar pairs with angular separations of $<8\arcsec$, from which our binary quasars were identified. We determine the projected correlation function of quasars ($\bar W_{\rm p}$) in four bins of proper transverse scale over the range $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc. The implied small-scale quasar clustering amplitude from the projected correlation function, integrated across our entire redshift range, is $A=24.1\pm3.6$ at $\sim 26.6 ~h^{-1}$\,kpc. Our sample is the first spectroscopically confirmed sample of quasar pairs that is sufficiently large to study how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ kpc. We find that empirical descriptions of how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ Mpc also adequately describe the evolution of quasar clustering at $\sim 25 ~h^{-1}$ kpc., Comment: 16 pages, 8 figures, 6 tables, Accepted for publication in MNRAS

We present the first results from an on-going survey to characterize the circumgalactic medium (CGM) of the massive high-redshift galaxies detected as submillimeter galaxies (SMGs). We constructed a parent sample of 163 SMG-QSO pairs with separations less than $\sim$36 arcsec by cross-matching far-infrared-selected galaxies from Herschel with spectroscopically confirmed QSOs. The Herschel sources were selected to match the properties of SMGs. We determined the sub-arcsecond positions of six Herschel sources with the Very Large Array and obtained secure redshift identification for three of those with near-infrared spectroscopy. The QSO sightlines probe transverse proper distances of 112, 157, and 198 kpc at foreground redshifts of 2.043, 2.515, and 2.184, respectively, which are comparable to the virial radius of the $\sim10^{13}$ Msun halos expected to host SMGs. High-quality absorption-line spectroscopy of the QSOs reveals systematically strong HI Lyman-alpha absorption around all three SMGs, with rest-frame equivalent widths of $\sim2-3$ \AA. However, none of the three absorbers exhibits compelling evidence for optically thick HI gas or metal absorption, in contrast to the dominance of strong neutral absorbers in the CGM of luminous $z \sim 2$ QSOs. The low covering factor of optically thick HI gas around SMGs tentatively indicates that SMGs may not have as prominent cool gas reservoirs in their halos as the co-eval QSOs and that they may inhabit less massive halos than previously thought., Comment: ApJ accepted after minor revision

The lifetime of quasars is fundamental for understanding the growth of supermassive black holes, and is an important ingredient in models of the reionization of the intergalactic medium. However, despite various attempts to determine quasar lifetimes, current estimates from a variety of methods are uncertain by orders of magnitude. This work combines cosmological hydrodynamical simulations and 1D radiative transfer to investigate the structure and evolution of the He II Ly$\alpha$ proximity zones around quasars at $z \simeq 3-4$. We show that the time evolution in the proximity zone can be described by a simple analytical model for the approach of the He II fraction $x_{\rm HeII}\left( t \right)$ to ionization equilibrium, and use this picture to illustrate how the transmission profile depends on the quasar lifetime, quasar UV luminosity, and the ionization state of helium in the ambient IGM (i.e. the average He II fraction, or equivalently the metagalactic He II ionizing background). A significant degeneracy exists between the lifetime and the average He II fraction, however the latter can be determined from measurements of the He II Ly$\alpha$ optical depth far from quasars, allowing the lifetime to be measured. We advocate stacking existing He II quasar spectra at $z\sim 3$, and show that the shape of this average proximity zone profile is sensitive to lifetimes as long as $\sim 30$ Myr. At higher redshift $z\sim 4$ where the He II fraction is poorly constrained, degeneracies will make it challenging to determine these parameters independently. Our analytical model for He II proximity zones should also provide a useful description of the properties of H I proximity zones around quasars at $z \simeq 6-7$., Comment: 26 pages, 18 figures, accepted to ApJ

The ionising continuum from active galactic nuclei (AGN) is fundamental for interpreting their broad emission lines and understanding their impact on the surrounding gas. Furthermore, it provides hints on how matter accretes onto supermassive black holes. Using HST's Wide Field Camera 3 we have constructed the first stacked ultraviolet (rest-frame wavelengths 600-2500\AA) spectrum of 53 luminous quasars at z=2.4, with a state-of-the-art correction for the intervening Lyman forest and Lyman continuum absorption. The continuum slope ($f_\nu \propto \nu^{\alpha_\nu}$) of the full sample shows a break at ~912\AA\ with spectral index $\alpha_\nu=-0.61\pm0.01$ at $\lambda>912$\AA\ and a softening at shorter wavelengths ($\alpha_\nu=-1.70 \pm 0.61$ at $\lambda\leq 912$\AA). Our analysis proves that a proper intergalactic medium absorption correction is required to establish the intrinsic continuum emission of quasars. We interpret our average ultraviolet spectrum in the context of photoionisation, accretion disk models, and quasar contribution to the ultraviolet background. We find that observed broad line ratios are consistent with those predicted assuming an ionising slope of $\alpha_\mathrm{ion}=$-2.0, similar to the observed ionising spectrum in the same wavelength range. The continuum break and softening are consistent with accretion disk plus X-ray corona models when black hole spin is taken into account. Our spectral energy distribution yields a 30% increase to previous estimates of the specific quasar emissivity, such that quasars may contribute significantly to the total specific Lyman limit emissivity estimated from the Ly$\alpha$ forest at z<3.2., Comment: 19 pagine, 13 figure, replaced with version that matches the one published in MNRAS

We present a first study of the effect of local photoionising radiation on gas cooling in smoothed particle hydrodynamics simulations of galaxy formation. We explore the combined effect of ionising radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm used for gravity, so it is computationally efficient and well resolved. The method foregoes calculating absorption and scattering in favour of a constant escape fraction for young stars to keep the calculation efficient enough to simulate the entire evolution of a galaxy in a cosmological context to the present day. This allows us to quantify the effect of the local photoionisation feedback through the whole history of a galaxy`s formation. The simulation of a Milky Way like galaxy using the local photoionisation model forms ~ 40 % less stars than a simulation that only includes a standard uniform background UV field. The local photoionisation model decreases star formation by increasing the cooling time of the gas in the halo and increasing the equilibrium temperature of dense gas in the disc. Coupling the local radiation field to gas cooling from the halo provides a preventive feedback mechanism which keeps the central disc light and produces slowly rising rotation curves without resorting to extreme feedback mechanisms. These preliminary results indicate that the effect of local photoionising sources is significant and should not be ignored in models of galaxy formation., Comment: Accepted for Publication in MNRAS, 13 pages, 13 figures

The fraction of AGN luminosity obscured by dust and re-emitted in the mid-IR is critical for understanding AGN evolution, unification, and parsec-scale AGN physics. For unobscured (Type-1) AGN, where we have a direct view of the accretion disk, the dust covering factor can be measured by computing the ratio of re-processed mid-IR emission to intrinsic nuclear bolometric luminosity. We use this technique to estimate the obscured AGN fraction as a function of luminosity and redshift for 513 Type-1 AGN from the XMM-COSMOS survey. The re-processed and intrinsic luminosities are computed by fitting the 18-band COSMOS photometry with a custom SED-fitting code, which jointly models emission from: hot-dust in the AGN torus, the accretion disk, and the host-galaxy. We find a relatively shallow decrease of the luminosity ratio as a function of Lbol, which we interpret as a corresponding decrease in the obscured fraction. In the context of the receding torus model, where dust sublimation reduces the covering factor of more luminous AGN, our measurements require a torus height which increases with luminosity as h ~ Lbol^{0.3-0.4}. Our obscured fraction-luminosity relation agrees with determinations from SDSS censuses of Type-1 and Type-2 quasars, and favors a torus optically thin to mid-IR radiation. We find a much weaker dependence of obscured fraction on 2-10 keV luminosity than previous determinations from X-ray surveys, and argue that X-ray surveys miss a significant population of highly obscured Compton-thick AGN. Our analysis shows no clear evidence for evolution of obscured fraction with redshift., Comment: 33 pages, 24 figures, ApJ accepted

We report the discovery of a multiply lensed Ly Alpha (Lya) emitter at z = 3.90 behind the massive galaxy cluster WARPS J1415.1+3612 at z = 1.026. Images taken by the Hubble Space Telescope(HST) using ACS reveal a complex lensing system that produces a prominent, highly magnified arc and a triplet of smaller arcs grouped tightly around a spectroscopically confirmed cluster member. Spectroscopic observations using FOCAS on Subaru confirm strong Lya emission in the source galaxy and provide redshifts for more than 21 cluster members, from which we obtain a velocity dispersion of 807+/-185 km/s. Assuming a singular isothermal sphere profile, the mass within the Einstein ring (7.13+/-0.38") corresponds to a central velocity dispersion of 686+15-19 km/s for the cluster, consistent with the value estimated from cluster member redshifts. Our mass profile estimate from combining strong lensing and dynamical analyses is in good agreement with both X-ray and weak lensing results., Comment: 5 pages, 2 figures, accepted by ApJL. See http://www.supernova.lbl.gov/ for additional information pertaining to the HST Cluster SN Survey. (This added URL is the only change in this version.)

Gravitational lensing is a powerful tool for the study of the distribution of dark matter in the Universe. The cold-dark-matter model of the formation of large-scale structures predicts the existence of quasars gravitationally lensed by concentrations of dark matter so massive that the quasar images would be split by over 7 arcsec. Numerous searches for large-separation lensed quasars have, however, been unsuccessful. All of the roughly 70 lensed quasars known, including the first lensed quasar discovered, have smaller separations that can be explained in terms of galaxy-scale concentrations of baryonic matter. Although gravitationally lensed galaxies with large separations are known, quasars are more useful cosmological probes because of the simplicity of the resulting lens systems. Here we report the discovery of a lensed quasar, SDSS J1004+4112, which has a maximum separation between the components of 14.62 arcsec. Such a large separation means that the lensing object must be dominated by dark matter. Our results are fully consistent with theoretical expectations based on the cold-dark-matter model., Comment: 10 pages, 3 figures, to appear in the 18th&25th Dec issue of Nature (Letters to Nature)

We study the recently discovered gravitational lens SDSS J1004+4112, the first quasar lensed by a cluster of galaxies. It consists of four images with a maximum separation of 14.62''. The system has been confirmed as a lensed quasar at z=1.734 on the basis of deep imaging and spectroscopic follow-up observations. We present color-magnitude relations for galaxies near the lens plus spectroscopy of three central cluster members, which unambiguously confirm that a cluster at z=0.68 is responsible for the large image separation. We find a wide range of lens models consistent with the data, but they suggest four general conclusions: (1) the brightest cluster galaxy and the center of the cluster potential well appear to be offset by several kpc; (2) the cluster mass distribution must be elongated in the North--South direction, which is consistent with the observed distribution of cluster galaxies; (3) the inference of a large tidal shear (~0.2) suggests significant substructure in the cluster; and (4) enormous uncertainty in the predicted time delays between the images means that measuring the delays would greatly improve constraints on the models. We also compute the probability of such large separation lensing in the SDSS quasar sample, on the basis of the CDM model. The lack of large separation lenses in previous surveys and the discovery of one in SDSS together imply a mass fluctuation normalization \sigma_8=1.0^{+0.4}_{-0.2} (95% CL), if cluster dark matter halos have an inner slope -1.5. Shallower profiles would require higher values of \sigma_8. Although the statistical conclusion might be somewhat dependent on the degree of the complexity of the lens potential, the discovery is consistent with the predictions of the abundance of cluster-scale halos in the CDM scenario. (Abridged), Comment: 21 pages, 24 figures, 5 tables, accepted for publication in ApJ

The number of mass clusters and their distribution in redshift are very sensitive to the density of matter Omega_m and the equation of state of dark energy w. Using weak lens gravitational tomography one can detect clusters of dark matter, weigh them, image their projected mass distribution, and determine their 3-D location. The degeneracy curve in the Omega_m - w plane is nearly orthogonal to that from CMB or SN measurements. Thus, a combination of CMB data with weak lens tomography of clusters can yield precision measurements of Omega_m and w, independently of the SN observations. The Large Synoptic Survey Telescope (LSST) will repeatedly survey 30,000 square degrees of the sky in multiple wavelengths. LSST will create a 3-D tomographic assay of mass overdensities back to half the age of the universe by measuring the shear and color-redshift of billions of high redshift galaxies. By simultaneously measuring several functions of cosmic shear and mass cluster abundance, LSST will provide a number of independent constraints on the dark energy density and the equation of state., Comment: 9 pages, 6 figures, for Proc. 5th International UCLA Symposium on Sources and Detection of Dark Matter, Marina del Rey, February 2002, ed. D. Cline