Search

We present NOrthern Extended Millimeter Array observations of CO(8-7), (9-8), and (10-9) lines, as well as the underlying continuum for two far-infrared luminous quasars: SDSS J2054-0005 at $\rm z=6.0389$ and SDSS J0129-0035 at $\rm z=5.7788$. Both quasars were previously detected in CO (2-1) and (6-5) transitions, making them candidates for studying the CO Spectral Line Energy Distribution (SLED) of quasars at $z \sim 6$. Utilizing the radiative transfer code CLOUDY, we fit the CO SLED with two heating mechanisms, including the photo-dissociation region (PDR) and X-ray-dominated region (XDR) for both objects. The CO SLEDs can be fitted by either a dense PDR component with an extremely strong far-ultraviolet radiation field (gas density $ n_{\rm H} \sim 10^6 \, \rm cm^{-3}$ and field strength $G_0 \gtrsim 10^6$) or a two-component model including a PDR and an XDR. However, the line ratios, including \tir and previous \cii and \ci measurements, argue against a very high PDR radiation field strength. Thus, the results prefer a PDR+XDR origin for the CO SLED. The excitation of the high-J CO lines in both objects is likely dominated by the central AGN. We then check the CO (9-8)-to-(6-5) line luminosity ratio $r_{96}$ for all $z \sim 6$ quasars with available CO SLEDs (seven in total) and find that there are no clear correlations between $r_{96}$ and both \fir and the AGN UV luminosities. This further demonstrates the complexity of the CO excitation powered by both the AGN and nuclear star formation in these young quasar host galaxies., Comment: 19 pages, 9 figures, accepted for publication in ApJ

ASPIRE (A SPectroscopic survey of bIased halos in the Reionization Era) is a quasar legacy survey primarily using JWST to target a sample of 25 $z>6$ quasars with NIRCam slitless spectroscopy and imaging. The first study in this series found evidence of a strong overdensity of galaxies around J0305-3150, a luminous quasar at $z=6.61$, within a single NIRCam pointing obtained in JWST Cycle 1. Here, we present the first results of a JWST Cycle 2 mosaic that covers 35 arcmin$^2$ with NIRCam imaging/WFSS of the same field to investigate the spatial extent of the putative protocluster. The F356W grism data targets [OIII]+H$\beta$ at $5.310$ cMpc., Comment: submitted to ApJ june 2024. 18 pages, 7 figures. companion paper II by Champagne+24b is also available on arXiv today

Understanding the rapid formation of supermassive black holes (SMBHs) in the early universe requires an understanding of how stellar mass grows in the host galaxies. Here, we perform an analysis of rest-frame optical spectra and imaging from JWST of two quasar host galaxies at z>6 which exhibit Balmer absorption lines. These features in the stellar continuum indicate a lack of young stars, similar to low-redshift post-starburst galaxies whose star formation was recently quenched. We find that the stellar mass (log(M_* / M_sun) > 10.6) of each quasar host grew in a starburst episode at redshift 7 or 8. One of the targets exhibits little ongoing star formation, as evidenced by the photometric signature of the Balmer break and a lack of spatially resolved H-alpha emission, placing it well below the star formation main sequence at z = 6. The other galaxy is transitioning to a quiescent phase; together, the two galaxies represent the most distant massive post-starburst galaxies known. The maturity of these two galaxies is further supported by the stellar velocity dispersions of their host galaxies, placing them slightly above the upper end of the local M_BH - sigma_* relation. The properties of our two post-starburst galaxies, each hosting an active SMBH with log(M_BH / M_sun) > 9, suggests that black holes played a major role in shaping the formation of the first massive galaxies in the Universe., Comment: 24 pages, 7 figures, submitted to a Nature journal

We present millimeter observations of the host galaxy of the most distant blazar known, VLASSJ041009.05-013919.88 (hereafter J0410-0139) at z=7, using ALMA and NOEMA observations. The ALMA data reveal a 2e42 erg/s [CII] 158um emission line at z=6.9964 with a [CII]-inferred star-formation rate of 58 Msun/yr. We estimate a dynamical mass of 4.6e9 Msun, implying a black hole mass to host a dynamical mass ratio of 0.15. The 238 GHz continuum (rest-frame IR) decreased by ~33% from the NOEMA to the ALMA observations taken ~10 months apart. The VLA 3-10 GHz radio flux densities showed a ~37% decrease in a similar time frame, suggesting a causal connection. At face value, J0410-0139 would have the lowest [CII]-to-IR luminosity ratio of a z>5.7 quasar reported to date (~1e-4). However, if only <20% of the measured IR luminosity were due to thermal emission from dust, the [CII]-to-IR luminosity ratio would be typical of (U)LIRGS, and the star formation rates derived from [CII] and IR luminosities would be consistent. These results provide further evidence that synchrotron emission significantly contributes to the observed rest-frame IR emission of J0410-0139, similar to what has been reported in some radio-loud AGN at z<1., Comment: Version after addressing referee report. See also companion paper: arXiv:2407.07236

Relativistic jets are thought to play a crucial role in the formation of massive galaxies and supermassive black holes. Here we report multi-wavelength and multi-epoch observations of the quasar VLASSJ0410-0139 at redshift z=7, powered by a 7e8 solar-mass black hole. Its radio variability, X-ray properties, and compact radio emission on parsec scales reveal that J0410-0139 is a blazar with a relativistic jet aligned with our line of sight. This blazar's existence implies that many more similar (unaligned) jetted sources must exist at z=7. One scenario is that we observe an intrinsically low-power radio jet, but we see it at high luminosity due to relativistic beaming effects. In this case, a large fraction (>80%) of the UV bright quasars must have a similar jet to match the number density expected from the UV quasar luminosity function. These jets can enhance the growth of supermassive black holes and substantially affect their host galaxies. However, the implications would be even more severe if the quasar belongs to the top 10% radio luminous quasars, as measured if the beaming enhancement is less than a factor of 10-15. In this scenario, there should be hundreds to thousands of radio-quiet quasars at z=7 with intrinsic properties similar to J0410-0139 -- in strong tension with the number density of bright quasars derived from their UV luminosity function. To reconcile these results, most black hole growth at z=7 must happen in an obscured phase, as some models predict. The existence of supermassive black holes in the epoch of reionization is facilitated by significant jet-enhanced or obscured super-Eddington accretion., Comment: Submitted

We present Atacama Large Millimetre/submillimetre Array (ALMA) observations of the [O$_{\rm III}$] 88 $\mu$m emission of a sample of thirteen galaxies at $z$ = 6 to 7.6 selected as [C$_{\rm II}$]-emitting companion sources of quasars. To disentangle the origins of the luminous Oxygen line in the $z$ > 6 Universe, we looked at emission-line galaxies that are selected through an excellent star-formation tracer [C$_{\rm II}$] with star-formation rates between 9 and 162 M$_{\odot}$/yr. Direct observations reveal [O$_{\rm III}$] emission in just a single galaxy (L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ = 2.3), and a stacked image shows no [O$_{\rm III}$] detection, providing deep upper limits on the L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ ratios in the $z > 6$ Universe (L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ < 1.2 at 3${\sigma}$). While the fidelity of this sample is high, no obvious optical/near-infrared counterpart is seen in the JWST imaging available for four galaxies. Additionally accounting for low-redshift CO emitters, line stacking shows that our sample-wide result remains robust: The enhanced L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ reported in the first billion years of the Universe is likely due to the selection towards bright, blue Lyman-break galaxies with high surface star-formation rates or young stellar populations. The deep upper limit on the rest-frame 90 $\mu$m continuum emission (< 141 $\mu$Jy at 3${\sigma}$), implies a low average dust temperature (T$_{\rm dust}$ < 30K) and high dust mass (M$_{\rm dust}$ ~ 10$^8$ M$_{\odot}$). As more normal galaxies are explored in the early Universe, synergy between JWST and ALMA is fundamental to further investigate the ISM properties of the a broad range of samples of high-$z$ galaxies., Comment: 20 pages; 13 figures; accepted for publication in MNRAS

We present optical and near-infrared (NIR) spectroscopic observations for a sample of $45$ quasars at $6.50 < z \leq 7.64$ with absolute magnitudes at $1450$ {\AA} in the range $-28.82 \leq M_{1450} \leq -24.13$ and their composite spectrum. The median redshift and $M_{1450}$ of the quasars in the sample are $z_{\rm{median}}=6.71$ and $M_{1450,\rm{median}} \simeq -26.1$, respectively. The NIR spectra are taken with echelle spectrographs, complemented with additional data from optical long slit instruments, and then reduced consistently using the open-source Python-based spectroscopic data reduction pipeline PypeIt. The median value of the mean signal-to-noise ratios of the spectra in J, H, and K band (median $\langle \rm{SNR}_{\lambda} \rangle$) is: median $\langle \rm{SNR}_{J} \rangle=9.7$, median $\langle \rm{SNR}_{H} \rangle=10.3$, and median $\langle \rm{SNR}_{K} \rangle=11.7$; demonstrating the good data quality. This work presents the largest medium/moderate-resolution sample of quasars at $z>6.5$ from ground-based instruments. Its homogeneity and reproducibility make it ideally suited for several scientific goals, i.e., the study of the quasar proximity zones and damping wings, the Ly$\alpha$ forest, the intergalactic medium's metal content, as well as other properties such as the distribution of SMBH masses and Eddington ratios. Our composite spectrum is compared to others at both high and low-$z$ from the literature, showing differences in the strengths of many emission lines, probably due to differences in luminosity among the samples, but a consistent continuum slope, which proves that the same spectral features are preserved in quasars at different redshift ranges., Comment: 18 pages, 10 figures, submitted to MNRAS

We present JWST/NIRSpec Integral Field Spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at $z=6.2342$ caught as its host galaxy interacts with companion galaxies. We detect spatially extended emission of several emission lines (H$\alpha$, H$\beta$, [OIII], [NII], [SII], HeII), which we use to study the properties of the ionized phase of the interstellar medium: the source and hardness of the photoionizing radiation field, metallicity, dust reddening, electron density and temperature, and star formation. We also marginally detect continuum starlight emission associated with the companion sources. We find that at least two independent satellite galaxies are part of the system. While the quasar host appears highly enriched and obscured, with AGN-like photoionization conditions, the western companion shows minimal dust extinction, low metallicity ($Z\sim0.4$ Z$_\odot$), and star-formation driven photoionization. The eastern companion shows higher extinction and metallicity ($Z\sim0.8$ Z$_\odot$) compared to the western companion, and it is at least partially photoionized by the nearby quasar. We do not find any indication of AGN in the companion sources. Our study shows that while the quasar host galaxy is already very massive ($M_{\rm dyn}>10^{11}$ M$_\odot$), it is still rapidly building up by accreting two relatively massive ($M_{\rm star}\sim 10^{10}$ M$_\odot$) companion sources. This dataset showcases the power of JWST in exposing the build-up of massive galaxies in the first Gyr of the Universe., Comment: 15 pages, 16 figures. Accepted for publication in A&A

The [CII] 158 $\mu$m emission line and the underlying far-infrared (FIR) dust continuum are important tracers for studying star formation and kinematic properties of early galaxies. We present a survey of the [CII] emission lines and FIR continua of 31 luminous quasars at $z>6.5$ using the Atacama Large Millimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) at sub-arcsec resolution. This survey more than doubles the number of quasars with [CII] and FIR observations at these redshifts and enables statistical studies of quasar host galaxies deep into the epoch of reionization. We detect [CII] emission in 27 quasar hosts with a luminosity range of $L_{\rm [CII]}=(0.3-5.5)\times10^9~L_\odot$ and detect the FIR continuum of 28 quasar hosts with a luminosity range of $L_{\rm FIR}=(0.5-13.0)\times10^{12}~L_\odot$. Both $L_{\rm [CII]}$ and $L_{\rm FIR}$ are correlated ($\rho\simeq0.4$) with the quasar bolometric luminosity, albeit with substantial scatter. The quasar hosts detected by ALMA are clearly resolved with a median diameter of $\sim$5 kpc. About 40% of the quasar host galaxies show a velocity gradient in [CII] emission, while the rest show either dispersion-dominated or disturbed kinematics. Basic estimates of the dynamical masses of the rotation-dominated host galaxies yield $M_{\rm dyn}=(0.1-7.5)\times10^{11}~M_\odot$. Considering our findings alongside those of literature studies, we found that the ratio between $M_{\rm BH}$ and $M_{\rm dyn}$ is about ten times higher than that of local $M_{\rm BH}-M_{\rm dyn}$ relation on average but with substantial scatter (the ratio difference ranging from $\sim$0.6 to 60) and large uncertainties., Comment: accepted for publication in ApJ

Theoretical models predict that z~6 quasars are hosted in the most massive halos of the underlying dark matter distribution and thus would be immersed in protoclusters of galaxies. However, observations report inconclusive results. We investigate the 1.1 pMpc2 environment of the z = 7.54 luminous quasar ULAS J1342+0928. We search for Lyman-break galaxy candidates (LBG) using deep imaging from the Hubble Space Telescope (HST) in the ACS/F814W, WFC3/F105W/F125W bands, and Spitzer/IRAC at 3.6 $\mu$m and 4.5 $\mu$m. We report a zphot = $7.69^{+0.33}_{-0.23}$ LBG with magF125W = 26.41 at 223 projected-pkpc from the quasar. We find no HST counterpart to one [CII]-emitter previously found with ALMA at 27 projected-pkpc and $z[CII]=7.5341\pm0.0009$ (Venemans et al. 2020). We estimate the completeness of our LBG candidates using results from CANDELS/GOODS deep blank field searches sharing a similar filter setup. We find that >50% of the z~7.5 LBGs with magF125W >25.5 are missed due to the absence of a filter redward of the Lyman-break in F105W, hindering the UV color accuracy of the candidates. We conduct a QSO-LBG clustering analysis revealing a low LBG excess of $0.46^{+1.52}_{-0.08}$ in this quasar field, consistent with an average or low-density field. Consequently, this result does not present strong evidence of an LBG overdensity around ULAS J1342+0928. Furthermore, we identify two LBG candidates with a zphot matching a confirmed z=6.84 absorber along the line-of-sight to the quasar. All these galaxy candidates are excellent targets for follow-up observations with JWST and/or ALMA to confirm their redshift and physical properties., Comment: Accepted to ApJ, 16 pages, 8 figures

We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm \mu m$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ per spectral element). Based on the H$\beta$ line from the broad line region, we obtain an estimate of the black hole mass $M_{\rm BH,H\beta}\sim 2.7\times 10^{9}\ \rm M_{\odot}$. This value is within a factor $\lesssim 1.5$ of the H$\alpha$-based black hole mass from the same spectrum ($M_{\rm BH, H\alpha}\sim 1.93\times 10^{9}\ \rm M_{\odot}$) and is consistent with a previous estimate relying on the MgII $\lambda 2799$ ($M_{\rm BH, MgII}\sim 2.65\times 10^{9}\ \rm M_{\odot}$). All these $M_{\rm BH}$ are within the $\sim 0.5$ dex intrinsic scatter of the adopted mass calibrations. The high Eddington ratio of PJ308-21 $\lambda_{\rm Edd,H\beta}\sim 0.67$ ($\lambda_{\rm Edd,H\alpha}\sim 0.96$) is in line with the overall quasar population at $z \gtrsim 6$. The relative strengths of the [OIII], FeII and H$\beta$ lines are consistent with the empirical "Eigenvector 1" correlations as observed for low redshift quasars. We find evidence for blueshifted [OIII] $\lambda 5007$ emission with a velocity offset $\Delta v_{\rm [OIII]}=-1922\pm 39$ km s$^{-1}$ from the systemic velocity and a $\rm FWHM([OIII])=2776^{+75}_{-74}$ km s$^{-1}$. This may be the signature of an outflow from the nuclear region, despite the true values of $\Delta v_{\rm [OIII]}$ and $\rm FWHM([OIII])$ are likely more uncertain due to the blending with H$\beta$ and FeII lines. Our study demonstrates the unique capabilities of NIRSpec in capturing quasar spectra at cosmic dawn and studying their properties in unprecedented detail., Comment: 17 pages, 9 figures, accepted for publication in A&A

Probing the molecular gas reservoirs of z>~6 quasar (QSO) host galaxies is fundamental to understanding the coevolution of star formation and black hole growth in these extreme systems. Yet, there is still an inhomogeneous coverage of molecular gas tracers. To measure the average excitation and mass of the molecular gas reservoirs in the brightest z>6.5 QSO hosts, we combined new observations of CO(2-1) emission with existing observations of CO(6-5), CO(7-6), [C I], [C II], and dust-continuum emission. We reduced and analysed the VLA observations of CO(2-1) in three z=6.5-6.9 QSO hosts -- the highest redshift observations of CO(2-1) to date. By combining these with the nine z=5.7-6.4 QSO hosts for which CO(2-1) has already been observed, we studied the spread in molecular gas masses and CO excitation. Two of our three QSOs, were undetected in CO(2-1), implying more highly excited CO than in the z=6.4 QSO J1148+5251. We detected CO(2-1) at $5.1\sigma$ for our highest-redshift target, J2348-3054, yielding a molecular gas mass of $(1.2\pm0.2)\times 10^{10}\, \mathrm{M}_\odot$. This molecular gas mass is equivalent to the lower limit on the dynamical mass measured from resolved [C II] observations, implying little mass in stars or neutral gas within the [C II]-emitting region. On average, these QSO hosts have far higher CO(6-5)-, CO(7-6)-, and [C II] vs CO(2-1) line ratios than local AGN hosts; with a mean CO(6-5)-to-CO(1-0) line luminosity ratio of $r_{6,1}=0.9\pm0.2$. Our new CO(2-1) observations show that even at 780 Myr after the Big Bang, QSO host galaxies can already have molecular gas masses of $10^{10}$ M$_\odot$, consistent with a picture in which these z>6 QSOs reside in massive starbursts that are coevolving with the supermassive black holes. Our results imply the presence of extremely dense and warm molecular gas reservoirs illuminated by strong interstellar radiation fields., Comment: published in A&A

Protoclusters, the progenitors of galaxy clusters, trace large scale structures in the early Universe and are important to our understanding of structure formation and galaxy evolution. To date, only a handful of protoclusters have been identified in the Epoch of Reionization (EoR). As one of the rarest populations in the early Universe, distant quasars that host active supermassive black holes are thought to reside in the most massive dark matter halos at that cosmic epoch, and could thus potentially pinpoint some of the earliest protoclusters. In this letter, we report the discovery of a massive protocluster around a luminous quasar at $z=6.63$. This protocluster is anchored by the quasar, and includes three [CII] emitters at $z\sim6.63$, 12 spectroscopically confirmed Ly$\alpha$ emitters (LAEs) at $6.5430$ in its central region (i.e., $R\lesssim2$ cMpc). We estimate that this protocluster will collapse into a galaxy cluster with a mass of $6.9^{+1.2}_{-1.4}\times10^{15}~M_\odot$ at the current epoch, more massive than the most massive clusters known in the local Universe such as Coma. In the quasar vicinity, we discover a double-peaked LAE which implies that the quasar has a UV lifetime greater than 0.8 Myrs and has already ionized its surrounding intergalactic medium., Comment: Accepted for publication in ApJL

We present rest-frame optical spectroscopy using JWST/NIRSpec IFU for the radio galaxy TN J1338-1942 at z=4.1, one of the most luminous galaxies in the early Universe with powerful extended radio jets. Previous observations showed evidence for strong, large-scale outflows on the basis of its large (~150 kpc) halo detected in Ly-alpha, and high velocity [O II] emission features detected in ground-based IFU data. Our NIRSpec/IFU observations spatially resolve the emission line properties across the host galaxy in great detail. We find at least five concentrations of line emission, coinciding with discrete continuum features previously detected in imaging from HST and JWST, over an extent of ~2'' (~15 kpc). The spectral diagnostics enabled by NIRSpec unambiguously trace the activity of the obscured AGN plus interaction between the interstellar medium and the radio jet as the dominant mechanisms for the ionization state and kinematics of the gas in the system. A secondary region of very high ionization lies at roughly 5 kpc distance from the nucleus, and within the context of an expanding cocoon enveloping the radio lobe, this may be explained by strong shock-ionization of the entrained gas. However, it could also signal the presence of a second obscured AGN, which may also offer an explanation for an intriguing outflow feature seen perpendicular to the radio axis. The presence of a dual SMBH system in this galaxy would support that large galaxies in the early Universe quickly accumulated their mass through the merging of smaller units (each with their own SMBH), at the centers of large overdensities. The inferred black hole mass to stellar mass ratio of 0.01-0.1 for TNJ1338 points to a more rapid assembly of black holes compared to the stellar mass of galaxies at high redshifts, consistent with other recent observations., Comment: 17 pages, 11 figures, submitted to MNRAS, comments welcome!

We present observations using the NOrthern Extended Millimetre Array (NOEMA) of CO and $\rm H_{2}O$ emission lines, and the underlying dust continuum in two quasars at $z \sim 6$, i.e., P215-16 at $z$ = 5.78 and J1429+5447 at $z$ = 6.18. Notably, among all published CO SLEDs of quasars at $z \sim 6$, the two systems reveal the highest and the lowest CO level of excitation, respectively. Our radiative transfer modeling of the CO SLED of P215-16 suggests that the molecular gas heated by AGN could be a plausible origin for the high CO excitation. For J1429+5447, we obtain the first well-sampled CO SLED (from transitions from 2-1 to 10-9) of a radio-loud quasar at $z\gtrsim 6$. Analysis of the CO SLED suggests that a single photo-dissociation region (PDR) component could explain the CO excitation in the radio-loud quasar J1429+5447. This work highlights the utility of the CO SLED in uncovering the ISM properties in these young quasar-starburst systems at the highest redshift. The diversity of the CO SLEDs reveals the complexities in gas conditions and excitation mechanisms at their early evolutionary stage., Comment: Accepted for publication in Apj

Blazars are radio-loud Active Galactic Nuclei (AGN) whose jets have a very small angle to our line of sight. Observationally, the radio emission are mostly compact or a compact-core with a 1-sided jet. With 2.5$^{\prime\prime}$ resolution at 3 GHz, the Very Large Array Sky Survey (VLASS) enables us to resolve the structure of some blazar candidates in the sky north of Decl. $-40$ deg. We introduce an algorithm to classify radio sources as either blazar-like or non-blazar-like based on their morphology in the VLASS images. We apply our algorithm to three existing catalogs, including one of known blazars (Roma-BzCAT) and two of blazar candidates identified by WISE colors and radio emission (WIBRaLS, KDEBLLACS). We show that in all three catalogs, there are objects with morphology inconsistent with being blazars. Considering all the catalogs, more than 12% of the candidates are unlikely to be blazars, based on this analysis. Notably, we show that 3% of the Roma-BzCAT "confirmed" blazars could be a misclassification based on their VLASS morphology. The resulting table with all sources and their radio morphological classification is available online., Comment: 21 pages, 11 figures, accepted for publication in The Astrophysical Journal

We present 0.075 (~400 pc) resolution ALMA observations of the [CII] and dust continuum emission from the host galaxy of the z=6.5406 quasar, P036+03. We find that the emission arises from a thin, rotating disk with an effective radius of 0.21" (1.1 kpc). The velocity dispersion of the disk is consistent with a constant value of 66.4+-1.0 km/s, yielding a scale height of 80+-30 pc. The [CII] velocity field reveals a distortion that we attribute to a warp in the disk. Modeling this warped disk yields an inclination estimate of 40.4+-1.3 degrees and a rotational velocity of 116+-3 km/s. The resulting dynamical mass estimate of (1.96+-0.10) x 10^10 Msun is lower than previous estimates, which strengthens the conclusion that the host galaxy is less massive than expected based on local scaling relations between the black hole mass and the host galaxy mass. Using archival MUSE Ly-alpha observations, we argue that counter-rotating halo gas could provide the torque needed to warp the disk. We further detect a region with excess (15-sigma) dust continuum emission, which is located 1.3 kpc northwest of the galaxy's center and is gravitationally unstable (Toomre-Q < 0.04). We posit this is a star-forming region whose formation was triggered by the warp, because the region is located within a part of the warped disk where gas can efficiently lose angular momentum. The combined ALMA and MUSE imaging provides a unique view of how gas interactions within the disk-halo interface can influence the growth of massive galaxies within the first billion years of the universe., Comment: Accepted for publication in the Astrophysical Journal

We present Atacama Large Millimeter/submillimeter Array (ALMA) [C II] 158 $\mu \rm{m}$ and dust continuum observations of the $z=6.79$ quasar J0109--3047 at a resolution of $0."045$ ($\sim$300 pc). The dust and [C II] emission are enclosed within a $\sim 500\, \rm{pc}$ radius, with the central beam ($r<144\ \rm{pc}$) accounting for $\sim$25\% (8\%) of the total continuum ([C II]) emission. The far--infrared luminosity density increases radially from $\sim$5 $\times 10^{11} L_\odot\ \rm{kpc}^{-2}$ to a central value of $\sim$70 $\times 10^{11} L_\odot\ \rm{kpc}^{-2}$ (SFRD $\sim$50-700 $M_\odot\ \rm{yr}^{-1}\ \rm{kpc}^{-2}$). The [C II] kinematics are dispersion-dominated with a constant velocity dispersion of $137 \pm 6 \,\rm{km\ s}^{-1}$. The constant dispersion implies that the underlying mass distribution is not centrally peaked, consistent with the expectations of a flat gas mass profile. The lack of an upturn in velocity dispersion within the central beam is inconsistent with a black hole mass greater than $M_{\rm{BH}}<6.5\times 10^{8}\ M_\odot\ (2\sigma$ level), unless highly fine-tuned changes in the ISM properties conspire to produce a decrease of the gas mass in the central beam comparable to the black hole mass. Our observations therefore imply either that a) the black hole is less massive than previously measured or b) the central peak of the far-infrared and [C II] emission are not tracing the location of the black hole, as suggested by the tentative offset between the near-infrared position of the quasar and the ALMA continuum emission., Comment: Accepted for publication in ApJ. 10 pages + appendices. Main results: Figures 1+4. Updated to match published version, references fixed

We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) galaxy redshift survey at 3-4 $\mu$m among JWST Cycle-1 programs and provide extensive legacy values for studying the formation of the earliest supermassive black holes (SMBHs), the assembly of galaxies, early metal enrichment, and cosmic reionization. In this first ASPIRE paper, we report the discovery of a filamentary structure traced by the luminous quasar J0305-3150 and ten [OIII] emitters at $z=6.6$. This structure has a 3D galaxy overdensity of $\delta_{\rm gal}=12.6$ over 637 cMpc$^3$, one of the most overdense structures known in the early universe, and could eventually evolve into a massive galaxy cluster. Together with existing VLT/MUSE and ALMA observations of this field, our JWST observations reveal that J0305-3150 traces a complex environment where both UV-bright and dusty galaxies are present, and indicate that the early evolution of galaxies around the quasar is not simultaneous. In addition, we discovered 31 [OIII] emitters in this field at other redshifts, $5.3

Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE)" program. Our JWST spectra cover the quasars' emission between rest frame $\sim$ 4100 and 5100 \r{A}. The profiles of these quasars' broad H$\beta$ emission lines span a FWHM from 3000 to 6000 $\rm{km~s^{-1}}$. The H$\beta$-based virial black hole (BH) masses, ranging from 0.6 to 2.1 billion solar masses, are generally consistent with their MgII-based BH masses. The new measurements based on the more reliable H$\beta$ tracer thus confirm the existence of billion solar-mass BHs in the reionization epoch. In the observed [OIII] $\lambda\lambda$4960,5008 doublets of these luminous quasars, broad components are more common than narrow core components ($\le~1200~\rm{km~s^{-1}}$), and only one quasar shows stronger narrow components than broad. Two quasars exhibit significantly broad and blueshifted [OIII] emission, thought to trace galactic-scale outflows, with median velocities of $-610~\rm{km~s^{-1}}$ and $-1430~\rm{km~s^{-1}}$ relative to the [CII] $158\,\mu$m line. All eight quasars show strong optical FeII emission, and follow the Eigenvector 1 relations defined by low-redshift quasars. The entire ASPIRE program will eventually cover 25 quasars and provide a statistical sample for the studies of the BHs and quasar spectral properties., Comment: 13 pages, 4 figures, accepted for publication in ApJL

The formation of the first supermassive black holes is expected to have occurred in some most pronounced matter and galaxy overdensities in the early universe. We have conducted a sub-mm wavelength continuum survey of 54 $z\sim6$ quasars using the Submillimeter Common-User Bolometre Array-2 (SCUBA2) on the James Clerk Maxwell Telescope (JCMT) to study the environments around $z \sim 6$ quasars. We identified 170 submillimeter galaxies (SMGs) with above 3.5$\sigma$ detections at 450 or 850 \um\, maps. Their FIR luminosities are 2.2 - 6.4 $\times$ 10$^{12} L_{\odot}$, and star formation rates are $\sim$ 400 - 1200 M$_{\odot}$ yr$^{-1}$. We also calculated the SMGs differential and cumulative number counts in a combined area of $\sim$ 620 arcmin$^2$. To a $4\sigma$ detection (at $\sim$ 5.5 mJy), SMGs overdensity is $0.68^{+0.21}_{-0.19}$($\pm0.19$), exceeding the blank field source counts by a factor of 1.68. We find that 13/54 quasars show overdensities (at $\sim$ 5.5 mJy) of $\delta_{SMG}\sim$ 1.5 - 5.4. The combined area of these 13 quasars exceeds the blank field counts with the overdensity to 5.5 mJy of \dsmg $\sim$ $2.46^{+0.64}_{-0.55}$($\pm0.25$) in the regions of $\sim$ 150 arcmin$^2$. However, the excess is insignificant on the bright end (e.g., 7.5 mJy). We also compare results with previous environmental studies of Lyman alpha emitters (LAEs) and Lyman-Break Galaxies (LBGs) on a similar scale. Our survey presents the first systematic study of the environment of quasars at $z\sim6$. The newly discovered SMGs provide essential candidates for follow-up spectroscopic observations to test whether they reside in the same large-scale structures as the quasars and search for protoclusters at an early epoch., Comment: 12 pages, 8 figures and 2 tables, accepted for publication in ApJ

Characterizing the physical conditions (density, temperature, ionization state, metallicity, etc) of the interstellar medium is critical to our understanding of the formation and evolution of galaxies. Here we present a multi-line study of the interstellar medium in the host galaxy of a quasar at z~6.4, i.e., when the universe was 840 Myr old. This galaxy is one of the most active and massive objects emerging from the dark ages, and therefore represents a benchmark for models of the early formation of massive galaxies. We used the Atacama Large Millimeter Array to target an ensemble of tracers of ionized, neutral, and molecular gas, namely the fine-structure lines: [OIII] 88$\mu$m, [NII] 122$\mu$m, [CII] 158$\mu$m, and [CI] 370$\mu$m and the rotational transitions of CO(7-6), CO(15-14), CO(16-15), and CO(19-18); OH 163.1$\mu$m and 163.4$\mu$m; and H$_2$O 3(0,3)-2(1,2), 3(3,1)-4(0,4), 3(3,1)-3(2,2), 4(0,4)-3(1,3), 4(3,2)-4(2,3). All the targeted fine-structure lines are detected, as are half of the targeted molecular transitions. By combining the associated line luminosities, the constraints on the dust temperature from the underlying continuum emission, and predictions from photoionization models of the interstellar medium, we find that the ionized phase accounts for about one third of the total gaseous mass budget, and is responsible for half of the total [CII] emission. It is characterized by high density (n~180 cm$^{-3}$), typical of HII regions. The spectral energy distribution of the photoionizing radiation is comparable to that emitted by B-type stars. Star formation also appears to drive the excitation of the molecular medium. We find marginal evidence for outflow-related shocks in the dense molecular phase, but not in other gas phases. This study showcases the power of multi-line investigations in unveiling the properties of the star-forming medium in galaxies at cosmic dawn., Comment: 16 pages, 14 figures. Accepted for publication in A&A

Feedback and outflows in galaxies that are associated with a quasar phase are expected to be pivotal in quenching the most massive galaxies. However, observations targeting the molecular outflow phase, which dominates both the mass and momentum and removes the immediate fuel for star formation, are limited in high-z QSO hosts. Massive quiescent galaxies found at z ~ 4 are predicted to have already quenched star formation by z ~ 5 and undergone their most intense growth at z > 6. Here, we present two ALMA detections of molecular outflows, traced by blue-shifted absorption of the OH 119 micron doublet, from a sample of three z > 6 infrared luminous QSO hosts: J2310+1855 and P183+05. OH 119 micron is also detected in emission in P183+05, and tentatively in the third source: P036+03. Using similar assumptions as for high-z Dusty Star-Forming Galaxy outflows, we find that our QSOs drive molecular outflows with comparable mass outflow rates, and that are comparably energetic except for J2310+1855's significantly lower outflow energy flux. We do not find evidence, nor require additional input from the central AGN to drive the molecular outflow in J2310+1855 but can not rule out an AGN contribution in P183+05 if a significant AGN contribution to L_FIR is assumed and/or if the outflow covering fraction is high (> 53%), which evidence from the literature suggests is unlikely in these sources. Differences observed in the blue-shifted absorption spectral properties may instead be caused by the QSO hosts' more compact dust continuum, limiting observations to lower altitude and more central regions of the outflow., Comment: Erratum: 4 pages, 3 figures, 1 table Original: 19 pages, 7 figures, 3 tables

The identification of bright quasars at z>6 enables detailed studies of supermassive black holes, massive galaxies, structure formation, and the state of the intergalactic medium within the first billion years after the Big Bang. We present the spectroscopic confirmation of 55 quasars at redshifts 5.6

We present the $z\!\approx\!6$ type-1 quasar luminosity function (QLF) based on the Pan-STARRS1 (PS1) quasar survey. The PS1 sample includes 125 quasars at $z\approx5.7-6.2$ with $-28\lesssim M_{1450}\lesssim-25$. Complemented by 48 fainter quasars from the SHELLQs survey, we evaluate the $z\approx6$ QLF over $-28\lesssim M_{1450}\lesssim-22$. Adopting a double power law with an exponential evolution of the quasar density ($\Phi(z)\propto10^{k(z-6)}$; $k=-0.7$), we use a maximum likelihood method to model our data. We find a break magnitude of $M^*=-26.38_{-0.60}^{+0.79}\,\text{mag}$, a faint end slope of $\alpha=-1.70_{-0.19}^{+0.29}$, and a steep bright end slope with $\beta=-3.84_{-1.21}^{+0.63}$. Based on our new QLF model we determine the quasar comoving spatial density at $z\!\approx\!6$ to be $n( M_{1450}<-26)=1.16_{-0.12}^{+0.13}\,\text{cGpc}^{-3}$. In comparison with the literature, we find the quasar density to evolve with a constant value of $k\approx-0.7$ from $z\approx7$ to $z\approx4$. Additionally, we derive an ionizing emissivity of $\epsilon_{912}(z=6) =7.23_{-1.02}^{+1.65}\times 10^{22}\,\text{erg}\,\text{s}^{-1}\text{Hz}^{-1}\text{cMpc}^{-3}$ based on the QLF measurement. Given standard assumptions and the recent measurement of the mean free path of Becker et al. (2021) at $z\approx6$ we calculate an HI photoionizing rate of $\Gamma_{\text{HI}}(z{=}6) \approx 6 \times10^{-16}\,\text{s}^{-1}$, strongly disfavoring a dominant role of quasars in hydrogen reionization., Comment: accepted by the Astrophysical Journal; see companion paper by Ba\~nados et al. (arXiv:2212.04452)

The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting their underlying, previously-undetected host galaxies. Here we report rest-frame optical images and spectroscopy of two HSC-SSP quasars at $z>6$ with JWST. Using NIRCam imaging at 3.6$\mu$m and 1.5$\mu$m and subtracting the light from the unresolved quasars, we find that the host galaxies are massive (stellar masses of $13\times$ and $3.4\times$ $10^{10}$ M$_{\odot}$, respectively), compact, and disk-like. NIRSpec medium-resolution spectroscopy shows stellar absorption lines in the more massive quasar, confirming the detection of the host. Velocity-broadened gas in the vicinity of these quasars enables measurements of their black hole masses ($1.4\times 10^9$ and $2.0\times$ $10^{8}$ M$_{\odot}$, respectively). Their location in the black hole mass - stellar mass plane is consistent with the distribution at low redshift, suggesting that the relation between black holes and their host galaxies was already in place less than a billion years after the Big Bang., Comment: Matched to the published Nature version of the article. 27 pages, 4 main figures, 1 table, 6 supplementary figures, 2 supplementary table

Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral at $z>7$, and largely ionized by $z\sim5.5$. However, most methods to measure the IGM neutral fraction are highly model-dependent, and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low ($\bar{x}_{\rm HI} \lesssim 10^{-3}$) or close to unity ($\bar{x}_{\rm HI}\sim 1$). In particular, the neutral fraction evolution of the IGM at the critical redshift range of $z=6-7$ is poorly constrained. We present new constraints on $\bar{x}_{\rm HI}$ at $z\sim5.1-6.8$, by analyzing deep optical spectra of $53$ quasars at $5.736.1$ are consistent with the redshift evolution of the neutral fraction of the IGM derived from the Planck 2018., Comment: 13 pages, 6 figures, 2 tables, accepted for publication in ApJ

We present measurements of black hole masses and Eddington ratios for a sample of 38 bright (M$_{1450}$ < -24.4 mag) quasars at 5.8 < z < 7.5, derived from VLT/X-shooter near-IR spectroscopy of their broad CIV and MgII emission lines. The black hole masses (on average M$_{BH}$ ~ 4.6 x 10$^9$ M$_{\odot}$) and accretion rates (with Eddington ratios ranging between 0.1 and 1.0) are broadly consistent with that of similarly luminous 0.3 < z < 2.3 quasars, but there is evidence for a mild increase in the median Eddington ratio going towards z > 6. Combined with deep ALMA observations of the [CII] 158 $\mu$m line from the quasar host galaxies and VLT/MUSE investigations of the extended Ly$\alpha$ halos, this study provides fundamental clues to models of the formation and growth of the first massive galaxies and black holes. Compared to local scaling relations, z > 5.7 black holes appear to be over-massive with respect to their host galaxies, and their accretion properties do not change with host galaxy morphology. Under the assumption that the kinematics of the T ~ 10$^4$ K gas, traced by the extended Ly$\alpha$ halos, are dominated by the gravitational potential of the dark matter halo, we report a similar relation between the black hole mass and circular velocity to the one reported for z ~ 0 galaxies. These results paint a picture where the first supermassive black holes reside in massive halos at z > 6 and lead the first stages of galaxy formation by rapidly growing in mass with a duty cycle of order unity. However, this duty cycle needs to drastically drop towards lower redshifts, while the host galaxies continue forming stars at a rate of hundreds of solar masses per year, sustained by the large reservoirs of cool gas surrounding them., Comment: Accepted for publication on AAS journals. 24 pages + appendices, 15 figures, 5 tables (including full list of z>5.7 quasars with MgII based black hole mass estimates). For access to the data and codes used in this work, please contact the authors

We present a comparison of the interstellar medium traced by [CII] (ALMA), and ionised halo gas traced by Lya (MUSE), in and around QSO host galaxies at z~6. To date, 18 QSOs at this redshift have been studied with both MUSE and high-resolution ALMA imaging; of these, 8 objects display a Lya halo. Using datacubes matched in velocity resolution, we compare and contrast the spatial and kinematic information of the Lya halos and the host galaxies' [CII] (and dust-continuum) emission. We find that the Lya halos extend typically 3-30 times beyond the interstellar medium of the host galaxies. The majority of the Lya halos do not show ordered motion in their velocity fields, whereas most of the [CII] velocity fields do. In those cases where a velocity gradient can be measured in Lya, the kinematics do not align with those derived from the [CII] emission. This implies that the Lya emission is not tracing the outskirts of a large rotating disk that is a simple extension of the central galaxy seen in [CII] emission. It rather suggests that the kinematics of the halo gas are decoupled from those of the central galaxy. Given the scattering nature of Lya, these results need to be confirmed with JWST IFU observations that can constrain the halo kinematics further using the non-resonant Ha line., Comment: Published in ApJ

There are only five radio-loud quasars currently known within 1 Gyr from the Big Bang ($z>6$) and the properties of their host galaxies have not been explored in detail. We present a NOrthern Extended Millimeter Array (NOEMA) survey of [CII] (158 $\mu$m) and underlying continuum emission of four $z>6$ radio-loud quasars, revealing their diverse properties. J0309+2717 ($z=6.10$) has a bright [CII] line and underlying continuum, implying a starburst with a star-formation rate SFR=$340-1200$ $M_\odot$ yr$^{-1}$. J1429+5447 ($z=6.18$) has a SFR=$520-870$ $M_{\odot}$yr$^{-1}$ and its [CII] profile is consistent with two Gaussians, which could be interpreted as a galaxy merger. J1427+3312 ($z=6.12$) has a moderate SFR=$30-90$ $M_\odot$ yr$^{-1}$. Notably, this is a broad absorption line quasar and we searched for the presence of high-velocity outflows in the host galaxy. Although the NOEMA data reveal a tentative broad component of the [CII] line as wide as $\sim$1400~km~s$^{-1}$, the sensitivity of our current data is not sufficient to confirm it. Finally, P172+18 ($z=6.82$) is undetected in both [CII] and continuum, implying a SFR$<22-40$ $M_{\odot}$yr$^{-1}$. The broad range of SFRs is similar to what is observed in radio-quiet quasars at similar redshifts. If radio jets do not significantly contribute to both [CII] and IR luminosities, this suggest no feedback from the jet on the star formation in the host galaxy., Comment: submitted to A&A

We present high spatial resolution (0.4", 2.2kpc) observations of the CO(6-5), CO(7-6) and [CI] lines and dust continuum emission from the interstellar medium in the host galaxy of the quasar J0305$-$3150 at $z=6.6$. These, together with archival [CII] data at comparable spatial resolution, enable studies of the spatial distribution and kinematics between the ISM in different phases. When comparing the radial profiles of CO, [CII] and the dust continuum, we find that the CO and dust continuum exhibit similar spatial distributions, both of which are less extended than the [CII], indicating that the CO and dust continuum are tracing the same gas component, while the [CII] is tracing a more extended one. In addition, we derive the radial profiles of the [CII]/CO, [CII]/far-infrared (FIR), CO/FIR, and dust continuum $S_{98.7 \rm GHz}/S_{258.1 \rm GHz}$ ratios. We find a decreasing $S_{98.7 \rm GHz}/S_{258.1 \rm GHz}$ ratio with radius, possibly indicating a decrease of dust optical depth with increasing radius. We also detect some of the ISM lines and continuum emission in the companion galaxies previously discovered in the field around J0305$-$3150. Through comparing the line-to-line and line-to-FIR ratios, we find no significant differences between the quasar and its companion galaxies., Comment: 12 pages, 6 figures, accepted for publication in ApJ

We investigate the molecular gas content of z~6 quasar host galaxies using the IRAM / Northern Extended Millimeter Array. We target the 3mm dust continuum, and the line emission from CO(6-5), CO(7-6), [CI]2-1 in 10 infra-red-luminous quasars that have been previously studied in their 1mm dust continuum and [CII] line emission. We detect CO(7-6) at various degrees of significance in all the targeted sources, thus doubling the number of such detections in z~6 quasars. The 3mm to 1mm flux density ratios are consistent with a modified black body spectrum with a dust temperature $T_{dust}$~47 K and an optical depth $\tau_{\nu}$=0.2 at the [CII] frequency. Our study provides us with four independent ways to estimate the molecular gas mass, $M_{H2}$, in the targeted quasars. This allows us to set constraints on various parameters used in the derivation of molecular gas mass estimates, such as the mass per luminosity ratios $\alpha_{CO}$ and $\alpha_{[CII]}$, the gas-to-dust ratio $\delta_{g/d}$, and the carbon abundance [C]/H2. Leveraging either on the dust, CO, [CI], or [CII] emission yields mass estimates of the entire sample in the range $M_{H2}$~$10^{10}$ to $10^{11}$ M$_{\odot}$. We compare the observed luminosities of dust, [CII], [CI], and CO(7-6) with predictions from photo-dissociation and X-ray dominated regions. We find that the former provide better model fits to our data, assuming that the bulk of the emission arises from dense ($n_H>10^4$ cm$^{-3}$) clouds with a column density $N_{H}$~$10^{23}$ cm$^{-2}$, exposed to a radiation field with intensity $G_0$~$10^3$ (in Habing units). Our analysis reiterates the presence of massive reservoirs of molecular gas fueling star formation and nuclear accretion in $z$~6 quasar host galaxies. It also highlights the power of combined 3mm and 1mm observations for quantitative studies of the dense gas content in massive galaxies at cosmic dawn., Comment: 13 pages, 7 figures. Accepted for publication in A&A

We quantify galaxy overdensities around three high-redshift quasars with known [CII] 158um companions: PJ231-20 (z=6.59), PJ308-21 (z=6.24) and J0305-3150 (z=6.61). Recent SCUBA2 imaging revealed the presence of 17 submillimeter galaxies (SMG) with sky separations $0.7'< \theta < 2.4'$ from these three quasars. We present ALMA Band 6 follow-up observations of these SCUBA2-selected SMGs to confirm their nature and redshift. We also search for continuum-undetected [CII] emitters in the ALMA pointings and make use of archival MUSE observations to search for Lyman-$\alpha$ Emitters (LAE) associated with the quasars. While most of the SCUBA2-selected sources are detected with ALMA in the continuum, no [CII] line emission could be detected, indicating that they are not at the quasar redshifts. Based on the serendipitous detection of CO 7-6 and [CI] emission lines, we find that four SMGs in the field of PJ231-20 are at z~2.4, which is coincident with the redshift of a Mg II absorber in the quasar rest-frame UV spectrum. We report the discovery of 2 LAEs within <0.6 cMpc of PJ231-20 at the same redshift, indicating a LAE overdensity around this quasar. Taken together, these observations provide new constraints on the large-scale excess of Lyman-$\alpha$- and [CII]-emitting galaxies around z>6 quasars and suggest that only wide-field observations, such as MUSE, ALMA or JWST mosaics, can reveal a comprehensive picture of large-scale structure around quasars in the first billion years of the Universe., Comment: Accepted for publication in The Astrophysical Journal

We present 0".035 resolution (~200 pc) imaging of the 158 um [CII] line and the underlying dust continuum of the z=6.9 quasar J234833.34-305410.0. The 18 h ALMA observations reveal extremely compact emission (diameter ~1 kpc) that is consistent with a simple, almost face-on, rotation-supported disk with a significant velocity dispersion of ~160 km/s. The gas mass in just the central 200 pc is ~4x10^9 M_sun, about a factor two higher than that of the central supermassive black hole. Consequently we do not resolve the black hole's sphere of influence, and find no kinematic signature of the central supermassive black hole. Kinematic modeling of the [CII] line shows that the dynamical mass at large radii is consistent with the gas mass, leaving little room for a significant mass contribution by stars and/or dark matter. The Toomre-Q parameter is less than unity throughout the disk, and thus is conducive to star formation, consistent with the high infrared luminosity of the system. The dust in the central region is optically thick, at a temperature >132 K. Using standard scaling relations of dust heating by star formation, this implies an unprecedented high star formation rate density of >10^4 M_sun / yr / kpc^2. Such a high number can still be explained with the Eddington limit for star formation under certain assumptions, but could also imply that the central supermassive black hole contributes to the heating of the dust in the central 110 pc., Comment: Accepted for publication in The Astrophysical Journal

We report X-ray observations of the most distant known gravitationally lensed quasar, J0439+1634 at $z=6.52$, which is also a broad absorption line (BAL) quasar, using the XMM-Newton Observatory. With a 130 ks exposure, the quasar is significantly detected as a point source at the optical position with a total of 358$^{+19}_{-19}$ net counts using the EPIC instrument. By fitting a power-law plus Galactic absorption model to the observed spectra, we obtain a spectral slope of $\Gamma=1.45^{+0.10}_{-0.09}$. The derived optical-to-X-ray spectral slope $\alpha_{\rm{ox}}$ is $-2.07^{+0.01}_{-0.01}$, suggesting that the X-ray emission of J0439+1634 is weaker by a factor of 18 than the expectation based on its 2500 Angstrom luminosity and the average $\alpha_{\rm{ox}}$ vs. luminosity relationship. This is the first time that an X-ray weak BAL quasar at $z>6$ has been observed spectroscopically. Its X-ray weakness is consistent with the properties of BAL quasars at lower redshift. By fitting a model including an intrinsic absorption component, we obtain intrinsic column densities of $N_{\rm{H}}=2.8^{+0.7}_{-0.6}\times10^{23}\,\rm{cm}^{-2}$ and $N_{\rm{H}}= 4.3^{+1.8}_{-1.5}\times10^{23}\,\rm{cm}^{-2}$, assuming a fixed $\Gamma$ of 1.9 and a free $\Gamma$, respectively. The intrinsic rest-frame 2--10 keV luminosity is derived as $(9.4-15.1)\times10^{43}\,\rm{erg\,s}^{-1}$, after correcting for lensing magnification ($\mu=51.3$). The absorbed power-law model fitting indicates that J0439+1634 is the highest redshift obscured quasar with a direct measurement of the absorbing column density. The intrinsic high column density absorption can reduce the X-ray luminosity by a factor of $3-7$, which also indicates that this quasar could be a candidate of intrinsically X-ray weak quasar., Comment: 8 pages, 3 figures, 1 table; Accepted for publication in ApJL

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

The interactions between radio jets and the interstellar medium play a defining role for the co-evolution of central supermassive black holes and their host galaxies, but observational constraints on these feedback processes are still very limited at redshifts $z > 2$. We investigate the radio-loud quasar PSO J352.4034-15.3373 at $z \sim 6$ at the edge of the Epoch of Reionization. This quasar is among the most powerful radio emitters and the first one with direct evidence of extended radio jets ($\sim$1.6 kpc) at these high redshifts. We analyze NOEMA and ALMA millimeter data targeting the CO (6-5) and [CII] far-infrared emission lines, respectively, and the underlying continuum. The broad $440\pm 80$ km s$^{-1}$ and marginally resolved [CII] emission line yields a systemic redshift of $z\!=\!5.832 \pm 0.001$. Additionally, we report a strong 215 MHz radio continuum detection, $88\pm 7$ mJy, using the GMRT. This measurement significantly improves the constraints at the low-frequency end of the spectral energy distribution of this quasar. In contrast to what is typically observed in high-redshift radio-quiet quasars, we show that cold dust emission alone cannot reproduce the millimeter continuum measurements. This is evidence that the strong synchrotron emission from the quasar contributes substantially to the emission even at millimeter (far-infrared in the rest-frame) wavelengths. This quasar is an ideal system to probe the effects of radio jets during the formation of a massive galaxy within the first Gyr of the Universe., Comment: 14 pages, 5 figures

We report ALMA observations of the dust continuum and {\cii} emission of the host galaxy of J0439+1634, a gravitationally lensed quasar at $z=6.5$. Gravitational lensing boosts the source-plane resolution to $\sim0\farcs15$ $(\sim0.8\text{ kpc})$. The lensing model derived from the ALMA data is consistent with the fiducial model in \citet{fan19} based on {\it HST} imaging. The host galaxy of J0439+1634 can be well-fitted by a S\'ersic profile consistent with an exponential disk, both in the far-infrared (FIR) continuum and the {\cii} emission. The overall magnification is $4.53\pm0.05$ for the continuum and $3.44\pm0.05$ for the {\cii} line. The host galaxy of J0439+1634 is a compact ultra-luminous infrared galaxy, with a total star formation rate (SFR) of $1.56\times10^{3}M_\odot/\text{year}$ after correcting for lensing and an effective radius of $0.74$ kpc. The resolved regions in J0439+1634 follow the ``{\cii} deficit," where the {\cii}-to-FIR ratio decreases with FIR surface brightness. The reconstructed velocity field of J0439+1634 appears to be rotation-like. The maximum line-of-sight rotation velocity of 130 km/s at a radius of 2 kpc. However, our data cannot be fit by an axisymmetric thin rotating disk, and the inclination of the rotation axis, $i$, remains unconstrained. We estimate the dynamical mass of the host galaxy to be $7.9\sin^{-2}(i)\times10^{9}M_\odot$. J0439+1634 is likely to have a high gas-mass fraction and an oversized SMBH compared to local relations. The SFR of J0439+1634 reaches the maximum possible values, and the SFR surface density is close to the highest value seen in any star-forming galaxy currently known in the universe., Comment: 14 pages, 7 figures. Accepted by ApJ

We present a method of selecting quasars up to redshift $\approx$ 6 with random forests, a supervised machine learning method, applied to Pan-STARRS1 and WISE data. We find that, thanks to the increasing set of known quasars we can assemble a training set that enables supervised machine learning algorithms to become a competitive alternative to other methods up to this redshift. We present a candidate set for the redshift range 4.8 to 6.3 which includes the region around z = 5.5 where quasars are difficult to select due to photometric similarity to red and brown dwarfs. We demonstrate that under our survey restrictions we can reach a high completeness ($66 \pm 7 \%$ below redshift 5.6 / $83^{+6}_{-9}\%$ above redshift 5.6) while maintaining a high selection efficiency ($78^{+10}_{-8}\%$ / $94^{+5}_{-8}\%$). Our selection efficiency is estimated via a novel method based on the different distributions of quasars and contaminants on the sky. The final catalog of 515 candidates includes 225 known quasars. We predict the candidate catalog to contain an additional $148^{+41}_{-33}$ new quasars below redshift 5.6 and $45^{+5}_{-8}$ above and make the catalog publicly available. Spectroscopic follow-up observations of 37 candidates lead us to discover 20 new high redshift quasars (18 at $4.6\le z\le5.5$, 2 $z\sim5.7$). These observations are consistent with our predictions on efficiency. We argue that random forests can lead to higher completeness because our candidate set contains a number of objects that would be rejected by common color cuts, including one of the newly discovered redshift 5.7 quasars., Comment: Accepted by AJ

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 <2%), but shows no optical or near-infrared counterpart. Further follow-up is required to establish whether these two sources are physically associated., Comment: submitted to ApJ on Nov 29, 2020; accepted on Jan 31, 2021. See the companion paper by Momjian et al

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., Comment: 33 pages. Kinematic fitting code is available at https://github.com/mneeleman/qubefit

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., Comment: 4 figures, 1 table; accepted for publication in ApJL

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 analyze relative abundances and ionization conditions in a strong absorption system at z=6.84, seen in the spectrum of the z=7.54 background quasar ULAS J134208.10+092838.61. Singly ionized C, Si, Fe, Mg, and Al measurements are consistent with a warm neutral medium that is metal-poor but not chemically pristine. Firm non-detections of C IV and Si IV imply that any warm ionized phase of the IGM or CGM has not yet been enriched past the ultra-metal-poor regime (<0.001Z_{solar}), unlike lower redshift DLAs where these lines are nearly ubiquitous. Relative abundances of the heavy elements 794 Myr after the Big Bang resemble those of metal-poor damped Lyman Alpha systems at intermediate redshift and Milky Way halo stars, and show no evidence of enhanced [alpha/Fe], [C/Fe] or other signatures of yields dominated by massive stars. A detection of the CII* fine structure line reveals local sources of excitation from heating, beyond the level of photo-excitation supplied by the CMB. We estimate the total and [CII] cooling rates, balancing against ISM heating sources to develop an heuristic two-phase model of the neutral medium. The implied heating requires a surface density of star formation slightly exceeding that of the Milky Way but not at the level of a strong starburst. For a typical (assumed) NHI=10^{20.6}, an abundance of [Fe/H]=-2.2 matches the columns of species in the neutral phase. To remain undetected in C IV, a warm ionized phase would either need much lower [C/H]<-4.2 over an absorption path of 1 kpc, or else a very small absorption path (a few pc). While still speculative, these results suggest a significant reduction in heavy element enrichment outside of neutral star forming regions of the ISM, as would be expected in early stages of galactic chemical evolution., Comment: (21 pages, 8 figures)

We present a study of the [CII] 158micron line and underlying far-infrared (FIR) continuum emission of 27 quasar host galaxies at z~6, traced by the Atacama Large Millimeter/submillimeter Array at a spatial resolution of ~1 physical kpc. The [CII] emission in the bright, central regions of the quasars have sizes of 1.0-4.8kpc. The dust continuum emission is typically more compact than [CII]. We find that 13/27 quasars (approximately one-half) have companion galaxies in the field, at projected separations of 3-90kpc. The position of dust emission and the Gaia-corrected positions of the central accreting black holes are cospatial (typical offsets <0.1"). This suggests that the central black holes are located at the bottom of the gravitational wells of the dark matter halos in which the z>6 quasar hosts reside. Some outliers with offsets of ~500pc can be linked to disturbed morphologies, most likely due to ongoing or recent mergers. We find no correlation between the central brightness of the FIR emission and the bolometric luminosity of the accreting black hole. The FIR-derived star-formation rate densities (SFRDs) in the host galaxies peak at the galaxies' centers, at typical values between 100 and 1000 M_sun/yr/kpc^2. These values are below the Eddington limit for star formation, but similar to those found in local ultraluminous infrared galaxies. The SFRDs drop toward larger radii by an order of magnitude. Likewise, the [CII]/FIR luminosity ratios of the quasar hosts are lowest in their centers (few x10^-4) and increase by a factor of a few toward the galaxies' outskirts, consistent with resolved studies of lower-redshift sources., Comment: 24 pages, 11 figures. Published in ApJ

We study the interstellar medium in a sample of 27 high-redshift quasar host galaxies at z>6, using the [CII] 158um emission line and the underlying dust continuum observed at ~1kpc resolution with ALMA. By performing uv-plane spectral stacking of both the high and low spatial resolution data, we investigate the spatial and velocity extent of gas, and the size of the dust-emitting regions. We find that the average surface brightness profile of both the [CII] and the dust continuum emission can be described by a steep component within a radius of 2kpc, and a shallower component with a scale length of 2kpc, detected up to ~10kpc. The surface brightness of the extended emission drops below ~1% of the peak at radius of ~5kpc, beyond which it constitutes 10-20% of the total measured flux density. Although the central component of the dust continuum emission is more compact than that of the [CII] emission, the extended components have equivalent profiles. The observed extended components are consistent with those predicted by hydrodynamical simulations of galaxies with similar infrared luminosities, where the dust emission is powered by star formation. The [CII] spectrum measured in the mean uv-plane stacked data can be described by a single Gaussian, with no observable [CII] broad-line emission (velocities in excess of >500km/s), that would be indicative of outflows. Our findings suggest that we are probing the interstellar medium and associated star formation in the quasar host galaxies up to radii of 10kpc, whereas we find no evidence for halos or outflows., Comment: 32 pages, 20 figures, 1 table, accepted for publication in ApJ

Ultra-luminous quasars ($M_{1450} \leq -29$) provide us with a rare view into the nature of the most massive and most rapidly accreting supermassive black holes (SMBHs). Following the discovery of two of these extreme sources, J0341${+}$1720 ($M_{1450}=-29.56$, $z=3.71$) and J2125${-}$1719 ($M_{1450}=-29.39$, $z=3.90$), in the Extremely Luminous Quasar Survey (ELQS) and its extension to the Pan-STARRS\,1 footprint (PS-ELQS), we herein present an analysis of their rest-frame UV to optical spectroscopy. Both quasars harbor very massive SMBHs with $M_{\rm{BH}}=6.73_{-0.83}^{+0.75}\times10^{9}\,M_{\odot}$ and $M_{\rm{BH}}=5.45_{-0.55}^{+0.60}\times10^{9}\,M_{\odot}$, respectively, showing evidence of accretion above the Eddington limit ($L_{\rm{bol}}/L_{\rm{Edd}}=2.74_{-0.27}^{+0.39}$ and $L_{\rm{bol}}/L_{\rm{Edd}}=3.01_{-0.30}^{+0.34}$). NOEMA 3 millimeter observations of J0341${+}$1720 reveal a highly star-forming ($\rm{SFR}\approx1500\,M_{\odot}\,\rm{yr}^{-1}$), ultra-luminous infrared galaxy ($L_{\rm{TIR}}\approx1.0\times10^{13}\,L_{\odot}$) host, which, based on an estimate of its dynamical mass, is only ${\sim}30$ times more massive than the SMBH it harbors at its center. As examples of luminous super-Eddington accretion, these two quasars provide support for theories, which explain the existence of billion solar mass SMBHs ${\sim}700$ million years after the Big Bang by moderate super-Eddington growth from standard SMBH seeds., Comment: Accepted for publication in the Astrophysical Journal (ApJ)

We present a study of strong intervening absorption systems in the near-IR spectra of 31 luminous quasars at $z>5.7$. The quasar spectra were obtained with {\it Gemini} GNIRS that provide continuous wavelength coverage from $\sim$0.9 to $\sim$2.5 $\mu$m. We detect 32 strong Mg II doublet absorbers with rest-frame equivalent width $W_r$ ($\lambda2796$) $>1.0$ \AA at $2.2 < z < 6.0$. Each Mg II absorber is confirmed by at least two associated Fe II absorption lines in the rest-frame wavelength range of $\sim 1600-2600$ \AA. We find that the comoving line density ($dN/dX$) of the strong Fe II-bearing Mg II absorbers decreases towards higher redshift at $z>3$, consistent with previous studies. Compared with strong Mg II absorbers detected in damped Ly$\alpha$ systems at 2 $

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